camagorra24
19 de Mayo del 2006

Quisiera pedir el gran favor, si alguien me puede facilitar el codigo fuente de alguna forma de solucion de eliminacion gaussiana en matlab para fines netamente educativos.
Gracias
Atte: Carlos Mario G

Benjamin
19 de Mayo del 2006

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

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 nn')
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 nn')
for i=1:n
fprintf('x(%d) = %d n',i,x(i)); % Presentación de resultados
end