Go to the documentation of this file.00001 #ifndef mbl_gamma_h
00002 #define mbl_gamma_h
00003
00004
00005
00006
00007
00008
00009
00010 double mbl_log_gamma(double a);
00011
00012
00013 double mbl_gamma_p(double a, double x);
00014
00015
00016 double mbl_gamma_q(double a, double x);
00017
00018
00019 double mbl_log_gamma_q(double a, double x);
00020
00021
00022
00023 inline double mbl_erf(double x)
00024 { return (x<0)?-mbl_gamma_p(0.5,x*x):mbl_gamma_p(0.5,x*x); }
00025
00026 #endif // mbl_gamma_h