Clase 4 2022 2 Newton Raphson

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 63

UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS

FACULTAD DE GEOLOGÍA, MINAS, METALÚRGICA,...


Escuela Académico-Profesional de Ingeniería Metalúrgica

Métodos Numéricos

Prof. Gladys Melgarejo Estremadoyro 2022-2


MÉTODO
DE
NEWTON RAPHSON
Es uno de los métodos numéricos más conocido y poderoso para encontrar
las raíces de 𝑓𝑓 𝑥𝑥 .

Supóngase 𝑓𝑓 ∈ 𝐶𝐶 2 [𝑎𝑎, 𝑏𝑏].

Sea 𝑥𝑥 ∈ [𝑎𝑎, 𝑏𝑏] una aproximación a “𝑝𝑝”; 𝑓𝑓 ′ 𝑥𝑥 ≠ 0 y ⃒𝑥𝑥 − 𝑝𝑝𝑝 es muy


pequeño.
Consideremos el polinomio de Taylor (De primer grado para 𝑓𝑓 𝑥𝑥 alrededor de 𝑥𝑥 ).

2
𝑥𝑥 − 𝑥𝑥
𝑓𝑓 𝑥𝑥 = 𝑓𝑓 𝑥𝑥 + 𝑥𝑥 − 𝑥𝑥 𝑓𝑓 ′ 𝑥𝑥 + 𝑓𝑓 ′′ 𝜉𝜉 𝑥𝑥
2!

ó
2
1
𝑥𝑥 − 𝑥𝑥 2
𝑓𝑓 𝑥𝑥 = 𝑓𝑓 𝑥𝑥 + 𝑥𝑥 − 𝑥𝑥 𝑓𝑓 𝑥𝑥 + 𝑓𝑓 𝜉𝜉 𝑥𝑥
2!

Donde 𝜉𝜉 esta entre 𝑥𝑥 y 𝑥𝑥. ( 𝜉𝜉 =Error).

Como 𝑝𝑝 es raíz de 𝑓𝑓 𝑥𝑥 es decir; 𝑓𝑓 𝑝𝑝 = 0 ,


Entonces:
2
1
𝑝𝑝 − 𝑥𝑥 2
𝑓𝑓 𝑝𝑝 = 𝑓𝑓 𝑥𝑥 + 𝑝𝑝 − 𝑥𝑥 𝑓𝑓 𝑥𝑥 + 𝑓𝑓 𝜉𝜉 𝑥𝑥
2!
2
1
𝑝𝑝 − 𝑥𝑥 2
0 = 𝑓𝑓 𝑥𝑥 + 𝑝𝑝 − 𝑥𝑥 𝑓𝑓 𝑥𝑥 + 𝑓𝑓 𝜉𝜉 𝑥𝑥
2!
Como ⃒𝑥𝑥 − 𝑝𝑝𝑝 es muy pequeño, entonces al derivar

⃒𝑥𝑥 − 𝑝𝑝 ⃒2 → 0 Es despreciable, luego:


1
0 ≅ 𝑓𝑓 𝑥𝑥 + 𝑝𝑝 − 𝑥𝑥 𝑓𝑓 𝑥𝑥
1
−𝑓𝑓 𝑥𝑥 ≅ 𝑝𝑝 − 𝑥𝑥 𝑓𝑓 𝑥𝑥

𝑓𝑓 𝑥𝑥
− 1
≅ 𝑝𝑝 − 𝑥𝑥
𝑓𝑓 𝑥𝑥

𝑓𝑓 𝑥𝑥 𝑓𝑓 𝑝𝑝 𝑛𝑛−1
∴ 𝑝𝑝 ≅ 𝑥𝑥 − 1
⇒ 𝑝𝑝 𝑛𝑛 ≅ 𝑝𝑝 𝑛𝑛−1 − 1 , 𝑛𝑛 ≥ 1
𝑓𝑓 𝑥𝑥 𝑓𝑓 𝑝𝑝 𝑛𝑛−1

La expresión anterior es una mejor aproximación a “𝑝𝑝” que 𝑥𝑥 mismo.


Entonces la fórmula de recurrencia es llamada formula de Newton en una forma
generalizada.
𝒇𝒇 𝒙𝒙𝒏𝒏−𝟏𝟏
𝒙𝒙 𝒏𝒏 ≅ 𝒙𝒙 𝒏𝒏−𝟏𝟏 − 𝟏𝟏 , 𝒏𝒏 ≥ 𝟏𝟏
𝒇𝒇 𝒙𝒙𝒏𝒏−𝟏𝟏

Donde la sucesión 𝑥𝑥𝑛𝑛 𝑛𝑛≥1 es una sucesión convergente hacia 𝑝𝑝

Tenemos que:

𝒇𝒇 𝒙𝒙𝒏𝒏
𝒙𝒙 𝒏𝒏+𝟏𝟏 = 𝒙𝒙 𝒏𝒏 − ′ , 𝒏𝒏 ≥ 𝟎𝟎
𝒇𝒇 𝒙𝒙𝒏𝒏

Nota

1) 𝑓𝑓 ′ 𝑥𝑥 ≠ 0 , pues la recta tangente que pasa por el punto (𝑥𝑥0 , 𝑓𝑓(𝑥𝑥0 )) , no debe ser paralelo al eje de
las X.
2) 𝐶𝐶 2 [𝑎𝑎, 𝑏𝑏].= Conjunto de todas las funciones que tienen segunda derivada, continua en x
ALGORITMO DE NEWTON RAPHSON

INGRESAR: 𝑥𝑥0 = Valor inicial


𝜉𝜉 = Error ó Tolerancia
𝑁𝑁 = Numero de iteraciones
𝑓𝑓 𝑥𝑥 = Proporcionar la función
𝑓𝑓′ 𝑥𝑥 = Proporcionar la derivada de la función
1) Hacer 𝑖𝑖 = 1
2) Mientras 𝑖𝑖 ≤ 𝑁𝑁 , seguir los paso
𝑓𝑓 𝑥𝑥𝑖𝑖
3) Hacer 𝑥𝑥𝑖𝑖+1 = 𝑥𝑥𝑖𝑖 −
𝑓𝑓′ 𝑥𝑥𝑖𝑖

4) Si ⃒𝑥𝑥𝑛𝑛+1 − 𝑥𝑥𝑛𝑛 ⃒ < 𝜉𝜉


5) Entonces :
6) Imprimir “Solución Aproximada ”
7) x = xi+1
8) Hacer i=i+1
9) Imprimir “Método no Converge a una raiz y terminar”
EJEMPLO

Se sabe que una aleación de níquel con cromo la densidad de dislocaciones


que puede presentar dicha estructura está determinada por la función
𝑓𝑓 𝑥𝑥 = 𝑥𝑥 3 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 + 7𝑥𝑥, determinar la densidad experimental por el Método de
Newton Raphson con un error estimado menor e igual a 0.01 y con 𝑥𝑥0 = 6.
(Considere 6 dígitos en la parte decimal)
• 𝑓𝑓 𝑥𝑥 = 𝑥𝑥 3 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 + 7𝑥𝑥
• 𝑓𝑓 ′ (𝑥𝑥) = 3𝑥𝑥 2 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 + 𝑥𝑥 3 cos 𝑥𝑥 + 7
𝑥𝑥0 =6

𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 1

𝑓𝑓 𝑥𝑥0 (6)3 𝑠𝑠𝑠𝑠𝑠𝑠 (6) + 7(6)


𝑥𝑥1 = 𝑥𝑥0 − ′ → 𝑥𝑥1 = 6 −
𝑓𝑓 𝑥𝑥0 3(6)2 𝑠𝑠𝑠𝑠𝑠𝑠 (6) + (6)3 cos(6) + 7

𝑥𝑥1 = 6.099630
S𝑒𝑒𝑒𝑒 𝑛𝑛 = 2

𝑓𝑓 𝑥𝑥1 (6.099630 )𝟑𝟑 𝒔𝒔𝒔𝒔𝒔𝒔 (6.099630 ) + 𝟕𝟕(6.099630 )


𝑥𝑥2 = 𝑥𝑥1 − ′ → 𝒙𝒙𝟐𝟐 = 6.099630 −
𝑓𝑓 𝑥𝑥1 𝟑𝟑(6.099630 )𝟐𝟐 𝒔𝒔𝒔𝒔𝒔𝒔 (6.099630 ) + (6.099630 )𝟑𝟑 𝐜𝐜𝐜𝐜𝐜𝐜(6.099630 ) + 𝟕𝟕

𝑥𝑥2 = 6.093552

𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 3

𝑓𝑓 𝑥𝑥2 (6.093552)3 𝑠𝑠𝑠𝑠𝑠𝑠 (6.093552) + 7(6.093552)


𝑥𝑥3 = 𝑥𝑥2 − ′ → 𝑥𝑥3 = 6.093552 −
𝑓𝑓 𝑥𝑥2 3(6.093552)2 𝑠𝑠𝑠𝑠𝑠𝑠 (6.093552) + (6.093552)3 cos(6.093552) + 7

𝑥𝑥3 = 6.093529
n 𝒙𝒙𝒏𝒏 𝝃𝝃𝒔𝒔 = ⃒𝒙𝒙𝒏𝒏+𝟏𝟏 − 𝒙𝒙𝒏𝒏 |

0 𝑥𝑥0 = 6
1 6.099630 0.099630
2 6.093552 0.006078
3 6.093529 0.000023
EJEMPLO

En la refinación electrolítica de cobre el barro se forma en los ánodos, este puede


perjudicar la pureza del producto. En cierta concentradora la cantidad de barro en
el ánodo (en Kg) sigue la función polinomial: 𝒇𝒇 𝒙𝒙 = 𝒙𝒙𝟐𝟐 𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄 − 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔 − 𝟐𝟐𝟐𝟐
donde x es el tiempo transcurrido en horas. Determine el cruce por cero de la
función, mediante el método de Newton Raphson.
Use 𝒙𝒙𝟎𝟎 = 𝟐𝟐𝟐𝟐.

(4 dígitos)
Sea

𝒇𝒇 𝒙𝒙 = 𝒙𝒙𝟐𝟐 𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄 − 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔 − 𝟐𝟐𝟐𝟐

𝒇𝒇′ 𝒙𝒙 = 2𝒙𝒙 𝒄𝒄𝒄𝒄𝒄𝒄𝒄𝒄 − 𝒙𝒙2 𝒔𝒔𝒔𝒔𝒔𝒔𝒔𝒔 − 6𝑳𝑳𝑳𝑳𝑳𝑳 − 6

𝑥𝑥0 = 23

𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 1

𝑓𝑓 𝑥𝑥0 (23)3 𝑠𝑠𝑠𝑠𝑠𝑠 (23) + 7(23)


𝑥𝑥1 = 𝑥𝑥0 − ′ → 𝑥𝑥1 = 23 −
𝑓𝑓 𝑥𝑥0 3(23)2 𝑠𝑠𝑠𝑠𝑠𝑠 (23) + (23)3 cos(23) + 7

𝑥𝑥1 = 24.8567
INTRODUCCIÓN

S𝑒𝑒𝑒𝑒 𝑛𝑛 = 2

𝑓𝑓 𝑥𝑥1 (24.8567)𝟑𝟑 𝒔𝒔𝒔𝒔𝒔𝒔 (24.8567) + 𝟕𝟕(24.8567)


𝑥𝑥2 = 𝑥𝑥1 − ′ → 𝒙𝒙𝟐𝟐 = 24.8567 −
𝑓𝑓 𝑥𝑥1 𝟑𝟑(24.8567)𝟐𝟐 𝒔𝒔𝒔𝒔𝒔𝒔 (24.8567) + (24.8567)𝟑𝟑 𝐜𝐜𝐜𝐜𝐜𝐜(24.8567) + 𝟕𝟕

𝑥𝑥2 = 24.3840

𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 3

𝑓𝑓 𝑥𝑥2 (24.3840)3 𝑠𝑠𝑠𝑠𝑠𝑠 (24.3840) + 7(24.3840)


𝑥𝑥3 = 𝑥𝑥2 − ′ → 𝑥𝑥3 = 24.3840 −
𝑓𝑓 𝑥𝑥2 3(24.3840)2 𝑠𝑠𝑠𝑠𝑠𝑠 (24.3840) + (24.3840)3 cos(24.3840) + 7

𝑥𝑥3 = 24.5206
n 𝒙𝒙𝒏𝒏 𝝃𝝃𝒔𝒔 = ⃒𝒙𝒙𝒏𝒏+𝟏𝟏 − 𝒙𝒙𝒏𝒏 |

0 𝑥𝑥0 = 23
1 24.8567 1.8567
2 24.3840 0.4726
3 24.5206 0.1366
4 24.5307 0.0101
5 24.5308 0.0001
6 24.5308 0.0000
EJEMPLO:

El nivel de 𝐶𝐶𝐶𝐶2 inicial en partes por millón (ppm) en el Observatorio de


Manua Loa tuvo un valor inicial de 8 , asumiendo que la cantidad de 𝐶𝐶𝐶𝐶2
𝑥𝑥
decrece mediante la función 𝑓𝑓 𝑥𝑥 = 3𝑥𝑥−5 − 𝑐𝑐𝑐𝑐𝑐𝑐 − 6 , determine la
2

concentración cuando el error aproximado sea menor e igual a 0.001

https://picazo.eltiempo.es/2015-el-ano-en-el-que-superamos-las-400-ppm-de-c02/
𝑥𝑥
𝑓𝑓 𝑥𝑥 = 3𝑥𝑥−5 − 𝑐𝑐𝑐𝑐𝑐𝑐 −6
2

′ 𝑥𝑥−5
1 𝑥𝑥
𝑓𝑓 𝑥𝑥 = 3 𝑙𝑙𝑙𝑙 3 + 𝑠𝑠𝑠𝑠𝑠𝑠
2 2

𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 1

𝑓𝑓 𝑥𝑥0 38−5 − cos 8/2 − 6


𝑥𝑥1 = 𝑥𝑥0 − ′ → 𝑥𝑥1 = 8 −
𝑓𝑓 𝑥𝑥0 1
38−5 𝑙𝑙𝑙𝑙(3) + sen 8/2
2

𝑥𝑥1 = 7.260567
S𝑒𝑒𝑒𝑒 𝑛𝑛 = 2

𝑓𝑓 𝑥𝑥1 37.260567−5 + 𝑐𝑐𝑐𝑐𝑐𝑐 7.260567/2 − 6


𝑥𝑥2 = 𝑥𝑥1 − ′ → 𝑥𝑥2 = 7.260567 −
𝑓𝑓 𝑥𝑥1 1
37.260567−5 𝑙𝑙𝑙𝑙(3) + 𝑠𝑠𝑠𝑠𝑠𝑠 7.260567/2
2

𝑥𝑥2 = 6.729556

𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 3

𝑓𝑓 𝑥𝑥2 36.729556−5 − 𝑐𝑐𝑐𝑐𝑐𝑐 6.729556/2 − 6


𝑥𝑥3 = 𝑥𝑥2 − ′ → 𝑥𝑥3 = 6.729556 −
𝑓𝑓 𝑥𝑥2 1
36.729556−5 𝑙𝑙𝑙𝑙(3) + 𝑠𝑠𝑠𝑠𝑠𝑠 6.729556/2
2

𝑥𝑥3 = 6.499873
n 𝒙𝒙𝒏𝒏 𝐸𝐸𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎
0 𝑥𝑥0 = 8
1 7.260567 0.101842
2 6.729556 0.078907
3 6.499873 0.035337
4 6.466344 0.005185
5 6.465731 0.000095 ≤ 𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎
𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨 𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨 − 𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨 𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨
𝑬𝑬𝒂𝒂𝒂𝒂𝒂𝒂𝒂𝒂𝒂𝒂 = 𝑬𝑬𝑬𝑬𝑬𝑬𝑬𝑬𝑬𝑬 𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨 = ⃒ ⃒
𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨 𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨
EJEMPLO:

Aplique el método de Newton Raphson para determinar el cruce por

cero de la función 𝑓𝑓 𝑥𝑥 = 𝑐𝑐𝑐𝑐𝑐𝑐 𝑥𝑥 − 3𝑥𝑥 , con 𝑥𝑥0 = 0.5 ,


convergencia con 4 cifras decimales
1) 𝑓𝑓 𝑥𝑥 = 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 − 3𝑥𝑥
2) 𝑓𝑓 ′ 𝑥𝑥 = −𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 − 3

𝑓𝑓 𝑥𝑥𝑛𝑛 𝑐𝑐𝑐𝑐𝑐𝑐 𝑥𝑥𝑛𝑛 −3(𝑥𝑥𝑛𝑛 )


3) 𝑥𝑥𝑛𝑛+1 = 𝑥𝑥𝑛𝑛 − ⟹ 𝑥𝑥𝑛𝑛+1 = 𝑥𝑥𝑛𝑛 −
𝑓𝑓′ 𝑥𝑥𝑛𝑛 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥𝑛𝑛 −3

Tenemos: x0 = 0.5

𝑐𝑐𝑐𝑐𝑐𝑐 𝑥𝑥0 −3 𝑥𝑥0


𝑛𝑛 = 0 ⇒ 𝑥𝑥1 = 𝑥𝑥0 − ⇒ 𝑥𝑥1 = 0.5
𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥0 −3

𝑐𝑐𝑐𝑐𝑐𝑐 0.5 − 3 0.5


𝑛𝑛 = 1 ⇒ 𝑥𝑥1 = 0.5 − ⇒ 𝑥𝑥1 = 0.3211
𝑠𝑠𝑠𝑠𝑠𝑠 0.5 − 3
n 𝒙𝒙𝒏𝒏

0 0.5000
1 0.3211
2 0.3168
3 0.3168
4 0.3168
5 0.3168
6 0.3168
n Xn
f(x)= 0 0.5000 -0.0109
f'(X)= -0.0607

1 0.3211 -0.0003
-0.0579

2 0.3168 0.0000
-0.0578

3 0.3168 0.0000
-0.0578

4 0.3168 0.0000
-0.0578

5 0.3168 0.0000
-0.0578

6 0.3168 0.0000
-0.0578
EJEMPLO:

2
Resolver 𝑓𝑓 𝑥𝑥 = 𝑥𝑥 − 2 − 𝑙𝑙𝑙𝑙 𝑥𝑥 utilizando el método de

Newton Raphson con un error estimado igual a cero, cuando el

valor inicial es 𝑥𝑥0 = 1.


𝟐𝟐
1) 𝒇𝒇 𝒙𝒙 = 𝒙𝒙 − 𝟐𝟐 − 𝒍𝒍𝒍𝒍 𝒙𝒙

′ 𝟏𝟏 𝒇𝒇 𝒙𝒙𝒏𝒏 (𝒙𝒙𝒏𝒏 −𝟐𝟐)𝟐𝟐 −𝒍𝒍𝒍𝒍𝒙𝒙𝒏𝒏


2) 𝒇𝒇 𝒙𝒙 = 𝟐𝟐 𝒙𝒙 − 𝟐𝟐 − ⟹ 𝒙𝒙𝒏𝒏+𝟏𝟏 = 𝒙𝒙𝒏𝒏 −
𝒙𝒙 𝒇𝒇′ 𝒙𝒙𝒏𝒏
⟹ 𝒙𝒙𝒏𝒏+𝟏𝟏 = 𝒙𝒙𝒏𝒏 − 𝟏𝟏
𝟐𝟐 𝒙𝒙𝒏𝒏 −𝟐𝟐 −
𝒙𝒙𝒏𝒏

Tenemos:

(𝒙𝒙𝟎𝟎 − 𝟐𝟐)𝟐𝟐 − 𝒍𝒍𝒍𝒍𝒙𝒙𝟎𝟎 (𝒙𝒙𝟎𝟎 − 𝟐𝟐)𝟐𝟐 − 𝒍𝒍𝒍𝒍𝒙𝒙𝟎𝟎


𝒏𝒏 = 𝟎𝟎 ⇒ 𝒙𝒙𝟏𝟏 = 𝒙𝒙𝟎𝟎 − ⇒ 𝒙𝒙𝟏𝟏 = 𝟏𝟏 −
𝟏𝟏 𝟏𝟏
𝟐𝟐 𝒙𝒙𝟎𝟎 − 𝟐𝟐 − 𝟐𝟐 𝒙𝒙𝟎𝟎 − 𝟐𝟐 −
𝒙𝒙𝟎𝟎 𝒙𝒙𝟎𝟎

⇒ 𝒙𝒙𝟏𝟏 = 1.33333333

|𝒙𝒙𝟏𝟏 − 𝒙𝒙𝟎𝟎 | = 0.33333333

(𝒙𝒙𝟏𝟏 −𝟐𝟐)𝟐𝟐 −𝒍𝒍𝒍𝒍𝒙𝒙𝟏𝟏


𝒏𝒏 = 𝟏𝟏 ⇒ 𝒙𝒙𝟐𝟐 = 𝒙𝒙𝟏𝟏 − 𝟏𝟏 ⇒ 𝒙𝒙𝟐𝟐 = 1.40857927
𝟐𝟐 𝒙𝒙𝟏𝟏 −𝟐𝟐 −
𝒙𝒙𝟏𝟏
i Xi 𝛏𝛏𝐬𝐬 =│Xi+1 - Xi │

0 1.00000000

1 1.33333333 0.33333333

2 1.40857927 0.07524594

3 1.41238156 0.00380229

4 1.41239117 0.00000961

5 1.41239117 0.00000000
n Xi
f (x) = 0 1.00000000 1.00000000
f '(x)= -3.00000000

1 1.33333333 0.15676237
-2.08333333

2 1.40857927 0.00719689
-1.89277665

3 1.41238156 0.00001809
-1.88326084

4 1.41239117 0.00000000
-1.88323681

5 1.41239117 -6.17219078
-1.88323681
EJEMPLO:

Sea 𝑓𝑓 𝑥𝑥 = −0.9𝑥𝑥 2 + 1.7𝑥𝑥 + 2.5 𝑐𝑐𝑐𝑐𝑐𝑐 𝑥𝑥0 = 5. Calcule la raíz


de 𝑓𝑓 𝑥𝑥 ,aplicando el método de Newton Raphson, con
𝜉𝜉𝑠𝑠 < 0.01%
(5 dígitos )
1) 𝒇𝒇 𝒙𝒙 = −𝟎𝟎. 𝟗𝟗𝒙𝒙𝟐𝟐 + 𝟏𝟏. 𝟕𝟕𝟕𝟕 + 𝟐𝟐. 𝟓𝟓
2) 𝒇𝒇′ 𝒙𝒙 = − 𝟏𝟏. 𝟖𝟖𝟖𝟖 + 𝟏𝟏. 𝟕𝟕
𝒇𝒇 𝒙𝒙𝒏𝒏 −𝟎𝟎.𝟗𝟗𝒙𝒙𝒏𝒏 𝟐𝟐 +𝟏𝟏.𝟕𝟕𝒙𝒙𝒏𝒏 +𝟐𝟐.𝟓𝟓
3) 𝒙𝒙𝒏𝒏+𝟏𝟏 = 𝒙𝒙𝒏𝒏 − 𝒇𝒇′ 𝒙𝒙𝒏𝒏
⟹ 𝒙𝒙𝒏𝒏+𝟏𝟏 = 𝒙𝒙𝒏𝒏 −
− 𝟏𝟏.𝟖𝟖𝒙𝒙𝒏𝒏 +𝟏𝟏𝟏𝟏

Tenemos:

−0.9𝑥𝑥0 2 +1.7𝑥𝑥0 +2.5 −0.9𝑥𝑥0 2 +1.7𝑥𝑥0 +2.5


𝒏𝒏 = 𝟎𝟎 ⇒ 𝒙𝒙𝟏𝟏 = 𝑥𝑥0 − ⇒ 𝑥𝑥1 = 5 −
− 1.8𝑥𝑥0 +17 − 1.8𝑥𝑥0 +17
⇒ 𝑥𝑥1 = 3.424657534.

|𝒙𝒙𝟏𝟏 − 𝒙𝒙𝟎𝟎 | = 𝟏𝟏. 𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓𝟓

−0.9𝑥𝑥1 2 +1.7𝑥𝑥1 +2.5


𝑛𝑛 = 1 ⇒ 𝑥𝑥2 = 𝑥𝑥1 − ⇒ 𝑥𝑥2 = 2.924356997
− 1.8𝑥𝑥1 +17
n 𝒙𝒙𝒏𝒏 𝝃𝝃𝒔𝒔 = ⃒𝒙𝒙𝒏𝒏+𝟏𝟏 − 𝒙𝒙𝒏𝒏 |

0 𝑥𝑥0 = 5
1 3.424657534 1.575342466
2 2.924356997 0.500300537
3 2.861146976 0.063210021
4 2.860104689 0.001042287
5 2.860104406 0.000000283 < 𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
EJEMPLO:

Sea 𝑓𝑓 𝑥𝑥 = 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 − 𝑥𝑥 + 1 , 𝑥𝑥 𝜖𝜖 1, 3 . Calcule la raíz de


𝑓𝑓 𝑥𝑥 , aplicando el método de Newton Raphson.
Sea 𝑓𝑓 𝑥𝑥 = 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 − 𝑥𝑥 + 1
𝑓𝑓 ′ 𝑥𝑥 = 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 − 1
𝑓𝑓 ′′ 𝑥𝑥 = −𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠
Escogiendo el intervalo:
𝑓𝑓 1 = 0.841470985 𝑓𝑓 3 = −1.858879992

𝑓𝑓 1 𝑓𝑓 3 < 0 ⇒ ∃ 𝑥𝑥 ∗ ∈ 1, 3

El mejor extremo es decir, el valor inicial

𝑓𝑓 1 = 0.841470985 𝑓𝑓 3 = −1.858879992

𝑓𝑓 ′′ 1 = −0.841470985 𝑓𝑓 ′′ 3 = −0.141120008

Tenemos que: 𝑥𝑥0 = 3


Entonces:
𝑓𝑓 𝑥𝑥𝑛𝑛 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥𝑛𝑛 − 𝑥𝑥𝑛𝑛 + 1
𝑥𝑥𝑛𝑛+1 = 𝑥𝑥𝑛𝑛 − ⟹ 𝑥𝑥𝑛𝑛+1 = 𝑥𝑥𝑛𝑛 −
𝑓𝑓 ′ 𝑥𝑥𝑛𝑛 cos 𝑥𝑥𝑛𝑛 − 1

Tenemos:

𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥0 − 𝑥𝑥0 + 1


𝑥𝑥1 = 𝑥𝑥0 −
cos 𝑥𝑥0 − 1

𝑠𝑠𝑠𝑠𝑠𝑠 3−3+1
⇒ 𝑥𝑥1 = 3 −
cos 3−1

⇒ 𝑥𝑥1 = 2.065885929
n Xn

0
3.000000000
1
2.065885929
2
1.939820941
3
1.934572691
4
1.934563211
5
1.934563211
OBSERVACION

1) Utilizando el método de Newton Raphson se obtiene resultados en pocas iteraciones.

2) El error 𝑅𝑅𝑛𝑛+1 es proporcional al cuadrado del error anterior.


Supongamos que el error en una iteración es 10−𝑛𝑛 el error siguiente es proporcional al
cuadrado del error anterior es decir 10−2𝑛𝑛 , el que sigue es 10−4𝑛𝑛 ……; esto nos
permite afirmar que cada iteración duplica aproximadamente el número de dígitos
correctos.

3) El método de Newton Raphson algunas veces no converge sino que oscila , esto ocurre si
no hay raíz real , si la raíz es un punto de inflexión ó está el valor inicial muy alejado de la
raíz buscada ó se entrapa.

4) Necesidad de conocer el valor de la derivada de 𝑓𝑓 𝑥𝑥 en cada aproximación lo cual es


una dificultad del método.
MÉTODO
DE
NEWTON RAPHSON
MODIFICADO
Sean los ejemplos:
Ejemplo 1.-

𝑓𝑓 𝑥𝑥 = 𝑥𝑥 3 − 5𝑥𝑥 2 + 7𝑥𝑥 − 3
𝑓𝑓 𝑥𝑥 = 𝑥𝑥 − 3 𝑥𝑥 − 1 𝑥𝑥 − 1

Ejemplo 2.-

𝑓𝑓 𝑥𝑥 = 𝑥𝑥 4 − 6𝑥𝑥 3 + 12𝑥𝑥 2 − 10𝑥𝑥 + 3


𝑓𝑓 𝑥𝑥 = 𝑥𝑥 − 3 𝑥𝑥 − 1 𝑥𝑥 − 1 𝑥𝑥 − 1
Según Ralston y Rabinowitz (1978), propusieron una alternativa para probar que
el método de Newton Raphson converge en forma cuadrática en vez de lineal,
cuando tiene raíces múltiples para esto se ha propuesto las siguientes
modificaciones:
𝒇𝒇 𝒙𝒙
𝒖𝒖 𝒙𝒙 = ′ (1)
𝒇𝒇 𝒙𝒙
La cual tiene raíces.
Si la ecuación de Newton Raphson es

𝑢𝑢 𝑥𝑥𝑖𝑖
𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 − 2
𝑢𝑢′ 𝑥𝑥𝑖𝑖

De (1) derivando tenemos:


𝑓𝑓′ 𝑥𝑥𝑖𝑖 𝑓𝑓′ 𝑥𝑥𝑖𝑖 − 𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓′′ 𝑥𝑥𝑖𝑖
𝑢𝑢′ 𝑥𝑥𝑖𝑖 = 3
𝑓𝑓′ 𝑥𝑥𝑖𝑖 2
Reemplazando (3) en (2), tenemos:

𝑓𝑓 𝑥𝑥𝑖𝑖
𝑓𝑓 ′ 𝑥𝑥𝑖𝑖
𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 − ′
𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓 ′ 𝑥𝑥𝑖𝑖 − 𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓 ′′ 𝑥𝑥𝑖𝑖
𝑓𝑓 ′ 𝑥𝑥𝑖𝑖 2

2
𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓′ 𝑥𝑥𝑖𝑖
𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 −
𝑓𝑓′ 𝑥𝑥𝑖𝑖 𝑓𝑓′ 𝑥𝑥𝑖𝑖 2 − 𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓′′ 𝑥𝑥𝑖𝑖 ]
Por lo tanto
𝒇𝒇 𝒙𝒙𝒊𝒊 𝒇𝒇′ 𝒙𝒙𝒊𝒊
𝒙𝒙 𝒊𝒊+𝟏𝟏 = 𝒙𝒙 𝒊𝒊 −
𝒇𝒇′ 𝒙𝒙𝒊𝒊 𝟐𝟐 − 𝒇𝒇 𝒙𝒙𝒊𝒊 𝒇𝒇′′ 𝒙𝒙𝒊𝒊
EJEMPLO:

El coeficiente de la fricción para el flujo turbulento en un tubo está


dado por la función 𝑓𝑓 𝑥𝑥 = 𝑥𝑥 2 − 2𝑥𝑥𝑒𝑒 −𝑥𝑥 + 𝑒𝑒 −2𝑥𝑥 , con 𝑥𝑥0 = 0.5 .
Emplee el método de Newton Raphson Modificado, e iterar hasta
que 𝜉𝜉𝑠𝑠 < 10−5 . (6 dígitos)
1) 𝒇𝒇 𝒙𝒙 = 𝒙𝒙𝟐𝟐 − 𝟐𝟐𝟐𝟐𝒆𝒆−𝒙𝒙 + 𝒆𝒆−𝟐𝟐𝟐𝟐
2) 𝒇𝒇′ 𝒙𝒙 = 𝟐𝟐𝟐𝟐 − 𝟐𝟐𝒆𝒆−𝒙𝒙 + 𝟐𝟐𝒙𝒙𝒙𝒙−𝒙𝒙 − 𝟐𝟐𝒆𝒆−𝟐𝟐𝟐𝟐
3) 𝒇𝒇′′ 𝒙𝒙 = 𝟐𝟐 + 𝟒𝟒𝒆𝒆−𝒙𝒙 − 𝟐𝟐𝒙𝒙𝒙𝒙−𝒙𝒙 + 𝟒𝟒𝒆𝒆−𝟐𝟐𝟐𝟐

𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓 ′ 𝑥𝑥𝑖𝑖
𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 − ′
𝑓𝑓 𝑥𝑥𝑖𝑖 2 − 𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓 ′′ 𝑥𝑥𝑖𝑖

𝟏𝟏𝟏𝟏−𝟓𝟓 = 0.0001
n 𝒙𝒙𝒏𝒏
⃒𝒙𝒙𝒏𝒏+𝟏𝟏 − 𝒙𝒙𝒏𝒏 |
0 0.500000000

1 0.568013734
0.068013734
2 0.567380913
0.000632821
3 0.567208171
0.000172742
4 0.567129245
0.000078926
5 0.567146333
0.000017087
6 0.567142632 0.000003701 < 𝟏𝟏𝟏𝟏−𝟓𝟓 = 𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
EJEMPLO:

Con la finalidad de evitar contaminación gaseosa, se realiza la lixiviación


por agitación de un mineral mediante la siguiente función:
𝑓𝑓 𝑥𝑥 = 3𝑥𝑥 3 + 5𝑥𝑥 2 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 − 8 .
La solución inicial de lixiviación es de 1.3 gpl de H+, halle la solución final
de dicho proceso si el error relativo estimado es cero, empleando el Método
de Newton Raphson Modificado (9 dígitos)
1) 𝑓𝑓 𝑥𝑥 = 3𝑥𝑥 3 + 5𝑥𝑥 2 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 − 8

2)𝑓𝑓𝑓 𝑥𝑥 = 9𝑥𝑥 2 + 10𝑥𝑥 − 𝑐𝑐𝑐𝑐𝑐𝑐 𝑥𝑥


𝟑𝟑) 𝑓𝑓 ′′ (𝑥𝑥) = 18𝑥𝑥 + 10 + 𝑠𝑠𝑠𝑠𝑠𝑠(𝑥𝑥)

𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓 ′ 𝑥𝑥𝑖𝑖
𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 − ′
𝑓𝑓 𝑥𝑥𝑖𝑖 2 − 𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓 ′′ 𝑥𝑥𝑖𝑖
𝑥𝑥0 =1.3

PRIMERA ITERACION

𝑓𝑓 1.3 = 3(1.3)3 + 5 1.3 2


− 𝑠𝑠𝑠𝑠𝑠𝑠 1.3 − 8 = 6.077441815

𝑓𝑓 ′ 1.3 = 9(1.3)2 +10 1.3 − 𝑐𝑐𝑐𝑐𝑐𝑐 1.3 = 27.94250117

𝑓𝑓 ′′ 1.3 = 18 1.3 + 10 + 𝑠𝑠𝑠𝑠𝑠𝑠 1.3 = 34.36355819

𝑓𝑓(𝑥𝑥0 ) 𝑓𝑓𝑓(𝑥𝑥0 )
𝑥𝑥1 = 𝑥𝑥0 − = 1.003083062
𝑓𝑓𝑓(𝑥𝑥0 ) 2 − 𝑓𝑓(𝑥𝑥0 ) 𝑓𝑓𝑓𝑓(𝑥𝑥0 )

𝒙𝒙𝟏𝟏 = 1.003083062

𝜀𝜀𝑠𝑠 = 1.003083062 − 1.300000000 = 29.600433878


SEGUNDA ITERACION:

𝑓𝑓 1.003083062 = 3(1.003083062)3 + 5 1.003083062 2


− 𝑠𝑠𝑠𝑠𝑠𝑠 1.003083062 − 8 = −0.784421436
2
𝑓𝑓 (1.003083062) = 9 1.003083062 + 10 1.003083062 − 𝑐𝑐𝑐𝑐𝑐𝑐 1.003083062 = 18.54870584
𝑓𝑓 ′′ 1.003083062 = 18 1.003083062 + 10 + 𝑠𝑠𝑠𝑠𝑠𝑠 1.003083062 = 28.89862788

𝑓𝑓(𝑥𝑥1 ) 𝑓𝑓𝑓(𝑥𝑥1 )
𝑥𝑥2 = 𝑥𝑥1 − = 1.042758770
𝑓𝑓𝑓(𝑥𝑥1 ) 2 − 𝑓𝑓(𝑥𝑥1 ) 𝑓𝑓𝑓𝑓(𝑥𝑥1 )

𝒙𝒙𝟐𝟐 = 1.042758770

𝜀𝜀𝑠𝑠 = 1.042758770 − 1.003083062 = 3.804878914


TERCERA ITERACION:

𝑓𝑓 1.042758770 = 3(1.003083062)3 + 5 1.003083062 2


− 𝑠𝑠𝑠𝑠𝑠𝑠 1.003083062 − 8 = −0.025549952

2
𝑓𝑓 ′(1.042758770) = 9 1.003083062 + 10 1.003083062 − 𝑐𝑐𝑐𝑐𝑐𝑐 1.003083062 = 19.70986122

𝑓𝑓 ′′ 1.042758770 = 18 1.003083062 + 10 + 𝑠𝑠𝑠𝑠𝑠𝑠 1.003083062 = 29.63345535

𝑓𝑓(𝑥𝑥2 ) 𝑓𝑓′(𝑥𝑥2 )
𝑥𝑥3 = 𝑥𝑥2 − = 1.044052552
𝑓𝑓′(𝑥𝑥2 ) 2 − 𝑓𝑓(𝑥𝑥2 ) 𝑓𝑓′′(𝑥𝑥2 )

𝒙𝒙𝟑𝟑 = 1.044052552

𝜺𝜺𝒔𝒔 = 1.044052552 − 1.042758770 = 0.123919190


n 𝒙𝒙𝒏𝒏
⃒𝒙𝒙𝒏𝒏+𝟏𝟏 − 𝒙𝒙𝒏𝒏 |
0
1.300000000
1
1.003083062 29.600433878
2
1.042758770 3.804878914
3
1.044052552 0.123919190
4
1.044053812 0.000120725
5
1.044053812 0.000000000
EJEMPLO:

Encuentre la raíz de 𝑡𝑡𝑡𝑡𝑡𝑡 𝑥𝑥 − 0.1x en 𝑥𝑥0 = 𝜋𝜋 , mediante


el método de Newton Raphson modificado y con una
tolerancia de 0.0001
1) 𝑓𝑓 𝑥𝑥 = 𝑡𝑡𝑡𝑡𝑡𝑡 𝑥𝑥 − 0.1𝑥𝑥
2) 𝑓𝑓 ′ 𝑥𝑥 = 𝑠𝑠𝑠𝑠𝑠𝑠 2 𝑥𝑥 − 0.1
3) 𝑓𝑓 ′′ 𝑥𝑥 = 2 𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥. 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡 = 2 𝑠𝑠𝑠𝑠𝑠𝑠 2 𝑥𝑥 𝑡𝑡𝑡𝑡𝑡𝑡 𝑥𝑥

𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓′ 𝑥𝑥𝑖𝑖


4) 𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 −
𝑓𝑓′ 𝑥𝑥𝑖𝑖 2 − 𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓′′ 𝑥𝑥𝑖𝑖

(𝑡𝑡𝑡𝑡𝑡𝑡 𝑥𝑥 − 0.1𝑥𝑥) (𝑠𝑠𝑠𝑠𝑠𝑠 2 𝑥𝑥 − 0.1)


𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 −
(𝑠𝑠𝑠𝑠𝑠𝑠 2 𝑥𝑥 − 0.1) 2 − 𝑡𝑡𝑡𝑡𝑡𝑡 𝑥𝑥 − 0.1𝑥𝑥 (2𝑠𝑠𝑠𝑠𝑠𝑠 2 𝑥𝑥𝑥𝑥𝑥𝑥𝑥𝑥𝑥𝑥)
n 𝒙𝒙𝒏𝒏
⃒𝒙𝒙𝒏𝒏+𝟏𝟏 − 𝒙𝒙𝒏𝒏 |

0 𝑥𝑥0 = 𝜋𝜋

1 3.490658504 0349065850

2 3.476054580 0.014603924

3 3.476140306 0.000085726 < 𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎𝟏𝟏


EJEMPLO:

La velocidad de producción de gas cloro en una celda electrolítica viene dada por

la siguiente función :

𝑥𝑥
𝑓𝑓 𝑥𝑥 = 2𝑐𝑐𝑐𝑐𝑐𝑐 𝑥𝑥 − 𝑡𝑡𝑡𝑡𝑡𝑡( ) donde utiliza Cloruro de Sodio (NaCl) acuoso
2

como electrólito y por la cual inicialmente pasa una corriente de 0.008727

Amperios (A). ¿A qué cantidad de corriente se obtendrá un error estimado menor

de 0.0040? . Emplee el método de Newton Raphson Modificado


𝑥𝑥
𝑓𝑓 𝑥𝑥 = 2𝑐𝑐𝑐𝑐𝑐𝑐 𝑥𝑥 − 𝑡𝑡𝑡𝑡𝑡𝑡
2
1 𝑥𝑥 2
𝑓𝑓 ′ 𝑥𝑥 = −2𝑠𝑠𝑠𝑠𝑠𝑠 𝑥𝑥 − sec
2 2
1 1 𝑥𝑥 𝑥𝑥 𝑥𝑥
𝑓𝑓 ′′ 𝑥𝑥 = −2 cos 𝑥𝑥 − ∗ ∗ 2𝑠𝑠𝑠𝑠𝑠𝑠 𝑠𝑠𝑠𝑠𝑠𝑠 𝑡𝑡𝑡𝑡𝑡𝑡
2 2 2 2 2

′′
1 𝑥𝑥 2 𝑥𝑥
𝑓𝑓 𝑥𝑥 = −2 cos 𝑥𝑥 − sec 𝑡𝑡𝑡𝑡𝑡𝑡
2 2 2
𝑥𝑥0 =0.008727
PRIMERA ITERACION

0.008727
𝑓𝑓 0.008727 = 2𝑐𝑐𝑐𝑐𝑐𝑐 0.008727 − 𝑡𝑡𝑡𝑡𝑡𝑡 =1.9956
2
2
1 0.008727
𝑓𝑓 ′ 0.008727 = −2𝑠𝑠𝑠𝑠𝑠𝑠 0.008727 − sec = −0.5175
2 2

2
′′
1 0.008727 0.008727
𝑓𝑓 0.008727 = −2 cos 0.008727 − sec 𝑡𝑡𝑡𝑡𝑡𝑡 = −2.0021
2 2 2

𝑓𝑓(𝑥𝑥0 ) 𝑓𝑓𝑓(𝑥𝑥0 )
𝑥𝑥1 = 𝑥𝑥0 − = 0.2510
𝑓𝑓𝑓(𝑥𝑥0 ) 2 − 𝑓𝑓(𝑥𝑥0 ) 𝑓𝑓𝑓𝑓(𝑥𝑥0 )

𝒙𝒙𝟏𝟏 = 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐

𝜀𝜀𝑠𝑠 = 0.2510 − 0.008727 = 0.2422


SEGUNDA ITERACION:

𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
𝑓𝑓 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = 2𝑐𝑐𝑐𝑐𝑐𝑐 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 − 𝑡𝑡𝑡𝑡𝑡𝑡 = 1.8112
2

2
1 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
𝑓𝑓 ′(𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐) = −2𝑠𝑠𝑠𝑠𝑠𝑠 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 − 𝑠𝑠𝑠𝑠𝑠𝑠 = −1.0046
2 2

1 𝟎𝟎.𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 2 𝟎𝟎.𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
′′
𝑓𝑓 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = −2 𝑐𝑐𝑐𝑐𝑐𝑐 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑡𝑡𝑡𝑡𝑡𝑡 = −2.0014
2 2 2

𝑓𝑓(𝑥𝑥1 ) 𝑓𝑓𝑓(𝑥𝑥1 )
𝑥𝑥2 = 𝑥𝑥1 − = 0.6437
𝑓𝑓𝑓(𝑥𝑥1 ) 2 − 𝑓𝑓(𝑥𝑥1 ) 𝑓𝑓𝑓𝑓(𝑥𝑥1 )

𝒙𝒙𝟐𝟐 = 0.6436

𝜀𝜀𝑠𝑠 = 0.6436 − 0.2510 = 0.3926


TERCERA ITERACION:

0.6437
𝑓𝑓 0.6437 = 2𝑐𝑐𝑐𝑐𝑐𝑐 0.6437 − 𝑡𝑡𝑡𝑡𝑡𝑡 = 1.2665
2

2
1 0.6437
𝑓𝑓 ′(0.6437) = −2𝑠𝑠𝑠𝑠𝑠𝑠 0.6437 − sec = −1.7557
2 2

2
′′
1 0.6437 0.6437
𝑓𝑓 0.6437 = −2 cos 0.6437 − sec 𝑡𝑡𝑡𝑡𝑡𝑡 = −1.7851
2 2 2

𝑓𝑓(𝑥𝑥2 ) 𝑓𝑓′(𝑥𝑥2 )
𝑥𝑥3 = 𝑥𝑥2 − 2
= 0.3461
𝑓𝑓′(𝑥𝑥2 ) − 𝑓𝑓(𝑥𝑥2 ) 𝑓𝑓′′(𝑥𝑥2 )

𝒙𝒙𝟑𝟑 = 1.0597

𝜺𝜺𝒔𝒔 = 𝟏𝟏. 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 − 𝟎𝟎. 𝟔𝟔𝟔𝟔𝟔𝟔𝟔𝟔 = 𝟎𝟎. 𝟒𝟒𝟒𝟒𝟒𝟒𝟒𝟒


n 𝒙𝒙𝒏𝒏
⃒𝒙𝒙𝒏𝒏+𝟏𝟏 − 𝒙𝒙𝒏𝒏 |

0 0.008727
1 0.2510 0.2422
2 0.6436 0.3926
3 1.0597 0.4161
4 1.2085 0.1487
5 1.2155 0.0071
6 1.2156 0 < 𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎
EJEMPLO:
Para un proceso minero metalúrgico se evalúa la etapa de conminución,
empleando el Método de Newton Raphson Modificado en la función

𝒇𝒇 𝒙𝒙 = 𝒙𝒙𝟑𝟑 − 𝟓𝟓𝒙𝒙𝟐𝟐 + 𝟕𝟕𝒙𝒙 − 𝟑𝟑, para ello se desea saber sobre el rendimiento de las
chancadoras primarias en función de su capacidad esperando obtener un
resultado con error relativo porcentual pequeño, con un valor inicial de 𝒙𝒙𝟎𝟎 = 𝟎𝟎 y el
valor verdadero es 1. (6 dígitos)

1) 𝒇𝒇 𝒙𝒙 = 𝒙𝒙𝟑𝟑 − 𝟓𝟓𝒙𝒙𝟐𝟐 + 𝟕𝟕𝒙𝒙 − 𝟑𝟑


2) 𝒇𝒇′ 𝒙𝒙 = 𝟑𝟑𝒙𝒙𝟐𝟐 − 𝟏𝟏𝟏𝟏𝒙𝒙 + 𝟕𝟕
3) 𝒇𝒇′′ 𝒙𝒙 = 𝟔𝟔𝟔𝟔 − 𝟏𝟏𝟏𝟏
𝒇𝒇 𝒙𝒙𝒊𝒊 𝒇𝒇′ 𝒙𝒙𝒊𝒊
𝒙𝒙 𝒊𝒊+𝟏𝟏 = 𝒙𝒙 𝒊𝒊 − ′
𝒇𝒇 𝒙𝒙𝒊𝒊 𝟐𝟐 − 𝒇𝒇 𝒙𝒙𝒊𝒊 𝒇𝒇′′ 𝒙𝒙𝒊𝒊
n 𝒙𝒙𝒏𝒏 𝑬𝑬𝑹𝑹%

0 𝑥𝑥0 = 0 100

1 1.105263158 10.5263158

2 1.003081664 0.3081664

3 1.000002402 0.0002395

Observamos que los valores convergen cuadráticamente


Como ejemplo, aplicar el Método de Newton Raphson al ejemplo anterior

n 𝒙𝒙𝒏𝒏 𝑬𝑬𝑹𝑹%
0 𝒙𝒙𝟎𝟎 = 0
1 0.42857 57.143
2 0.68571 31.429
3 0.83286 16.714
4 0.91332 8.668
5 0.95577 4.423
6 0.97764 2.236
7 0.98875 1.125

También podría gustarte