#include double gauss(double x){ double val,sigma=1; val = 1./(sqrt(2*M_PI)*sigma)*exp(-pow(x/sigma,2)/2); return val; }