]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 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 |