exm020603.m

来自「matlab电子书籍」· M 代码 · 共 8 行

A=sym([1 1/2 1/2 -1;1 1 -1 1;1 -1/4 -1 1;-8 -1 1 1]);
b=sym([0;10;0;1])
X1=A\b	

S=solve('d+n/2+p/2-q','d+n-p+q-10','d-n/4-p+q','-8*d-n+p+q-1');
disp(['  d','  n','  p','  q'])
disp([S.d,S.n,S.p,S.q])