program main

!*****************************************************************************80
!
!! legendre_polynomial_test() tests legendre_shifted_polynomial().
!
!  Licensing:
!
!    This code is distributed under the MIT license.
!
!  Modified:
!
!    14 March 2016
!
!  Author:
!
!    John Burkardt
!
  implicit none

  call timestamp ( )
  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'legendre_shifted_polynomial_test():'
  write ( *, '(a)' ) '  FORTRAN90 version.'
  write ( *, '(a)' ) '  Test legendre_shifted_polynomial().'

  call p01_polynomial_value_test ( )
!
!  Terminate.
!
  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'LEGENDRE_SHIFTED_POLYNOMIAL_TEST():'
  write ( *, '(a)' ) '  Normal end of execution.'
  write ( *, '(a)' ) ' '
  call timestamp ( )

  stop 0
end
subroutine p01_polynomial_value_test ( )

!*****************************************************************************80
!
!! P01_POLYNOMIAL_VALUE_TEST tests P01_POLYNOMIAL_VALUE.
!
!  Licensing:
!
!    This code is distributed under the MIT license.
!
!  Modified:
!
!    15 March 2016
!
!  Author:
!
!    John Burkardt
!
  implicit none

  integer, parameter :: rk = kind ( 1.0D+00 )

  integer n_data
  real ( kind = rk ) e
  real ( kind = rk ) fx1
  real ( kind = rk ) fx2
  integer, parameter :: m = 1
  integer n
  real ( kind = rk ), allocatable :: v(:,:)
  real ( kind = rk ) x
  real ( kind = rk ) x_vec(1)

  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) 'P01_POLYNOMIAL_VALUE_TEST:'
  write ( *, '(a)' ) '  P01_POLYNOMIAL_VALUE evaluates the shifted Legendre polynomial P01(n,x).'
  write ( *, '(a)' ) ' '
  write ( *, '(a)' ) '                        Tabulated                 Computed'
  write ( *, '(a)' ) '     N        X         P01(N,X)                  P01(N,X)                     Error'
  write ( *, '(a)' ) ' '

  n_data = 0

  do

    call p01_polynomial_values ( n_data, n, x, fx1 )

    x_vec(1) = x

    if ( n_data == 0 ) then
      exit
    end if

    allocate ( v(m,n+1) )
    call p01_polynomial_value ( m, n, x_vec, v )
    fx2 = v(1,n+1)
    deallocate ( v )

    e = fx1 - fx2

    write ( *, '(2x,i4,2x,f12.6,2x,g24.16,2x,g24.16,2x,g8.2)' ) &
      n, x, fx1, fx2, e

  end do

  return
end