@karlnapf 
the methods adapted from quadgk.m is written by @votjakovr, which is used in CLogitLikelihood::get_log_zeroth_moments().
I guess the implementation is a well-know algorithm. see cf: http://en.wikipedia.org/wiki/Gauss%E2%80%93Kronrod_quadrature_formula