TAJD

function genS_jl(I)
    s0 = 600.0
    r = 0.02
    sigma = 2.0
    T = 1.0
    M = 100
    dt = T/M
    a = (r - 0.5*sigma^2)*dt
    b = sigma*sqrt(dt)