function q=fxy(p)
x=p(1);y=p(2);
q(1)=2*x-y-exp(-x);
q(2)=-x+2*y-exp(-y);
end
>> x=fsolve(@fxy,[-5,-5]',optimset('Display','off'))
x =
0.5671
0.5671
fsloⅴe要求自定义函数的格式可以为
F=[2*x(1)-x(2)-exp(-x(1));-x(1)+2*x(2)-exp(-x(2))]
也可以为
F1=2*x(1)-x(2)-exp(-x(1));
F2=-x(1)+2*x(2)-exp(-x(2));
F=[F1 F2]
也可以为
fun=@(x)[2*x(1)-x(2)-exp(-x(1));-x(1)+2*x(2)-exp(-x(2))];
>>x0=[-5;-5]
>>[x,h]=fsolve(fun,x0)