c
c
c
      function rlgdre(n,x)
C     CALCULATES LEGENDRE POLYNOMIALS USING RECURRENCE  ABRAM./STEGUN
C     PAGE 782
      implicit real*8 (a-h,o-z)
      rlgdre=1
      r2=rlgdre
      if(n.eq.0)return
      rlgdre=x
      r1=rlgdre
      if(n.eq.1)return
      do 1 i=2,n
      rxr1=x*r1
      rlgdre=rxr1+rxr1-r2+(r2-rxr1)/i
      r2=r1
 1    r1=rlgdre
      return
      end