matlab 符号方程组求解, 提取符号方程组变成矩阵,循环迭代增长矩阵

clc
clear

%% 符号方程组求解
% syms a b c
% y1=a*b+2==c;
% y2=a+b==c;
% [a,b]=solve(y1,y2,a,b)%输出符号解
% c=2;
% subs(a)%输出c为2时候的数值解
% subs(b)%输出c为2时候的数值解

% syms x y; % 定义符号变量  
% eq1 = x^2 + y^2 == 4; % 定义第一个方程  
% eq2 = x^2 - 2*y == 0; % 定义第二个方程  
% % 求解方程组  
% sol = solve(eq1, eq2, x, y);
% 
% disp(sol.x); % 显示x的解  
% disp(sol.y); % 显示y的解

%% 提取符号方程组变成矩阵
syms x y z  a b c

% a=x*x;
% b=y*y;
% c=z*z;
eqns =[x*x+y*y - 2*z*z == 0,
    x*x + y*y + z*z == 1,
    2*y*y - z*z == -5];
eqns2=subs(eqns,[x*x  y*y  z*z],[a b c] )   %%将符号进行替换
eqns1=eqns
vars = [ a b c ];
[A,b] = equationsToMatrix(eqns2,vars)

%%循环迭代增长矩阵
C=[];
B=[];
A=[
    1  2
    3  4];
for i=1:4
%     C=[C;
%           A];
     
      if(i==1)
          B=A;
      else
          B=[B;
                A];
            out=B;  
      end
      
end