求某期权获利的概率

 1 from scipy import stats
 2 import numpy as np
 3 
 4 # np.log = ln 求自然对数
 5 
 6 # stats.norm.cdf(d1, 0.0, 1.0) = N 累积正态分布
 7 
 8 v = 0.195
 9 t = 62 / 365
10 Vt = v * np.sqrt(t)
11 print(Vt)
12 
13 q = 3.5168 # 盈亏平衡点
14 p = 3.254 # 标的物市价
15 P = stats.norm.cdf(np.log(q / p) / Vt, 0.0, 1.0)
16 
17 print('在62天后低于{}的可能性是{}'.format(q, P))

感觉离成功越来越近！！