Método de la regla falsa

Las primeras dos iteraciones de regula falsi. La curva roja muestra la función f; las líneas azules, las secantes.

Como en el método de bisección, se parte de un intervalo inicial [a₀,b₀] con f(a₀) y f(b₀) de signos opuestos, lo que garantiza que en su interior hay al menos una raíz (véase Teorema de Bolzano). El algoritmo va obteniendo sucesivamente en cada paso un intervalo más pequeño [a_k, b_k] que sigue incluyendo una raíz de la función f.

A partir de un intervalo [a_k, b_k] se calcula un punto interior c_k:

${\displaystyle c_{k}={\frac {f(b_{k})a_{k}-f(a_{k})b_{k}}{f(b_{k})-f(a_{k})}}}$

Dicho punto es la intersección de la recta que pasa por (a,f(a_k)) y (b,f(b_k)) con el eje de abscisas (igual a como se hace en el método de la secante).

Se evalúa entonces f(c_k). Si es suficientemente pequeño, c_k es la raíz buscada. Si no, el próximo intervalo [a_k+1, b_k+1] será:

• [a_k, c_k] si f(a_k) y f(c_k) tienen signos opuestos;
• [c_k, b_k] en caso contrario.

Se puede demostrar que bajo ciertas condiciones el método de la falsa posición tiene orden de convergencia lineal, por lo que suele converger más lentamente a la solución de la ecuación que el método de la secante, aunque el método de la falsa posición siempre converge a una solución de la ecuación.

El algoritmo tiene el inconveniente de que si la función es convexa o cóncava cerca de la solución, el extremo del intervalo más alejado de la solución queda fijo variando únicamente el más cercano, convergiendo muy lentamente.

Un ejemplo de este fenómeno se da en la función:

${\displaystyle f(x)=2x^{3}-4x^{2}+3x\,}$

comenzando con [−1,1]. El extremo izquierdo del intervalo, −1, nunca cambia; el extremo derecho se aproxima a 0 linealmente.

La situación en que el método falla es fácil de detectar (el mismo extremo del intervalo se elige dos veces seguidas) y fácil de corregir eligiendo un c_k diferente, como:

${\displaystyle c_{k}={\frac {{\frac {1}{2}}f(b_{k})a_{k}-f(a_{k})b_{k}}{{\frac {1}{2}}f(b_{k})-f(a_{k})}}}$

o

${\displaystyle c_{k}={\frac {f(b_{k})a_{k}-{\frac {1}{2}}f(a_{k})b_{k}}{f(b_{k})-{\frac {1}{2}}f(a_{k})}}}$

restándole peso a uno de los extremos del intervalo para obligar a que el próximo c_k ocurra de ese lado de la función.

El factor 2 usado arriba, garantiza una convergencia superlineal (asintóticamente, el algoritmo ejecuta dos pasos normales por cada paso modificado). Hay otras formas que dan incluso mejores tasas de convergencia. El ajuste mencionado arriba, y otras modificaciones similares se conocen como Algoritmo Illinois. Ford[1] resume y analiza las variantes superlineales del método regula falsi modificado. A juzgar por la bibliografía, estos métodos eran bien conocidos en los años 1970 pero han sido olvidados en los textos actuales.