syms x y
eqns = [x+2*y==2;
        3*x+y==1,];
[A,b] = equationsToMatrix(eqns)
a=[A,b];
disp(a);
q=input('value');
w=input('value 2');
e=input(' value 3 ');
fprintf("step -1\n");
for i=1:3
a(1,i)=a(1,i)*(q);
end
disp(a);
fprintf("step -2\n");
for i=1:3
a(2,i)=a(2,i)-(a(1,i)*w);
end
disp(a);
fprintf("step -3\n");
for i=1:3
a(1,i)=a(1,i)+a(2,i);
end
disp(a);
fprintf("step -4\n");
for i=1:3
a(2,i)=a(2,i)/e;
end
disp(a);
p=a(1,3);
o=a(2,3);
fprintf("X = %d\n" ,p);
fprintf("Y = %d\n",o);