Документ взят из кэша поисковой машины. Адрес оригинального документа : http://xmm.vilspa.esa.es/sas/7.0.0/doc/ssclib/node69.html
Дата изменения: Wed Jun 28 23:12:27 2006
Дата индексирования: Sat Dec 22 08:35:04 2007
Кодировка:

XMM-Newton Science Analysis System


ssclib (ssclib-3.29) [xmmsas_20060628_1801-7.0.0]

Meta Index / Home Page / Poissonian statistics and source

Poisson probability

This returns the probability $P(i)$ of the occurrence of an integer $i$ according to the Poisson distribution

$$P(i) = \frac{a^i exp(-a)}{i!}$$

where $a$ is the average or expectation value of $i$.

Note that the argument may also be a real number. In this case what is returned is

$$p(r) = \frac{a^r exp(-a)}{\Gamma(r+1)}$$

The value $p$ is not quite a probability density: it would need to be normalized by

$$\int_{0}^{\inf} dr \; \frac{a^r exp(-a)}{\Gamma(r+1)}$$

in order for this to be true. However it does have the property that, if $r=i$, $p(r)=P(i)$.

  interface poissonProb
    real function poissonProbSingle(av, realI)
      real(single), intent(in) :: av, realI
    end function poissonProbSingle

    real function poissonProbInt32(av, i)
      real(single), intent(in) :: av
      integer(int32), intent(in) :: i
    end function poissonProbInt32
  end interface