xref: /phasta/shapeFunction/src/Fn.c (revision 16223cb9c3f88b34f2cb94151b5cf5ffc1aac5e2) 
1 /*-------------------------------------------------------------------------
2   Scientific Computation Research Center, RPI, Troy NY
3   (C) Copyright 1995, RPI-SCOREC
4 
5   Project  : shapeFuntions
6   Author(s): Saikat Dey
7   Creation : Oct., 95
8   Modifi.  :
9   Function :
10              mode shape for a triangular face. i,j are highest order of
11              r,s in the polynomial, maple generated code.
12 -------------------------------------------------------------------------*/
13 
14 #ifdef __cplusplus
15 extern "C" {
16 #endif
17 
18 double Fn(int i, int j, double r, double s) {
19    double f = 0.0;
20    double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18;
21    double t19,t20,t21,t22,t23,t24,t25,t26,t27,t28,t29,t30,t31,t32,t33,t34;
22    double t35,t36,t37,t38,t39,t40,t41,t42,t43,t44,t45,t47,t50,t51;
23 
24    /* p=3 0 0 */
25    if( i==0 && j==0 ) {
26       f = 1.0;
27    /* p=4 0 1 */
28    } else if( i==0 && j==1 ) {
29       f = s-1.0/3.0;
30    /* p=4 1 0 */
31    } else if( i==1 && j==0 ) {
32       f = r-1.0/3.0;
33    /* p=5 0 2 */
34    } else if( i==0 && j==2 ) {
35       t1 = s*s;
36       f = t1-3.0/4.0*s+3.0/28.0;
37    /* p=5 1 1 */
38    } else if( i==1 && j==1 ) {
39       f = r*s-r/4-s/4+1.0/14.0;
40    /* p=5 2 0 */
41    } else if( i==2 && j==0 ) {
42       t1 = r*r;
43       f = t1-3.0/4.0*r+3.0/28.0;
44    /* p=6 0 3 */
45    } else if( i==0 && j==3 ) {
46       t1 = s*s;
47       f = t1*s-6.0/5.0*t1+2.0/5.0*s-1.0/30.0;
48    /* p=6 1 2 */
49    } else if( i==1 && j==2 ) {
50       t1 = s*s;
51       f = r*t1-3.0/5.0*r*s+r/15-t1/5+2.0/15.0*s-1.0/60.0;
52    /* p=6 2 1 */
53    } else if( i==2 && j==1 ) {
54       t1 = r*r;
55       f = t1*s-t1/5-3.0/5.0*r*s+2.0/15.0*r+s/15-1.0/60.0;
56    /* p=6 3 0 */
57    } else if( i==3 && j==0 ) {
58       t1 = r*r;
59       f = r*t1-6.0/5.0*t1+2.0/5.0*r-1.0/30.0;
60    /* p=7 0 4 */
61    } else if( i==0 && j==4 ) {
62       t1 = s*s;
63       t2 = t1*t1;
64       f = t2-5.0/3.0*t1*s+10.0/11.0*t1-2.0/11.0*s+1.0/99.0;
65    /* p=7 1 3 */
66    } else if( i==1 && j==3 ) {
67       t1 = s*s;
68       t2 = t1*s;
69       f = r*t2-r*t1+3.0/11.0*r*s-r/55-t2/6+2.0/11.0*t1-3.0/55.0*s+2.0/495.0;
70    /* p=7 2 2 */
71    } else if( i==2 && j==2 ) {
72       t1 = r*r;
73       t2 = s*s;
74       f = t1*t2-t1*s/2+t1/22-r*t2/2+3.0/11.0*r*s-3.0/110.0*r+t2/22-3.0/110.0*s
75 +1.0/330.0;
76    /* p=7 3 1 */
77    } else if( i==3 && j==1 ) {
78       t1 = r*r;
79       t2 = r*t1;
80       f = t2*s-t2/6-t1*s+2.0/11.0*t1+3.0/11.0*r*s-3.0/55.0*r-s/55+2.0/495.0;
81    /* p=7 4 0 */
82    } else if( i==4 && j==0 ) {
83       t1 = r*r;
84       t2 = t1*t1;
85       f = t2-5.0/3.0*r*t1+10.0/11.0*t1-2.0/11.0*r+1.0/99.0;
86    /* p=8 0 5 */
87    } else if( i==0 && j==5 ) {
88       t1 = s*s;
89       t2 = t1*t1;
90       f = t2*s-15.0/7.0*t2+150.0/91.0*t1*s-50.0/91.0*t1+75.0/1001.0*s-3.0/
91 1001.0;
92    /* p=8 1 4 */
93    } else if( i==1 && j==4 ) {
94       t1 = s*s;
95       t2 = t1*t1;
96       t4 = t1*s;
97       f = r*t2-10.0/7.0*r*t4+60.0/91.0*r*t1-10.0/91.0*r*s+5.0/1001.0*r-t2/7+
98 20.0/91.0*t4-10.0/91.0*t1+20.0/1001.0*s-1.0/1001.0;
99    /* p=8 2 3 */
100    } else if( i==2 && j==3 ) {
101       t1 = r*r;
102       t2 = s*s;
103       t3 = t2*s;
104       f = t1*t3-6.0/7.0*t1*t2+18.0/91.0*t1*s-t1/91-3.0/7.0*r*t3+36.0/91.0*r*
105 t2-9.0/91.0*r*s+6.0/1001.0*r+3.0/91.0*t3-3.0/91.0*t2+9.0/1001.0*s-3.0/5005.0;
106    /* p=8 3 2 */
107    } else if( i==3 && j==2 ) {
108       t1 = r*r;
109       t2 = r*t1;
110       t3 = s*s;
111       f = t2*t3-3.0/7.0*t2*s+3.0/91.0*t2-6.0/7.0*t1*t3+36.0/91.0*t1*s-3.0/
112 91.0*t1+18.0/91.0*r*t3-9.0/91.0*r*s+9.0/1001.0*r-t3/91+6.0/1001.0*s-3.0/5005.0;
113    /* p=8 4 1 */
114    } else if( i==4 && j==1 ) {
115       t1 = r*r;
116       t2 = t1*t1;
117       t4 = r*t1;
118       f = t2*s-t2/7-10.0/7.0*t4*s+20.0/91.0*t4+60.0/91.0*t1*s-10.0/91.0*t1
119 -10.0/91.0*r*s+20.0/1001.0*r+5.0/1001.0*s-1.0/1001.0;
120    /* p=8 5 0 */
121    } else if( i==5 && j==0 ) {
122       t1 = r*r;
123       t2 = t1*t1;
124       f = r*t2-15.0/7.0*t2+150.0/91.0*r*t1-50.0/91.0*t1+75.0/1001.0*r-3.0/
125 1001.0;
126    /* p=9 0 6 */
127    } else if( i==0 && j==6 ) {
128       t1 = s*s;
129       t2 = t1*t1;
130       f = t1*t2-21.0/8.0*t2*s+21.0/8.0*t2-5.0/4.0*t1*s+15.0/52.0*t1-3.0/104.0*
131 s+1.0/1144.0;
132    /* p=9 1 5 */
133    } else if( i==1 && j==5 ) {
134       t1 = s*s;
135       t2 = t1*t1;
136       t3 = t2*s;
137       t6 = t1*s;
138       f = r*t3-15.0/8.0*r*t2+5.0/4.0*r*t6-5.0/14.0*r*t1+15.0/364.0*r*s-r/728-
139 t3/8+t2/4-5.0/28.0*t6+5.0/91.0*t1-5.0/728.0*s+1.0/4004.0;
140    /* p=9 2 4 */
141    } else if( i==2 && j==4 ) {
142       t1 = r*r;
143       t2 = s*s;
144       t3 = t2*t2;
145       t5 = t2*s;
146       f = t1*t3-5.0/4.0*t1*t5+t1*t2/2-t1*s/14+t1/364-3.0/8.0*r*t3+r*t5/2-3.0/
147 14.0*r*t2+3.0/91.0*r*s-r/728+t3/40-t5/28+3.0/182.0*t2-s/364+1.0/8008.0;
148    /* p=9 3 3 */
149    } else if( i==3 && j==3 ) {
150       t1 = r*r;
151       t2 = r*t1;
152       t3 = s*s;
153       t4 = t3*s;
154       f = t2*t4-3.0/4.0*t2*t3+3.0/20.0*t2*s-t2/140-3.0/4.0*t1*t4+3.0/5.0*t1*
155 t3-9.0/70.0*t1*s+3.0/455.0*t1+3.0/20.0*r*t4-9.0/70.0*r*t3+27.0/910.0*r*s-3.0/
156 1820.0*r-t4/140+3.0/455.0*t3-3.0/1820.0*s+1.0/10010.0;
157    /* p=9 4 2 */
158    } else if( i==4 && j==2 ) {
159       t1 = r*r;
160       t2 = t1*t1;
161       t3 = s*s;
162       t6 = r*t1;
163       f = t2*t3-3.0/8.0*t2*s+t2/40-5.0/4.0*t6*t3+t6*s/2-t6/28+t1*t3/2-3.0/
164 14.0*t1*s+3.0/182.0*t1-r*t3/14+3.0/91.0*r*s-r/364+t3/364-s/728+1.0/8008.0;
165    /* p=9 5 1 */
166    } else if( i==5 && j==1 ) {
167       t1 = r*r;
168       t2 = t1*t1;
169       t3 = r*t2;
170       t6 = r*t1;
171       f = t3*s-t3/8-15.0/8.0*t2*s+t2/4+5.0/4.0*t6*s-5.0/28.0*t6-5.0/14.0*t1*s
172 +5.0/91.0*t1+15.0/364.0*r*s-5.0/728.0*r-s/728+1.0/4004.0;
173    /* p=9 6 0 */
174    } else if( i==6 && j==0 ) {
175       t1 = r*r;
176       t2 = t1*t1;
177       f = t1*t2-21.0/8.0*r*t2+21.0/8.0*t2-5.0/4.0*r*t1+15.0/52.0*t1-3.0/104.0*
178 r+1.0/1144.0;
179    /* p=10 0 7 */
180    } else if( i==0 && j==7 ) {
181       t1 = s*s;
182       t2 = t1*s;
183       t3 = t1*t1;
184       f = t2*t3-28.0/9.0*t1*t3+196.0/51.0*t3*s-245.0/102.0*t3+245.0/306.0*t2
185 -7.0/51.0*t1+7.0/663.0*s-1.0/3978.0;
186    /* p=10 1 6 */
187    } else if( i==1 && j==6 ) {
188       t1 = s*s;
189       t2 = t1*t1;
190       t3 = t1*t2;
191       t5 = t2*s;
192       t8 = t1*s;
193       f = r*t3-7.0/3.0*r*t5+35.0/17.0*r*t2-175.0/204.0*r*t8+35.0/204.0*r*t1-r
194 *s/68+r/2652-t3/9+14.0/51.0*t5-35.0/136.0*t2+35.0/306.0*t8-5.0/204.0*t1+s/442
195 -1.0/15912.0;
196    /* p=10 2 5 */
197    } else if( i==2 && j==5 ) {
198       t1 = r*r;
199       t2 = s*s;
200       t3 = t2*t2;
201       t4 = t3*s;
202       t7 = t2*s;
203       f = t1*t4-5.0/3.0*t1*t3+50.0/51.0*t1*t7-25.0/102.0*t1*t2+5.0/204.0*t1*s
204 -t1/1428-r*t4/3+10.0/17.0*r*t3-25.0/68.0*r*t7+5.0/51.0*r*t2-5.0/476.0*r*s+r
205 /3094+t4/51-5.0/136.0*t3+5.0/204.0*t7-5.0/714.0*t2+5.0/6188.0*s-1.0/37128.0;
206    /* p=10 3 4 */
207    } else if( i==3 && j==4 ) {
208       t1 = r*r;
209       t2 = t1*r;
210       t3 = s*s;
211       t4 = t3*t3;
212       t6 = t3*s;
213       f = t2*t4-10.0/9.0*t2*t6+20.0/51.0*t2*t3-5.0/102.0*t2*s+t2/612-2.0/3.0*
214 t1*t4+40.0/51.0*t1*t6-5.0/17.0*t1*t3+2.0/51.0*t1*s-t1/714+2.0/17.0*r*t4-5.0/
215 34.0*r*t6+r*t3/17-r*s/119+r/3094-t4/204+t6/153-t3/357+2.0/4641.0*s-1.0/55692.0;
216    /* p=10 4 3 */
217    } else if( i==4 && j==3 ) {
218       t1 = r*r;
219       t2 = t1*t1;
220       t3 = s*s;
221       t4 = t3*s;
222       t8 = t1*r;
223       f = t2*t4-2.0/3.0*t2*t3+2.0/17.0*t2*s-t2/204-10.0/9.0*t8*t4+40.0/51.0*
224 t8*t3-5.0/34.0*t8*s+t8/153+20.0/51.0*t1*t4-5.0/17.0*t1*t3+t1*s/17-t1/357-5.0/
225 102.0*r*t4+2.0/51.0*r*t3-r*s/119+2.0/4641.0*r+t4/612-t3/714+s/3094-1.0/55692.0;
226    /* p=10 5 2 */
227    } else if( i==5 && j==2 ) {
228       t1 = r*r;
229       t2 = t1*t1;
230       t3 = r*t2;
231       t4 = s*s;
232       t9 = t1*r;
233       f = t3*t4-t3*s/3+t3/51-5.0/3.0*t2*t4+10.0/17.0*t2*s-5.0/136.0*t2+50.0/
234 51.0*t9*t4-25.0/68.0*t9*s+5.0/204.0*t9-25.0/102.0*t1*t4+5.0/51.0*t1*s-5.0/714.0
235 *t1+5.0/204.0*r*t4-5.0/476.0*r*s+5.0/6188.0*r-t4/1428+s/3094-1.0/37128.0;
236    /* p=10 6 1 */
237    } else if( i==6 && j==1 ) {
238       t1 = r*r;
239       t2 = t1*t1;
240       t3 = t1*t2;
241       t5 = r*t2;
242       t8 = t1*r;
243       f = t3*s-t3/9-7.0/3.0*t5*s+14.0/51.0*t5+35.0/17.0*t2*s-35.0/136.0*t2
244 -175.0/204.0*t8*s+35.0/306.0*t8+35.0/204.0*t1*s-5.0/204.0*t1-r*s/68+r/442+s
245 /2652-1.0/15912.0;
246    /* p=10 7 0 */
247    } else if( i==7 && j==0 ) {
248       t1 = r*r;
249       t2 = t1*r;
250       t3 = t1*t1;
251       f = t2*t3-28.0/9.0*t1*t3+196.0/51.0*r*t3-245.0/102.0*t3+245.0/306.0*t2
252 -7.0/51.0*t1+7.0/663.0*r-1.0/3978.0;
253    /* p=11 0 8 */
254    } else if( i==0 && j==8 ) {
255       t1 = s*s;
256       t2 = t1*t1;
257       t3 = t2*t2;
258       t4 = t1*s;
259       f = t3-18.0/5.0*t2*t4+504.0/95.0*t1*t2-392.0/95.0*t2*s+588.0/323.0*t2
260 -147.0/323.0*t4+98.0/1615.0*t1-6.0/1615.0*s+3.0/41990.0;
261    /* p=11 1 7 */
262    } else if( i==1 && j==7 ) {
263       t1 = s*s;
264       t2 = t1*s;
265       t3 = t1*t1;
266       t4 = t2*t3;
267       t6 = t1*t3;
268       t8 = t3*s;
269       f = r*t4-14.0/5.0*r*t6+294.0/95.0*r*t8-98.0/57.0*r*t3+490.0/969.0*r*t2
270 -49.0/646.0*t1*r+49.0/9690.0*r*s-r/9690-t4/10+28.0/95.0*t6-98.0/285.0*t8+196.0/
271 969.0*t3-245.0/3876.0*t2+49.0/4845.0*t1-7.0/9690.0*s+1.0/62985.0;
272    /* p=11 2 6 */
273    } else if( i==2 && j==6 ) {
274       t1 = r*r;
275       t2 = s*s;
276       t3 = t2*t2;
277       t4 = t2*t3;
278       t6 = t3*s;
279       t9 = t2*s;
280       t20 = 70.0/323.0*r*t9-105.0/2584.0*r*t2+21.0/6460.0*r*s-r/12920+3.0/190.0
281 *t4-7.0/190.0*t6+21.0/646.0*t3-35.0/2584.0*t9+7.0/2584.0*t2-3.0/12920.0*s+1.0/
282 167960.0;
283       f = t1*t4-21.0/10.0*t1*t6+63.0/38.0*t1*t3-35.0/57.0*t1*t9+35.0/323.0*t1
284 *t2-21.0/2584.0*t1*s+7.0/38760.0*t1-3.0/10.0*r*t4+63.0/95.0*r*t6-21.0/38.0*r*t3
285 +t20;
286    /* p=11 3 5 */
287    } else if( i==3 && j==5 ) {
288       t1 = r*r;
289       t2 = t1*r;
290       t3 = s*s;
291       t4 = t3*t3;
292       t5 = t4*s;
293       t8 = t3*s;
294       t17 = t2*t5-3.0/2.0*t2*t4+15.0/19.0*t2*t8-10.0/57.0*t2*t3+5.0/323.0*t2*s-
295 t2/2584-3.0/5.0*t1*t5+18.0/19.0*t1*t4-10.0/19.0*t1*t8+40.0/323.0*t1*t3-15.0/
296 1292.0*t1*s+t1/3230;
297       t23 = 9.0/95.0*r*t5-3.0/19.0*r*t4+30.0/323.0*r*t8-15.0/646.0*r*t3+3.0/
298 1292.0*r*s-3.0/45220.0*r-t5/285+2.0/323.0*t4-5.0/1292.0*t8+t3/969-s/9044+1.0/
299 293930.0;
300       f = t17+t23;
301    /* p=11 4 4 */
302    } else if( i==4 && j==4 ) {
303       t1 = r*r;
304       t2 = t1*t1;
305       t3 = s*s;
306       t4 = t3*s;
307       t6 = t1*r;
308       t7 = t3*t3;
309       t12 = 1.0/352716.0-t2*t4+t1/1938-5.0/3876.0*t6+t3/1938-5.0/3876.0*t4+t7
310 /969-s/13566-2.0/57.0*r*t7-20.0/57.0*t6*t3+2.0/969.0*r*s+20.0/19.0*t6*t4;
311       t24 = -5.0/323.0*r*t3+40.0/323.0*t1*t3-t6*t7+40.0/969.0*t6*s+6.0/19.0*t1*
312 t7-20.0/57.0*t1*t4+t2*t7-5.0/323.0*t1*s+6.0/19.0*t2*t3-2.0/57.0*t2*s+40.0/969.0
313 *r*t4+t2/969-r/13566;
314       f = t12+t24;
315    /* p=11 5 3 */
316    } else if( i==5 && j==3 ) {
317       t1 = r*r;
318       t2 = t1*t1;
319       t3 = r*t2;
320       t4 = s*s;
321       t5 = t4*s;
322       t12 = t1*r;
323       t16 = t3*t5-3.0/5.0*t3*t4+9.0/95.0*t3*s-t3/285-3.0/2.0*t2*t5+18.0/19.0*t2
324 *t4-3.0/19.0*t2*s+2.0/323.0*t2+15.0/19.0*t12*t5-10.0/19.0*t12*t4+30.0/323.0*t12
325 *s-5.0/1292.0*t12;
326       t23 = -10.0/57.0*t1*t5+40.0/323.0*t1*t4-15.0/646.0*t1*s+t1/969+5.0/323.0*
327 r*t5-15.0/1292.0*r*t4+3.0/1292.0*r*s-r/9044-t5/2584+t4/3230-3.0/45220.0*s+1.0/
328 293930.0;
329       f = t16+t23;
330    /* p=11 6 2 */
331    } else if( i==6 && j==2 ) {
332       t1 = r*r;
333       t2 = t1*t1;
334       t3 = t2*t1;
335       t4 = s*s;
336       t7 = r*t2;
337       t12 = t1*r;
338       t20 = 70.0/323.0*t12*s-35.0/2584.0*t12+35.0/323.0*t1*t4-105.0/2584.0*t1*s
339 +7.0/2584.0*t1-21.0/2584.0*r*t4+21.0/6460.0*r*s-3.0/12920.0*r+7.0/38760.0*t4-s
340 /12920+1.0/167960.0;
341       f = t3*t4-3.0/10.0*t3*s+3.0/190.0*t3-21.0/10.0*t7*t4+63.0/95.0*t7*s-7.0
342 /190.0*t7+63.0/38.0*t2*t4-21.0/38.0*t2*s+21.0/646.0*t2-35.0/57.0*t12*t4+t20;
343    /* p=11 7 1 */
344    } else if( i==7 && j==1 ) {
345       t1 = r*r;
346       t2 = t1*r;
347       t3 = t1*t1;
348       t4 = t2*t3;
349       t6 = t1*t3;
350       t8 = r*t3;
351       f = t4*s-t4/10-14.0/5.0*t6*s+28.0/95.0*t6+294.0/95.0*t8*s-98.0/285.0*t8
352 -98.0/57.0*t3*s+196.0/969.0*t3+490.0/969.0*t2*s-245.0/3876.0*t2-49.0/646.0*t1*s
353 +49.0/4845.0*t1+49.0/9690.0*r*s-7.0/9690.0*r-s/9690+1.0/62985.0;
354    /* p=11 8 0 */
355    } else if( i==8 && j==0 ) {
356       t1 = r*r;
357       t2 = t1*t1;
358       t3 = t2*t2;
359       t4 = t1*r;
360       f = t3-18.0/5.0*t2*t4+504.0/95.0*t2*t1-392.0/95.0*r*t2+588.0/323.0*t2
361 -147.0/323.0*t4+98.0/1615.0*t1-6.0/1615.0*r+3.0/41990.0;
362    /* p=12 0 9 */
363    } else if( i==0 && j==9 ) {
364       t1 = s*s;
365       t2 = t1*t1;
366       t3 = t2*t2;
367       t5 = t1*s;
368       f = t3*s-45.0/11.0*t3+540.0/77.0*t2*t5-72.0/11.0*t2*t1+756.0/209.0*t2*s
369 -252.0/209.0*t2+840.0/3553.0*t5-90.0/3553.0*t1+9.0/7106.0*s-1.0/49742.0;
370    /* p=12 1 8 */
371    } else if( i==1 && j==8 ) {
372       t1 = s*s;
373       t2 = t1*t1;
374       t3 = t2*t2;
375       t5 = t1*s;
376       t6 = t2*t5;
377       t8 = t2*t1;
378       t10 = t2*s;
379       f = r*t3-36.0/11.0*r*t6+48.0/11.0*r*t8-168.0/55.0*r*t10+252.0/209.0*r*
380 t2-56.0/209.0*r*t5+112.0/3553.0*t1*r-6.0/3553.0*r*s+r/35530-t3/11+24.0/77.0*t6
381 -24.0/55.0*t8+336.0/1045.0*t10-28.0/209.0*t2+112.0/3553.0*t5-14.0/3553.0*t1+4.0
382 /17765.0*s-1.0/248710.0;
383    /* p=12 2 7 */
384    } else if( i==2 && j==7 ) {
385       t1 = r*r;
386       t2 = s*s;
387       t3 = t2*s;
388       t4 = t2*t2;
389       t5 = t3*t4;
390       t7 = t2*t4;
391       t9 = t4*s;
392       t19 = t1*t5-28.0/11.0*t1*t7+28.0/11.0*t1*t9-14.0/11.0*t1*t4+70.0/209.0*t1
393 *t3-28.0/627.0*t2*t1+28.0/10659.0*t1*s-t1/21318-3.0/11.0*r*t5+8.0/11.0*r*t7
394 -42.0/55.0*r*t9+84.0/209.0*r*t4;
395       t23 = -70.0/627.0*r*t3+56.0/3553.0*r*t2-7.0/7106.0*r*s+r/53295+t5/77-2.0/
396 55.0*t7+42.0/1045.0*t9-14.0/627.0*t4+70.0/10659.0*t3-7.0/7106.0*t2+7.0/106590.0
397 *s-1.0/746130.0;
398       f = t19+t23;
399    /* p=12 3 6 */
400    } else if( i==3 && j==6 ) {
401       t1 = s*s;
402       t2 = t1*t1;
403       t3 = t2*t1;
404       t5 = r*r;
405       t6 = r*t5;
406       t7 = t1*s;
407       t10 = t2*s;
408       t21 = 6.0/77.0*r*t3-5.0/11.0*t6*t7+27.0/209.0*r*t2-9.0/55.0*r*t10+15.0/
409 11.0*t6*t2-21.0/11.0*t6*t10+t6*t3-1.0/1492260.0-10.0/209.0*r*t7+15.0/209.0*t6*
410 t1+12.0/3553.0*t5*s-6.0/11.0*t5*t3+12.0/11.0*t5*t10-9.0/11.0*t5*t2;
411       t27 = 60.0/209.0*t5*t7-t6*s/209-t5/14212-10.0/209.0*t1*t5-5.0/14212.0*t1+
412 30.0/3553.0*t1*r-9.0/14212.0*r*s+6.0/1045.0*t10+20.0/10659.0*t7-t2/209+s/35530+
413 t6/10659+r/71060-t3/385;
414       f = t21+t27;
415    /* p=12 4 5 */
416    } else if( i==4 && j==5 ) {
417       t1 = s*s;
418       t2 = t1*t1;
419       t4 = t2*s;
420       t6 = r*r;
421       t7 = r*t6;
422       t13 = t1*s;
423       t15 = t6*t6;
424       t22 = 60.0/1463.0*r*t2-2.0/77.0*r*t4+100.0/77.0*t7*t2-25.0/49742.0*r*s+
425 400.0/74613.0*t1*r-10.0/11.0*t7*t4-200.0/4389.0*t6*t1-1.0/2089164.0-100.0/
426 4389.0*r*t13+15.0/1463.0*t15*s-30.0/77.0*t6*t2-10.0/77.0*t15*t1+50.0/77.0*t15*
427 t13+20.0/77.0*t6*t4+300.0/1463.0*t6*t13;
428       t29 = t15*t4-15.0/11.0*t15*t2+100.0/24871.0*t6*s-50.0/4389.0*t7*s+200.0/
429 1463.0*t1*t7-t15/4389-5.0/49742.0*t6-25.0/149226.0*t1-5.0/4389.0*t2+t4/1463+5.0
430 /298452.0*s-50.0/77.0*t7*t13+50.0/74613.0*t13+r/74613+20.0/74613.0*t7;
431       f = t22+t29;
432    /* p=12 5 4 */
433    } else if( i==5 && j==4 ) {
434       t1 = r*r;
435       t2 = t1*t1;
436       t3 = r*t2;
437       t4 = s*s;
438       t5 = t4*s;
439       t10 = t4*t4;
440       t12 = t1*r;
441       t21 = -10.0/11.0*t5*t3+20.0/77.0*t3*t4+200.0/1463.0*t1*t5-5.0/4389.0*t2
442 -25.0/149226.0*t1-2.0/77.0*t3*s+15.0/1463.0*r*t10+50.0/77.0*t12*t10-25.0/
443 49742.0*r*s+100.0/24871.0*r*t4-200.0/4389.0*t1*t4-30.0/77.0*t2*t4+100.0/77.0*t5
444 *t2-10.0/77.0*t1*t10-15.0/11.0*t2*t10;
445       t29 = 400.0/74613.0*t1*s+60.0/1463.0*t2*s-100.0/4389.0*t12*s+t3*t10+300.0
446 /1463.0*t12*t4-50.0/4389.0*r*t5-50.0/77.0*t12*t5+50.0/74613.0*t12+20.0/74613.0*
447 t5-t10/4389+t3/1463-5.0/49742.0*t4+s/74613+5.0/298452.0*r-1.0/2089164.0;
448       f = t21+t29;
449    /* p=12 6 3 */
450    } else if( i==6 && j==3 ) {
451       t1 = s*s;
452       t2 = r*r;
453       t3 = t2*t2;
454       t4 = t3*t2;
455       t6 = t1*s;
456       t9 = r*t2;
457       t12 = t3*r;
458       t15 = -t1/14212-1.0/1492260.0-6.0/11.0*t1*t4+t4*t6+6.0/77.0*t4*s+20.0/
459 10659.0*t9-5.0/11.0*t9*t6-9.0/14212.0*r*s+6.0/1045.0*t12+12.0/11.0*t12*t1+t6
460 /10659-9.0/55.0*t12*s+s/71060+r/35530;
461       t27 = -t4/385+27.0/209.0*t3*s-10.0/209.0*t9*s+30.0/3553.0*t2*s+15.0/209.0
462 *t2*t6+60.0/209.0*t9*t1-r*t6/209+15.0/11.0*t3*t6-t3/209-5.0/14212.0*t2-10.0/
463 209.0*t2*t1-21.0/11.0*t12*t6+12.0/3553.0*t1*r-9.0/11.0*t3*t1;
464       f = t15+t27;
465    /* p=12 7 2 */
466    } else if( i==7 && j==2 ) {
467       t1 = r*r;
468       t2 = t1*r;
469       t3 = t1*t1;
470       t4 = t3*t2;
471       t5 = s*s;
472       t8 = t3*t1;
473       t11 = t3*r;
474       t16 = t4*t5-3.0/11.0*t4*s+t4/77-28.0/11.0*t8*t5+8.0/11.0*t8*s-2.0/55.0*t8
475 +28.0/11.0*t11*t5-42.0/55.0*t11*s+42.0/1045.0*t11-14.0/11.0*t5*t3+84.0/209.0*t3
476 *s-14.0/627.0*t3;
477       t23 = 70.0/209.0*t5*t2-70.0/627.0*t2*s+70.0/10659.0*t2-28.0/627.0*t1*t5+
478 56.0/3553.0*t1*s-7.0/7106.0*t1+28.0/10659.0*r*t5-7.0/7106.0*r*s+7.0/106590.0*r-
479 t5/21318+s/53295-1.0/746130.0;
480       f = t16+t23;
481    /* p=12 8 1 */
482    } else if( i==8 && j==1 ) {
483       t1 = r*r;
484       t2 = t1*t1;
485       t3 = t2*t2;
486       t5 = t1*r;
487       t6 = t5*t2;
488       t8 = t2*t1;
489       t10 = r*t2;
490       f = t3*s-t3/11-36.0/11.0*t6*s+24.0/77.0*t6+48.0/11.0*t8*s-24.0/55.0*t8
491 -168.0/55.0*t10*s+336.0/1045.0*t10+252.0/209.0*t2*s-28.0/209.0*t2-56.0/209.0*t5
492 *s+112.0/3553.0*t5+112.0/3553.0*t1*s-14.0/3553.0*t1-6.0/3553.0*r*s+4.0/17765.0*
493 r+s/35530-1.0/248710.0;
494    /* p=12 9 0 */
495    } else if( i==9 && j==0 ) {
496       t1 = r*r;
497       t2 = t1*t1;
498       t3 = t2*t2;
499       t5 = t1*r;
500       t9 = t3*r-45.0/11.0*t3+540.0/77.0*t5*t2-72.0/11.0*t2*t1+756.0/209.0*r*t2
501 -252.0/209.0*t2+840.0/3553.0*t5-90.0/3553.0*t1+9.0/7106.0*r-1.0/49742.0;
502    /* p=13 0 10 */
503    } else if( i==0 && j==10 ) {
504       t1 = s*s;
505       t2 = t1*t1;
506       t3 = t2*t2;
507       t6 = t1*s;
508       f = t3*t1-55.0/12.0*t3*s+825.0/92.0*t3-225.0/23.0*t2*t6+150.0/23.0*t2*
509 t1-63.0/23.0*t2*s+315.0/437.0*t2-50.0/437.0*t6+75.0/7429.0*t1-25.0/59432.0*s+
510 1.0/178296.0;
511    /* p=13 1 9 */
512    } else if( i==1 && j==9 ) {
513       t1 = s*s;
514       t2 = t1*t1;
515       t3 = t2*t2;
516       t4 = t3*s;
517       t7 = t1*s;
518       t8 = t2*t7;
519       t10 = t2*t1;
520       t12 = t2*s;
521       f = r*t4-15.0/4.0*r*t3+135.0/23.0*r*t8-1260.0/253.0*r*t10+630.0/253.0*r
522 *t12-189.0/253.0*r*t2+630.0/4807.0*r*t7-60.0/4807.0*r*t1+45.0/81719.0*r*s-5.0/
523 653752.0*r-t4/12+15.0/46.0*t3-135.0/253.0*t8+120.0/253.0*t10-63.0/253.0*t12+
524 378.0/4807.0*t2-70.0/4807.0*t7+120.0/81719.0*t1-45.0/653752.0*s+1.0/980628.0;
525    /* p=13 2 8 */
526    } else if( i==2 && j==8 ) {
527       t1 = r*r;
528       t2 = s*s;
529       t3 = t2*t2;
530       t4 = t3*s;
531       t8 = t3*t3;
532       t10 = t2*s;
533       t11 = t3*t10;
534       t17 = t3*t2;
535       t19 = -588.0/253.0*t1*t4-28.0/4807.0*r*t2+24.0/81719.0*r*s+t1*t8-3.0*t1*
536 t11-4.0/4807.0*t1*s-63.0/253.0*r*t3-42.0/253.0*t1*t10+210.0/253.0*t1*t3+84.0/
537 23.0*t1*t17-42.0/1265.0*t4+28.0/81719.0*t2+t1/81719;
538       t26 = -3.0/653752.0*r-14.0/4807.0*t10-3.0/163438.0*s+1.0/3268760.0+18.0/
539 23.0*r*t11-r*t8/4+252.0/4807.0*r*t10-252.0/253.0*r*t17+168.0/253.0*r*t4+84.0/
540 4807.0*t2*t1+63.0/4807.0*t3-9.0/253.0*t11+t8/92+12.0/253.0*t17;
541       f = t19+t26;
542    /* p=13 3 7 */
543    } else if( i==3 && j==7 ) {
544       t1 = r*r;
545       t2 = r*t1;
546       t3 = s*s;
547       t4 = t3*t3;
548       t5 = t3*s;
549       t6 = t4*t5;
550       t7 = t4*s;
551       t9 = t4*t3;
552       t14 = 1.0/7354710.0-t2/43263+14.0/4807.0*t4-t6/506-35.0/43263.0*t5-7.0/
553 980628.0*s+28.0/245157.0*t3-294.0/253.0*t1*t7+4.0/245157.0*t1-r/326876-7.0/
554 1265.0*t7+4.0/759.0*t9+28.0/23.0*t1*t9-14.0/4807.0*r*t3+14.0/81719.0*r*s-t1*t6
555 /2;
556       t31 = -14.0/14421.0*t1*s-21.0/253.0*r*t4+84.0/4807.0*t1*t3-35.0/253.0*t1*
557 t5+140.0/253.0*t1*t4+105.0/4807.0*r*t5+49.0/23.0*t2*t7-7.0/3.0*t2*t9+7.0/4807.0
558 *t2*s-7.0/253.0*t3*t2+175.0/759.0*t2*t5-245.0/253.0*t2*t4+3.0/46.0*r*t6+t2*t6
559 -42.0/253.0*r*t9+42.0/253.0*r*t7;
560       f = t14+t31;
561    /* p=13 4 6 */
562    } else if( i==4 && j==6 ) {
563       t1 = s*s;
564       t2 = t1*s;
565       t4 = r*r;
566       t5 = t4*t4;
567       t7 = t1*t1;
568       t8 = t7*s;
569       t9 = t7*t1;
570       t15 = 1.0/11767536.0+50.0/4807.0*r*t2+t5/19228+75.0/4807.0*t4*t1+15.0/
571 19228.0*t7-5.0/1961256.0*r-t8/1012+5.0/10626.0*t9-25.0/86526.0*t2-5.0/1307504.0
572 *s+25.0/490314.0*t1+5.0/245157.0*t4-105.0/253.0*t4*t8+5.0/23.0*t4*t9-25.0/
573 14421.0*r*t1+10.0/81719.0*r*s-5.0/4807.0*t4*s;
574       t21 = t4*r;
575       t34 = -15.0/506.0*r*t7-25.0/253.0*t4*t2+75.0/253.0*t4*t7+105.0/92.0*t5*t7
576 -175.0/506.0*t5*t2+15.0/4807.0*t21*s-5.0/86526.0*t21+t5*t9-7.0/4.0*t5*t8+25.0/
577 506.0*t5*t1-3.0/1012.0*t5*s+35.0/23.0*t21*t8-5.0/6.0*t21*t9-25.0/506.0*t21*t1+
578 250.0/759.0*t21*t2-525.0/506.0*t21*t7-5.0/253.0*r*t9+10.0/253.0*r*t8;
579       f = t15+t34;
580    /* p=13 5 5 */
581    } else if( i==5 && j==5 ) {
582       t1 = s*s;
583       t2 = t1*t1;
584       t4 = r*r;
585       t7 = t1*s;
586       t10 = t4*t4;
587       t11 = t10*r;
588       t14 = t2*s;
589       t19 = t4*r;
590       t26 = 1.0/13728792.0-75.0/7084.0*r*t2-125.0/100947.0*t4*s+500.0/33649.0*
591 t4*t1-125.0/1771.0*t4*t7+250.0/1771.0*t4*t2-5.0/4.0*t11*t2-375.0/506.0*t10*t7
592 -5.0/4.0*t10*t14+375.0/67298.0*r*t7+250.0/1771.0*t10*t1-75.0/7084.0*t10*s+375.0
593 /67298.0*t19*s+625.0/1771.0*t19*t7-375.0/506.0*t19*t2+25.0/3542.0*r*t14+t11*t14
594 +25.0/3542.0*t11*s;
595       t35 = -t14/7084-25.0/253.0*t4*t14-25.0/201894.0*t7+50.0/1716099.0*t4-25.0
596 /253.0*t11*t1+50.0/1716099.0*t1-25.0/9152528.0*r-t11/7084-125.0/100947.0*r*t1+
597 125.0/1144066.0*r*s+25.0/46.0*t11*t7-25.0/9152528.0*s-125.0/1771.0*t19*t1-25.0/
598 201894.0*t19+25.0/46.0*t19*t14+15.0/67298.0*t10+75.0/46.0*t10*t2+15.0/67298.0*
599 t2;
600       f = t26+t35;
601    /* p=13 6 4 */
602    } else if( i==6 && j==4 ) {
603       t1 = r*r;
604       t2 = t1*t1;
605       t3 = t2*t1;
606       t4 = s*s;
607       t5 = t4*t4;
608       t7 = r*t2;
609       t9 = t4*s;
610       t25 = t3*t5-105.0/253.0*t4*t7-5.0/6.0*t3*t9-5.0/4807.0*t4*r+10.0/253.0*t7
611 *s+10.0/81719.0*r*s+35.0/23.0*t7*t9-25.0/14421.0*t1*s-3.0/1012.0*r*t5+75.0/
612 4807.0*t4*t1-25.0/506.0*t1*t9+25.0/506.0*t5*t1-7.0/4.0*t5*t7+105.0/92.0*t5*t2
613 -525.0/506.0*t2*t9+15.0/4807.0*r*t9+75.0/253.0*t4*t2;
614       t28 = r*t1;
615       t34 = -15.0/506.0*t2*s-5.0/253.0*t3*s+50.0/4807.0*t28*s-25.0/253.0*t28*t4
616 +5.0/23.0*t3*t4+250.0/759.0*t28*t9-175.0/506.0*t28*t5+15.0/19228.0*t2+25.0/
617 490314.0*t1-5.0/86526.0*t9+5.0/10626.0*t3-5.0/1307504.0*r+t5/19228-5.0/
618 1961256.0*s-t7/1012+5.0/245157.0*t4-25.0/86526.0*t28+1.0/11767536.0;
619       f = t25+t34;
620    /* p=13 7 3 */
621    } else if( i==7 && j==3 ) {
622       t1 = s*s;
623       t3 = r*r;
624       t4 = t3*t3;
625       t5 = t4*t3;
626       t6 = t1*s;
627       t10 = t3*r;
628       t11 = t4*t10;
629       t14 = t4*r;
630       t16 = 1.0/7354710.0-s/326876-14.0/14421.0*r*t1+4.0/245157.0*t1-7.0/3.0*t5
631 *t6+4.0/759.0*t5-7.0/980628.0*r+14.0/81719.0*r*s+140.0/253.0*t4*t1+t11*t6+14.0/
632 4807.0*t4-t11/506-t11*t1/2-t6/43263-7.0/1265.0*t14-294.0/253.0*t14*t1;
633       t31 = 3.0/46.0*t11*s-35.0/43263.0*t10+175.0/759.0*t10*t6-21.0/253.0*t4*s
634 -245.0/253.0*t4*t6+49.0/23.0*t14*t6-35.0/253.0*t10*t1+7.0/4807.0*r*t6+28.0/
635 245157.0*t3-7.0/253.0*t3*t6+42.0/253.0*t14*s+84.0/4807.0*t3*t1-14.0/4807.0*t3*s
636 -42.0/253.0*t5*s+105.0/4807.0*t10*s+28.0/23.0*t5*t1;
637       f = t16+t31;
638    /* p=13 8 2 */
639    } else if( i==8 && j==2 ) {
640       t1 = r*r;
641       t2 = r*t1;
642       t4 = s*s;
643       t7 = t1*t1;
644       t8 = t7*t1;
645       t11 = t7*r;
646       t14 = t7*t7;
647       t17 = 1.0/3268760.0+28.0/81719.0*t1+252.0/4807.0*t2*s+84.0/4807.0*t4*t1
648 -28.0/4807.0*t1*s-252.0/253.0*t8*s+84.0/23.0*t8*t4+168.0/253.0*t11*s-42.0/253.0
649 *t2*t4+t14/92-t14*s/4-14.0/4807.0*t2+t14*t4;
650       t18 = t7*t2;
651       t26 = -9.0/253.0*t18+63.0/4807.0*t7+24.0/81719.0*r*s+210.0/253.0*t7*t4
652 -63.0/253.0*t7*s+18.0/23.0*t18*s-3.0*t18*t4-42.0/1265.0*t11-3.0/653752.0*s-4.0/
653 4807.0*t4*r+t4/81719+12.0/253.0*t8-3.0/163438.0*r-588.0/253.0*t11*t4;
654       f = t17+t26;
655    /* p=13 9 1 */
656    } else if( i==9 && j==1 ) {
657       t1 = r*r;
658       t2 = t1*t1;
659       t3 = t2*t2;
660       t4 = t3*r;
661       t7 = r*t1;
662       t8 = t7*t2;
663       t10 = t2*t1;
664       t12 = t2*r;
665       f = t4*s-t4/12-15.0/4.0*t3*s+15.0/46.0*t3+135.0/23.0*t8*s-135.0/253.0*
666 t8-1260.0/253.0*t10*s+120.0/253.0*t10+630.0/253.0*t12*s-63.0/253.0*t12-189.0/
667 253.0*t2*s+378.0/4807.0*t2+630.0/4807.0*t7*s-70.0/4807.0*t7-60.0/4807.0*t1*s+
668 120.0/81719.0*t1+45.0/81719.0*r*s-45.0/653752.0*r-5.0/653752.0*s+1.0/980628.0;
669    /* p=13 10 0 */
670    } else if( i==10 && j==0 ) {
671       t1 = r*r;
672       t2 = t1*t1;
673       t3 = t2*t2;
674       t6 = r*t1;
675       f = t3*t1-55.0/12.0*t3*r+825.0/92.0*t3-225.0/23.0*t2*t6+150.0/23.0*t2*
676 t1-63.0/23.0*t2*r+315.0/437.0*t2-50.0/437.0*t6+75.0/7429.0*t1-25.0/59432.0*r+
677 1.0/178296.0;
678    /* p=14 0 11 */
679    } else if( i==0 && j==11 ) {
680       t1 = s*s;
681       t2 = t1*s;
682       t3 = t1*t1;
683       t4 = t3*t3;
684       f = t4*t2-66.0/13.0*t4*t1+726.0/65.0*t4*s-363.0/26.0*t4+3267.0/299.0*t3
685 *t2-8316.0/1495.0*t3*t1+2772.0/1495.0*t3*s-594.0/1495.0*t3+297.0/5681.0*t2-22.0
686 /5681.0*t1+66.0/482885.0*s-3.0/1931540.0;
687    /* p=14 1 10 */
688    } else if( i==1 && j==10 ) {
689       t1 = s*s;
690       t2 = t1*t1;
691       t3 = t2*t2;
692       t4 = t3*t1;
693       t6 = t3*s;
694       t9 = t1*s;
695       t10 = t2*t9;
696       t12 = t2*t1;
697       t14 = t2*s;
698       t20 = r*t4-55.0/13.0*r*t6+99.0/13.0*r*t3-99.0/13.0*r*t10+1386.0/299.0*r*
699 t12-2646.0/1495.0*r*t14+126.0/299.0*r*t2-18.0/299.0*r*t9+27.0/5681.0*r*t1-r*s
700 /5681+r/482885;
701       t21 = -t4/13+22.0/65.0*t6-33.0/52.0*t3+198.0/299.0*t10-126.0/299.0*t12+
702 252.0/1495.0*t14-63.0/1495.0*t2+36.0/5681.0*t9-3.0/5681.0*t1+2.0/96577.0*s-1.0/
703 3863080.0;
704       f = t20+t21;
705    /* p=14 2 9 */
706    } else if( i==2 && j==9 ) {
707       t1 = s*s;
708       t2 = t1*t1;
709       t3 = t2*t2;
710       t6 = r*r;
711       t7 = t2*t1;
712       t9 = t1*s;
713       t10 = t2*t9;
714       t12 = t3*s;
715       t14 = t2*s;
716       t22 = 54.0/65.0*r*t3+648.0/312455.0*r*t1-252.0/65.0*t6*t7+324.0/65.0*t6*
717 t10+t6*t12-3.0/42493880.0+378.0/16445.0*t14+2646.0/1495.0*t6*t14-567.0/82225.0*
718 t2-7938.0/16445.0*t6*t2+81.0/312455.0*t6*s-45.0/13.0*t6*t3-378.0/16445.0*r*t9
719 -108.0/16445.0*t6*t1+252.0/3289.0*t6*t9;
720       t29 = -7938.0/16445.0*r*t14+1512.0/1495.0*r*t7+2268.0/16445.0*r*t2-756.0/
721 16445.0*t7-9.0/260.0*t3-t6/312455+81.0/1495.0*t10-81.0/65.0*r*t10-27.0/312455.0
722 *r*s-36.0/312455.0*t1+6.0/5311735.0*r+3.0/325.0*t12+378.0/312455.0*t9+27.0/
723 5311735.0*s-3.0/13.0*r*t12;
724       f = t22+t29;
725    /* p=14 3 8 */
726    } else if( i==3 && j==8 ) {
727       t1 = r*r;
728       t2 = r*t1;
729       t3 = s*s;
730       t4 = t3*s;
731       t7 = t3*t3;
732       t8 = t7*t3;
733       t10 = t7*s;
734       t13 = t7*t7;
735       t15 = t7*t4;
736       t25 = -1764.0/16445.0*t4*t2+168.0/16445.0*t3*t2+1008.0/325.0*t2*t8-588.0/
737 325.0*t2*t10+882.0/1495.0*t2*t7+t2*t13-3.0/106234700.0+36.0/7475.0*t15-6.0/
738 1562275.0*t1-72.0/1562275.0*r*s-54.0/325.0*r*t15+18.0/325.0*r*t13+1512.0/
739 1562275.0*r*t3-504.0/325.0*t1*t8+432.0/325.0*t1*t15-6.0/13.0*t1*t13-756.0/
740 82225.0*r*t4+432.0/1562275.0*t1*s;
741       t35 = -504.0/82225.0*t3*t1+1008.0/16445.0*t4*t1-5292.0/16445.0*t1*t7+
742 7056.0/7475.0*t1*t10-10584.0/82225.0*r*t10+1512.0/7475.0*r*t8-36.0/13.0*t2*t15+
743 756.0/16445.0*r*t7-126.0/82225.0*t7-t13/650-504.0/82225.0*t8+336.0/82225.0*t10
744 -36.0/82225.0*t2*s+504.0/1562275.0*t4+48.0/26558675.0*s-56.0/1562275.0*t3+9.0/
745 1562275.0*t2+18.0/26558675.0*r;
746       f = t25+t35;
747    /* p=14 4 7 */
748    } else if( i==4 && j==7 ) {
749       t1 = r*r;
750       t2 = r*t1;
751       t3 = s*s;
752       t4 = t3*s;
753       t7 = t3*t3;
754       t8 = t7*t3;
755       t12 = t7*t4;
756       t14 = t1*t1;
757       t22 = t7*s;
758       t24 = 4.0/312455.0*t2+8.0/15935205.0*r-14.0/937365.0*t3+14.0/15935205.0*s
759 +7.0/62491.0*t4-490.0/3289.0*t4*t2+56.0/3289.0*t3*t2+112.0/65.0*t2*t8+196.0/
760 299.0*t2*t7-28.0/937365.0*r*s-r*t12/65-1.0/63740820.0-294.0/16445.0*t14*t3+56.0
761 /3289.0*r*t7+14.0/16445.0*t14*s-10.0/13.0*t2*t12+49.0/299.0*t14*t4-49.0/65.0*
762 t14*t7+56.0/1495.0*r*t8+588.0/1495.0*t1*t22;
763       t39 = -14.0/3289.0*r*t4+84.0/312455.0*t1*s-84.0/16445.0*t3*t1+140.0/
764 3289.0*t4*t1-588.0/3289.0*t1*t7-588.0/16445.0*r*t22-4.0/937365.0*t1-t14/82225+
765 588.0/325.0*t14*t22+t14*t12-28.0/13.0*t14*t8+168.0/312455.0*r*t3-28.0/65.0*t1*
766 t8+12.0/65.0*t1*t12+t12/2990-7.0/16445.0*t7+14.0/16445.0*t22-98.0/65.0*t2*t22
767 -14.0/16445.0*t8-14.0/16445.0*t2*s;
768       f = t24+t39;
769    /* p=14 5 6 */
770    } else if( i==5 && j==6 ) {
771       t1 = r*r;
772       t2 = r*t1;
773       t3 = s*s;
774       t4 = t3*t3;
775       t5 = t4*t3;
776       t6 = t4*s;
777       t8 = t3*s;
778       t16 = t1*t1;
779       t21 = t16*r;
780       t24 = 75.0/3289.0*t3*t2+6.0/13.0*t2*t5-21.0/26.0*t2*t6+315.0/598.0*t4*t2
781 -3.0/124982.0*r*s+126.0/65.0*t16*t6-15.0/13.0*t16*t5+45.0/124982.0*r*t3-t1*t5
782 /13+21.0/598.0*t21*t3-15.0/6578.0*r*t8;
783       t44 = t21*t5-21.0/13.0*t21*t6+105.0/299.0*t16*t8-63.0/32890.0*t21*s-63.0/
784 52.0*t16*t4+9.0/3289.0*t16*s+42.0/299.0*t1*t6-t1/187473+t21/32890-3.0/65780.0*
785 t16-1.0/84987760.0;
786       f = 3.0/124982.0*t2-3.0/32890.0*t5+3.0/16445.0*t6-9.0/65780.0*t4-9.0/
787 6578.0*t2*s+r/2124694+3.0/62491.0*t8+3.0/5311735.0*s-t3/124982-525.0/3289.0*t2*
788 t8+t24+18.0/62491.0*t1*s-15.0/3289.0*t3*t1+100.0/3289.0*t1*t8-315.0/3289.0*t4*
789 t1-7.0/26.0*t21*t8-63.0/6578.0*r*t6+3.0/598.0*r*t5+63.0/65.0*t21*t4+45.0/6578.0
790 *r*t4-315.0/6578.0*t16*t3+t44;
791    /* p=14 6 5 */
792    } else if( i==6 && j==5 ) {
793       t1 = r*r;
794       t2 = s*s;
795       t3 = t2*t2;
796       t6 = t1*t1;
797       t7 = r*t6;
798       t9 = t3*s;
799       t12 = t6*t1;
800       t13 = t2*s;
801       t24 = t1*r;
802       t27 = 3.0/124982.0*t13+t12*t9-315.0/3289.0*t6*t2-1.0/84987760.0-3.0/
803 124982.0*r*s+3.0/5311735.0*r+s/2124694-9.0/6578.0*r*t13+105.0/299.0*t24*t3+3.0/
804 62491.0*t24-7.0/26.0*t24*t9;
805       t44 = 75.0/3289.0*t1*t13-3.0/65780.0*t3-63.0/52.0*t6*t3-21.0/26.0*t7*t13
806 -3.0/32890.0*t12-t12*t2/13+3.0/16445.0*t7+42.0/299.0*t7*t2-15.0/6578.0*t24*s-t1
807 /124982+21.0/598.0*t1*t9;
808       f = -315.0/6578.0*t3*t1-15.0/3289.0*t1*t2-63.0/6578.0*t7*s-63.0/32890.0
809 *r*t9+9.0/3289.0*r*t3+6.0/13.0*t12*t13+3.0/598.0*t12*s-21.0/13.0*t7*t9+45.0/
810 124982.0*t1*s+126.0/65.0*t7*t3+t27+t9/32890+45.0/6578.0*t6*s-9.0/65780.0*t6+
811 18.0/62491.0*r*t2-t2/187473+315.0/598.0*t6*t13-525.0/3289.0*t24*t13+100.0/
812 3289.0*t24*t2+63.0/65.0*t6*t9-15.0/13.0*t12*t3+t44;
813    /* p=14 7 4 */
814    } else if( i==7 && j==4 ) {
815       t1 = r*r;
816       t2 = t1*r;
817       t3 = t1*t1;
818       t4 = t3*t2;
819       t5 = s*s;
820       t6 = t5*s;
821       t9 = r*t3;
822       t12 = t5*t5;
823       t20 = t3*t1;
824       t24 = -1.0/63740820.0-14.0/937365.0*t1-10.0/13.0*t4*t6+12.0/65.0*t4*t5
825 -98.0/65.0*t6*t9+56.0/3289.0*t6*t1+t4*t12-588.0/3289.0*t3*t5-294.0/16445.0*t1*
826 t12-84.0/16445.0*t1*t5+588.0/325.0*t12*t9-588.0/16445.0*t9*s+14.0/16445.0*r*t12
827 +56.0/1495.0*t20*s+t4/2990-14.0/16445.0*t20-t4*s/65+14.0/16445.0*t9+140.0/
828 3289.0*t2*t5-t12/82225;
829       t39 = -49.0/65.0*t12*t3-7.0/16445.0*t3+7.0/62491.0*t2-28.0/937365.0*r*s
830 -4.0/937365.0*t5+84.0/312455.0*r*t5+14.0/15935205.0*r+4.0/312455.0*t6+8.0/
831 15935205.0*s-14.0/16445.0*r*t6+56.0/3289.0*t3*s+49.0/299.0*t12*t2-490.0/3289.0*
832 t6*t2+196.0/299.0*t6*t3+588.0/1495.0*t9*t5-14.0/3289.0*t2*s-28.0/65.0*t20*t5
833 -28.0/13.0*t20*t12+112.0/65.0*t20*t6+168.0/312455.0*t1*s;
834       f = t24+t39;
835    /* p=14 8 3 */
836    } else if( i==8 && j==3 ) {
837       t2 = s*s;
838       t4 = t2*s;
839       t6 = r*r;
840       t7 = t6*r;
841       t8 = t6*t6;
842       t10 = t8*t8;
843       t16 = t8*t6;
844       t18 = t8*t7;
845       t22 = 48.0/26558675.0*r-72.0/1562275.0*r*s+432.0/1562275.0*r*t2-36.0/
846 82225.0*r*t4+504.0/1562275.0*t7+756.0/16445.0*t8*s-126.0/82225.0*t8-6.0/13.0*
847 t10*t2+1008.0/16445.0*t7*t2+882.0/1495.0*t8*t4+18.0/325.0*t10*s-1764.0/16445.0*
848 t7*t4-504.0/325.0*t16*t2-504.0/82225.0*t16-36.0/13.0*t18*t4+36.0/7475.0*t18
849 -54.0/325.0*t18*s-756.0/82225.0*t7*s;
850       t23 = t8*r;
851       t35 = 336.0/82225.0*t23+7056.0/7475.0*t23*t2-504.0/82225.0*t6*t2-5292.0/
852 16445.0*t8*t2+168.0/16445.0*t6*t4-56.0/1562275.0*t6-588.0/325.0*t23*t4+t10*t4+
853 432.0/325.0*t18*t2+1512.0/7475.0*t16*s+1008.0/325.0*t16*t4-10584.0/82225.0*t23*
854 s+1512.0/1562275.0*t6*s+9.0/1562275.0*t4+18.0/26558675.0*s-3.0/106234700.0-6.0/
855 1562275.0*t2-t10/650;
856       f = t22+t35;
857    /* p=14 9 2 */
858    } else if( i==9 && j==2 ) {
859       t1 = s*s;
860       t4 = r*r;
861       t5 = t4*t4;
862       t6 = t5*t5;
863       t9 = t6*r;
864       t12 = r*t4;
865       t15 = t5*r;
866       t18 = t5*t12;
867       t20 = t5*t4;
868       t24 = 81.0/312455.0*r*t1-27.0/312455.0*r*s+54.0/65.0*t6*s+2268.0/16445.0*
869 t5*s+t9*t1-3.0/13.0*t9*s+252.0/3289.0*t12*t1-45.0/13.0*t6*t1+2646.0/1495.0*t15*
870 t1-378.0/16445.0*t12*s-81.0/65.0*t18*s-3.0/42493880.0+1512.0/1495.0*t20*s+648.0
871 /312455.0*t4*s-7938.0/16445.0*t15*s;
872       t29 = -36.0/312455.0*t4+324.0/65.0*t18*t1-7938.0/16445.0*t5*t1-108.0/
873 16445.0*t4*t1-252.0/65.0*t20*t1+81.0/1495.0*t18-756.0/16445.0*t20+378.0/16445.0
874 *t15+3.0/325.0*t9+378.0/312455.0*t12-567.0/82225.0*t5-9.0/260.0*t6+6.0/
875 5311735.0*s+27.0/5311735.0*r-t1/312455;
876       f = t24+t29;
877    /* p=14 10 1 */
878    } else if( i==10 && j==1 ) {
879       t1 = r*r;
880       t2 = t1*t1;
881       t3 = t2*t2;
882       t4 = t3*t1;
883       t6 = t3*r;
884       t9 = r*t1;
885       t10 = t2*t9;
886       t12 = t2*t1;
887       t14 = r*t2;
888       t16 = t4*s-t4/13-55.0/13.0*t6*s+22.0/65.0*t6+99.0/13.0*t3*s-33.0/52.0*t3
889 -99.0/13.0*t10*s+198.0/299.0*t10+1386.0/299.0*t12*s-126.0/299.0*t12-2646.0/
890 1495.0*t14*s;
891       t21 = 252.0/1495.0*t14+126.0/299.0*t2*s-63.0/1495.0*t2-18.0/299.0*t9*s+
892 36.0/5681.0*t9+27.0/5681.0*t1*s-3.0/5681.0*t1-r*s/5681+2.0/96577.0*r+s/482885
893 -1.0/3863080.0;
894       f = t16+t21;
895    /* p=14 11 0 */
896    } else if( i==11 && j==0 ) {
897       t1 = r*r;
898       t2 = r*t1;
899       t3 = t1*t1;
900       t4 = t3*t3;
901       f = t4*t2-66.0/13.0*t4*t1+726.0/65.0*r*t4-363.0/26.0*t4+3267.0/299.0*t3
902 *t2-8316.0/1495.0*t3*t1+2772.0/1495.0*t3*r-594.0/1495.0*t3+297.0/5681.0*t2-22.0
903 /5681.0*t1+66.0/482885.0*r-3.0/1931540.0;
904    /* p=15 0 12 */
905    } else if( i==0 && j==12 ) {
906       t1 = s*s;
907       t2 = t1*t1;
908       t3 = t2*t2;
909       t5 = t1*s;
910       f = t3*t2-39.0/7.0*t3*t5+286.0/21.0*t3*t1-1210.0/63.0*t3*s+121.0/7.0*t3
911 -363.0/35.0*t2*t5+484.0/115.0*t2*t1-132.0/115.0*t2*s+33.0/161.0*t2-11.0/483.0*
912 t5+22.0/15295.0*t1-2.0/45885.0*s+1.0/2340135.0;
913    /* p=15 1 11 */
914    } else if( i==1 && j==11 ) {
915       t1 = s*s;
916       t2 = t1*s;
917       t3 = t1*t1;
918       t4 = t3*t3;
919       t5 = t4*t2;
920       t7 = t4*t1;
921       t9 = t4*s;
922       t12 = t3*t2;
923       t14 = t3*t1;
924       t16 = t3*s;
925       t22 = r*t5-33.0/7.0*r*t7+605.0/63.0*r*t9-3025.0/273.0*r*t4+726.0/91.0*r*
926 t12-242.0/65.0*r*t14+1694.0/1495.0*r*t16-66.0/299.0*t3*r+55.0/2093.0*r*t2-11.0/
927 6279.0*r*t1+11.0/198835.0*r*s-r/1789515;
928       t23 = -t5/14+22.0/63.0*t7-605.0/819.0*t9+242.0/273.0*t4-121.0/182.0*t12+
929 484.0/1495.0*t14-154.0/1495.0*t16+44.0/2093.0*t3-11.0/4186.0*t2+22.0/119301.0*
930 t1-11.0/1789515.0*s+2.0/30421755.0;
931       f = t22+t23;
932    /* p=15 2 10 */
933    } else if( i==2 && j==10 ) {
934       t1 = s*s;
935       t2 = t1*t1;
936       t3 = r*r;
937       t5 = t1*s;
938       t6 = t2*t5;
939       t9 = t2*t2;
940       t13 = t2*s;
941       t15 = t2*t1;
942       t17 = t9*t1;
943       t19 = t9*s;
944       t25 = t2/299-t3*s/12558-550.0/91.0*t3*t6+20.0/2093.0*r*t5+275.0/42.0*t3*
945 t9+5.0/2093.0*t3*t1-10.0/299.0*t3*t5-77.0/65.0*t3*t13+44.0/13.0*t3*t15+t3*t17
946 -55.0/14.0*t3*t19-275.0/182.0*r*t9-21.0/299.0*r*t2+55.0/63.0*r*t19+462.0/1495.0
947 *r*t13+t3/1193010;
948       t32 = -11.0/182.0*t6+3.0/79534.0*t1-3.0/4186.0*r*t1-s/715806-t5/2093-r
949 /3579030+11.0/299.0*t15-11.0/13.0*r*t15+r*s/39767+132.0/91.0*r*t6+1.0/
950 60843510.0-3.0/14.0*r*t17+11.0/182.0*t9-55.0/1638.0*t19+t17/126-21.0/1495.0*t13
951 +77.0/299.0*t3*t2;
952       f = t25+t32;
953    /* p=15 3 9 */
954    } else if( i==3 && j==9 ) {
955       t1 = r*r;
956       t2 = s*s;
957       t4 = t2*s;
958       t6 = t2*t2;
959       t8 = t6*s;
960       t11 = t1*r;
961       t12 = t6*t2;
962       t15 = t6*t6;
963       t16 = t15*s;
964       t21 = t6*t4;
965       t28 = 48.0/23023.0*t2*t1-84.0/3289.0*t4*t1+252.0/1495.0*t1*t6-42.0/65.0*
966 t1*t8+27.0/2187185.0*r*s+2.0/334639305.0-40.0/13.0*t11*t12-21.0/65.0*t11*t6+t11
967 *t16-45.0/14.0*t11*t15+3.0/23023.0*t11*s-36.0/115115.0*r*t2-180.0/91.0*t1*t21+
968 12.0/3289.0*r*t4+10.0/7.0*t1*t15-9.0/115115.0*t1*s-r/6561555+8.0/1495.0*t12+
969 14.0/299.0*t11*t4+96.0/65.0*t1*t12;
970       t39 = -3.0/7.0*t1*t16-15.0/91.0*r*t15-378.0/16445.0*r*t6+r*t16/21+126.0/
971 1495.0*r*t8+2.0/2187185.0*t1+108.0/455.0*r*t21+24.0/2187185.0*t2-12.0/3289.0*
972 t11*t2-42.0/16445.0*t8+12.0/16445.0*t6-t16/819+2.0/455.0*t15+84.0/65.0*t11*t8-s
973 /2187185-2.0/16445.0*t4-t11/690690-12.0/65.0*r*t12-3.0/455.0*t21+30.0/7.0*t11*
974 t21;
975       f = t28+t39;
976    /* p=15 4 8 */
977    } else if( i==4 && j==8 ) {
978       t1 = s*s;
979       t2 = t1*t1;
980       t3 = t2*s;
981       t4 = r*r;
982       t5 = r*t4;
983       t6 = t2*t1;
984       t10 = t2*t2;
985       t14 = t1*s;
986       t17 = 4.0/937365.0*t1+2.0/9867.0*t2-28.0/49335.0*t3-80.0/39.0*t5*t6+224.0
987 /195.0*t5*t3-14.0/39.0*t5*t2-5.0/7.0*t5*t10+16.0/69069.0*t5*s-56.0/9867.0*t5*t1
988 +56.0/897.0*t5*t14-8.0/49335.0*t1*r;
989       t18 = t2*t14;
990       t20 = t4*t4;
991       t31 = 40.0/21.0*t5*t18-14.0/195.0*t20*t14+16.0/455.0*r*t18-40.0/91.0*t4*
992 t18-8.0/115115.0*t4*s+16.0/9867.0*t4*t1-56.0/195.0*t4*t3+32.0/65.0*t6*t4-10.0/
993 819.0*r*t10+28.0/65.0*t20*t2-4.0/16445.0*t20*s;
994       t37 = 112.0/4485.0*t3*r-56.0/39.0*t20*t3-t5/345345-t18/1365-18.0/7.0*t20*
995 t18+t20/345345+t10/4095+16.0/2187185.0*r*s-2.0/49335.0*t14+4.0/4485.0*t6-4.0/
996 19684665.0*s;
997       t47 = -8.0/195.0*r*t6-2.0/19684665.0*r+t20*t10+1.0/334639305.0+28.0/
998 4485.0*t20*t1+8.0/3.0*t20*t6-28.0/3289.0*r*t2+28.0/299.0*t4*t2-56.0/3289.0*t4*
999 t14+10.0/63.0*t4*t10+16.0/9867.0*r*t14+2.0/2187185.0*t4;
1000       f = t17+t31+t37+t47;
1001    /* p=15 5 7 */
1002    } else if( i==5 && j==7 ) {
1003       t1 = s*s;
1004       t2 = t1*t1;
1005       t3 = t2*t1;
1006       t6 = r*r;
1007       t7 = t6*t6;
1008       t8 = t7*r;
1009       t9 = t2*s;
1010       t11 = t1*s;
1011       t12 = t11*t2;
1012       t14 = t6*r;
1013       t22 = 2.0/1003917915.0-t3*r/117+r*s/187473+14.0/9.0*t8*t9-15.0/14.0*t7*
1014 t12-100.0/117.0*t14*t3+28.0/39.0*t14*t9+10.0/29601.0*t14*s-70.0/9867.0*t14*t1+
1015 175.0/2691.0*t14*t11-r*t1/9867+25.0/63.0*t14*t12;
1016       t33 = -35.0/234.0*t7*t11+r*t12/273-50.0/819.0*t6*t12-2.0/29601.0*t6*s+
1017 40.0/29601.0*t1*t6-2.0*t8*t3-350.0/29601.0*t11*t6+16.0/117.0*t6*t3+28.0/39.0*t7
1018 *t2-35.0/9867.0*r*t2-t8/148005-7.0/49335.0*t9;
1019       t38 = -t11/59202-s/8436285-r/11810799+2.0/13455.0*t3+140.0/2691.0*t6*t2+
1020 25.0/29601.0*r*t11+4.0/3936933.0*t6+2.0/937365.0*t1+2.0/29601.0*t2-35.0/117.0*
1021 t2*t14+2.0/207207.0*t7-t12/16380;
1022       t50 = -t14/207207+7.0/13455.0*t8*s+t8*t12+14.0/897.0*t7*t1-7.0/585.0*t1*
1023 t8+20.0/9.0*t7*t3-70.0/117.0*t8*t2+14.0/117.0*t11*t8-7.0/9867.0*t7*s+7.0/897.0*
1024 r*t9-14.0/117.0*t6*t9-70.0/39.0*t7*t9;
1025       f = t22+t33+t38+t50;
1026    /* p=15 6 6 */
1027    } else if( i==6 && j==6 ) {
1028       t1 = r*r;
1029       t2 = r*t1;
1030       t3 = s*s;
1031       t4 = t3*t3;
1032       t5 = t4*t3;
1033       t7 = t4*s;
1034       t10 = t3*s;
1035       t12 = t1*t1;
1036       t14 = t12*t1;
1037       t19 = r*t12;
1038       t21 = -25.0/117.0*t2*t5+14.0/39.0*t7*t2-35.0/156.0*t4*t2+1.0/573667380.0
1039 -175.0/19734.0*t1*t10+t1/749892-35.0/156.0*t12*t10-t14*s/780+5.0/9867.0*r*t10-
1040 t1*s/13156+25.0/19734.0*t3*t1-3.0/2.0*t19*t5;
1041       t34 = 35.0/1196.0*t4*t1-7.0/156.0*t7*t1+t1*t5/39+21.0/26.0*t12*t4-21.0/
1042 13156.0*r*t4-21.0/13156.0*t12*s+7.0/2990.0*t7*r-35.0/26.0*t7*t12-7.0/156.0*t19*
1043 t3+14.0/39.0*t19*t10+t14*t5+t14*t3/39;
1044       t42 = 5.0/9867.0*t2*s-175.0/19734.0*t3*t2+175.0/2691.0*t2*t10-t3*r/13156+
1045 t3/749892-t2/118404+t12/39468-3.0/2.0*t14*t7-7.0/197340.0*t19-7.0/197340.0*t7+
1046 3.0/624910.0*r*s-t10/118404;
1047       t51 = 7.0/3.0*t19*t7+t4/39468+t14/53820-r*t5/780-s/11248380+t5/53820-r
1048 /11248380+5.0/6.0*t14*t4+35.0/1196.0*t12*t3-25.0/117.0*t14*t10+5.0/6.0*t12*t5
1049 -35.0/26.0*t19*t4+7.0/2990.0*t19*s;
1050       f = t21+t34+t42+t51;
1051    /* p=15 7 5 */
1052    } else if( i==7 && j==5 ) {
1053       t1 = r*r;
1054       t2 = t1*t1;
1055       t3 = t2*t1;
1056       t4 = s*s;
1057       t5 = t4*t4;
1058       t7 = t1*r;
1059       t8 = t7*t2;
1060       t11 = t5*s;
1061       t13 = t4*s;
1062       t14 = t2*r;
1063       t20 = 20.0/9.0*t3*t5+2.0/1003917915.0-50.0/819.0*t8*t4+t8*s/273+t11*t8-
1064 t13/207207-7.0/49335.0*t14+r*s/187473+14.0/9.0*t14*t11-2.0*t3*t11+14.0/117.0*t7
1065 *t11-35.0/234.0*t7*t5;
1066       t33 = 16.0/117.0*t4*t3+25.0/29601.0*t7*s-350.0/29601.0*t7*t4+175.0/2691.0
1067 *t7*t13-2.0/29601.0*r*t4-35.0/117.0*t2*t13+25.0/63.0*t8*t13+10.0/29601.0*r*t13-
1068 t1*s/9867+40.0/29601.0*t4*t1-70.0/9867.0*t1*t13+14.0/897.0*t1*t5;
1069       t47 = -7.0/585.0*t1*t11+28.0/39.0*t2*t5-7.0/9867.0*r*t5-35.0/9867.0*t2*s+
1070 7.0/13455.0*r*t11-70.0/117.0*t2*t11-14.0/117.0*t14*t4+28.0/39.0*t14*t13-100.0/
1071 117.0*t3*t13-15.0/14.0*t8*t5-70.0/39.0*t14*t5+140.0/2691.0*t4*t2;
1072       t50 = 7.0/897.0*t14*s+2.0/13455.0*t3+2.0/207207.0*t5+4.0/3936933.0*t4-s
1073 /11810799-t11/148005-t8/16380-r/8436285+2.0/29601.0*t2-t7/59202+2.0/937365.0*t1
1074 -t3*s/117;
1075       f = t20+t33+t47+t50;
1076    /* p=15 8 4 */
1077    } else if( i==8 && j==4 ) {
1078       t1 = r*r;
1079       t2 = t1*t1;
1080       t3 = t2*t2;
1081       t4 = s*s;
1082       t5 = t4*s;
1083       t7 = t1*r;
1084       t9 = t4*t4;
1085       t10 = t2*t1;
1086       t12 = t7*t2;
1087       t13 = 1.0/334639305.0-5.0/7.0*t3*t5-2.0/49335.0*t7+2.0/9867.0*t2+16.0/
1088 2187185.0*r*s+t3/4095-t5/345345+t9/345345+8.0/3.0*t10*t9-t12/1365+4.0/4485.0*
1089 t10;
1090       t22 = -10.0/819.0*t3*s-2.0/19684665.0*s-4.0/19684665.0*r+2.0/2187185.0*t4
1091 -40.0/91.0*t12*t4+16.0/455.0*t12*s+10.0/63.0*t4*t3-14.0/195.0*t7*t9+32.0/65.0*
1092 t10*t4+16.0/9867.0*t7*s-56.0/3289.0*t7*t4;
1093       t27 = t2*r;
1094       t34 = -8.0/115115.0*r*t4-14.0/39.0*t2*t5+28.0/299.0*t4*t2-56.0/195.0*t27*
1095 t4-28.0/49335.0*t27+224.0/195.0*t27*t5-28.0/3289.0*t2*s+28.0/4485.0*t1*t9+4.0/
1096 937365.0*t1-8.0/195.0*t10*s+40.0/21.0*t12*t5;
1097       t47 = 16.0/69069.0*r*t5-8.0/49335.0*t1*s+16.0/9867.0*t4*t1-56.0/9867.0*t1
1098 *t5+28.0/65.0*t2*t9-4.0/16445.0*r*t9-80.0/39.0*t10*t5-18.0/7.0*t12*t9-56.0/39.0
1099 *t27*t9+112.0/4485.0*t27*s+t3*t9+56.0/897.0*t7*t5;
1100       f = t13+t22+t34+t47;
1101    /* p=15 9 3 */
1102    } else if( i==9 && j==3 ) {
1103       t1 = s*s;
1104       t3 = r*r;
1105       t4 = t3*t3;
1106       t6 = t4*t3;
1107       t9 = t3*r;
1108       t10 = t4*t9;
1109       t12 = t4*t4;
1110       t16 = t12*r;
1111       t17 = t1*s;
1112       t27 = r*t4;
1113       t30 = -9.0/115115.0*r*t1+252.0/1495.0*t4*t1+96.0/65.0*t6*t1-12.0/65.0*t6*
1114 s+2.0/334639305.0+108.0/455.0*t10*s-15.0/91.0*t12*s+27.0/2187185.0*r*s+10.0/7.0
1115 *t12*t1+t16*t17+12.0/3289.0*t9*s-84.0/3289.0*t1*t9+14.0/299.0*t9*t17-21.0/65.0*
1116 t4*t17+30.0/7.0*t10*t17+3.0/23023.0*r*t17-12.0/3289.0*t3*t17-378.0/16445.0*t4*s
1117 -42.0/65.0*t27*t1+84.0/65.0*t27*t17;
1118       t39 = -40.0/13.0*t6*t17+t16*s/21-3.0/7.0*t16*t1+126.0/1495.0*t27*s+8.0/
1119 1495.0*t6+12.0/16445.0*t4-2.0/16445.0*t9-t16/819-s/6561555+2.0/455.0*t12-t17
1120 /690690-180.0/91.0*t10*t1+2.0/2187185.0*t1-3.0/455.0*t10-r/2187185-42.0/16445.0
1121 *t27+48.0/23023.0*t3*t1+24.0/2187185.0*t3-36.0/115115.0*t3*s-45.0/14.0*t12*t17;
1122       f = t30+t39;
1123    /* p=15 10 2 */
1124    } else if( i==10 && j==2 ) {
1125       t1 = r*r;
1126       t2 = t1*t1;
1127       t3 = t2*t2;
1128       t4 = t3*r;
1129       t5 = t3*t1;
1130       t6 = s*s;
1131       t8 = t2*t1;
1132       t10 = r*t1;
1133       t11 = t2*t10;
1134       t15 = 1.0/60843510.0-55.0/1638.0*t4+t5*t6+44.0/13.0*t8*t6-t10/2093+11.0/
1135 299.0*t8+t2/299+t5/126+11.0/182.0*t3-s/3579030-550.0/91.0*t11*t6-11.0/182.0*t11
1136 -r/715806+t6/1193010+132.0/91.0*t11*s+r*s/39767;
1137       t18 = t2*r;
1138       t32 = 275.0/42.0*t3*t6-3.0/14.0*t5*s-77.0/65.0*t18*t6-55.0/14.0*t4*t6+
1139 77.0/299.0*t2*t6+462.0/1495.0*t18*s+55.0/63.0*t4*s-21.0/1495.0*t18+3.0/79534.0*
1140 t1-3.0/4186.0*t1*s+20.0/2093.0*t10*s-10.0/299.0*t10*t6-r*t6/12558-11.0/13.0*t8*
1141 s+5.0/2093.0*t6*t1-21.0/299.0*t2*s-275.0/182.0*t3*s;
1142       f = t15+t32;
1143    /* p=15 11 1 */
1144    } else if( i==11 && j==1 ) {
1145       t1 = r*r;
1146       t2 = r*t1;
1147       t3 = t1*t1;
1148       t4 = t3*t3;
1149       t5 = t4*t2;
1150       t7 = t4*t1;
1151       t9 = r*t4;
1152       t12 = t3*t2;
1153       t14 = t3*t1;
1154       t16 = t5*s-t5/14-33.0/7.0*t7*s+22.0/63.0*t7+605.0/63.0*t9*s-605.0/819.0*
1155 t9-3025.0/273.0*t4*s+242.0/273.0*t4+726.0/91.0*t12*s-121.0/182.0*t12-242.0/65.0
1156 *t14*s+484.0/1495.0*t14;
1157       t17 = t3*r;
1158       t23 = 1694.0/1495.0*t17*s-154.0/1495.0*t17-66.0/299.0*t3*s+44.0/2093.0*t3
1159 +55.0/2093.0*t2*s-11.0/4186.0*t2-11.0/6279.0*t1*s+22.0/119301.0*t1+11.0/
1160 198835.0*r*s-11.0/1789515.0*r-s/1789515+2.0/30421755.0;
1161       f = t16+t23;
1162    /* p=15 12 0 */
1163    } else if( i==12 && j==0 ) {
1164       t1 = r*r;
1165       t2 = t1*t1;
1166       t3 = t2*t2;
1167       t5 = r*t1;
1168       f = t3*t2-39.0/7.0*t3*t5+286.0/21.0*t3*t1-1210.0/63.0*t3*r+121.0/7.0*t3
1169 -363.0/35.0*t2*t5+484.0/115.0*t2*t1-132.0/115.0*t2*r+33.0/161.0*t2-11.0/483.0*
1170 t5+22.0/15295.0*t1-2.0/45885.0*r+1.0/2340135.0;
1171     }
1172     return f;
1173 }
1174 #ifdef __cplusplus
1175 }
1176 #endif
1177