function St = g(sigma,sigma_0,lambda,c,del_t)

St = sigma_0 + exp(-lambda*del_t)*(sigma-sigma_0) + ...
          sqrt(c*(1-exp(-2*lambda*del_t))) * randn;