Puede encontrar las soluciones numéricas para una amplia gama de funciones utilizando el método de Newton, siempre que conozca la derivada de su función. Si no, o la función no tiene una derivada en ese rango, puede usar el método Secante.
En el caso de las raíces, el método de Newton está bien y trataré de explicarlo aquí:
Imagine que desea calcular [math] \ sqrt {2} [/ math]. Entonces, en realidad estás tratando de resolver [matemáticas] x = \ sqrt {2} \ rightarrow x ^ 2 = 2 [/ matemáticas].
[matemáticas] x ^ 2 = 2 \ flecha derecha x ^ 2-2 = 0 [/ matemáticas]
Así que vamos a definir [matemáticas] f (x) = x ^ 2-2 [/ matemáticas] e intentamos encontrar su raíz con el método de Newton.
El método de Newton es el siguiente :
1. Adivina la raíz y llámalo [math] x_0 [/ math]
2. Siga el siguiente proceso hasta que obtenga su respuesta con la precisión que desea:
[matemáticas] x_ {n + 1} = x_n- {f (x_n) \ over f ‘(x_n)} [/ matemáticas]
Entonces, para [matemáticas] f (x) = x ^ 2-2 [/ matemáticas]:
[matemáticas] f ‘(x) = 2x [/ matemáticas]
El método de Newton en este caso sería [matemáticas] x_ {n + 1} = x_n- {x_n ^ 2-2 \ over 2x_n} [/ matemáticas]
Supongamos que nuestra primera suposición es 1, entonces [math] x_0 = 1 [/ math] y:
[matemáticas] x_1 = 1- {1 ^ 2-2 \ más de 2 \ veces 1} = 1.5 [/ matemáticas]
[matemáticas] x_2 = 1.5- {1.5 ^ 2-2 \ más de 2 \ veces 1.5} = 1.41666 [/ matemáticas]
[matemática] x_3 = 1.41666- {1.41666 ^ 2-2 \ over 2 \ times 1.41666} = 1.414215 [/ math]
Con solo 3 pasos hemos alcanzado una precisión de 5 puntos decimales.
- Un ejemplo para la raíz cúbica :
[matemáticas] x = \ sqrt [3] {2} \ flecha derecha x ^ 3 = 2 [/ matemáticas]
[matemáticas] f (x) = x ^ 3-2 [/ matemáticas]
[matemáticas] f ‘(x) = 3x ^ 2 [/ matemáticas]
[matemáticas] x_ {n + 1} = x_n- {x_n ^ 3-2 \ over 3x ^ 2_n} [/ matemáticas]
Deje [math] x_0 = 1 [/ math]
[matemáticas] x_1 = 1.33333333333333 [/ matemáticas]
[matemáticas] x_2 = 1.26388888888889 [/ matemáticas]
[matemáticas] x_3 = 1.25993349344998 [/ matemáticas]
[matemáticas] x_4 = 1.25992105001777 [/ matemáticas]
Con solo 4 pasos hemos alcanzado una precisión de 7 puntos decimales (Gracias Angie Segura por señalar el error que cometí en el proceso).
- Un ejemplo un poco más complicado para encontrar la solución para ([math] \ sin_x + \ cos_x = \ sqrt {x} [/ math])
[matemáticas] \ sin_x + \ cos_x = \ sqrt {x} \ rightarrow \ sin_x + \ cos_x- \ sqrt {x} = 0 [/ matemáticas]
[matemáticas] f (x) = \ sin_x + \ cos_x- \ sqrt {x} [/ matemáticas]
[matemáticas] f ‘(x) = \ cos_x- \ sin_x- {1 \ over 2 \ sqrt {x}} [/ matemáticas]
[matemáticas] x_ {n + 1} = x_n – {\ sin_x + cos_x- \ sqrt {x} \ over \ cos_x- \ sin_x- {1 \ over 2 \ sqrt {x}}} [/ math]
Deje [math] x_0 = 1 [/ math]
[matemáticas] x_1 = 1.47652048902019 [/ matemáticas]
[matemáticas] x_2 = 1.38098724894177 [/ matemáticas]
[matemáticas] x_3 = 1.37734243869716 [/ matemáticas]
[matemáticas] x_4 = 1.37733687713611 [/ matemáticas]
Con solo 4 pasos hemos alcanzado una precisión de 10 puntos decimales.
Para mayor precisión, puede continuar el proceso.
*** *** ***
Hay otra forma de calcular la enésima raíz. Como sabe, exp (x) e ln (x) están predefinidos en la mayoría de los compiladores e intérpretes que utilizan expansiones de la serie Taylor. Por lo tanto, puede usarlos para encontrar la raíz enésima para N real distinto de cero (no necesariamente entero) de la siguiente manera:
[matemáticas] \ sqrt [N] {x} = x ^ {{1 \ over N}} = e ^ {\ ln {x ^ {1 \ over N}}} = e ^ {{1 \ over N} \ ln x} = e ^ {\ ln x \ sobre N} [/ math]
Es obvio que este es un caso especial que solo funciona para calcular la enésima raíz. Para otras ecuaciones no lineales como mi tercer ejemplo, me temo que está atascado con el método de Newton y necesita elegir sus conjeturas iniciales para evitar los problemas mencionados (mínimos y máximos locales) por caso, o idear un algoritmo que si f ‘(x) = 0 en el proceso, comience el proceso desde arriba con otra aproximación inicial.
*** *** ***
Solo por mencionarlo, el método Secante utiliza 2 conjeturas iniciales y no necesita la derivada para su cálculo, pero converge más lentamente que el método de Newton.
El método es como el de Newton:
1. Tome 2 conjeturas iniciales.
2. Repita el siguiente proceso hasta llegar a la precisión requerida:
[matemáticas] x_n = x_ {n-1} -f (x_ {n-1}) {x_ {n-1} -x_ {n-2} \ over f (x_ {n-1}) – f (x_ {n-2})} [/ matemáticas]