PDF/spdf/wate32.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:31:58  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:30:07  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13       SUBROUTINE WATE32
  14 C  32 POINT GAUSSIAN QUADRATURE ROUTINE
  15 #include "pdf/expdp.inc"
  16      +       X(16),W(16)
  17 #include "pdf/expdp.inc"
  18      +               XI(32),WI(32),XX(33)
  19       COMMON/W5051I9/XI,WI,XX,NTERMS
  20       NTERMS=32
  21       X(1)=0.048307665687738316235D0
  22       X(2)=0.144471961582796493485D0
  23       X(3)=0.239287362252137074545D0
  24       X(4)=0.331868602282127649780D0
  25       X(5)=0.421351276130635345364D0
  26       X(6)=0.506899908932229390024D0
  27       X(7)=0.587715757240762329041D0
  28       X(8)=0.663044266930215200975D0
  29       X(9)=0.732182118740289680387D0
  30       X(10)=0.794483795967942406963D0
  31       X(11)=0.849367613732569970134D0
  32       X(12)=0.896321155766052123965D0
  33       X(13)=0.934906075937739689171D0
  34       X(14)=0.964762255587506430774D0
  35       X(15)=0.985611511545268335400D0
  36       X(16)=0.997263861849481563545D0
  37       W(1)=0.096540088514727800567D0
  38       W(2)=0.095638720079274859419D0
  39       W(3)=0.093844399080804565639D0
  40       W(4)=0.091173878695763884713D0
  41       W(5)=0.087652093004403811143D0
  42       W(6)=0.083311924226946755222D0
  43       W(7)=0.078193895787070306472D0
  44       W(8)=0.072345794108848506225D0
  45       W(9)=0.065822222776361846838D0
  46       W(10)=0.058684093478535547145D0
  47       W(11)=0.050998059262376176196D0
  48       W(12)=0.042835898022226680657D0
  49       W(13)=0.034273862913021433103D0
  50       W(14)=0.025392065309262059456D0
  51       W(15)=0.016274394730905670605D0
  52       W(16)=0.007018610009470096600D0
  53       NTERMH = NTERMS/2
  54       DO 1 I=1,NTERMH
  55       XI(I)=-X(17-I)
  56       WI(I)=W(17-I)
  57       XI(I+16)=X(I)
  58       WI(I+16)=W(I)
  59     1 CONTINUE
  60       DO 2 I=1,NTERMS
  61     2 XX(I)=0.5D0*(XI(I)+1.0D0)
  62       XX(33)=1.0D0
  63       RETURN
  64       END