Removed 'Written by' line
[pspp-builds.git] / lib / gsl-extras / hypergeometric.c
index 2e9cf69267db0f2354e0320daf44f42afd74ea04..67b757b2f93d6a0389d13b004c0df4b9dbc7a31e 100644 (file)
@@ -1,7 +1,7 @@
 /* cdf/hypergeometric.c
  *
  * Copyright (C) 2004 Free Software Foundation, Inc.
- * Written by Jason H. Stover.
+ * 
  *
  * This program is free software; you can redistribute it and/or modify
  * it under the terms of the GNU General Public License as published by
@@ -15,7 +15,7 @@
  *
  * You should have received a copy of the GNU General Public License
  * along with this program; if not, write to the Free Software
- * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307, USA.
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
  */
 
 /*
@@ -29,7 +29,7 @@
  *
  * References:
  *
- * T. Wu. An accurate computation of the hypergeometric distribution 
+ * T. Wu. An accurate computation of the hypergeometric distribution
  * function. ACM Transactions on Mathematical Software. Volume 19, number 1,
  * March 1993.
  *  This algorithm is not used, since it requires factoring the
@@ -38,9 +38,8 @@
  *  time than the algorithm used here.
  *
  * W. Feller. An Introduction to Probability Theory and Its Applications,
- * third edition. 1968. Chapter 2, section 6. 
+ * third edition. 1968. Chapter 2, section 6.
  */
-#include <config.h>
 #include <math.h>
 #include <gsl/gsl_math.h>
 #include <gsl/gsl_errno.h>
@@ -52,7 +51,7 @@
  * Pr (X <= k)
  */
 double
-gslextras_cdf_hypergeometric_P (const unsigned int k, 
+gslextras_cdf_hypergeometric_P (const unsigned int k,
                                 const unsigned int n0,
                                 const unsigned int n1,
                                 const unsigned int t)
@@ -126,7 +125,7 @@ gslextras_cdf_hypergeometric_P (const unsigned int k,
  * Pr (X > k)
  */
 double
-gslextras_cdf_hypergeometric_Q (const unsigned int k, 
+gslextras_cdf_hypergeometric_Q (const unsigned int k,
                                 const unsigned int n0,
                                 const unsigned int n1,
                                 const unsigned int t)
@@ -154,7 +153,7 @@ gslextras_cdf_hypergeometric_Q (const unsigned int k,
       P = 0.0;
       mode = (int) t*n0 / (n0+n1);
       relerr = 1.0;
-      
+
       if(k < mode)
        {
          i = mode;