fid = fopen('class2.output','w');
a=[2 1 0;6 2 6;-4 -3 9];
%
%
%
e1=eye(3);e1(2,1)=-3;e1(3,1)=2;
a1=e1*a;
disp(a1);
%
%
%
e2=eye(3);e2(3,2)=-1;
a2=e2*a1;
disp(a2);
%
%
%
u=a2;
l=inv(e1)*inv(e2);
disp(l);
disp(l*u);
%
%
%
disp(l*u);
%
%  Use of the diagonalization
%
b=[1;5;3];
bsolve=inv(a)*b;
disp(bsolve);
%
% Multiplying the b matrix by the inverse of l
%
b1=e1*b;
disp([a1 b1]);
%
%
%
b2=e2*b1;
disp([a2 b2]);
%
%
%
bLU=b2;
bLU(3,1)=b2(3,1)/3;
bLU(2,1)=-(b2(2,1)-6*bLU(3,1));
bLU(1,1)=(b2(1,1)-bLU(2,1))/2;
disp(bLU);
%
%
%
a=[7 3 1; 3 4 0; 1 0 2];
[v,d]=eig(a);
disp(d);
%
%
%
e1=eye(3);e1(2,1)=-3/7;e1(3,1)=-1/7;
a1=e1*a;
disp(a1);
%
%
%
e2=eye(3);e2(3,2)=-a1(3,2)/a1(2,2);
a2=e2*a1;
disp(a2);
%
%
%
u=a2;
l=inv(e1)*inv(e2);
disp(l);
disp(l*u);
%
%
%
t=zeros(3,3);t(1,1)=1/u(1,1);t(2,2)=1/u(2,2);t(3,3)=1/u(3,3);
ut=t*u;
disp(ut);