* * $Id$ * * $Log$ * Revision 1.4 2000/09/18 10:02:37 hristov * Makefile added to PDF8 * * Revision 1.1.1.2 1996/10/30 08:31:58 cernlib * Version 7.04 * * Revision 1.1.1.1 1996/04/12 15:30:07 plothow * Version 7.01 * * #include "pdf/pilot.h" SUBROUTINE WATE32 C 32 POINT GAUSSIAN QUADRATURE ROUTINE #if defined(CERNLIB_DOUBLE) DOUBLE PRECISION #endif #if defined(CERNLIB_SINGLE) REAL #endif + X(16),W(16) #if defined(CERNLIB_DOUBLE) DOUBLE PRECISION #endif #if defined(CERNLIB_SINGLE) REAL #endif + XI(32),WI(32),XX(33) COMMON/W5051I9/XI,WI,XX,NTERMS NTERMS=32 X(1)=0.048307665687738316235D0 X(2)=0.144471961582796493485D0 X(3)=0.239287362252137074545D0 X(4)=0.331868602282127649780D0 X(5)=0.421351276130635345364D0 X(6)=0.506899908932229390024D0 X(7)=0.587715757240762329041D0 X(8)=0.663044266930215200975D0 X(9)=0.732182118740289680387D0 X(10)=0.794483795967942406963D0 X(11)=0.849367613732569970134D0 X(12)=0.896321155766052123965D0 X(13)=0.934906075937739689171D0 X(14)=0.964762255587506430774D0 X(15)=0.985611511545268335400D0 X(16)=0.997263861849481563545D0 W(1)=0.096540088514727800567D0 W(2)=0.095638720079274859419D0 W(3)=0.093844399080804565639D0 W(4)=0.091173878695763884713D0 W(5)=0.087652093004403811143D0 W(6)=0.083311924226946755222D0 W(7)=0.078193895787070306472D0 W(8)=0.072345794108848506225D0 W(9)=0.065822222776361846838D0 W(10)=0.058684093478535547145D0 W(11)=0.050998059262376176196D0 W(12)=0.042835898022226680657D0 W(13)=0.034273862913021433103D0 W(14)=0.025392065309262059456D0 W(15)=0.016274394730905670605D0 W(16)=0.007018610009470096600D0 NTERMH = NTERMS/2 DO 1 I=1,NTERMH XI(I)=-X(17-I) WI(I)=W(17-I) XI(I+16)=X(I) WI(I+16)=W(I) 1 CONTINUE DO 2 I=1,NTERMS 2 XX(I)=0.5D0*(XI(I)+1.0D0) XX(33)=1.0D0 RETURN END