Sistema Dee Cuac I Ones
Sistema Dee Cuac I Ones
Sistema Dee Cuac I Ones
que la próxima diferencia entre valores sucesivos de x y o z sea menor que 0.02.
10𝑥 + 2𝑦 + 6𝑧 = 28
𝑥 + 10𝑦 + 9𝑧 = 7
−2𝑥 + 7𝑦 + 10𝑧 = 17
Termine los cálculos si no se verifica dicho creiterio de detención luego de haber efectuado 5
iteraciones.
Iteracion 1
X=2.8 ea=100
Y=0.42 ea=100
Z=1.9660 ea=100
--------------------------------------------
Iteracion 2
X=1.5364 ea=82.2442%
Y=-1.223 ea=134.3407%
Z=2.8634 ea=31.3406%
--------------------------------------------
Iteracion 3
X=1.3266 ea=15.8181
Y=-2.0097 ea=39.1439
Z=3.3721 ea=15.0858
--------------------------------------------
Iteracion 4
X=1.1787 ea=12.5472%
Y=-2.4528 ea=18.0633%
Z=3.6527 ea=7.6809%
--------------------------------------------
Iteracion 5
X=1.0989 ea=7.2546%
Y=-2.6973 ea=9.0657%
Z=3.8079 ea=4.0764%
--------------------------------------------
3)un ingeniero supervisa la producción de tres tipos de automóviles. Se requieren tres clases
de materiales (metal plástico y caucho) para la producción. La cantidad necesaria para producir
cada modelo de auto es de:
Automóvil Metal plastico Caucho
Modelo1 1500 25 100
Modelo2 1700 33 120
Modelo3 1900 42 160
Si se dispone de un total de 106 toneladas de metal, 2.17 toneladas de plástico y 8.2 toneladas
de caucho ¿Cuántos automóviles se pueden producir por dia de cada modelo?.
106000
B= 2170
8200
10
A-1xB =20
30
--------------------------------------------
Iteracion 3
X=3 ea=0.31580
Y=-2.5 ea=0.0145
Z=7 ea=0.0042
--------------------------------------------
Iteracion 4
X=3 ea=0.0011
Y=-2.5 ea=0.0005
Z=7 ea=0.0000
--------------------------------------------
5)Use el método de Gauss jordan para obtener la solución al siguiente sistema de ecuaciones
lineales.
--> rref(A)
ans =
1. 0. 0. 3.
0. 1. 0. -2.5
0. 0. 1. 7.
Emplee el método de Gauss-Jordán y muestre todas las etapas del cálculo necesario ea, eb ,ec
voltaje en voltios.
13 | -3 |-4
15 | 47 |-12
-3 | -3 | 9
13 | -3 | -4 :120
15 | 47 | -12 :0
-3 | -3 | 9 :15
15 | 47 | -12 :0
13 | -3 | -4:120
-3 | -3 | 9 : 15
15 | 47 | -12 :0
-3 | -3 |9 :15
15 | 47 | -12 :0
0 | 32/5 | 33/5: 15
15 | 47 | -12 :0
0 | -656/15 | 17/5 :120
0 |0 | 291/41 :1335/41
15 | 47 | -12 :0
0 |0 |1 :445/97
15 | 47 | -12 :0
0 | -656/15 | 0 :10127/97
0 |0 | 1 :445/97
15 | 47 | 0 :5340/97
0 | -656/15 | 0 :10127/97
0 |0 | 1 :445/97
15 | 47 | 0 :5340/97
0 | 1 | 0: -3705/1552
0 | 0 | 1 :445/97
15 | 0 | 0 :167.251933
0 | 1 | 0: -3705/1552
0 | 0 | 1 :445/97
1 | 0 | 0 :10.072816
0 | 1 |0: -3705/1552
0 | 0 | 1 :445/97
X=10.072816;
Y=-2.1116505;
Z=4.3203883
function gauss_seidel
A=input('ingrese la matriz de coeficientes:');
B=input('ingrese las constantes:');
es=input('ingrese el %de error:');
x=[0 0 0];
xs=[0 0 0];
i=input('ingrese las iteraciones maximas: ');
for k=1:i
//calculando las nuevas aproximaciones
x(1)=(B(1)-A(1,2)*x(2)-A(1,3)*x(3))/A(1,1);
x(2)=(B(2)-A(2,1)*x(1)-A(2,3)*x(3))/A(2,2);
x(3)=(B(3)-A(3,1)*x(1)-A(3,2)*x(2))/A(3,3);
//calculando los porcentajes de error
error1=abs((xs(1)-x(1))/x(1))*100;
error2=abs((xs(2)-x(2))/x(2))*100;
error3=abs((xs(3)-x(3))/x(3))*100;
//muestra los resultados en la pantalla
printf('Iteracion%2d\n',k);
printf('\t X\t Y\t Z\n');
printf('\t%1.4f\t%1.4f\t%1.4f\n',x(1),x(2),x(3));
printf('\terror1\terror2\terror3\n');
printf('\t%1.4f\t%1.4f\t%1.4f\n',error1,error2,error3);
disp('--------------------------------------------');
//termina si se llega al error
if error1<=es & error2<=es & error3<=es
break
end
xs=x;
end
endfunction