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/double/gdtr.c | 130 ----------------------------------------------------- 1 file changed, 130 deletions(-) delete mode 100644 libm/double/gdtr.c (limited to 'libm/double/gdtr.c') diff --git a/libm/double/gdtr.c b/libm/double/gdtr.c deleted file mode 100644 index 6b27d9abb..000000000 --- a/libm/double/gdtr.c +++ /dev/null @@ -1,130 +0,0 @@ -/* gdtr.c - * - * Gamma distribution function - * - * - * - * SYNOPSIS: - * - * double a, b, x, y, gdtr(); - * - * y = gdtr( 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 - * gdtr domain x < 0 0.0 - * - */ - /* gdtrc.c - * - * Complemented gamma distribution function - * - * - * - * SYNOPSIS: - * - * double a, b, x, y, gdtrc(); - * - * y = gdtrc( 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 - * gdtrc domain x < 0 0.0 - * - */ - -/* gdtr() */ - - -/* -Cephes Math Library Release 2.8: June, 2000 -Copyright 1984, 1987, 1995, 2000 by Stephen L. Moshier -*/ - -#include -#ifdef ANSIPROT -extern double igam ( double, double ); -extern double igamc ( double, double ); -#else -double igam(), igamc(); -#endif - -double gdtr( a, b, x ) -double a, b, x; -{ - -if( x < 0.0 ) - { - mtherr( "gdtr", DOMAIN ); - return( 0.0 ); - } -return( igam( b, a * x ) ); -} - - - -double gdtrc( a, b, x ) -double a, b, x; -{ - -if( x < 0.0 ) - { - mtherr( "gdtrc", DOMAIN ); - return( 0.0 ); - } -return( igamc( b, a * x ) ); -} -- cgit v1.2.3