* add function lbinom that returns the logarithm of a binom
* calculate with doubles not __float128

This is a huge performance increase, but it will make the program crash for
high values of lmax. I will fix that later.