titulo



























FUNCIONES
Funcion#1
[ co,X ] = montante6( N,A );
if co
for i=1:N
sprintf('La composicion en el plato %f es de %d \n',i,X(i))
end
elsedisp('no converge')
end

Funcion#2
gauss jordan
function [co,X] = jordan6( N,A );
X=NaN;
for k=1:N
L=k;
P=A(k,k);
while P==0
L=L+1;
if L>N
co=false;
return
else
forj=k:N+1
T=A(k,j);
A(k,j)=A(L,j);
A(L,j)=T;
end
P=A(k,k);
end
end
for j=k:N+1
A(k,j)=A(k,j)/P;
end
for i=1:N
if i~=k
C=A(i,k);
for j=k:N+1A(i,j)=A(i,j)-(A(k,j)*C);
end
end
end
end
for i=1:N
X(i)=A(i,N+1);
end
co=true;
end

Funcion#3
montante
function [ co,X ] = montante6( N,A )
PA=1;
X=NaN;
for k=1:N
L=k;P=A(k,k);
while P==0
L=L+1;
if L>N
co=false;
return
else
for j=k:N+1
T=A(k,j);
A(k,j)=A(L,j);
A(l,j)=T;
end
end
P=A(k,k);
end
for i=1:N
if i~=kc=A(i,k);
for j=k:N+1
A(i,j)=((A(i,j)*P)-(c*A(k,j)))/PA;
end
end
end
PA=P;
end
for i=1:N
X(i)=A(i,N+1)/P;
end
co=true;
end

ans =
La composicion en el plato 1.000000 es de1.800000e-002

ans =
La composicion en el plato 2.000000 es de 4.319999e-004

ans =
La composicion en el plato 3.000000 es de 1.036786e-005

ans =
La composicion en el plato 4.000000 es de2.486887e-007

ans =
La composicion en el plato 5.000000 es de 5.828641e-009
es el mismo resultado en cualquiera de los 2 metodos



PROGRAMA PRINCIPAL
Programa principal
%programa principal columnaplatos
clc;
N=5;
A(1,1)=-512;
A(1,2)=500;
A(1,3)=0;
A(1,4)=0;
A(1,5)=0;
A(1,6)=-9;
A(2,1)=12;
A(2,2)=-512;
A(2,3)=500;
A(2,4)=0;
A(2,5)=0;
A(2,6)=0;
A(3,1)=0;
A(3,2)=12;
A(3,3)=-512;A(3,4)=500;
A(3,5)=0;
A(3,6)=0;
A(4,1)=0;
A(4,2)=0;
A(4,3)=12;
A(4,4)=-512;
A(4,5)=500;
A(4,6)=0;
A(5,1)=0;
A(5,2)=0;
A(5,3)=0;
A(5,4)=12;
A(5,5)=-512;
A(5,6)=0;
[co,X] = jordan6(...