function[vT,sigmaV]=fun1(D,E,sigmaE,r,sigmaV)
d2=1;
d3=0;
d1=0;
vT=0;
while(1)
sigmaV=E*sigmaE/(D*normcdf(d1)*exp(-r)+E);
d1=d2+sigmaV;
vT=E*sigmaE/(normcdf(d1)*sigmaV);
d3=(log(vT/D)+(r+0.5*sigmaV*sigmaV))/sigmaV-sigmaV;
if(abs(d3-d2)<0.0001)
break;
end
d2=d3;
end
这是源代码，当我调用该代码时第一次我随便带了几个值可以得到结果，可第二次再调用时就没有结果了，具体情况如下截图所示。 求指点，不甚感激。期末论文没几天就要交了，急！！！