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