Descripción
Método de Newton-Raphson para la búsqueda de un cero de una función univariable dada la función y su función derivada.
El método de Newton no garantiza encontrar una raiz de la función desde un punto inicial, y en caso de existir varias solo encuentra una dependiendo del punto inicial.
Esta función es idéntica a la función newton (que calcula numéricamente la deriva), pero más rápida y precisa
Cadena de entrada
sol_newton
Cadena de salida
sol_newton
Uso
sol_newton(<función>,<derivada>,<varname>,<punto_inicial>,<maxiter>[,<error>])
Parámetros | |||
---|---|---|---|
# | Parámetro | Descripción | Valor por defecto |
1 | función | Función de una variable a encontrar raiz | |
2 | derivada | Función derivada precalculada por el usuario (por supuesto, si la derivada es incorrecta el resultado será 'probablemente' incorrecto) | |
3 | varname | Variable de la que depende la función | |
4 | punto_inicial | Punto de inicio de la búsqueda (este valor influye para encontrar la solución y cual de ellas) | |
5 | maxiter | Máximo de iteraciones | |
6 | error | Margen de error | 1E-15 |
Ejemplos
sol_newton(x*ln(x)-1,ln(x)+1,x,1,100) ln(x)=x-1
RealDoble: 1.7632228343518965
sol_newton(x*ln(x)-1,ln(x)+1,x,0,100) ln(x)=x-1
FuncionException: <<<FuncionException>>> en funcion "sol_newton": solucion no encontrada --> sol_newton(Vector:[x*ln(x)-1,ln(x)+1,x,0,100])
sol_newton(x^2-1,2x,x,0.1,100) x2=1
RealDoble: 1
sol_newton(x^2-1,2x,x,-0.1,100) x2=1
RealDoble: -1
Véase también…
método de Newton, método de la bisección, método de la secante