examp6_6.m

来自「《高等应用数学问题MATLAB求解》的PDF版及习题解答」· M 代码 · 共 7 行

syms x y;
[x,y]=solve('x^2/2+x+3/2+2/y+5/(2*y^2)+3/x^3=0',...
    'y/2+3/(2*x)+1/x^4+5*y^4','x,y');
size(x)

err=[x.^2/2+x+3/2+2./y+5./(2*y.^2)+3./x.^3,y/2+3./(2*x)+1./x.^4+5*y.^4];
norm(double(eval(err)))