FACULTAD DE CIENCIAS

ESTADISTICA E INFORMATICA

DOCENTE:
JAVIER ROSAS ORLANDO JAVIER
CURSO:
METODOS NUMERICOS
TEMA:
METODO DE NEWTON PARA SISTEMAS DE ECUACIONES NO LINEALES

QUIÑONES CORTABRAZO JERICO JEAN PAUL

UNASAM
2020-I
 METODO DE NEWTON PARA SISTEMAS DE ECUACIONES NO LINEALES

Definición:

Consiste en elegir las coordenadas de un punto ( x 1, y 1) como aproximación del

punto de intersección de las funciones u(x,y) y v(x,y) que hacen que estas se anulen.

Obtener los valores de las funciones u(x ,y), v(x,y) evaluadas con las coordenadas (
x 1, y 1) y localizar los cuatro puntos u( x 1 ,y), v( x 1,y ),u(x, y 1) y v(x, y 1).
Trazar una recta tangente paralela a la secante que une los puntos u( x 1,y) y u(x, y 1) y
otra recta tangente paralela a la secante que une los puntos v( x 1,y) y v(x, y 1).

El punto de intersección de estas dos tangentes constituye una segunda

aproximación ( x 2 , y 2 ¿ del punto de intersección de las dos funciones.
 El proceso se repite n veces hasta que las coordenadas de punto de intersección (
x n , y n ¿ coincida prácticamente con el valor exacto de la intersección entre las dos
curvas.

Este procedimiento corresponde, analíticamente, a extender el uso de la derivada, ahora

para calcular la intersección entre dos funciones no lineales.

Al igual que para una sola ecuación, el cálculo de se basa en la expansión de la serie de
Taylor de primer orden, ahora de múltiples variables, para considerar la contribución de
más de una variable independiente en la determinación de la raíz.

Para dos variables, la serie de Taylor de primer orden se escribe, para cada ecuación no
lineal

dU i dU
U i +1=U i+( X i +1−X i ) +( y i+1 − y i) i
dx dy

dv i dv
vi +1=vi +( X i +1−X i ) +( y i+1 − y i) i
dx dy
.Pero U i +1=V i+1=0:
dU i dU i dU i dU i
U i+ x i+1−¿ x i+ y i+1 − y i=0
dx dx dy dy
dv i dv i dv i dv i
V i+ x i+1−¿ x i+ y i+1 − y i=0
dx dx dy dy

Que escribiendo en el orden conveniente:

dU i dU i dU dU i
x i+1 + y i+1 =−U i + i x i+ yi
dx dy dx dy
dv i dv i dv dv
x i+1 + y i+1 =−V i + i x i+ i y i
dx dy dx dy

Solución del sistema por determinantes:

dU i dU i dU dU i
x i+1 + y i+1 =[−U i + i x i+ y i] c1
dx dy dx dy
dv i dv i dv dv
x i+1 + y i+1 =[−V i + i x i+ i y i] c2
dx dy dx dy

dU i dU i
dx dy

Jacobiano
dv i dv i dv dU i dU i dv i
D= dx dy
= dyi dx
- dy dx

dU i
C1
dy
dv dv dv i dU i dU i dv i du dv dU i dU i dv i
Dx= C 2 dyi =−ui dyi +¿ dy dx
x i +¿
dy dy
y i+ vi i −¿ i
dy dx dy
x i−¿
dy dy
yi
 dv i dU i dU i dv i dv i dui
xi [ dy dx
−
dy dx ]
− ui
dy
+ vi
dy
X i +1=¿
dv i dU i dU i dv i
−
dy dx dy dx

dvi dui
ui −v i
dy dy
X i +1=¿ X i -
J

La solución del sistema es:

dvi dui
ui −v i
dy dy
X i +1=¿ X i -
J
Donde “J” es el
determinante Jacobiano
del sistema
dui dv i
vi −ui
dx dx
Y i+1 =¿ Y i-
J
Ejemplo:
Hallar las raices de :
x 2+xy-10=0
y+3x y 2-57=0
con X i =1.5 ,Y i =3.5
Solución
ITERACION 1:

U= x +xy-10=0 v= y+3x y 2-57=0

2

dui dui
=2 x+ y =2 ( 1.5 )+3.5=6.5 =x=1.5
dx dy
dv i dv i
=3 y =3(3.5)2=36.75
2
=1+ 6 xy=1+6 ( 1.5 ) ( 3.5 )=32.5
dx dy
dv i dui dui dv i
J= – =( 32.5 ) ( 6.5 )− (1.5 )( 36.75 )=156.125
dy dx dy dx
Entonces, reemplazamos x y y en u y v:
ui=(1.5)2+ (1.5 )( 3.5 ) −10=−2.5

vi =( 3.5 ) +3 ( 1.5 ) (3.5)2±57=1.625

dvi dui
ui −v i
dy dy
−2.5 ( 32.5 )−1.625(1.5)
X i +1=¿ X i - =1.5- 156.125 =2.0360
J

dui dv i
vi −ui
dx dx
 1.625(6.5)−(−2.5) ( 36.75 )
Y i+1 =¿ Y i- =3.5 - =2.8438
156.125
J
Tabla de iteración 1:
Iteración Xi Yi Ui Vi du/dx du/dy dv/dx dv/dy Jacobiano
1.62
1 1.5 3.5 -2.5 5 6.5 1.5 36.75 32.5 156.125
2 2.036 2.8438              

Y i=2.8438 X i=2.036

ITERACION 2
Y i=2.8438 X i=2.036
dui dui
=2 x+ y =2 ( 2.036 )+ 2.8438=6. 9158 =x=2.036
dx dy
dv i
=3 y =3(2.8438)2=24.2616
2
dx
dv i
=1+ 6 xy=1+6 ( 2.036 ) ( 2.8438 )=35.7399
dy
dv i dui dui dv i
J= – =( 35.7399 ) ( 6.9158 )−( 2.0360 )( 24.2616 ) =197.7734
dy dx dy dx
Entonces, reemplazamos x y y en u y v:
ui=(2.0360)2+ ( 2.036 ) (2.8438 )−10=−0.0647

vi =( 2.8438 ) +3 ( 2.036 ) (2.8438)2±57=−4.7596

dvi dui
ui −v i
dy dy
−0.0647 ( 35.7399 )−(−4.7596)(2.036)
X i +1=¿ X i - =2.036- 197.7734 =1.9987
J

dui dv i
vi −ui
dx dx
−4.7596(6.9158)−(−0.0647) ( 24.2616 )
Y i+1 =¿ Y i- =2.8438 - =3.0023
197.7734
J
Tabla de iteración 2:
Jacobian
Iteración Xi Yi Ui Vi du/dx du/dy dv/dx dv/dy o

1 1.5 3.5
-2.5 1.625 6.5 1.5 36.75 32.5 156.125
- -
2 2.036 2.8438 0.0647 4.7596 6.9158 2.036 24.2616 35.7399 197.77
3 1.9987 3.0023              

Y i=3.0023 X i =1.9987

ITERACION 3
Y i=3.0023 X i =1.9987
dui dui
=2 x+ y =2 ( 1.9987 )+ 3.0023=6.9997 =x=1.9987
dx dy
dv i
=3 y =3(3.0023)2=21.0414
2
dx
dv i
=1+ 6 xy=1+6 ( 1.9987 ) ( 3.0023 )=35. 0042
dy
dv i dui dui dv i
J= – =( 37.0042 ) ( 6.9997 )−( 1.9987 ) ( 27.0414 ) =204.9707
dy dx dy dx
Entonces, reemplazamos x y y en u y v:
ui=(1.9987)2+ ( 1.9987 ) (3.0023 )−10=−0.0045

vi =( 3.0023 ) +3 ( 1.9987 ) (3.0023)2−57=0.0499

dvi dui
ui −v i
dy dy
−0.0045 ( 375.0042) −(0.0499)(1.9987)
X i +1=¿ X i - =1.9987- 204.9707 =2
J

dui dv i
vi −ui
dx dx
 0.0499(6.9997)−(−0 .0045) ( 21.0414 )
Y i+1 =¿ Y i- =3.0023- =3
204.9707
J
ITERACION 4

Y i=3 X i=2

ui=(2)2+ (2 )( 3 ) −10=0

vi =( 3 ) +3 ( 2 ) (3)2−57=0

Jacobian
Iteración Xi Yi Ui Vi du/dx du/dy dv/dx dv/dy o

1 1.5 3.5
-2.5 1.625 6.5 1.5 36.75 32.5 156.125
- -
2 2.036 2.8438 0.0647 4.7596 6.9158 2.036 24.2616 35.7399 197.77
-
3 1.9987 3.0023 0.0045 0.0499 6.9997 1.9987 27.0414 37.0042 204.9707
4 2 3 0 0 0 0 0 0 0