La figura 3.
4 muestra el código en Matlab para solución de sistemas de ecuaciones por el método de Gauss simple
(eliminación Gaussiana). Al ejecutar este código el programa pide inicialmente que se indique el número de variables en el
sistema, posteriormente se introducen los coeficientes del sistema.
clear all
close all
clc
fprintf('Método de Gauss Simple\n\n')
x=input('Defina el no. de variables del sistema: ');
r=x;
c=r+1;
fprintf('Se tiene la matriz aumentada de %i x %i\n',r,c)
for i=1:r
for j=1:c
fprintf('Valor de a(%i,%i):',i,j)
v=input(' ');
A(i,j)=v;
end
end
fprintf('\nMatriz Aumentada del Sistema:\n')
A
fprintf('Procedimiento:\n')
pause
X=zeros(1,c-1);
for i=1:c-1
A(i,:)=A(i,:)/A(i,i)
for j=i+1:r
A(j,:)=A(j,:)-A(i,:)*A(j,i)
end
end
for i=r:-1:1
s=0;
for j=c-1:-1:1
s=s+A(i,j)*X(1,j);
end
X(1,i)=A(i,c)-s;
end
fprintf('Resultados:\n');
for i=r:-1:1
fprintf('x%i = %f\n',i,X(1,i))
end
Figura 3.4 | Código en Matlab para el método de Gauss simple.
Al implementar dicho código se tiene lo siguiente:
Método de Gauss Simple
Defina el no. de variables del sistema: 3
Se tiene la matriz aumentada de 3 x 4
Valor de a(1,1): 3
Valor de a(1,2): 9
Valor de a(1,3): 3
Valor de a(1,4): 3
Valor de a(2,1): 2
Valor de a(2,2): -1
Valor de a(2,3): 4
Valor de a(2,4): 22
Valor de a(3,1): 5
Valor de a(3,2): 3
Valor de a(3,3): -1
Valor de a(3,4): 11
Matriz Aumentada del Sistema:
A=
3 9 3 3
2 -1 4 22
5 3 -1 11
Procedimiento:
A=
1 3 1 1
2 -1 4 22
5 3 -1 11
A=
1 3 1 1
0 -7 2 20
5 3 -1 11
A=
1 3 1 1
0 -7 2 20
0 -12 -6 6
A=
1.0000 3.0000 1.0000 1.0000
0 1.0000 -0.2857 -2.8571
0 -12.0000 -6.0000 6.0000
A=
1.0000 3.0000 1.0000 1.0000
0 1.0000 -0.2857 -2.8571
0 0 -9.4286 -28.2857
A=
1.0000 3.0000 1.0000 1.0000
0 1.0000 -0.2857 -2.8571
0 0 1.0000 3.0000
Resultados:
x3 = 3.000000
x2 = -2.000000
x1 = 4.000000