This is a R port of the original work available for Matlab.
It is an implementation of the Heston option pricing function based on the work by Christian Kahl, Peter Jäckel and Roger Lord. The standard Heston model is characterized by high numerical instability of the integral which has been reduced in this implementation.
TODO