Método de la secante
# Método de la secante

En análisis numérico el método de la secante es un método para encontrar los ceros de una función de forma iterativa.

Es una variación del método de Newton-Raphson donde en vez de calcular la derivada de la función en el punto de estudio, teniendo en mente la definición de derivada, se aproxima la pendiente a la recta que une la función evaluada en el punto de estudio y en el punto de la iteración anterior. Este método es de especial interés cuando el coste computacional de derivar la función de estudio y evaluarla es demasiado elevado, por lo que el método de Newton no resulta atractivo.

En otras palabras, el método de la secante es un algoritmo de la raíz de investigación que utiliza una serie de raíces de las líneas secantes para aproximar mejor la raíz de una función f. El método de la secante se puede considerar como una aproximación en diferencias finitas del método de Newton-Raphson. Sin embargo, este método fue desarrollado independientemente de este último.

## El método

El método se define por la relación de recurrencia:

    $x_{n+1}=x_{n}-{\frac {x_{n}-x_{n-1}}{f(x_{n})-f(x_{n-1})}}f(x_{n}).$

Como se puede ver, este método necesitará dos aproximaciones iniciales de la raíz para poder inducir una pendiente inicial.

## Derivación del método

El método se basa en obtener la ecuación de la recta que pasa por los puntos $(x_{n−1},$ $f(x_{n−1}))$ y $(x_n,$ $f(x_n)).$ A dicha recta se le llama secante por cortar la gráfica de la función. En la imagen de arriba a la derecha se toman los puntos iniciales $x_0$ y $x_1,$ se construye una línea por los puntos $(x_0,$ $f(x_0))$ y $(x_1,$ $f(x_1)).$ En forma punto-pendiente, esta línea tiene la ecuación mostrada anteriormente. Posteriormente se escoge como siguiente elemento de la relación de recurrencia, $x_{n+1},$ la intersección de la recta secante con el eje de abscisas obteniendo la fórmula, y un nuevo valor. Seguimos este proceso, hasta llegar a un nivel suficientemente alto de precisión (una diferencia lo suficientemente pequeñas entre $x_n$ y $x_{n-1}).$

## Convergencia

El orden de convergencia de este método, en un punto cercano a la solución, es $\varphi$ donde

    $\varphi ={\frac {1+{\sqrt {5}}}{2}}\approx 1,618$

es el número áureo, por lo que se trata de una convergencia superlineal inferior a la del método de Newton-Raphson. En caso de que la aproximación inicial sea demasiado lejana o la raíz no sea simple, este método no asegura la convergencia y tiene un comportamiento similar al de Newton-Raphson.

## Comparación con otros métodos de búsqueda de raíces

El método de bisección necesita de muchas iteraciones comparado con el método de la secante, ya que el proceso que este sigue es mucho más preciso que el de bisección, el cual solo divide por mitades sucesivamente hasta dar con un valor aproximado al real y por consecuente conlleva un número significativamente mayor de iteraciones.

El método de la regla falsa utiliza la misma fórmula que el método de la secante. Sin embargo, no se aplica la fórmula en $x_{n−1}$ y $x_n,$ como el método de la secante, pero en $x_n$ y en la última iteración $x_k$ tal que $f(x_k)$ y $f(x_n)$ tiene un signo diferente. Esto significa que el método de regla falsa siempre converge.

La fórmula de recurrencia del método de la secante se puede derivar de la fórmula para el método de Newton-Raphson:

    $x_{n}=x_{n-1}-{\frac {f(x_{n-1})}{f'(x_{n-1})}}.$

utilizando la aproximación de diferencias finitas:

    $f'(x_{n-1})\approx {\frac {f(x_{n-1})-f(x_{n-2})}{x_{n-1}-x_{n-2}}}.$

Si comparamos el método de Newton-Raphson con el método de la secante, vemos que el método de Newton-Raphson converge más rápido (para 2 en contra α ≈ 1,6). Sin embargo, el método de Newton-Raphson requiere la evaluación de ambos f y su derivada en cada paso, mientras que el método de la secante sólo requiere la evaluación de f. Por lo tanto, el método de la secante puede muy bien ser más rápido en la práctica.

## Ejercicio de ejemplo

Utilice el método de la secante para encontrar una raíz real de la ecuación polinomial: $F(x)=x^3+2x^2+10x-20=0.$

Utilizando la ecuación:

    $x_{n+1}=x_{n}-{\frac {x_{n}-x_{n-1}}{f(x_{n})-f(x_{n-1})}}f(x_{n}).$

Obtenemos:

    $x_{n+1}=x_{n}-{\frac {(x_{n}-x_{n-1})(x_{n}^{3}+2x_{n}^{2}+10x_{n}-20)}{(x_{n}^{3}+2x_{n}^{2}+10x_{n}-20)-((x_{n-1}^{3}+2x_{n-1}^{2}+10x_{n-1}-20))}}$

Y mediante $x_0=0$ y $x_1=1$ se calcula $x_2$

    $x_{2}=1-{\frac {(1-0)(1^{3}+2(1)^{2}+10(1)-20)}{(1^{3}+2(1)^{2}+10(1)-20)-(0^{3}+2(0)^{2}+10(0)-20)}}=1,53846$

Los valores posteriores son los siguientes:

| n | $x_n$ | $|x_{n+1}-x_{n}|$ |
| - | - | - | - | - |
| 0 | 0.00000 | |
| 1 | 1.00000 | 1.00000 |
| 2 | 1.53846 | 0.53846 |
| 3 | 1.35031 | 0.18815 |
| 4 | 1.36792 | 0.01761 |
| 5 | 1.36881 | 0.00090 |

Ahí tenemos el resultado, cuando $|x_{n+1}-x_{n}|$ ≤ $\epsilon =10^{-3}$

Comprobando el resultado graficando la función utilizando software obtenemos:

Si bien no se converge a la raíz tan rápido como resolviéndolo utilizando el método Newton-Raphson, la velocidad de convergencia no es tan lenta como resolviéndolo por el método de punto fijo; entonces se tiene para este ejemplo una velocidad de convergencia intermedia.

## Código en Java

Programa escrito en Java correspondiente al ejemplo f(x) = exp( x ) - ( 2 * x^3 )

## Código en Fortran 90

Programa escrito en Fortran 90 correspondiente al ejemplo f(x) = $x^3$ + $2x^2$ + 10x - 20

Para compilar en GNU/Linux con compilador de GNU, se escribe en una terminal:

## Código en Matlab

Programa escrito en Matlab para ejecutar el método de la secante.

Alternativa al código anterior.

## Código en Julia

Definición de la función:

function secant(f,a,b,rtol,maxIters)

iter = 0

while abs(b-a) > rtol*abs(b) && iter < maxIters

c,a = a,b

b = b + (b-c)/(f(c)/f(b)-1)

iter = iter + 1

end

return b

end

Utilización de la función aplicada al ejercicio propuesto:

secant( x-> x^3+ 2*x^2 +10*x -20, 1, 2, 1e-15, 10 )

1.3688081078213725

