mat= [10,2,3;4,11,6;7,8,16];
b = [10,20,40]';
%{
r = input('Enter rows of Matrix\n');
c = input('Enter Column of Matrix\n');
disp('Enter matrix element row major order\n');
mat = zeros(r,c);
for i=1:r
    for j=1:c
        mat(i,j) = input('');
    end
end
disp(mat);
b = zeros(r:1);
disp('Enter The right hand side matrix\n');
for i=1:r
    b(i) = input('');
end
%}
r = input('Enter rows of Matrix\n');
c = input('Enter Column of Matrix\n');
disp('Enter matrix element row major order\n');
for i = 1:r
    fprintf('Enter %d row\n',i);
    mat(i,:) = input('\');
end
disp('Enter The right hand side matrix\n');
for i=1:r
    b(i) = input('');
end

display('Ans from LU Decomposition Method');
LU_Decomposition_1505014( mat,b )
display('Ans from Gauss Siedel Method');
Gauss_Seidel_1505014(mat,b,0.01,50);