Basically, I now how to create the equation of gbm(in Matlab for example). But what if I want to perform Monte-Carlo simulation using GBM with stochastic volatility represented by GARCH(1,1). How GBM/GARCH(1,1) SDE looks and what are peculiarities of GBM/GARCH(1,1) modelling ?