#求助 fsolve指令使用遇到问题

先谈一些我认识的fsolve指令

近期在使用fsolve指令求解一个非线性代数方程式，总共有2002个未知数，2002个方程，我有自己去了解一些关于fsolve的用法，首先他主要有三种算法，trust-region-dogleg（信任区间狗腿法，默认方法）、trust-region（信任区间法）、levenberg-marquardt（简称LM），其他的帮助会有在其他博主的文章中找到，再次不赘述

讲述一下自己的问题

我将自己的函数放在b1的matlab中，

function F=b1(w)
global dt a nx y0 g S0
F(1:nx,1)=1/dt*a*w(1:nx)+1/2*w(nx+1:2*nx)-y0(1:nx)-w(2*nx+1)*ones(nx,1);
F(nx+1:2*nx,1)=1/dt*a*w(nx+1:2*nx)+1/2*w(nx+1:2*nx).^2./w(1:nx)+g/4*w(1:nx).^2...
 -1/2*a*g*w(1:nx).*S0-y0(nx+1:2*nx)-w(2*nx+2)*ones(nx,1);
F(2*nx+1,1)=w(nx+1)-4.42;
F(2*nx+2,1)=w(nx)-2;

global的指令是全域变数，即变数dt、a、nx、y0、g、S0在b1和case1.m文件中同时使用。
在case1.m中的指令如下

options = optimoptions( 'fsolve' , 'Display' , 'it