program main

c*********************************************************************72
c
cc MAIN is the main program for TOMS351_PRB.
c
c  Discussion:
c
c    TOMS351_PRB tests the TOMS351 library.
c
c  Modified:
c
c    05 January 2006
c
c  Author:
c
c    John Burkardt
c
      implicit none

      call timestamp ( )
      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'TOMS351_PRB'
      write ( *, '(a)' ) '  FORTRAN77 version'
      write ( *, '(a)' ) '  Test the TOMS351 library.'

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

      stop
      end
      subroutine test01 ( )

c*********************************************************************72
c
cc TEST01 tests ROMINT.
c
c  Modified:
c
c    05 January 2006
c
c  Author:
c
c    John Burkardt
c
      implicit none

      real a
      real b
      real eps
      real err
      external f01
      integer maxe
      integer maxe_in
      integer n
      real pi
      real val

      eps = 0.000001E+00
      pi = 3.14159265E+00

      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'TEST01'
      write ( *, '(a)' ) '  Test ROMINT, for Romberg integration.'
      write ( *, '(a)' ) ' '
      write ( *, '(a,g14.6)' ) '  Error tolerance EPS = ', eps
      write ( *, '(a)' ) ' '

      write ( *, '(a)' )
     &  '      A         B   MAXE  MAXE     N      ERR             VAL'
      write ( *, '(a)' )
     &  '                    in    out'
      write ( *, '(a)' ) ' '

      a = 0.0E+00
      b = pi / 2.0E+00

      do 10 maxe_in = 1, 5

        maxe = maxe_in

        call romint ( val, err, eps, a, b, n, maxe, f01 )

        write ( *,
     &    '(2x,f8.4,2x,f8.4,2x,i2,4x,i2,2x,i4,2x,g14.6,2x,g14.6)' )
     &    a, b, maxe_in, maxe, n, err, val

10    continue

      return
      end
      subroutine test02 ( )

c*********************************************************************72
c
cc TEST02 tests ROMINT.
c
c  Modified:
c
c    05 January 2006
c
c  Author:
c
c    John Burkardt
c
      implicit none

      real a
      real b
      real eps
      real err
      external f02
      integer maxe
      integer maxe_in
      integer n
      real val

      eps = 0.000001E+00

      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'TEST02'
      write ( *, '(a)' ) '  Test ROMINT, for Romberg integration.'
      write ( *, '(a)' ) ' '
      write ( *, '(a,g14.6)' ) '  Error tolerance EPS = ', eps
      write ( *, '(a)' ) ' '

      write ( *, '(a)' )
     &  '      A         B   MAXE  MAXE     N      ERR             VAL'
      write ( *, '(a)' )
     &  '                    in    out'
      write ( *, '(a)' ) ' '

      a = 0.0E+00
      b = 4.3E+00

      do 10 maxe_in = 1, 5

        maxe = maxe_in

        call romint ( val, err, eps, a, b, n, maxe, f02 )

        write ( *,
     &    '(2x,f8.4,2x,f8.4,2x,i2,4x,i2,2x,i4,2x,g14.6,2x,g14.6)' )
     &    a, b, maxe_in, maxe, n, err, val

10    continue

      return
      end
      function f01 ( x )

c*********************************************************************72
c
cc F01 evaluates test function 1.
c
      implicit none

      real f01
      real x

      f01 = cos ( x )

      return
      end
      function f02 ( x )

c*********************************************************************72
c
cc F02 evaluates test function 2.
c
      implicit none

      real f02
      real x

      f02 = exp ( - x * x )

      return
      end
      subroutine timestamp ( )

c*********************************************************************72
c
cc TIMESTAMP prints out the current YMDHMS date as a timestamp.
c
c  Discussion:
c
c    This FORTRAN77 version is made available for cases where the
c    FORTRAN90 version cannot be used.
c
c  Modified:
c
c    16 September 2005
c
c  Author:
c
c    John Burkardt
c
c  Parameters:
c
c    None
c
      implicit none

      character ( len = 8 ) date
      character ( len = 10 ) time

      call date_and_time ( date, time )

      write ( *, '(a8,2x,a10)' ) date, time

      return
      end