#ifndef _utl_NormalDistribution_h_
#define _utl_NormalDistribution_h_

// $Id$

#include <utl/Math.h>

namespace utl {

  typedef unsigned long ulong;

  inline double LogPoissonPDF(const ulong k, const double x)
  { return -x - LogGamma(k+1) + k * log(x); }

  inline double PoissonPDF(const ulong k, const double x)
  { return exp(LogPoissonPDF(k,x)); }

}

#endif

// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// End: