The first numerical experiments are performed using the BSM model [1]}, [2]}. The stock dynamics driven by the BSM model (a geometric Brownian process) are given by
\(S_T=S_0e^{(r-q-{1\over 2}\sigma ^2)T+\sigma W_T},\) 
