# gaussian(x,mean,sigma) is Gaussian function;

gaussian(0,0,1):= 1/sqrt(2)/sqrt(pi);