定义包含非线性方程组的匿名函数，然后用fsolve求解匿名函数中的非线性方程组

fun=@(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];%设定初值
x=fsolve(fun,x0)%调用fsolve函数，求解方程组

计算结果：

x =

    0.3532    0.6061

该方式的优点是不用单独建立方程组函数，且便于开展for循环计算