Since I wrote the post about Gaussian quadrature rules I’ve been using this integration procedure quite a lot. At last I’ve decided to make a Haskell package providing a non-adaptive Gaussian quadrature for numeric integration. I did it primary just for myself, so I could use it wherever and whenever I need. However, maybe someone else will find it useful.
As usual, if you find some bugs or want to request a new feature or just have questions, I’ll be glad to help.