program main

c*********************************************************************72
c
cc ASA032_PRB tests ASA032.
c
c  Modified:
c
c    04 January 2008
c
c  Author:
c
c    John Burkardt
c
      implicit none

      call timestamp ( )
      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'ASA032_PRB:'
      write ( *, '(a)' ) '  FORTRAN77 version'
      write ( *, '(a,i1,a,i1,a,i1)' ) '  GCC version ',
     &   __GNUC__ , '.',
     &   __GNUC_MINOR__, '.',
     &   __GNUC_PATCHLEVEL__ 
      write ( *, '(a)' ) '  Test the ASA032 library.'

      call test01 ( )
c
c  Terminate.
c
      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'ASA032_PRB:'
      write ( *, '(a)' ) '  Normal end of execution.'
      write ( *, '(a)' ) ' '
      call timestamp ( )

      stop
      end
      subroutine test01 ( )

c*********************************************************************72
c
cc TEST01 demonstrates the use of GAMAIN.
c
c  Modified:
c
c    17 January 2006
c
c  Author:
c
c    John Burkardt
c
      implicit none

      double precision a
      double precision fx
      double precision fx2
      double precision gamain
      integer ifault
      integer n_data
      double precision x

      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'TEST01:'
      write ( *, '(a)' ) '  GAMAIN computes the incomplete Gamma '
      write ( *, '(a)' ) '  function.  Compare to tabulated values.'
      write ( *, '(a)' ) ' '
      write ( *, '(a,a)' )
     &  '       A            X          ',
     &  '   FX                        FX                      DIFF'
      write ( *, '(a,a)' )
     &  '                               ',
     &  ' (tabulated)                (GAMAIN)'
      write ( *, '(a)' ) ' '

      n_data = 0

10    continue

        call gamma_inc_values ( n_data, a, x, fx )

        if ( n_data .eq. 0 ) then
          go to 20
        end if

        fx2 = gamain ( x, a, ifault )

        write ( *,
     &  '(2x,f12.8,2x,f12.8,2x,g24.16,2x,g24.16,2x,g10.4)' )
     &  a, x, fx, fx2, dabs ( fx - fx2 )

      go to 10

20    continue

      return
      end