From 7ce331c01ce6eb7b3f5c715a38a24359da9c6ee2 Mon Sep 17 00:00:00 2001 From: Eric Andersen Date: Thu, 22 Nov 2001 14:04:29 +0000 Subject: Totally rework the math library, this time based on the MacOs X math library (which is itself based on the math lib from FreeBSD). -Erik --- libm/ldouble/gdtrl.c | 130 --------------------------------------------------- 1 file changed, 130 deletions(-) delete mode 100644 libm/ldouble/gdtrl.c (limited to 'libm/ldouble/gdtrl.c') diff --git a/libm/ldouble/gdtrl.c b/libm/ldouble/gdtrl.c deleted file mode 100644 index 9a41790cb..000000000 --- a/libm/ldouble/gdtrl.c +++ /dev/null @@ -1,130 +0,0 @@ -/* gdtrl.c - * - * Gamma distribution function - * - * - * - * SYNOPSIS: - * - * long double a, b, x, y, gdtrl(); - * - * y = gdtrl( a, b, x ); - * - * - * - * DESCRIPTION: - * - * Returns the integral from zero to x of the gamma probability - * density function: - * - * - * x - * b - - * a | | b-1 -at - * y = ----- | t e dt - * - | | - * | (b) - - * 0 - * - * The incomplete gamma integral is used, according to the - * relation - * - * y = igam( b, ax ). - * - * - * ACCURACY: - * - * See igam(). - * - * ERROR MESSAGES: - * - * message condition value returned - * gdtrl domain x < 0 0.0 - * - */ - /* gdtrcl.c - * - * Complemented gamma distribution function - * - * - * - * SYNOPSIS: - * - * long double a, b, x, y, gdtrcl(); - * - * y = gdtrcl( a, b, x ); - * - * - * - * DESCRIPTION: - * - * Returns the integral from x to infinity of the gamma - * probability density function: - * - * - * inf. - * b - - * a | | b-1 -at - * y = ----- | t e dt - * - | | - * | (b) - - * x - * - * The incomplete gamma integral is used, according to the - * relation - * - * y = igamc( b, ax ). - * - * - * ACCURACY: - * - * See igamc(). - * - * ERROR MESSAGES: - * - * message condition value returned - * gdtrcl domain x < 0 0.0 - * - */ - -/* gdtrl() */ - - -/* -Cephes Math Library Release 2.3: March, 1995 -Copyright 1984, 1995 by Stephen L. Moshier -*/ - -#include -#ifdef ANSIPROT -extern long double igaml ( long double, long double ); -extern long double igamcl ( long double, long double ); -#else -long double igaml(), igamcl(); -#endif - -long double gdtrl( a, b, x ) -long double a, b, x; -{ - -if( x < 0.0L ) - { - mtherr( "gdtrl", DOMAIN ); - return( 0.0L ); - } -return( igaml( b, a * x ) ); -} - - - -long double gdtrcl( a, b, x ) -long double a, b, x; -{ - -if( x < 0.0L ) - { - mtherr( "gdtrcl", DOMAIN ); - return( 0.0L ); - } -return( igamcl( b, a * x ) ); -} -- cgit v1.2.3