December 29th, 2011, 1:45 pm
QuoteOriginally posted by: ehYou are probably right. I would argue, however, that the characteristic function approach is more straightforward to implement than the "closed-form" VG formula (loads of horrible special functions to evaluate).Another approach (in C++) is to use modified Bessel of 2nd kind and gamma functions in the awesome formula for VG option price. What is then left is the degenerate hypergeometric function for which there is nothing in Boost but it is a combination of gamma and a benign integral; the latter can be computed with a rough-and-ready quadrature method. A rough guess is it would take < 2 hours and debug to get it working.Unfortunately, there is no Quadrature library in Boost at this moment. //SSRN down, so can't check NIG closed solution, but its pdf is easy to compute.
Last edited by
Cuchulainn on December 28th, 2011, 11:00 pm, edited 1 time in total.