00001 00008 #include <math.h> 00009 00011 #define ITMAX 100000 00013 #define EPS 3.0e-7 00015 #define FPMIN 1.0e-30 00016 00017 //void nrerror(const char *msg); 00018 00019 double gammln(double xx); 00020 00021 double gammp(double a, double x); 00022 00023 double gammq(double a, double x); 00024 00025 void gser(double *gamser, double a, double x, double *gln); 00026 00027 void gcf(double *gammcf, double a, double x, double *gln); 00028 00029 double erf(double x);