This is an excellent work. However, after verifying this code using a group of test cases, I found that the results for most the test cases are not correct when compared with those derived by Maple. Does anyone has the same problem with me?
It's quite helpful to me. Thx
I also hope to know how to extend the code using other quadrature rules like Gauss-Hermite and Gauss-laguerre quadratures.
