Guaramy <Guaramy <at> hotmail.com> writes: > [SNIP] > i think the problem is the gamma function, does anyone know how to compute > gamma with imaginary numbers? Try lngamma_complex() in the gsl package? Ben Bolker