double gaussian_antideriv ( double mu, double sigma, double x );
double gaussian_value ( int n, double mu, double sigma, double x );
double hermite_polynomial_value ( int n, double x );
double r8_erf ( double x );