fixed Gauss-Lobatto integration for zero function (thanks to Peter
Caspers)
added wrapper for boost cumulative normal distribution function