ADAS Subroutine a8gcf

       subroutine a8gcf( gammcf, a    , x    , gln )
c---------------------------------------------------------------------
c
c  ********************* fortran77 subroutine a8gcf ******************
c
c  purpose: to evaluate the continued fraction expansion for the 
c           incomplete gamma function gamma(a)*q(a,x) - based on 
c           numerical recipes  
c
c  calling program:  a8gamg.for
c
c  input:
c            (r*8)  a      = parameter of q(a,x)
c            (r*8)  x      = paramete of q(a,x)
c  output:
c            (r*8)  gammcf = incomplete gamma function gamma(a)*q(a,x)
c            (r*8)  gln    = ln(gamma(a))
c
c
c  routines: 
c             none
c             a8gaml    adas    obtains log(gamma(a))
c
c  author:   Hugh Summers, University of Strathclyde ext.4196
c
c
c  version 1.1                             date:     25/06/99
c  modified: Hugh Summers
c		- first release
c
c-----------------------------------------------------------------------
      REAL*8              A,           GAMMCF,      GLN,         X