I'm looking to simulate the stochastic price and volatility process (Heston model) using some form of Euler method for Monte Carlo approximation of option prices. The results that I get are acceptable for deep in the money options and at the money options but not very satisfying at all for deep out of the money options. I want to reduce the variance for faster convergence and the importance sampling method seems suitable but the problem doesn't seem to be trivial at all. Does anyone have an idea