MATLAB数值求解带参数方程

MATLAB里面数值求解方程的函数为fsolve，MATLAB自带的帮助文档里有如何直接使用的例子，但没有如何解带参数方程的例子，其实只需要对待求解方程稍作修改即可。比如我们求如下关于$x_1$和$x_2$的方程：
$$\begin{gathered} e^{-e^{-(x_1+x_2)}}-x_2(1+x_1^2)-c_1=0 \\ {x}_1\cos (x_2)+x_2\sin (x_1)-c_2=0 \end{gathered}$$

其中$c_1$和$c_2$为外部设定的参数，那么首先自定义一个函数：

function myfun1 = root2d(x,c1,c2)

myfun1(1) = exp(-exp(-(x(1)+x(2)))) - x(2)*(1+x(1)^2)-c1;
myfun1(2) = x(1)*cos(x(2)) + x(2)*sin(x(1)) -c2;

然后在主程序文件中给定$c_1$和$c_2$的值，以及$x_1$和$x_2$的迭代初值，利用@符号将待求解变量进行区分，即可调用fsolve进行求解。

c1=0;
c2=0.5;
x0 = [0,0];
fun = @(x)root2d(x,c1,c2);
x = fsolve(fun,x0)