Gauss-jordan Sencillo Matlab 2x5c

%** %** %** %**

Metodo de Gauss Jordan por pasos Ing. Daniel Castro Garrido Materia: Computacion Aplicada IPN

** ** ** **

clear; clc; fprintf('RESOLUCION DE ECUACIONES LINEALES POR EL METODO GAUSS-JORDAN \n\n'); m=input('Cuantas ecuaciones son: '); n=input('Cuantas incognitas tienen las ecuaciones: ');

%** En los siguentes for anidados se da entrada a los %** a los elementos de la matriz 1 a 1

**

for i=1:m for j=1:n fprintf('El elemento A[%d %d]',i,j); A(i,j)= input (' ') end end A pause %En seguida se agregan los resultados de las ecuaciones for j=1:m fprintf('Agregue la columna de resultados R[%d]: ',j) R(j)= input(' ') end B=[A R'] %********************************************************* %** En seguida se normalizan los pivotes y se hacen cero** %** todos los numeros por debajo de ellos ** %*********************************************************

for i=1:m B(i,:)=B(i,:)/B(i,i); for j=i+1:m B(j,:)=B(j,:)-B(i,:)*B(j,i); j=j+1; B pause end i=i+1; B pause end fprintf('HOLA HAREMOS CERO LOS ELEMENTOS INFERIORES') %******************************************************

%** En la siguiente seccion se hacen cero los numeros** %** que estan por encima de la diagonal principal ** %******************************************************

for i=m:-1:2 for j=i-1:-1:1 B(j,:)=B(j,:)-B(i,:)*B(j,i); j=j-1; B pause end i=i-1; B pause end %Presentando las variables for i=1:m for j=n+1

fprintf('\nLa variable X[%d] es:\n',i) B(i,j) end end