13f708241SDavid Schultz 23a8617a8SJordan K. Hubbard /* 33a8617a8SJordan K. Hubbard * ==================================================== 43a8617a8SJordan K. Hubbard * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 53a8617a8SJordan K. Hubbard * 63f708241SDavid Schultz * Developed at SunSoft, a Sun Microsystems, Inc. business. 73a8617a8SJordan K. Hubbard * Permission to use, copy, modify, and distribute this 83a8617a8SJordan K. Hubbard * software is freely granted, provided that this notice 93a8617a8SJordan K. Hubbard * is preserved. 103a8617a8SJordan K. Hubbard * ==================================================== 113a8617a8SJordan K. Hubbard * 123a8617a8SJordan K. Hubbard */ 133a8617a8SJordan K. Hubbard 14*99843eb8SSteve Kargl /* gamma(x) 153a8617a8SJordan K. Hubbard * Return the logarithm of the Gamma function of x. 163a8617a8SJordan K. Hubbard * 17*99843eb8SSteve Kargl * Method: call gamma_r 183a8617a8SJordan K. Hubbard */ 193a8617a8SJordan K. Hubbard 203a8617a8SJordan K. Hubbard #include "math.h" 213a8617a8SJordan K. Hubbard #include "math_private.h" 223a8617a8SJordan K. Hubbard 233a8617a8SJordan K. Hubbard extern int signgam; 243a8617a8SJordan K. Hubbard 25a82bbc73SAlfred Perlstein double gamma(double x)26*99843eb8SSteve Karglgamma(double x) 273a8617a8SJordan K. Hubbard { 28*99843eb8SSteve Kargl return gamma_r(x,&signgam); 293a8617a8SJordan K. Hubbard } 30