program DRDPOI
c>> 2001-05-25 DRDPOI Krogh Minor change for making .f90 version.
c>> 1996-05-28 DRDPOI Krogh  Changed Fortran 90 code.
c>> 1994-10-19 DRDPOI Krogh  Changes to use M77CON
c>> 1994-07-06 DRDPOI WVS set up for chgtyp
c
c     Evaluate the Probability Function Q(n,lamda) of the Poisson
c     distribution by using DCDPOI.
c
c--D replaces "?": DR?POI, ?cdpoi, ?gamik
      double precision LAMDAS(4), P, Q(4)
      integer N, IERR, J
      data LAMDAS /0.5d0, 1.0d0, 1.5d0, 2.0d0/
c
10    format ('     Poisson Probability Function Q(lamda | n)'/
     1        '  N    LAMDA =',4f14.7/)
30    format (i3,15x,1p,4g14.7)
      call dgamik (0.0d0, 0.0d0, 0.0d0, 0)
      print 10, lamdas
      do 50 n = 1, 10
         do 40 j = 1, 4
            call dcdpoi (n, lamdas(j), p, q(j), ierr)
40       continue
         print 30, n, q
50    continue
      stop
      end