double pn( double _n)
{
/* =~ -log P(n) for gaussian n */

  return _n*_n;
}
