# 基于Black - Scholes 公式的期权定价公式

def bsm_call_value(s,,k,T,r,sigma):

from math import log,sqrt,exp

from scipy import stats

s=float(s)

d1=(log(s/k)+(r+0.5*sigma**2)*T)/(sqrt(T))

d2=(log(s/k)+(r-0.5*sigma**2)*T)/(sqrt(T))

value=(s*stats/norm.cdf(d1,0.0,1.0))-K*exp(-r*T)*stats.norm.cdf(d2,0.0,1.0)

return value

当前价 s

行权价 k

上面的函数只是计算隐含波动率所需的基本函数。当然，还需要单独的期权报价，也就是基于VSTOXX指数的欧式看涨期权，以及生成隐含波动率的代码，下面就用IPython来完成。

比如把参考日期t=0，

指数的收盘价为v=7，6639

对于无风险短期利率，假定年率r=0,01

这样就可以利用

print '期权价格 : %.4f' % bsm_call_value(s,,k,T,r,sigma):