z = 6;
w = 5;
syms x y
eqn1 = x + y + z + w == 13;
eqn2 = 2*x + 3*y - w == -1;
A = [6 5; zeros(10,2)]; % create matrix with x and y values as
L = 10;
k = 2;
% using while loop
while L > 0
sol = solve([eqn1, eqn2], [x,y]);
xSol = (sol.x);
ySol = (sol.y);
A(k,:) = [double(xSol) double(ySol)]; % use double
Z = double(xSol);
W = double(ySol);
syms z w x y;
eqn1 = x + y + z + w == 13;
eqn2 = 2*x + 3*y - w == -1;
eqn1 = subs(eqn1,[z,w],[Z,W]); % substitute new values for z and w in eqn1
eqn2 = subs(eqn2,[z,w],[Z,W]); % substitute new values for z and w in eqn2
k = k+1;
L = L-1;
end
A % final matrix with appended values
