clear all
clc
fprintf('Resolución del sistema Ax = b por \n');
fprintf('Eliminación Gaussiana con sustitución hacia atrás \n\n')
n=input('Ingrese el orden del sistema = ');
fprintf('\n');
fprintf('Matriz coeficientes del sistema \n\n')

for i=1:n
for j=1:n
fprintf('coeficiente A(%d,%d) = ', i,j ) % Definición de la matriz
c(i,j)=input(' '); % coeficientes del sistema
end
end
fprintf('\nVector términos independientes \n\n')
for i=1:n
fprintf('término b(%d) = ', i ) % Definición del vector de
c(i,n+1)=input(' '); % términos independientes
end

for i=1:n-1
for j=i+1:n
L=c(j,i)/c(i,i); % Cálculo de los coeficientes L
c(j,:)=c(j,:)-L*c(i,:); % Eliminación
end
end

x(n)=c(n,n+1)/c(n,n); % Cálculo de la última incógnita

for i=n-1:-1:1
Suma=0;
for j=i+1:n
Suma=Suma+c(i,j)*x(j);
end
x(i)=(c(i,n+1)-Suma)/c(i,i); % Cálculo del resto de las incógnitas
end
fprintf('\nSolución \n\n')
for i=1:n
fprintf('x(%d) = %d \n',i,x(i)); % Presentación de resultados
end

Responder