import matplotlib.pyplot as pltimport numpy as np#rect=[0.1,5.0,0.1,0.1]fig=plt.figure()#time spanT=2#drift factor飘移率mu=0.1#volatility波动率sigma=0.04#t=0初试价S0=20#length of stepsdt=0.01N=round(T/dt)t=np.linspace(0,T,N)#布朗运动W=np.random.standard_normal(size=N)W=np.cumsum(W)*np.sqrt(dt)X=(mu-0.5*sigma**2)*t+sigma*WS=S0*np.exp(X)plt.plot(t,S,lw=2)plt.show()