program HermiteExample
!
!       ===========================================================
!       Purpose : This program computes the zeros of Hermite 
!                 polynomial Ln(x) in the interval [-ì,ì] and the
!                 corresponding weighting coefficients for Gauss-
!                 Hermite integration using subroutine HERZO
!       Input :   n    --- Order of the Hermite polynomial
!                 X(n) --- Zeros of the Hermite polynomial
!                 W(n) --- Corresponding weighting coefficients
!       ===========================================================
!
            use hermite
            implicit none
            integer, parameter:: dble = kind(1d0)
            
            integer, parameter:: N = 5
            integer j
            real(dble) X(N), W(N)
            
            call HERZO(N,X,W)
            print *, "  Nodes and weights for Gauss-Hermite integration"
            print *, " " 
            print *, "  i             xi                      Wi"
            print *, " -----------------------------------------------------"
            do j=1,N
              print *, j,X(j),W(j)
            end do
 
        end program HermiteExample