SUBROUTINE Laguer_fun ( y, n, h, m )
USE Numeric_Kinds_Module
USE FileUnits_Module

!
!--------------------- Laguer ------------------------------------------
!
! Laguer -- computes Laguere polynomials
!
! --- on input the user supplies --
!
!     y  -- the dimensionless argument of the Laguere polynomial
!     n  -- the order of the polynomial (n=0 ==> h0=1)
!     m  -- n+1 (required), also the length of the space in h
!
! --- on exit the routine returns --
!
!     h  -- vector of length m, containing h(1)=h0 h(2)=h1, etc.
!     Laguer -- function value returned is h(m) = nth Laguere poly.
!
!-----------------------------------------------------------------------
!
INTEGER(KIND=IW_Kind) :: n, m, i
REAL(KIND=WP_Kind)    :: h(0:n), y
!
!-------------------- execution begins here ----------------------------
!
h(0) = One
IF( n .eq. 0 ) RETURN
!
h(1) = One - y
IF( n .eq. 1 ) RETURN
!
DO i  = 1, n-1
   h(i+1) = ((2*i+1-y)*h(i) - i*h(i-1))/(i+1)
ENDDO

RETURN
END