Gaussian_Elimination

clc;
c=[1 2 3 4;1 3 5 7;-1 -2 5 6;2 3 7 8];
b=[30 50 34 61];
b=transpose(b);
a=[c b]
a(2,:)=a(2,:)-a(1,:)/a(1,1)*a(2,1);
a(3,:)=a(3,:)-a(1,:)/a(1,1)*a(3,1);
a(4,:)=a(4,:)-a(1,:)/a(1,1)*a(4,1);
a(3,:)=a(3,:)-a(2,:)/a(2,2)*a(3,2);
a(4,:)=a(4,:)-a(2,:)/a(2,2)*a(4,2);
a(4,:)=a(4,:)-a(3,:)/a(3,3)*a(4,3);
x(4)=a(4,5)/a(4,4)
x(3)=(a(3,5)-a(3,4)*x(4))/a(3,3)
x(2)=(a(2,5)-a(2,4)*x(4)-a(2,3)*x(3))/a(2,2)
x(1)=(a(1,5)-a(1,4)*x(4)-a(1,3)*x(3)-a(1,2)*x(2))/a(1,1)