double f(double x) { return exp(-pow(x,2.0)); }