PDF/spdf/wate32.F

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