#include <math.h>

double phi(double x)
{
   double pi=3.14159265;
   double p=exp(-0.5*x*x)/sqrt(2*pi);
   return p;
}