Esope source files

   1 : C LEGEND    SOURCE    CHAT      05/01/13    01:14:36     5004           
   2 :       SUBROUTINE legend(n,x,w)
   3 : C     ------------------------------------------------------------------
   4 : C     This routine generates the Gauss-Lobatto normalized points for the 
   5 : C     Legendre polynomials in the monodimensional interval (A,B)=(-1,+1)
   6 : c     and the correspondent weights
   7 : C     A,B: end points of the interval                                   
   8 : C     X: vector of Gauss-Lobatto nodes (dimension NP)                  
   9 : C     W: vector of Gauss-Lobatto weights (dimension NP)                  
  10 : C     N: The degree of polynomial
  11 : C     (NP: number of points on the interval  = N + 1)
  12 : C     ------------------------------------------------------------------
  13 :       IMPLICIT INTEGER(I-N)
  14 :       IMPLICIT REAL*8 (A-H,O-Z)
  15 : c
  16 :       DIMENSION X(*),W(*)
  17 : c
  18 :       NP = N + 1
  19 :       DO 2 I=1,NP                                               
  20 :       W(I)=0.D0                                                
  21 : 2     X(I)=0.D0                                               
  21 :       CALL ZWGLJD(X,W,NP,0.D0,0.D0)                          
  22 : c
  23 :       RETURN                               
  24 :       END                                 
  25 :  
  26 :