clc
input('This is the GAUSS program: a way to find the value of different variables in a matrix')
clc
disp('Hello! /n Welcome to my cyberspace... /n Today is: ') 
date
clc
disp('Please wait while we are processing the current program') 
pause(4)
clc
disp('Now we are back')
pause(2)
clc
disp('This is the desired program you requested')
disp('   The Gaussian elimination method: ')
a=[1,1,1,1,-2;1,2,3,4,-10;2,3,4,1,4;3,4,1,2,-10]          %This is a sample initiation for the matrix a
[s,t]=size(a);   %Rows and Columns will be saved in these two parameters
% Converting to an up-triangular matrix
for w=1:s-1
    for i=w+1:s
        if a(w,w)==0 & w<s
            for n=1:s+1
                 var = a(w+1,n);
                 a(w+1,n)=a(w,n);
                 a(w,n)=var;
             end
         end
% Calculating m(i,j) in the matrix
   m(i,w)=a(i,w)/a(w,w);
      for j=1:s+1
         a(i,j)=a(i,j)-m(i,w)*a(w,j);
     end
   end
end   %End of the first 'for' loop
% Considering the variables 
r(s)=a(s,s+1)/a(s,s);
 for w=s-1:-1:1
     var=0 ;
     for j=w+1:s
         var=var+a(w,j)*r(j);
     end
 r(w)=(a(w,s+1)-var)/a(w,w);
 end
 disp 'The final answer to your matrix is as follow: ';disp(r)
 disp 'This is the end of our program.   Have a nice time'