Abstract:
SIMGBM returns a vector of a sample trajectory of GBM on the time interval [0,N]: dX(t) = MU*X(t)*dt + SIGMA*X(t)*dW(t), given starting value of the process X0, drift MU, volatility SIGMA, time step size DELTA, array of normally distributed pseudorandom numbers NO (array NO is simulated if not provided as an input variable) and method (direct integration, Euler scheme, Milstein scheme, 2nd order Milstein scheme).