Metodos_Iterativos

%GaussSeidel
clc
clear
a11=input('ingrese un valor para a11: ');
a21=input('ingrese un valor para a21: ');
a31=input('ingrese un valor para a31: ');
a12=input('ingrese un valor para a12: ');a22=input('ingrese un valor para a22: ');
a32=input('ingrese un valor para a32: ');
a13=input('ingrese un valor para a13: ');
a23=input('ingrese un valor para a23: ');
a33=input('ingrese un valor paraa33: ');
b1=input('ingrese un valor para b1: ');
b2=input('ingrese un valor para b2: ');
b3=input('ingrese un 5valor para b3: ');
a11p=abs(a21+a31);
a22p=abs(a12+a32);
a33p=abs(a13+a23);

ifabs(a11)>=a11p && abs(a22)>=a22p && abs(a33)>=a33p
x1=0;
x2=0;
x3=0;
e=1;
n=1;
tol=0.01;
while e>tol
xn1=(b1 - a31*x3 - a21*x2)/a11;
e= max([abs(xn1-x1)])
x1=xn1
xn2=(b2 - a12*x1 -a32*x3)/a22;
e= max([abs(xn1-x1) abs(xn2-x2)])
x2=xn2
xn3=(b3 - a23*x2 - a13*x1)/a33;
e= max([abs(xn1-x1) abs(xn2-x2) abs(xn3-x3)])
x3=xn3
n=n+1
end
iteraciones=n
else
fprintf('elsistema no es diagonalmente dominate')

end

%Jacobi
clc
clear
a11=input('ingrese un valor para a11: ');
a21=input('ingrese un valor para a21: ');
a31=input('ingrese un valor para a31: ');a12=input('ingrese un valor para a12: ');
a22=input('ingrese un valor para a22: ');
a32=input('ingrese un valor para a32: ');
a13=input('ingrese un valor para a13: ');
a23=input('ingrese un valor para a23:');
a33=input('ingrese un valor para a33: ');
b1=input('ingrese un valor para b1: ');
b2=input('ingrese un valor para b2: ');
b3=input('ingrese un 5valor para b3: ');

a11p=abs(a21+a31);a22p=abs(a12+a32);
a33p=abs(a13+a23);

if abs(a11)>=a11p && abs(a22)>=a22p && abs(a33)>=a33p

x1=0;
x2=0;
x3=0;
tol=0.01
e=1;
while e > tol
xn1=(b1 - a31*x3 - a21*x2)/a11;xn2=(b2 - a12*x1 - a32*x3)/a22;
xn3=(b3 - a23*x2 - a13*x1)/a33;
e= max([abs(xn1-x1) abs(xn2-x2) abs(xn3-x3)]);
x1=xn1
x2=xn2
x3=xn3
end
else
fprintf('el sistema no es...