2018-05-27 py4fi（1）--欧式看涨期权

通过蒙特卡洛模拟方法估计欧式看涨期权的价值

s0 =100 # 初始股票指数水平

K =105  # 欧式看涨期权的行权价格

T =1.0  # 到期时间1年

r =0.05  # 固定无风险短期利率5%

sigma =0.2  # 固定波动率 20%

# 估值算法

I =100000

z = np.random.standard_normal(I)

从标准正态分布中取得I个（伪）随机数z（i）

ST = s0 * np.exp((r-0.5*sigma**2)*T+sigma*np.sqrt(T)*z)

BSM模型

hT = np.maximum(ST-K,0)

计算到期时期权的所有内在价值

C0 = np.exp(-r*T)*sum(hT)/I

欧式期权的蒙特卡洛估算函数

print "Value of the European Call Option %5.3f"% C0

打印结果