Matthias Schabel writes: > return Gamma(a)*x^a*gamma_star(a,x);
This looks complex-valued for x<0 and non-integer a. The function gamma_star would be ok though. -- regards, Brian Gough Network Theory Ltd, Publishing Free Software Manuals --- http://www.network-theory.co.uk/ _______________________________________________ Bug-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-gsl
