Serving the Quantitative Finance Community

 
User avatar
babyboy
Topic Author
Posts: 1
Joined: March 6th, 2004, 12:58 pm

Compare 2 Heston

March 29th, 2004, 6:17 pm

Please,I'v 2 pricers for pricing vanilla options with Heston model :and the 2 models gives me differents results.Wich one is the good? or the more robust?? And what does that mean : in Heston93.xls : Integral ( fct(k),k=0..64 ) with 8 subintervals of length = 8 Thanks by advance
Attachments
Heston93_pdf.zip.zip
(95 KiB) Downloaded 70 times
Heston m.zip
(34.86 KiB) Downloaded 61 times
 
User avatar
mucki
Posts: 2
Joined: July 29th, 2002, 6:47 pm

Compare 2 Heston

March 30th, 2004, 12:32 pm

Hi,I compraed the heston93_pdf to my algorithm: If rho = 0 it seems that my algorithm returns the same results, for rho different from 0 I obtain different prices. However, your worksheet seems to apply Gauss Legendre algorithm to compute the integral while I am using Gauss Laguerre (150 nodes). Don't know wthere my prices are better or yours. Pls let me know when you learn more about it.
 
User avatar
AVt
Posts: 90
Joined: December 29th, 2001, 8:23 pm

Compare 2 Heston

March 30th, 2004, 6:05 pm

Hallo Mucki,it is the same code as in free code for Heston. Would you mind to give your figures and look at examples there thus having all at one place?
 
User avatar
mucki
Posts: 2
Joined: July 29th, 2002, 6:47 pm

Compare 2 Heston

March 31st, 2004, 8:14 am

Here's my excel sheet. The DLL file contains the code for the heston formula and must be copied into the same directory.
 
User avatar
mucki
Posts: 2
Joined: July 29th, 2002, 6:47 pm

Compare 2 Heston

March 31st, 2004, 8:15 am

Sorry this hasn't worked out the first time
Attachments
Heston.zip
(157.66 KiB) Downloaded 60 times
 
User avatar
AVt
Posts: 90
Joined: December 29th, 2001, 8:23 pm

Compare 2 Heston

March 31st, 2004, 11:33 am

Up to discounting it is the same (which is missing as "mib" pointed out in the other thread).The only thing might be numerical problems with high order Laguerre: one comes throughthe oscillating integral (see below) and the other through significant digits and roots forthe Laguerre polynomial: for order 64 the 'last' root is around 250 (i do not know it for a degree of 150) where it islikely that you get rounding errors _even if you have the exact root_ (which itself may bea problem). Otherwise stated: high order does not mean high exactness, it might result inthe converse. As i have only 15 digits in Excel i dicided to take a lower order.For highly oscillating integrands your nodes may be placed not that well w r t periodics,and for the following input the DLL gives a price of 29,9999999832 - below 30:Stock 100Strike 70Current Variance 0,04Long term Variance 0,04Mean reversion speed 1,15Volvol 0,39Interest Rate 0Time to expiry 0,03Correlation -0,64Market Price of Vol Risk 0Practical that is not a deep problem, but i think you will have exacter results takinglower order (which can adopt to the periodics) and do repeated finite integrationsand if you like Gauss Laguerre do one at the end.