#include <math.h>
#define PI 3.14159265


double log2, lgpi, hlg2pi, pi, pii;
void init()
{
  log2=log(2);
  pi = PI;
  pii = 1/ pi;
  lgpi=log(pi);
  hlg2pi=(log2+lgpi)/2;
}