我想用fsolve在MatLab中解一个非线性方程组，但我也要对两个变量的函数求导。以下是步骤中的问题：
步骤1：定义非线性函数系统F：

function F = root2d(x)
 
F(1) = (exp(-x(1)+2)/(1+exp(-x(1)+2)+exp(-x(2)+1)));
F(2) = (exp(-x(2)+1)/(1+exp(-x(1)+2)+exp(-x(2)+1)));

end

步骤2：我想找出F(1)关于x(1)的导数和F(2)关于x(2)的导数，使得：

function F = root2d(x)
 
F(1) = (exp(-x(1)+2)/(1+exp(-x(1)+2)+exp(-x(2)+1)));
F(2) = (exp(-x(2)+1)/(1+exp(-x(1)+2)+exp(-x(2)+1)));

f1_d=diff(F(1),x(1))
f1_d=diff(F(2),x(2))

end

步骤3：我希望我的函数是原始函数加上导数：

function F = root2d(x)
 
F(1) = (exp(-x(1)+2)/(1+exp(-x(1)+2)+exp(-x(2)+1)));
F(2) = (exp(-x(2)+1)/(1+exp(-x(1)+2)+exp(-x(2)+1)));

f1_d=diff(F(1),x(1));
f1_d=diff(F(2),x(2));

F(1)=F(1)+f1_d;
F(2)=F(2)+f2_d;

end

步骤4：在主文件中，我将在fsolve中使用此函数来求解x(1)和x(2)的非线性方程组：

syms x
fun = @root2d;
x0 = [0,0];
x = fsolve(fun,x0);
root2d(x)

function F = root2d(x)
 
F(1) = (exp(-x(1)+2)/(1+exp(-x(1)+2)+exp(-x(2)+1)));
F(2) = (exp(-x(2)+1)/(1+exp(-x(1)+2)+exp(-x(2)+1)));

f1_d=diff(F(1),x(1));
f1_d=diff(F(2),x(2));

F(1)=F(1)+f1_d;
F(2)=F(2)+f2_d;
end

你能帮帮我吗，怎么才能重写我的代码，以便解导数，然后再计算和解出系统呢？手工计算导数是不可行的，因为稍后我必须开发这个程序，我将拥有大约100个方程，而不是2个。