(42S) GAUSS LEGENDRE 16 POINT INTEGRATION

[Gerald H]

"Func name" should be the name of a programme taking one argument from the X register & returning a real value to the X register.

Func name → 00
Lower lim → 01
Upper lim → 02
# division → 03

Code:

0.    { 379-Byte Prgm }
1.    LBL “GL16”
2.    CLX
3.    STO 04
4.    RCL 02
5.    RCL 01
6.    STO 07
7.    –
8.    RCL 03
9.    STO 05
10.    STO+ ST X
11.    /
12.    STO 06
13.    LBL 17
14.    STO+ 07
15.    16
16.    STO 08
17.    LBL 18
18.    RCL 07
19.    XEQ IND 08
20.    RCL* 06
21.    STO 09
22.    –
23.    XEQ IND 00
24.    STO 10
25.    RCL 07
26.    RCL+ 09
27.    XEQ IND 00
28.    RCL+ 10
29.    DSE 08
30.    XEQ IND 08
31.    *
32.    STO+ 04
33.    DSE 08
34.    GTO 18
35.    RCL 06
36.    STO+ 07
37.    DSE 05
38.    GTO 17
39.    STO* 04
40.    RCL 04
41.    RTN
42.    LBL 01
43.    271524594118E-13
44.    RTN
45.    LBL 02 
46.    989400934992E-12
47.    RTN
48.    LBL 03
49.    622535239386E-13
50.    RTN
51.    LBL 04
52.    944575023073E-12
53.    RTN
54.    LBL 05
55.    951585116825E-13
56.    RTN
57.    LBL 06
58.    865631202388E-12
59.    RTN
60.    LBL 07
61.    124628971256E-12
62.    RTN
63.    LBL 08
64.    755404408355E-12
65.    RTN
66.    LBL 09
67.    149595988817E-12
68.    RTN
69.    LBL 10
70.    617876244403E-12
71.    RTN
72.    LBL 11
73.    169156519395E-12
74.    RTN
75.    LBL 12
76.    458016777657E-12
77.    RTN
78.    LBL 13
79.    182603415045E-12
80.    RTN
81.    LBL 14
82.    281603550779E-12
83.    RTN
84.    LBL 15
85.    189450610455E-12
86.    RTN
87.    LBL 16
88.    950125098376E-13
89.    END