用fsolve计算的一个非线性方程组运算好久没结果

这个是我编的函数：

function F=myfun(x,e,d,r,v)

F=[x(1)-(log2(x(4)/d)+(r+0.5*x(3).^2))/(x(3));

     x(2)-x(1)-x(3);v-x(4)*x(3)*normcdf(x(1),0,1)/e;

e-normcdf(x(1),0,1)*x(4)-normcdf(x(2),0,1)*d*exp(-1*r);];

end

这个是在命令窗口我输入的：

for i=1:817;

e=E(i);r=r(i);d=D(i);v=Evreturnbxlqpricee(i);

x0 = [1,1,1,1];

x(i)=fsolve(@(x)myfun(x,e,d,r,v),x0,optimset('display','iter','MaxFunEvals',10e+1000,'MaxIter',10e+1000));

end
运算好久都没有结果，求大神们指点啊

写了一个小程序，可以解目标方程组。小程序在zip里，略大，以后有空再压缩。已打包成exe文件，不需要搭建环境，即点即用。结果可以代入原方程组检验。

function F=myfun(x,ee,dd,rr,vv)

F=[x(1)-(log2(x(4)/dd)+(rr+0.5*x(3).^2))/(x(3));
      x(2)-x(1)+x(3);vv-x(4)*x(3)*normcdf(x(1),0,1)/ee;
ee-normcdf(x(1),0,1)*x(4)+normcdf(x(2),0,1)*dd*exp(-1*rr);];
for i=1:1:length(e);vp=[1;v];ee=e(i);dd=d(i);rr=r(i);vv=vp(i);x0=[1,1,1,3000000];
options = optimoptions('fsolve','Display','off','Diagnostics','off','TolFun',1e-6,'TolX',1e-6,'MaxFunEvals',1e+100,'MaxIter',1e+100,'PlotFcn',@optimplotfirstorderopt);
[x,fval,exitflag] = fsolve(@(x) myfun(x,ee,dd,rr,vv),x0,options);
x1(i)=x(1);x2(i)=x(2);x3(i)=x(3);x4(i)=x(4);
end

有可能是初值给的不合理，建议改一下初值实验一下。