Codificacion matlab

%Método de Gauss 7 variables
clear
clc
A=[2 3 4 -5 6 8 7 -6;
-3 2 3 6 -4 3 1 3;
4 -5 8 -7 4 1 -1 9;
7 4 6 6 5 6 -5 7;
5 2 -7 5 3 8 2 -8;
4 3 1 4 -1 9 4 2;
9 8 4 -3 1 6 3 -12]

%Eliminación
for C=1:6
for F=C + 1:7
A(F,:)=(A(C,:)/A(C,C))*(-A(F,C))+A(F,:);
end
end
%Sustitución
for n=7:-1:1
s=0;for j= n+1:7
s=s+A(n,j)*x(j);
end
x(n)=(A(n,8)-s)/A(n,n);
end



%Método de Gauss 7 variables
clear
clc
A=[2 3 4 -5 6 8 7 -6;
-3 2 3 6 -4 3 1 3 ;
4 -5 8 -7 4 1 -1 9;
7 4 6 6 5 6 -5 7;
5 2 -7 5 3 8 2 -8;
4 3 1 4 -1 9 4 2;
9 8 4 -3 1 6 3 -12]
%Eliminación
for C=1:6
for F=C + 1:7
A(F,:)=(A(C,:)/A(C,C))*(-A(F,C))+A(F,:);
end
end%Sustitución

x=zeros(F,1);
x(F)=A(F,C)/A(F,F);
for i= F-1:-1:1
disp('Los valores son :')
x(i)=(A(i,C)-A(i,i+1:F)*x(i+1:F))/A(i,i)

end


%metodo de gausss para 2 variables
clc
a=[2 5
4 -2]
b=[10;6]
A=[a,b]
%Paso 1:eliminacion hacia adelante
A(2,:)=A(1,:)*(-A(2,1)/A(1,1))+A(2,:)
%Paso 2:sustitucion hacia atras
y=A(2,3)/A(2,2)
x=(A(1,3)-A(1,2)*y)/A(1,1)%metodo de gausss para 3 variables
clc
a=[2 1 5
4 3 7
6 8 1]
b=[13;27;32]
A=[a,b]
%Paso 1:eliminacion hacia adelante
A(2,:)=A(1,:)*(-A(2,1)/A(1,1))+A(2,:)
A(3,:)=A(1,:)*(-A(2,1)/A(1,1))+A(3,:)
A(3,:)=A(2,:)*(-A(3,2)/A(2,2))+A(3,:)



function z=Gauss2v(a,b)
%metodo de gausss para 2 variables
A=[a,b]
%Paso 1:eliminacion hacia adelanteA(2,:)=A(1,:)*(-A(2,1)/A(1,1))+A(2,:);
%Paso 2:sustitucion hacia atras
y=A(2,3)/A(2,2);
x=(A(1,3)-A(1,2)*y)/A(1,1);
z=[x,y];



%metodo de pivoteo
function z=Pivo2v(a,b)
%metodo de gausss para 2 variables
A=[a,b];
%Paso 1:eliminacion hacia adelante
if(abs(A(1,1))>abs(A(2,1)))
A(2,:)=A(1,:)*(-A(2,1)/A(1,1))+A(2,:);
else
z=A(1,:);
A(1,:)=A(2,:);
A(2,:)=z;A(2,:)=A(1,:)*(-A(2,1)/A(1,1))+A(2,:);
end
%Paso 2:sustitucion hacia atras
y=A(2,3)/A(2,2);
x=(A(1,3)-A(1,2)*y)/A(1,1);
z=[x,y];


%Gráfica de 2 funciones mulivariable

x=-5:0.01:5;

y=(5*x-4)/2;
plot (x,y)
grid on
hold on

x=-5:0.01:5;
y=(-3*x+10)/4;
plot (x,y,'r')
grid on
hold off

clc
%METODO DE JACOBI
x=0;y=0;z=0;error=10;tol=0.001
disp(' x1 y1 z1 ') %ponemos la etiquetawhile(error>tol)
x1=(10-2*y+3*z)/8; %ponemos funciones para hallar nuevo valor
y1=(-4+2*y+z)/5;
z1=(6-x-2*y)/-10;
error=sqrt((x1-x)^2+(y1-y)^2+(z1-z)^2);
%para cada valor reporte los errores
disp([x1,y1,z1,error])
x=x1;y=y1;z=z1;
end


clc
%METODO DE JACOBI
%con 3 digitos significativos
x=0;y=0;z=0;error=10;tol=0.5*10^(-3);
disp(' x1 y1 z1 ') %ponemos la etiquetawhile(error>tol)
x1=(10-2*y+3*z)/8; %ponemos funciones para hallar nuevo valor
y1=(-4+2*y+z)/5;
z1=(6-x-2*y)/-10;

er1=abs((x1-x)/x1);
er2=abs((y1-y)/y1);
er3=abs((z1-1)/z1);

%de estos errores buscamos el mayor
error=max([er1 er2 er3]);

error=sqrt((x1-x)^2+(y1-y)^2+(z1-z)^2);
%para cada valor reporte los errores
disp([x1,y1,z1,error])

x=x1;y=y1;z=z1;
end


clc
%metodo dejacobi-matricial
error=20;tol=0.001;
X0=[0 0 0]'
M=[ 0 -1/4 3/8
2/5 0 1/5
1/10 1/5 0 ]
C=[10/8 -4/5 -3/5]'
disp(' X Y Z error')
while(error>tol)
X1=M*X0+C;
error=norm(X1-X0);
disp([X1',error])
X0=X1;
end


clc
%metodo de gauss-seidel
x=0;y=0;z=0;error=10;tol=0.001
disp(' x1 y1 z1 ') %ponemos la etiqueta
while(error>tol)x1=(10-2*y+3*z)/8; %ponemos funciones para hallar nuevo valor
y1=(-4+2*x1+z)/5;
z1=(6-x1-2*y1)/-10;
error=sqrt((x1-x)^2+(y1-y)^2+(z1-z)^2);
%para cada valor reporte los errores
disp([x1,y1,z1,error])
x=x1;y=y1;z=z1;
end


clc
%para graficar sistemas de funciones no lineales-explicita
x=-0.5:0.01:10;
y=sqrt(10*x-x.^2-8);
plot(x,y)
hold on %para activarlo
plot(x,-y)
%grafica de...