%N gammq %F $MIR/src/subs/gamma.for
%D Find incomplete Gamma function Q(a,x) = 1 - P(a,x)
%P nebk
%: gamma function, utilities
%B
      real function gammq (a, x)

      implicit none
      real a, x

 Find the incomplete Gamma function Q(a,x) = 1 - P(a,x)
 Algorithm from Numerical Recipes p 162

  Input
    a,x     r     Arguments
