Mtodos Numricos y Programacin
Gua de Laboratorio Nro: 06
Tema: Sistemas No Lineales
Docente: Braulio Gutirrez Pari
Ejercicio 0.1 En algn lenguaje de su preferencia, implemente el mtodo de Newton para
sistema de ecuaciones no lineales
Una Implementacin Bsica:
function.....= NEWTON_NL(x, prec)
iter =....
while norm(Fn(x))>....
iter = ...............
x =.................
if ...........
error(parece que newton no converge);
end
end
Ejercicio 0.2 Del grfico adjunto
1. Plantear el sistema de ecuaciones no lineales
2. Obtener aproximaciones iniciales del grfico para las raices.
3. Resuelva el sistema de ecuaciones planteado en la parte (1) con el mtodo de Newton,
pruebe con los puntos iniciales
"
12
2
"
y.
106
12
2
Itere hasta que x(k) x(k1) <
Ejercicio 0.3 Se quiere resolver el sistema de ecuaciones no lineales
x + 2y 4 = 0
(x 6)2 y + 2 = 0
1. Hallar el sistema de ecuaciones y el Jacobiano
F (x) =
Fn0 (xk ) = J(x) =
"
"
df
dx1
dg
dx1
x1 + 2x2 4
(x1 6)2 x2 + 2
df
dx2
dg
dx2
"
2. Ejecutar el programa NEWTON_NL, con un punto inicial
"
1
2
con un parmetro de pre-
cisin 0,000001, observamos:
"
4._______________
4._______________
(Verifquelo!)
3. Ahora ejecutando NEWTON_NL, con otro punto inicial
"
9
3
te-nemos:
que converge a otra raz
"
8._______________
6._______________
(Verifquelo!).
Ejercicio 0.4 Haga sus respectivas grficas del ejercicio anterior y comente, porqu el mismo
sistema de ecuaciones no lineales, para puntos distintos converge a dos raices distintas
Ejercicio 0.5 Se quiere resolver el sistema de ecuaciones no lineales
(x 3)2 y + 4 = 0
x + 2y 16 = 0
2
Al ejecutar NEWTON_NL, verifique con los siguientes puntos iniciales
"
1
2
"
6
10
Ejercicio 0.6 Consideremos el siguiente sistema de ecuaciones no li-neales de 3 incgnitas y
3 ecuaciones:
7x1 x2 + 5x2 x23 sen x1 12 = 0
x41 + cos2 x2 + 2x33 8 = 0
(1)
6x1 + 2x2 x3 + 34 = 0
Observe que
7x1 x2 + 5x2 x23 sen x1 12
F (x) = x41 + cos2 x2 + 2x33 8
6x1 + 2x2 x3 + 34
7x2 x23 cos x1
7x1 + 5
2x3 sen x1
J(x) =
4x31
2 cos x2 sen x2
6x23
6
2
1
10
Si usamos como punto inicial x0 = 20 y como parmetro de precisin = 106 , una
50
implementacin computacional bsica nos otorga:
x13
4, 23134959407946
= 1, 56752981158965
5, 47684281234392
(verifquelo!)
Ejercicio 0.7 Se quiere resolver el sistema de ecuaciones no lineales
5x21 + 6x1 x2 + 5x22 4x1 + 4x2 4 = 0
x21 + x22 1 = 0
Al ejecutar NEWTON_NL con punto inicial
"
1
1
y con una precisin de prec=0.000001
Ejercicio 0.8 Se quiere resolver el sistema de ecuaciones no lineales
x3 + y 3 2xy = 0
x2 + y 2 1 = 0
3
Al ejecutar Newton, verifique con los siguientes puntos iniciales
"
1
2
Ejercicio 0.9 Se quiere resolver el sistema de ecuaciones no lineales
2x4 2x2 y + y 2 2y 3 + y 4 = 0
3x2 2xy 6y 2 + 3 = 0
Ejercicio 0.10 Resolver el sistema de ecuaciones no lineales
x3 + x2 y xz + 6 = 0
ex + ey z = 0
y 2 2xz 4 = 0
Use el mtodo de Newton, con el punto inicial 2 y con parmetro de precisin
1
6
= 10
Ejercicio 0.11 Resolver el sistema de ecuaciones no lineales
sen(xy) + z = 1
2 sen(xz) + y = 2
3 sen(yz) + x = 3
0
pruebe con un punto inicial 0 y con parmetro de precisin = 106 .y tambin anote su
0
Jacobiano
J(x) =
iter
x1
x2
1
2
..
.
x3
