求解矩阵方程组
 

a1=-0.4;b1=1;c1=0.84;d1=0.4;
syms a b c d;
A1=[a b;c d];
A1=solve(kron(eye(2),a1)+kron([A1*inv(eye(2)-d1*A1)],(b1*c1)))
[COLOR="Red"]??? Error using ==> sym.maple
at offset 185, `;` unexpected
Error in ==> sym.transpose at 18
B = maple('transpose',A);
Error in ==> sym.findsym at 26
sc = char(S(:).');
Error in ==> solve at 99
vars = ['[' findsym(sym(eqns),neqns) ']'];
Error in ==> sym.solve at 49
[varargout{1:max(1,nargout)}] = solve(S{:});[/COLOR]
A1是个2*2的矩阵
其实还有这样类似的多个方程 后面的需要用到前面的解