program main

c*********************************************************************72
c
cc MAIN is the main program for ELLIPSE_GRID_PRB.
c
c  Discussion:
c
c    ELLIPSE_GRID_PRB tests the ELLIPSE_GRID lilbrary.
c
c  Licensing:
c
c    This code is distributed under the GNU LGPL license.
c
c  Modified:
c
c    11 November 2011
c
c  Author:
c
c    John Burkardt
c
      implicit none

      call timestamp ( )
      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'ELLIPSE_GRID_TEST:'
      write ( *, '(a)' ) '  FORTRAN77 version'
      write ( *, '(a)' ) '  Test the ELLIPSE_GRID library.'

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

      return
      end
      subroutine ellipse_grid_test01 ( )

c*********************************************************************72
c
cc ELLIPSE_GRID_TEST01 tests ELLIPSE_GRID.
c
c  Licensing:
c
c    This code is distributed under the GNU LGPL license.
c
c  Modified:
c
c    11 November 2011
c
c  Author:
c
c    John Burkardt
c
      implicit none

      double precision c(2)
      character * ( 80 ) filename
      integer n
      integer ng
      double precision r(2)
      double precision xy(2,500)

      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 'TEST01:'
      write ( *, '(a)' ) '  ELLIPSE_GRID can define a grid of points'
      write ( *, '(a)' ) '  with N+1 points on the minor half axis,'
      write ( *, '(a)' ) '  based on any ellipse.'

      n = 8
      r(1) = 2.0D+00
      r(2) = 1.0D+00
      c(1) = 1.0D+00
      c(2) = 2.0D+00

      write ( *, '(a)' ) ' '
      write ( *, '(a,i8)' ) '  We use N = ', n
      write ( *, '(a,g14.6,a,g14.6,a)' ) 
     &  '  Radius R = (', r(1), ',', r(2), ')'
      write ( *, '(a,g14.6,a,g14.6,a)' ) 
     &  '  Center C = (', c(1), ',', c(2), ')'

      call ellipse_grid_count ( n, r, c, ng )

      write ( *, '(a)' ) ' '
      write ( *, '(a,i8)' ) '  Number of grid points will be ', ng

      call ellipse_grid ( n, r, c, ng, xy )

      call r82vec_print_part ( ng, xy, 20, 
     &  '  Part of the grid point array:' )

      filename = 'ellipse_grid_test01.xy'

      call r8mat_write ( filename, 2, ng, xy )

      write ( *, '(a)' ) ' '
      write ( *, '(a)' ) 
     &  '  Data written to the file "' // trim ( filename ) // '".'

      return
      end