%Please don't mess around with my code
clc;
R=input('How many variables are there in your system of linear equations? ');
clc;
disp ('Please enter the elements of the Coefficient matrix "A" and the Constant/Known matrix "B" starting from first equation');
A=zeros(R,R);
B=zeros(R,1);
for H=1:R
for i=1:R
fprintf ('Coefficient Matrix Column \b'), disp(i), fprintf ('\b\b Row \b'), disp(H);
I=input(' ');
A(H,i)=A(H,i)+I;
end
fprintf ('Known Matrix Row \b'), disp(H);
J=input(' ');
B(H,1)=B(H,1)+J;
end
clc;
[L,U] = lu(A);
Y = L\B;
X = U\Y;
fprintf ('The Coefficient Matrix "A" is \n'), disp(A);
fprintf ('The Constant/Known Matrix "B" is \n'), disp(B);
fprintf ('The Solution Matrix "X" is \n'), disp(X);
