图这么小，大家可以下载桌面放大看，谢谢．

是这样的：

N=8.67*10^10,M=0.3*10^10,r=0.0414%,T=2,Os=0.5499,q=0.5,S=17.32,K=19.68;

d1=((ln(V/(N*K))+r*T+0.5*T*o^2)/(o*sqrt(T));
d2=d1-o*sqrt(T);
-------------------------------------------------------------------------------------------

V＝S*N+V*((q*M)/(N+M*q))*NormalCDF(d1)-N*((q*M)/(N+M*q))*exp(-r*T)*K*NormalCDF(d2);
o=S*Os*N/(V-(V*q*M/(N+M*q))*Normal(d1));

求V,o---

然后我编代码：

Constant N=8.67*10^10,M=0.3*10^10,r=0.0414%,T=2,Os=0.5499,q=0.5,S=17.32,K=19.68;
Conststr d1=((ln(V/(N*K))+r*T+0.5*T*o^2)/(o*sqrt(T));d2=d1-o*sqrt(T);
Function
V-S*N-V*((q*M)/(N+M*q))*NormalCDF(d1)+N*((q*M)/(N+M*q))*exp(-r*T)*K*NormalCDF(d2);
o-S*Os*N/(V-(V*q*M/(N+M*q))*Normal(d1));


对不对？？不对的话，这么编？？

改动如下：

Constant N=8.67*10^10,M=0.3*10^10,r=0.0414%,T=2,Os=0.5499,q=0.5,S=17.32,K=19.68;
Conststr d1=((ln(V/(N*K))+r*T+0.5*T*o^2))/(o*sqrt(T)),d2=d1-o*sqrt(T);
Function V-S*N-V*((q*M)/(N+M*q))*NormalCDF(d1)+N*((q*M)/(N+M*q))*exp(-r*T)*K*NormalCDF(d2);
                 o-S*Os*N/(V-(V*q*M/(N+M*q))*NormalCDF(d1));

3.0下运行没问题。结果：

v: 1508382505829.07
o: 0.553015355345716