library(compoisson)


data(insurance)
fit = com.fit(Lemaire)
z = com.compute.z(fit$lambda, fit$nu)
com.log.density(0, fit$lambda, fit$nu, fit$z)
a = exp(com.log.difference(log(100), log(20))); # a = 80
b = exp(com.log.sum(log(100), log(20))); # b = 120
c = exp(com.log.factorial(4)); # c = 24