UNIVERSIDAD NACIONAL DE TRUJILLO
FACULTAD DE INGENIERÍA
ESCUELA PROFESIONAL DE INGENIERÍA DE MATERIALES
“EXAMEN PRÁCTICO 02”
INTEGRANTES
DOCENTE
Alberto Geldres
CURSO
Métodos Numéricos
TRUJILLO – PERÚ
2023
MÉTODO MULTIVARIABLE DE NEWTON RAPHSON
DIAGRAMA DE FLUJO
PRUEBA DE ESCRITORIO
Hallar las soluciones de siguientes funciones, considerando valores iniciales
de x = 0.5, y = 0.2 con un error permitido de 0.2
2 2
u ( x , y )= x + y −5
V ( x , y )=xy −2
Derivamos las funciones respecto a x e y
Para Iteración 0: x = 0.5, y = 0.2
2 2
U ( i )=( 0.5) +(0.2) −5=−4.71
V ( i )=0.5 ( 0.2 )−2=−1.9
Reemplazamos los valores iniciales en las derivadas
𝑈´(𝑥) = 2𝑥 = 2(0. 5) = 1
𝑈´(𝑦) = 2𝑦 = 2(0. 2) = 0. 4
𝑉´(𝑦) = 𝑥 =0. 5
𝑉(𝑥) = 𝑦 = y = 0.2
−4.71 ×0.5−(−1.9 ×0.4 )
x i+1=0.5− =4.297619048
1 ×0.5−0.4 × 0.2
−1.9 ×1−(−4.71 ×0.2 )
y i+ 1=0.2− =2. 48095238
1× 0.5−0.4 ×0.2
−|4.297619048−0.5|
E ax= × 100=100 %
4.297619048
|2.48095238−0.2|
Eay= ×100=91. 93857965 %
2.48095238
-Para la iteración 1:
x=4.297619048 , y =2.4809523
2 2
U (i)=(4.297619048) +(2.48095238) −5=19. 6246542
V (i)=4. 297619048(2. 48095238)−2=8. 66218821
Reemplazamos los valores iniciales en las derivadas
U ´ ( x )=2 x=2 ( 4. 297619048 )=8. 595238095
U ´ ( y )=2 y=2 ( 2. 48095238 )=4. 961904762
V ´ ( y )=x=4.297619048
V ( x )= y= y =2. 48095238
19.6246542 ×4.297619048−8.66218821 × 4.961904762
X i +1=4.297619048− =2.618352548
8.595238095 ×4.297619048−4.96190476 × 2.48095238
−1.9 ×1−−4.71 ×0.2
Y i+ 1=2. 48095238− =1.43478986
1× 0.5−0.4 ×0.2
|2.618352548−4.2976 19048|
Eax= × 100=64.13446884 %
2.618352548
|1.43478986−2.48095238|
Eay= × 100=72.91398935 %
1.43478986
- Para la Iteración 2
x=2. 618352548 , y=1. 43478986
U (i)=¿
V (i)=2. 618352548(1. 43478986)−2=1. 75678568
Reemplazamos los valores iniciales en las derivadas
U ´ (x )=2 x=2(2. 618352548)=5. 236705097
U ´ ( y)=2 y=2(1. 43478986)=2. 869579714
V ´ ( y )=x=2. 618352548
V (x )= y=1. 43478986
X i +1=2.075527772
Y i+ 1=1.0612931
|2.075527772−2.618352548|
Eax= ×100=26.15357808 %
2.075527772
|1.0612931−1.43478986|
Eay= ×100=35.19261192 %
1.0612931
x=2. 075527772 , y=1. 0612931
- Para Iteración 3 :
2 2
𝑈(𝑖) = (2. 075527772) + (1. 0612931) − 5 = 0. 434158568
𝑉(𝑖) = 2. 075527772(1. 0612931) − 2 = 0. 2027433
𝑈´(𝑥) = 2𝑥 = 2(2. 075527772) = 4. 151055544
Reemplazamos los valores iniciales en las derivadas
𝑈´(𝑦) = 2𝑦 = 2(1. 0612931) = 2. 122586193
𝑉´(𝑦) = 𝑥 = 2. 075527772
𝑉(𝑥) = 𝑦 = 1. 0612931
𝑋 = 2. 001541898
𝑖+1
●
𝑌 = 1. 00144201
𝑖+1
●
𝐸𝑎𝑥 = × 100 = 3. 696443925%
|2.001541898−2.075527772|
2.0015418
●
98
𝐸𝑎𝑦 = × 100 = 5. 976490812%
|1.00144201−1.0612931|
1.001442
●
01
x=2. 001541898 , y=1. 00144201
- Para Iteración 4 :
𝑈(𝑖) = (2. 001541898) + (1. 00144201) − 5 = 0. 009056062
2 2
𝑉(𝑖) = 2. 001541898(1. 00144201) − 2 = 0. 00442814
Reemplazamos los valores iniciales en las derivadas
𝑈´(𝑥) = 2𝑥 = 2(2. 001541898) = 4. 003083796
𝑈´(𝑦) = 2𝑦 = 2(1. 00144201) = 2. 002884014
𝑉´(𝑦) = 𝑥 =2. 001541898
𝑉(𝑥) = 𝑦 =1. 00144201
● 𝑋 = 2. 000000744
𝑖+1
𝑌 = 1. 00000074
𝑖+1
●
𝐸𝑎𝑥 = × 100 = 0. 077057683%
|2.000000744−2.001541898|
2.0000007
●
44
𝐸𝑎𝑦 = × 100 = 0. 144126708%
|1.00000074−1.00144201|
1.000000
●
74
x=2. 000000744 , y=1. 00000074
- Para Iteración 5 :
2 2
𝑈(𝑖) = (2. 000000744) + (1. 00000074) − 5 = 4. 45241𝐸 − 06
𝑉(𝑖) = 2. 000000744(1. 00000074) − 2 = 2. 2212𝐸 − 06
𝑈´(𝑥) = 2𝑥 = 2(2. 000000744) = 4. 000001487
Reemplazamos los valores iniciales en las derivadas
𝑈´(𝑦) = 2𝑦 = 2(1. 00000074) = 2. 00000074
𝑉´(𝑦) = 𝑥 =2. 000000744
𝑉(𝑥) = 𝑦 =y=1. 00000074
𝑋 = 2
𝑖+1
●
𝑌 = 1
𝑖+1
●
● 𝐸𝑎𝑥 = |2−2.000000744|
2
× 100 = 3. 71865𝐸 − 05
● 𝐸𝑎𝑦 = |1−1.00000074|
1
× 100 = 7. 38742𝐸 − 05
Para la iteración N°5 se consigue encontrar las soluciones de las ecuaciones para
error permitido de 0.2 . Las raíces convergen a x=2 ; y=1
CORRIDA EN MATLAB
⇒ FUNCIÓN RESUELTA
Link del video: https://www.youtube.com/watch?
v=rHzSj7LwcC8