1 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   2 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   3 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   4 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   5 | c
   6 | c.......................................................................
   7 | c
   8 | c  cumulative of the standard normal
   9 | c
  10 | c  in :  * $x$                   ---> xstcd
  11 | c
  12 | c  out : * $\Gamma (x)$
  13 | c                                ---> cdstcd
  14 | c
  15 | c.......................................................................
  16 |       subroutine stcd(xstcd,cdstcd)
  17 | c.......................................................................
  18 |       implicit double precision (a-h,o-z)
  19 | c.......................................................................
  20 |       real stcd2, erfc
  21 | c.......................................................................
  22 |       stcd1=-xstcd/dsqrt(2.d0)
  23 |       stcd2=real(stcd1)
  24 |       stcd2=erfc(stcd2)
  25 |       cdstcd=dble(stcd2)/2.d0
  26 | c.......................................................................
  27 |       return
  28 |       end