typedef double real;
#include <math.h>
extern double  probab(const int& ndof, const double& chisq);