定義包含非線性方程組的匿名函數，然後用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循環計算