主要你是第一个自定义函数设置有问题，才出现结果不一样。应这样改
syms x
F=@(x)[exp(-exp(-x(1)+x(2))) - x(2)*(1+x(1)^2);x(1)*cos(x(2)) + x(2)*sin(x(1)) - 0.5];
x0=[0;0];
fsolve(F,x0)
说明：第一个结果是由上述代码得到。第二个结果是用F = root2d(x)函数得到。
