Numerical solution with higher order accuracy for option pricing with stochastic volatility using a geometric transformation