program main

c*********************************************************************72
c
cc ellipsoid_grid_test() tests ellipsoid_grid().
c
c  Licensing:
c
c    This code is distributed under the MIT 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)' ) 'ellipsoid_grid_test():'
      write ( *, '(a)' ) '  FORTRAN77 version'
      write ( *, '(a)' ) '  Test ellipsoid_grid().'

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

      return
      end
      subroutine ellipsoid_grid_test01 ( )

c*********************************************************************72
c
cc ELLIPSOID_GRID_TEST01 tests ELLIPSOID_GRID.
c
c  Licensing:
c
c    This code is distributed under the MIT license.
c
c  Modified:
c
c    11 November 2011
c
c  Author:
c
c    John Burkardt
c
      implicit none

      double precision c(3)
      character * ( 80 ) filename
      integer n
      integer ng
      double precision r(3)
      double precision xyz(3,1200)

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

      n = 4
      r(1) = 2.0D+00
      r(2) = 1.0D+00
      r(3) = 1.5D+00
      c(1) = 1.0D+00
      c(2) = 2.0D+00
      c(3) = 1.5D+00

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

      call ellipsoid_grid_count ( n, r, c, ng )

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

      call ellipsoid_grid ( n, r, c, ng, xyz )

      call r83vec_print_part ( ng, xyz, 20, 
     &  '  Part of the grid point array:' )

      filename = 'ellipsoid_grid_test01.xyz'

      call r8mat_write ( filename, 3, ng, xyz )

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

      return
      end