Future derivative traders of Quantopian, please see the attached implementation of Heston's (1993) model for pricing options under stochastic volatility as described on Wilmott. The prices I am getting seem to be a bit off, which I think is due to the numerical integration routine or the complex logarithm routine I am using in Python. I would be interested in insights on how to improve the code.