%N legdr %F $MIR/src/subs/lspoly.for
%D Legendre polynomial utility
%P mchw
%: math, polynomials
%B
      subroutine legdr(x,p,np)

      implicit none
      integer np
      real x,p(np)

  Will put values of first NP Legendre polynomials of X in array P.
  Max NP is 9 (i.e., up to 8th order). P(N) is (N-1)th order
  BUG:  p(np) dangerous if compiled with run time array checking (-C)
