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 return the "body" mode for a tetrahedron. i,j,k are the highest
11 monomial orders of r,s,t in the resulting polynomial order
12 NOTE: code generated by Maple.
13 -------------------------------------------------------------------------*/
14
15 #ifdef __cplusplus
16 extern "C" {
17 #endif
18
Bn(int i,int j,int k,double r,double s,double t)19 double Bn(int i, int j, int k, double r, double s, double t) {
20 double f = 0.0;
21 double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18;
22 double t19,t20,t21,t22,t23,t24,t25,t26,t27,t28,t29,t30,t31,t32,t33,t34;
23 double t35,t36,t37,t38,t39,t40,t41,t42,t43,t44,t45;
24
25 /* p=4 0 0 0 */
26 if( i==0 && j==0 && k==0 ) {
27 f = 1.0;
28 /* p=5 0 0 1 */
29 } else if( i==0 && j==0 && k==1 ) {
30 f = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;
31 /* p=5 0 1 0 */
32 } else if( i==0 && j==1 && k==0 ) {
33 f = s+2.0/5.0*r-2.0/5.0;
34 /* p=5 1 0 0 */
35 } else if( i==1 && j==0 && k==0 ) {
36 f = r-2.0/7.0;
37 /* p=6 0 0 2 */
38 } else if( i==0 && j==0 && k==2 ) {
39 t1 = t*t;
40 t2 = s+r-1.0;
41 t4 = t2*t2;
42 f = t1+6.0/5.0*t*t2+3.0/10.0*t4;
43 /* p=6 0 1 1 */
44 } else if( i==0 && j==1 && k==1 ) {
45 f = (s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
46 /* p=6 0 2 0 */
47 } else if( i==0 && j==2 && k==0 ) {
48 t1 = s*s;
49 t2 = r-1.0;
50 t4 = t2*t2;
51 f = t1+6.0/7.0*t2*s+t4/7;
52 /* p=6 1 0 1 */
53 } else if( i==1 && j==0 && k==1 ) {
54 f = (r-2.0/9.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
55 /* p=6 1 1 0 */
56 } else if( i==1 && j==1 && k==0 ) {
57 f = (r-2.0/9.0)*(s+2.0/5.0*r-2.0/5.0);
58 /* p=6 2 0 0 */
59 } else if( i==2 && j==0 && k==0 ) {
60 t1 = r*r;
61 f = t1-2.0/3.0*r+1.0/12.0;
62 /* p=7 0 0 3 */
63 } else if( i==0 && j==0 && k==3 ) {
64 t1 = t*t;
65 t3 = s+r-1.0;
66 t5 = t3*t3;
67 f = t1*t+12.0/7.0*t3*t1+6.0/7.0*t*t5+4.0/35.0*t3*t5;
68 /* p=7 0 1 2 */
69 } else if( i==0 && j==1 && k==2 ) {
70 t2 = t*t;
71 t3 = s+r-1.0;
72 t5 = t3*t3;
73 f = (s+2.0/9.0*r-2.0/9.0)*(t2+6.0/5.0*t*t3+3.0/10.0*t5);
74 /* p=7 0 2 1 */
75 } else if( i==0 && j==2 && k==1 ) {
76 t1 = s*s;
77 t2 = r-1.0;
78 t4 = t2*t2;
79 f = (t1+2.0/3.0*s*t2+t4/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
80 /* p=7 0 3 0 */
81 } else if( i==0 && j==3 && k==0 ) {
82 t1 = s*s;
83 t3 = r-1.0;
84 t5 = t3*t3;
85 f = t1*s+4.0/3.0*t1*t3+s*t5/2+t5*t3/21;
86 /* p=7 1 0 2 */
87 } else if( i==1 && j==0 && k==2 ) {
88 t2 = t*t;
89 t3 = s+r-1.0;
90 t5 = t3*t3;
91 f = (r-2.0/11.0)*(t2+6.0/5.0*t*t3+3.0/10.0*t5);
92 /* p=7 1 1 1 */
93 } else if( i==1 && j==1 && k==1 ) {
94 f = (r-2.0/11.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
95 /* p=7 1 2 0 */
96 } else if( i==1 && j==2 && k==0 ) {
97 t2 = s*s;
98 t3 = r-1.0;
99 t5 = t3*t3;
100 f = (r-2.0/11.0)*(t2+6.0/7.0*s*t3+t5/7);
101 /* p=7 2 0 1 */
102 } else if( i==2 && j==0 && k==1 ) {
103 t1 = r*r;
104 f = (t1-6.0/11.0*r+3.0/55.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
105 /* p=7 2 1 0 */
106 } else if( i==2 && j==1 && k==0 ) {
107 t1 = r*r;
108 f = (t1-6.0/11.0*r+3.0/55.0)*(s+2.0/5.0*r-2.0/5.0);
109 /* p=7 3 0 0 */
110 } else if( i==3 && j==0 && k==0 ) {
111 t1 = r*r;
112 f = t1*r-12.0/11.0*t1+18.0/55.0*r-4.0/165.0;
113 /* p=8 0 0 4 */
114 } else if( i==0 && j==0 && k==4 ) {
115 t1 = t*t;
116 t2 = t1*t1;
117 t4 = s+r-1.0;
118 t6 = t4*t4;
119 t10 = t6*t6;
120 f = t2+20.0/9.0*t1*t*t4+5.0/3.0*t1*t6+10.0/21.0*t*t6*t4+5.0/126.0*t10;
121 /* p=8 0 1 3 */
122 } else if( i==0 && j==1 && k==3 ) {
123 t2 = t*t;
124 t4 = s+r-1.0;
125 t6 = t4*t4;
126 f = (s+2.0/11.0*r-2.0/11.0)*(t2*t+12.0/7.0*t2*t4+6.0/7.0*t*t6+4.0/35.0*
127 t6*t4);
128 /* p=8 0 2 2 */
129 } else if( i==0 && j==2 && k==2 ) {
130 t1 = s*s;
131 t2 = r-1.0;
132 t4 = t2*t2;
133 t6 = t*t;
134 t7 = s+r-1.0;
135 t9 = t7*t7;
136 f = (t1+6.0/11.0*s*t2+3.0/55.0*t4)*(t6+6.0/5.0*t*t7+3.0/10.0*t9);
137 /* p=8 0 3 1 */
138 } else if( i==0 && j==3 && k==1 ) {
139 t1 = s*s;
140 t3 = r-1.0;
141 t5 = t3*t3;
142 f = (t1*s+12.0/11.0*t1*t3+18.0/55.0*s*t5+4.0/165.0*t5*t3)*(t+2.0/3.0*s+
143 2.0/3.0*r-2.0/3.0);
144 /* p=8 0 4 0 */
145 } else if( i==0 && j==4 && k==0 ) {
146 t1 = s*s;
147 t2 = t1*t1;
148 t4 = r-1.0;
149 t6 = t4*t4;
150 t10 = t6*t6;
151 f = t2+20.0/11.0*t1*s*t4+12.0/11.0*t1*t6+8.0/33.0*s*t6*t4+t10/66;
152 /* p=8 1 0 3 */
153 } else if( i==1 && j==0 && k==3 ) {
154 t2 = t*t;
155 t4 = s+r-1.0;
156 t6 = t4*t4;
157 f = (r-2.0/13.0)*(t2*t+12.0/7.0*t2*t4+6.0/7.0*t*t6+4.0/35.0*t6*t4);
158 /* p=8 1 1 2 */
159 } else if( i==1 && j==1 && k==2 ) {
160 t4 = t*t;
161 t5 = s+r-1.0;
162 t7 = t5*t5;
163 f = (r-2.0/13.0)*(s+2.0/9.0*r-2.0/9.0)*(t4+6.0/5.0*t*t5+3.0/10.0*t7);
164 /* p=8 1 2 1 */
165 } else if( i==1 && j==2 && k==1 ) {
166 t2 = s*s;
167 t3 = r-1.0;
168 t5 = t3*t3;
169 f = (r-2.0/13.0)*(t2+2.0/3.0*s*t3+t5/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
170 /* p=8 1 3 0 */
171 } else if( i==1 && j==3 && k==0 ) {
172 t2 = s*s;
173 t4 = r-1.0;
174 t6 = t4*t4;
175 f = (r-2.0/13.0)*(t2*s+4.0/3.0*t2*t4+s*t6/2+t6*t4/21);
176 /* p=8 2 0 2 */
177 } else if( i==2 && j==0 && k==2 ) {
178 t1 = r*r;
179 t3 = t*t;
180 t4 = s+r-1.0;
181 t6 = t4*t4;
182 f = (t1-6.0/13.0*r+1.0/26.0)*(t3+6.0/5.0*t*t4+3.0/10.0*t6);
183 /* p=8 2 1 1 */
184 } else if( i==2 && j==1 && k==1 ) {
185 t1 = r*r;
186 f = (t1-6.0/13.0*r+1.0/26.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*
187 r-2.0/3.0);
188 /* p=8 2 2 0 */
189 } else if( i==2 && j==2 && k==0 ) {
190 t1 = r*r;
191 t3 = s*s;
192 t4 = r-1.0;
193 t6 = t4*t4;
194 f = (t1-6.0/13.0*r+1.0/26.0)*(t3+6.0/7.0*s*t4+t6/7);
195 /* p=8 3 0 1 */
196 } else if( i==3 && j==0 && k==1 ) {
197 t1 = r*r;
198 f = (t1*r-12.0/13.0*t1+3.0/13.0*r-2.0/143.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/
199 3.0);
200 /* p=8 3 1 0 */
201 } else if( i==3 && j==1 && k==0 ) {
202 t1 = r*r;
203 f = (t1*r-12.0/13.0*t1+3.0/13.0*r-2.0/143.0)*(s+2.0/5.0*r-2.0/5.0);
204 /* p=8 4 0 0 */
205 } else if( i==4 && j==0 && k==0 ) {
206 t1 = r*r;
207 t2 = t1*t1;
208 f = t2-20.0/13.0*t1*r+10.0/13.0*t1-20.0/143.0*r+1.0/143.0;
209 /* p=9 0 0 5 */
210 } else if( i==0 && j==0 && k==5 ) {
211 t1 = t*t;
212 t2 = t1*t1;
213 t4 = s+r-1.0;
214 t7 = t4*t4;
215 t11 = t7*t7;
216 f = t2*t+30.0/11.0*t2*t4+30.0/11.0*t1*t*t7+40.0/33.0*t1*t7*t4+5.0/22.0*
217 t*t11+t11*t4/77;
218 /* p=9 0 1 4 */
219 } else if( i==0 && j==1 && k==4 ) {
220 t2 = t*t;
221 t3 = t2*t2;
222 t5 = s+r-1.0;
223 t7 = t5*t5;
224 t11 = t7*t7;
225 f = (s+2.0/13.0*r-2.0/13.0)*(t3+20.0/9.0*t2*t*t5+5.0/3.0*t2*t7+10.0/
226 21.0*t*t7*t5+5.0/126.0*t11);
227 /* p=9 0 2 3 */
228 } else if( i==0 && j==2 && k==3 ) {
229 t1 = s*s;
230 t2 = r-1.0;
231 t4 = t2*t2;
232 t6 = t*t;
233 t8 = s+r-1.0;
234 t10 = t8*t8;
235 f = (t1+6.0/13.0*t2*s+t4/26)*(t6*t+12.0/7.0*t6*t8+6.0/7.0*t*t10+4.0/
236 35.0*t10*t8);
237 /* p=9 0 3 2 */
238 } else if( i==0 && j==3 && k==2 ) {
239 t1 = s*s;
240 t3 = r-1.0;
241 t5 = t3*t3;
242 t9 = t*t;
243 t10 = s+r-1.0;
244 t12 = t10*t10;
245 f = (t1*s+12.0/13.0*t1*t3+3.0/13.0*s*t5+2.0/143.0*t5*t3)*(t9+6.0/5.0*t*
246 t10+3.0/10.0*t12);
247 /* p=9 0 4 1 */
248 } else if( i==0 && j==4 && k==1 ) {
249 t1 = s*s;
250 t2 = t1*t1;
251 t4 = r-1.0;
252 t6 = t4*t4;
253 t10 = t6*t6;
254 f = (t2+20.0/13.0*t1*s*t4+10.0/13.0*t1*t6+20.0/143.0*s*t6*t4+t10/143)*(
255 t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
256 /* p=9 0 5 0 */
257 } else if( i==0 && j==5 && k==0 ) {
258 t1 = s*s;
259 t2 = t1*t1;
260 t4 = r-1.0;
261 t7 = t4*t4;
262 t11 = t7*t7;
263 f = t2*s+30.0/13.0*t2*t4+25.0/13.0*t1*s*t7+100.0/143.0*t1*t7*t4+15.0/
264 143.0*s*t11+2.0/429.0*t11*t4;
265 /* p=9 1 0 4 */
266 } else if( i==1 && j==0 && k==4 ) {
267 t2 = t*t;
268 t3 = t2*t2;
269 t5 = s+r-1.0;
270 t7 = t5*t5;
271 t11 = t7*t7;
272 f = (r-2.0/15.0)*(t3+20.0/9.0*t2*t*t5+5.0/3.0*t2*t7+10.0/21.0*t*t7*t5+
273 5.0/126.0*t11);
274 /* p=9 1 1 3 */
275 } else if( i==1 && j==1 && k==3 ) {
276 t4 = t*t;
277 t6 = s+r-1.0;
278 t8 = t6*t6;
279 f = (r-2.0/15.0)*(s+2.0/11.0*r-2.0/11.0)*(t4*t+12.0/7.0*t4*t6+6.0/7.0*t
280 *t8+4.0/35.0*t8*t6);
281 /* p=9 1 2 2 */
282 } else if( i==1 && j==2 && k==2 ) {
283 t2 = s*s;
284 t3 = r-1.0;
285 t5 = t3*t3;
286 t8 = t*t;
287 t9 = s+r-1.0;
288 t11 = t9*t9;
289 f = (r-2.0/15.0)*(t2+6.0/11.0*s*t3+3.0/55.0*t5)*(t8+6.0/5.0*t*t9+3.0/
290 10.0*t11);
291 /* p=9 1 3 1 */
292 } else if( i==1 && j==3 && k==1 ) {
293 t2 = s*s;
294 t4 = r-1.0;
295 t6 = t4*t4;
296 f = (r-2.0/15.0)*(t2*s+12.0/11.0*t2*t4+18.0/55.0*s*t6+4.0/165.0*t4*t6)*
297 (t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
298 /* p=9 1 4 0 */
299 } else if( i==1 && j==4 && k==0 ) {
300 t2 = s*s;
301 t3 = t2*t2;
302 t5 = r-1.0;
303 t7 = t5*t5;
304 t11 = t7*t7;
305 f = (r-2.0/15.0)*(t3+20.0/11.0*t2*s*t5+12.0/11.0*t2*t7+8.0/33.0*s*t7*t5
306 +t11/66);
307 /* p=9 2 0 3 */
308 } else if( i==2 && j==0 && k==3 ) {
309 t1 = r*r;
310 t3 = t*t;
311 t5 = s+r-1.0;
312 t7 = t5*t5;
313 f = (t1-2.0/5.0*r+1.0/35.0)*(t3*t+12.0/7.0*t3*t5+6.0/7.0*t*t7+4.0/35.0*
314 t7*t5);
315 /* p=9 2 1 2 */
316 } else if( i==2 && j==1 && k==2 ) {
317 t1 = r*r;
318 t5 = t*t;
319 t6 = s+r-1.0;
320 t8 = t6*t6;
321 f = (t1-2.0/5.0*r+1.0/35.0)*(s+2.0/9.0*r-2.0/9.0)*(t5+6.0/5.0*t*t6+3.0/
322 10.0*t8);
323 /* p=9 2 2 1 */
324 } else if( i==2 && j==2 && k==1 ) {
325 t1 = r*r;
326 t3 = s*s;
327 t4 = r-1.0;
328 t6 = t4*t4;
329 f = (t1-2.0/5.0*r+1.0/35.0)*(t3+2.0/3.0*s*t4+t6/12)*(t+2.0/3.0*s+2.0/
330 3.0*r-2.0/3.0);
331 /* p=9 2 3 0 */
332 } else if( i==2 && j==3 && k==0 ) {
333 t1 = r*r;
334 t3 = s*s;
335 t5 = r-1.0;
336 t7 = t5*t5;
337 f = (t1-2.0/5.0*r+1.0/35.0)*(s*t3+4.0/3.0*t3*t5+s*t7/2+t7*t5/21);
338 /* p=9 3 0 2 */
339 } else if( i==3 && j==0 && k==2 ) {
340 t1 = r*r;
341 t4 = t*t;
342 t5 = s+r-1.0;
343 t7 = t5*t5;
344 f = (t1*r-4.0/5.0*t1+6.0/35.0*r-4.0/455.0)*(t4+6.0/5.0*t*t5+3.0/10.0*t7)
345 ;
346 /* p=9 3 1 1 */
347 } else if( i==3 && j==1 && k==1 ) {
348 t1 = r*r;
349 f = (t1*r-4.0/5.0*t1+6.0/35.0*r-4.0/455.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/
350 3.0*s+2.0/3.0*r-2.0/3.0);
351 /* p=9 3 2 0 */
352 } else if( i==3 && j==2 && k==0 ) {
353 t1 = r*r;
354 t4 = s*s;
355 t5 = r-1.0;
356 t7 = t5*t5;
357 f = (t1*r-4.0/5.0*t1+6.0/35.0*r-4.0/455.0)*(t4+6.0/7.0*s*t5+t7/7);
358 /* p=9 4 0 1 */
359 } else if( i==4 && j==0 && k==1 ) {
360 t1 = r*r;
361 t2 = t1*t1;
362 f = (t2-4.0/3.0*t1*r+4.0/7.0*t1-8.0/91.0*r+1.0/273.0)*(t+2.0/3.0*s+2.0/
363 3.0*r-2.0/3.0);
364 /* p=9 4 1 0 */
365 } else if( i==4 && j==1 && k==0 ) {
366 t1 = r*r;
367 t2 = t1*t1;
368 f = (t2-4.0/3.0*t1*r+4.0/7.0*t1-8.0/91.0*r+1.0/273.0)*(s+2.0/5.0*r-2.0/
369 5.0);
370 /* p=9 5 0 0 */
371 } else if( i==5 && j==0 && k==0 ) {
372 t1 = r*r;
373 t2 = t1*t1;
374 f = t2*r-2.0*t2+10.0/7.0*t1*r-40.0/91.0*t1+5.0/91.0*r-2.0/1001.0;
375 /* p=10 0 0 6 */
376 } else if( i==0 && j==0 && k==6 ) {
377 t1 = t*t;
378 t2 = t1*t1;
379 t5 = s+r-1.0;
380 t7 = t5*t5;
381 t12 = t7*t7;
382 f = t2*t1+42.0/13.0*t2*t*t5+105.0/26.0*t2*t7+350.0/143.0*t1*t*t7*t5+
383 105.0/143.0*t1*t12+14.0/143.0*t*t12*t5+7.0/1716.0*t12*t7;
384 /* p=10 0 1 5 */
385 } else if( i==0 && j==1 && k==5 ) {
386 t2 = t*t;
387 t3 = t2*t2;
388 t5 = s+r-1.0;
389 t8 = t5*t5;
390 t12 = t8*t8;
391 f = (s+2.0/15.0*r-2.0/15.0)*(t3*t+30.0/11.0*t3*t5+30.0/11.0*t2*t*t8+
392 40.0/33.0*t2*t8*t5+5.0/22.0*t*t12+t12*t5/77);
393 /* p=10 0 2 4 */
394 } else if( i==0 && j==2 && k==4 ) {
395 t1 = s*s;
396 t2 = r-1.0;
397 t4 = t2*t2;
398 t6 = t*t;
399 t7 = t6*t6;
400 t9 = s+r-1.0;
401 t11 = t9*t9;
402 t15 = t11*t11;
403 f = (t1+2.0/5.0*s*t2+t4/35)*(t7+20.0/9.0*t6*t*t9+5.0/3.0*t6*t11+10.0/
404 21.0*t*t11*t9+5.0/126.0*t15);
405 /* p=10 0 3 3 */
406 } else if( i==0 && j==3 && k==3 ) {
407 t1 = s*s;
408 t3 = r-1.0;
409 t5 = t3*t3;
410 t9 = t*t;
411 t11 = s+r-1.0;
412 t13 = t11*t11;
413 f = (t1*s+4.0/5.0*t1*t3+6.0/35.0*s*t5+4.0/455.0*t3*t5)*(t9*t+12.0/7.0*
414 t11*t9+6.0/7.0*t*t13+4.0/35.0*t13*t11);
415 /* p=10 0 4 2 */
416 } else if( i==0 && j==4 && k==2 ) {
417 t1 = s*s;
418 t2 = t1*t1;
419 t4 = r-1.0;
420 t6 = t4*t4;
421 t10 = t6*t6;
422 t12 = t*t;
423 t13 = s+r-1.0;
424 t15 = t13*t13;
425 f = (t2+4.0/3.0*t1*s*t4+4.0/7.0*t1*t6+8.0/91.0*s*t6*t4+t10/273)*(t12+
426 6.0/5.0*t*t13+3.0/10.0*t15);
427 /* p=10 0 5 1 */
428 } else if( i==0 && j==5 && k==1 ) {
429 t1 = s*s;
430 t2 = t1*t1;
431 t4 = r-1.0;
432 t7 = t4*t4;
433 t11 = t7*t7;
434 f = (s*t2+2.0*t2*t4+10.0/7.0*t1*s*t7+40.0/91.0*t1*t7*t4+5.0/91.0*s*t11+
435 2.0/1001.0*t11*t4)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
436 /* p=10 0 6 0 */
437 } else if( i==0 && j==6 && k==0 ) {
438 t1 = s*s;
439 t2 = t1*t1;
440 t5 = r-1.0;
441 t7 = t5*t5;
442 t12 = t7*t7;
443 f = t2*t1+14.0/5.0*t2*s*t5+3.0*t2*t7+20.0/13.0*t1*s*t7*t5+5.0/13.0*t1*
444 t12+6.0/143.0*s*t12*t5+t12*t7/715;
445 /* p=10 1 0 5 */
446 } else if( i==1 && j==0 && k==5 ) {
447 t2 = t*t;
448 t3 = t2*t2;
449 t5 = s+r-1.0;
450 t8 = t5*t5;
451 t12 = t8*t8;
452 f = (r-2.0/17.0)*(t3*t+30.0/11.0*t3*t5+30.0/11.0*t2*t*t8+40.0/33.0*t2*
453 t8*t5+5.0/22.0*t*t12+t12*t5/77);
454 /* p=10 1 1 4 */
455 } else if( i==1 && j==1 && k==4 ) {
456 t4 = t*t;
457 t5 = t4*t4;
458 t7 = s+r-1.0;
459 t9 = t7*t7;
460 t13 = t9*t9;
461 f = (r-2.0/17.0)*(s+2.0/13.0*r-2.0/13.0)*(t5+20.0/9.0*t4*t*t7+5.0/3.0*
462 t4*t9+10.0/21.0*t*t9*t7+5.0/126.0*t13);
463 /* p=10 1 2 3 */
464 } else if( i==1 && j==2 && k==3 ) {
465 t2 = s*s;
466 t3 = r-1.0;
467 t5 = t3*t3;
468 t8 = t*t;
469 t10 = s+r-1.0;
470 t12 = t10*t10;
471 f = (r-2.0/17.0)*(t2+6.0/13.0*s*t3+t5/26)*(t8*t+12.0/7.0*t8*t10+6.0/7.0
472 *t*t12+4.0/35.0*t12*t10);
473 /* p=10 1 3 2 */
474 } else if( i==1 && j==3 && k==2 ) {
475 t2 = s*s;
476 t4 = r-1.0;
477 t6 = t4*t4;
478 t11 = t*t;
479 t12 = s+r-1.0;
480 t14 = t12*t12;
481 f = (r-2.0/17.0)*(t2*s+12.0/13.0*t2*t4+3.0/13.0*s*t6+2.0/143.0*t6*t4)*(
482 t11+6.0/5.0*t*t12+3.0/10.0*t14);
483 /* p=10 1 4 1 */
484 } else if( i==1 && j==4 && k==1 ) {
485 t2 = s*s;
486 t3 = t2*t2;
487 t5 = r-1.0;
488 t7 = t5*t5;
489 t11 = t7*t7;
490 f = (r-2.0/17.0)*(t3+20.0/13.0*t2*s*t5+10.0/13.0*t2*t7+20.0/143.0*s*t7*
491 t5+t11/143)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
492 /* p=10 1 5 0 */
493 } else if( i==1 && j==5 && k==0 ) {
494 t2 = s*s;
495 t3 = t2*t2;
496 t5 = r-1.0;
497 t8 = t5*t5;
498 t12 = t8*t8;
499 f = (r-2.0/17.0)*(s*t3+30.0/13.0*t3*t5+25.0/13.0*t2*s*t8+100.0/143.0*t2
500 *t8*t5+15.0/143.0*s*t12+2.0/429.0*t12*t5);
501 /* p=10 2 0 4 */
502 } else if( i==2 && j==0 && k==4 ) {
503 t1 = r*r;
504 t3 = t*t;
505 t4 = t3*t3;
506 t6 = s+r-1.0;
507 t8 = t6*t6;
508 t12 = t8*t8;
509 f = (t1-6.0/17.0*r+3.0/136.0)*(t4+20.0/9.0*t3*t*t6+5.0/3.0*t3*t8+10.0/
510 21.0*t*t8*t6+5.0/126.0*t12);
511 /* p=10 2 1 3 */
512 } else if( i==2 && j==1 && k==3 ) {
513 t1 = r*r;
514 t5 = t*t;
515 t7 = s+r-1.0;
516 t9 = t7*t7;
517 f = (t1-6.0/17.0*r+3.0/136.0)*(s+2.0/11.0*r-2.0/11.0)*(t5*t+12.0/7.0*t7
518 *t5+6.0/7.0*t*t9+4.0/35.0*t9*t7);
519 /* p=10 2 2 2 */
520 } else if( i==2 && j==2 && k==2 ) {
521 t1 = r*r;
522 t3 = s*s;
523 t4 = r-1.0;
524 t6 = t4*t4;
525 t9 = t*t;
526 t10 = s+r-1.0;
527 t12 = t10*t10;
528 f = (t1-6.0/17.0*r+3.0/136.0)*(t3+6.0/11.0*s*t4+3.0/55.0*t6)*(t9+6.0/
529 5.0*t*t10+3.0/10.0*t12);
530 /* p=10 2 3 1 */
531 } else if( i==2 && j==3 && k==1 ) {
532 t1 = r*r;
533 t3 = s*s;
534 t5 = r-1.0;
535 t7 = t5*t5;
536 f = (t1-6.0/17.0*r+3.0/136.0)*(s*t3+12.0/11.0*t3*t5+18.0/55.0*s*t7+4.0/
537 165.0*t7*t5)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
538 /* p=10 2 4 0 */
539 } else if( i==2 && j==4 && k==0 ) {
540 t1 = r*r;
541 t3 = s*s;
542 t4 = t3*t3;
543 t6 = r-1.0;
544 t8 = t6*t6;
545 t12 = t8*t8;
546 f = (t1-6.0/17.0*r+3.0/136.0)*(t4+20.0/11.0*t3*s*t6+12.0/11.0*t3*t8+8.0
547 /33.0*s*t8*t6+t12/66);
548 /* p=10 3 0 3 */
549 } else if( i==3 && j==0 && k==3 ) {
550 t1 = r*r;
551 t4 = t*t;
552 t6 = s+r-1.0;
553 t8 = t6*t6;
554 f = (t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0)*(t4*t+12.0/7.0*t4*t6+6.0/
555 7.0*t*t8+4.0/35.0*t8*t6);
556 /* p=10 3 1 2 */
557 } else if( i==3 && j==1 && k==2 ) {
558 t1 = r*r;
559 t6 = t*t;
560 t7 = s+r-1.0;
561 t9 = t7*t7;
562 f = (t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0)*(s+2.0/9.0*r-2.0/9.0)*(t6+
563 6.0/5.0*t*t7+3.0/10.0*t9);
564 /* p=10 3 2 1 */
565 } else if( i==3 && j==2 && k==1 ) {
566 t1 = r*r;
567 t4 = s*s;
568 t5 = r-1.0;
569 t7 = t5*t5;
570 f = (t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0)*(t4+2.0/3.0*s*t5+t7/12)*(t
571 +2.0/3.0*s+2.0/3.0*r-2.0/3.0);
572 /* p=10 3 3 0 */
573 } else if( i==3 && j==3 && k==0 ) {
574 t1 = r*r;
575 t4 = s*s;
576 t6 = r-1.0;
577 t8 = t6*t6;
578 f = (t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0)*(t4*s+4.0/3.0*t4*t6+s*t8/2
579 +t8*t6/21);
580 /* p=10 4 0 2 */
581 } else if( i==4 && j==0 && k==2 ) {
582 t1 = r*r;
583 t2 = t1*t1;
584 t5 = t*t;
585 t6 = s+r-1.0;
586 t8 = t6*t6;
587 f = (t2-20.0/17.0*t1*r+15.0/34.0*t1-r/17+1.0/476.0)*(t5+6.0/5.0*t*t6+
588 3.0/10.0*t8);
589 /* p=10 4 1 1 */
590 } else if( i==4 && j==1 && k==1 ) {
591 t1 = r*r;
592 t2 = t1*t1;
593 f = (t2-20.0/17.0*t1*r+15.0/34.0*t1-r/17+1.0/476.0)*(s+2.0/7.0*r-2.0/7.0
594 )*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
595 /* p=10 4 2 0 */
596 } else if( i==4 && j==2 && k==0 ) {
597 t1 = r*r;
598 t2 = t1*t1;
599 t5 = s*s;
600 t6 = r-1.0;
601 t8 = t6*t6;
602 f = (t2-20.0/17.0*t1*r+15.0/34.0*t1-r/17+1.0/476.0)*(t5+6.0/7.0*s*t6+t8
603 /7);
604 /* p=10 5 0 1 */
605 } else if( i==5 && j==0 && k==1 ) {
606 t1 = r*r;
607 t2 = t1*t1;
608 f = (t2*r-30.0/17.0*t2+75.0/68.0*t1*r-5.0/17.0*t1+15.0/476.0*r-3.0/
609 3094.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
610 /* p=10 5 1 0 */
611 } else if( i==5 && j==1 && k==0 ) {
612 t1 = r*r;
613 t2 = t1*t1;
614 f = (t2*r-30.0/17.0*t2+75.0/68.0*t1*r-5.0/17.0*t1+15.0/476.0*r-3.0/
615 3094.0)*(s+2.0/5.0*r-2.0/5.0);
616 /* p=10 6 0 0 */
617 } else if( i==6 && j==0 && k==0 ) {
618 t1 = r*r;
619 t2 = t1*t1;
620 f = t2*t1-42.0/17.0*t2*r+315.0/136.0*t2-35.0/34.0*t1*r+15.0/68.0*t1-9.0/
621 442.0*r+1.0/1768.0;
622 /* p=11 0 0 7 */
623 } else if( i==0 && j==0 && k==7 ) {
624 t1 = t*t;
625 t2 = t1*t;
626 t3 = t1*t1;
627 t6 = s+r-1.0;
628 t9 = t6*t6;
629 t11 = t9*t6;
630 t13 = t9*t9;
631 f = t3*t2+56.0/15.0*t3*t1*t6+28.0/5.0*t3*t*t9+56.0/13.0*t3*t11+70.0/
632 39.0*t2*t13+56.0/143.0*t1*t13*t6+28.0/715.0*t*t13*t9+8.0/6435.0*t13*t11;
633 /* p=11 0 1 6 */
634 } else if( i==0 && j==1 && k==6 ) {
635 t2 = t*t;
636 t3 = t2*t2;
637 t6 = s+r-1.0;
638 t8 = t6*t6;
639 t13 = t8*t8;
640 f = (s+2.0/17.0*r-2.0/17.0)*(t3*t2+42.0/13.0*t3*t*t6+105.0/26.0*t3*t8+
641 350.0/143.0*t2*t*t8*t6+105.0/143.0*t2*t13+14.0/143.0*t*t13*t6+7.0/1716.0*t13*t8
642 );
643 /* p=11 0 2 5 */
644 } else if( i==0 && j==2 && k==5 ) {
645 t1 = s*s;
646 t2 = r-1.0;
647 t4 = t2*t2;
648 t6 = t*t;
649 t7 = t6*t6;
650 t9 = s+r-1.0;
651 t12 = t9*t9;
652 t16 = t12*t12;
653 f = (t1+6.0/17.0*s*t2+3.0/136.0*t4)*(t7*t+30.0/11.0*t7*t9+30.0/11.0*t6*
654 t*t12+40.0/33.0*t6*t12*t9+5.0/22.0*t*t16+t16*t9/77);
655 /* p=11 0 3 4 */
656 } else if( i==0 && j==3 && k==4 ) {
657 t1 = s*s;
658 t3 = r-1.0;
659 t5 = t3*t3;
660 t9 = t*t;
661 t10 = t9*t9;
662 t12 = s+r-1.0;
663 t14 = t12*t12;
664 t18 = t14*t14;
665 f = (t1*s+12.0/17.0*t3*t1+9.0/68.0*s*t5+t5*t3/170)*(t10+20.0/9.0*t9*t*
666 t12+5.0/3.0*t9*t14+10.0/21.0*t*t14*t12+5.0/126.0*t18);
667 /* p=11 0 4 3 */
668 } else if( i==0 && j==4 && k==3 ) {
669 t1 = s*s;
670 t2 = t1*t1;
671 t4 = r-1.0;
672 t6 = t4*t4;
673 t10 = t6*t6;
674 t12 = t*t;
675 t14 = s+r-1.0;
676 t16 = t14*t14;
677 f = (t2+20.0/17.0*t1*s*t4+15.0/34.0*t1*t6+s*t6*t4/17+t10/476)*(t12*t+
678 12.0/7.0*t14*t12+6.0/7.0*t*t16+4.0/35.0*t16*t14);
679 /* p=11 0 5 2 */
680 } else if( i==0 && j==5 && k==2 ) {
681 t1 = s*s;
682 t2 = t1*t1;
683 t4 = r-1.0;
684 t7 = t4*t4;
685 t11 = t7*t7;
686 t15 = t*t;
687 t16 = s+r-1.0;
688 t18 = t16*t16;
689 f = (s*t2+30.0/17.0*t2*t4+75.0/68.0*t1*s*t7+5.0/17.0*t1*t7*t4+15.0/
690 476.0*s*t11+3.0/3094.0*t11*t4)*(t15+6.0/5.0*t*t16+3.0/10.0*t18);
691 /* p=11 0 6 1 */
692 } else if( i==0 && j==6 && k==1 ) {
693 t1 = s*s;
694 t2 = t1*t1;
695 t5 = r-1.0;
696 t7 = t5*t5;
697 t12 = t7*t7;
698 f = (t2*t1+42.0/17.0*t2*s*t5+315.0/136.0*t2*t7+35.0/34.0*t1*s*t7*t5+
699 15.0/68.0*t1*t12+9.0/442.0*s*t12*t5+t12*t7/1768)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0
700 );
701 /* p=11 0 7 0 */
702 } else if( i==0 && j==7 && k==0 ) {
703 t1 = s*s;
704 t2 = t1*s;
705 t3 = t1*t1;
706 t6 = r-1.0;
707 t9 = t6*t6;
708 t11 = t9*t6;
709 t13 = t9*t9;
710 f = t3*t2+56.0/17.0*t3*t1*t6+147.0/34.0*t3*s*t9+49.0/17.0*t3*t11+35.0/
711 34.0*t2*t13+42.0/221.0*t1*t13*t6+7.0/442.0*s*t13*t9+t13*t11/2431;
712 /* p=11 1 0 6 */
713 } else if( i==1 && j==0 && k==6 ) {
714 t2 = t*t;
715 t3 = t2*t2;
716 t6 = s+r-1.0;
717 t8 = t6*t6;
718 t13 = t8*t8;
719 f = (r-2.0/19.0)*(t3*t2+42.0/13.0*t3*t*t6+105.0/26.0*t3*t8+350.0/143.0*
720 t2*t*t8*t6+105.0/143.0*t2*t13+14.0/143.0*t*t13*t6+7.0/1716.0*t13*t8);
721 /* p=11 1 1 5 */
722 } else if( i==1 && j==1 && k==5 ) {
723 t4 = t*t;
724 t5 = t4*t4;
725 t7 = s+r-1.0;
726 t10 = t7*t7;
727 t14 = t10*t10;
728 f = (r-2.0/19.0)*(s+2.0/15.0*r-2.0/15.0)*(t5*t+30.0/11.0*t5*t7+30.0/
729 11.0*t4*t*t10+40.0/33.0*t4*t10*t7+5.0/22.0*t*t14+t14*t7/77);
730 /* p=11 1 2 4 */
731 } else if( i==1 && j==2 && k==4 ) {
732 t2 = s*s;
733 t3 = r-1.0;
734 t5 = t3*t3;
735 t8 = t*t;
736 t9 = t8*t8;
737 t11 = s+r-1.0;
738 t13 = t11*t11;
739 t17 = t13*t13;
740 f = (r-2.0/19.0)*(t2+2.0/5.0*s*t3+t5/35)*(t9+20.0/9.0*t8*t*t11+5.0/3.0*
741 t13*t8+10.0/21.0*t*t13*t11+5.0/126.0*t17);
742 /* p=11 1 3 3 */
743 } else if( i==1 && j==3 && k==3 ) {
744 t2 = s*s;
745 t4 = r-1.0;
746 t6 = t4*t4;
747 t11 = t*t;
748 t13 = s+r-1.0;
749 t15 = t13*t13;
750 f = (r-2.0/19.0)*(t2*s+4.0/5.0*t2*t4+6.0/35.0*s*t6+4.0/455.0*t6*t4)*(
751 t11*t+12.0/7.0*t13*t11+6.0/7.0*t*t15+4.0/35.0*t15*t13);
752 /* p=11 1 4 2 */
753 } else if( i==1 && j==4 && k==2 ) {
754 t2 = s*s;
755 t3 = t2*t2;
756 t5 = r-1.0;
757 t7 = t5*t5;
758 t11 = t7*t7;
759 t14 = t*t;
760 t15 = s+r-1.0;
761 t17 = t15*t15;
762 f = (r-2.0/19.0)*(t3+4.0/3.0*t2*s*t5+4.0/7.0*t2*t7+8.0/91.0*s*t5*t7+t11
763 /273)*(t14+6.0/5.0*t*t15+3.0/10.0*t17);
764 /* p=11 1 5 1 */
765 } else if( i==1 && j==5 && k==1 ) {
766 t2 = s*s;
767 t3 = t2*t2;
768 t5 = r-1.0;
769 t8 = t5*t5;
770 t12 = t8*t8;
771 f = (r-2.0/19.0)*(s*t3+2.0*t3*t5+10.0/7.0*t2*s*t8+40.0/91.0*t2*t8*t5+
772 5.0/91.0*s*t12+2.0/1001.0*t12*t5)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
773 /* p=11 1 6 0 */
774 } else if( i==1 && j==6 && k==0 ) {
775 t2 = s*s;
776 t3 = t2*t2;
777 t6 = r-1.0;
778 t8 = t6*t6;
779 t13 = t8*t8;
780 f = (r-2.0/19.0)*(t3*t2+14.0/5.0*t3*s*t6+3.0*t3*t8+20.0/13.0*t2*s*t8*t6
781 +5.0/13.0*t2*t13+6.0/143.0*s*t13*t6+t13*t8/715);
782 /* p=11 2 0 5 */
783 } else if( i==2 && j==0 && k==5 ) {
784 t1 = r*r;
785 t3 = t*t;
786 t4 = t3*t3;
787 t6 = s+r-1.0;
788 t9 = t6*t6;
789 t13 = t9*t9;
790 f = (t1-6.0/19.0*r+1.0/57.0)*(t4*t+30.0/11.0*t6*t4+30.0/11.0*t3*t*t9+
791 40.0/33.0*t3*t9*t6+5.0/22.0*t*t13+t13*t6/77);
792 /* p=11 2 1 4 */
793 } else if( i==2 && j==1 && k==4 ) {
794 t1 = r*r;
795 t5 = t*t;
796 t6 = t5*t5;
797 t8 = s+r-1.0;
798 t10 = t8*t8;
799 t14 = t10*t10;
800 f = (s+2.0/13.0*r-2.0/13.0)*(t1-6.0/19.0*r+1.0/57.0)*(t6+20.0/9.0*t5*t*
801 t8+5.0/3.0*t5*t10+10.0/21.0*t*t10*t8+5.0/126.0*t14);
802 /* p=11 2 2 3 */
803 } else if( i==2 && j==2 && k==3 ) {
804 t1 = r*r;
805 t3 = s*s;
806 t4 = r-1.0;
807 t6 = t4*t4;
808 t9 = t*t;
809 t11 = s+r-1.0;
810 t13 = t11*t11;
811 f = (t1-6.0/19.0*r+1.0/57.0)*(t3+6.0/13.0*s*t4+t6/26)*(t9*t+12.0/7.0*t9
812 *t11+6.0/7.0*t*t13+4.0/35.0*t13*t11);
813 /* p=11 2 3 2 */
814 } else if( i==2 && j==3 && k==2 ) {
815 t1 = r*r;
816 t3 = s*s;
817 t5 = r-1.0;
818 t7 = t5*t5;
819 t12 = t*t;
820 t13 = s+r-1.0;
821 t15 = t13*t13;
822 f = (t1-6.0/19.0*r+1.0/57.0)*(t3*s+12.0/13.0*t3*t5+3.0/13.0*s*t7+2.0/
823 143.0*t7*t5)*(t12+6.0/5.0*t*t13+3.0/10.0*t15);
824 /* p=11 2 4 1 */
825 } else if( i==2 && j==4 && k==1 ) {
826 t1 = r*r;
827 t3 = s*s;
828 t4 = t3*t3;
829 t6 = r-1.0;
830 t8 = t6*t6;
831 t12 = t8*t8;
832 f = (t1-6.0/19.0*r+1.0/57.0)*(t4+20.0/13.0*t3*s*t6+10.0/13.0*t3*t8+20.0
833 /143.0*s*t8*t6+t12/143)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
834 /* p=11 2 5 0 */
835 } else if( i==2 && j==5 && k==0 ) {
836 t1 = r*r;
837 t3 = s*s;
838 t4 = t3*t3;
839 t6 = r-1.0;
840 t9 = t6*t6;
841 t13 = t9*t9;
842 f = (t1-6.0/19.0*r+1.0/57.0)*(s*t4+30.0/13.0*t4*t6+25.0/13.0*t3*s*t9+
843 100.0/143.0*t3*t9*t6+15.0/143.0*s*t13+2.0/429.0*t13*t6);
844 /* p=11 3 0 4 */
845 } else if( i==3 && j==0 && k==4 ) {
846 t1 = r*r;
847 t4 = t*t;
848 t5 = t4*t4;
849 t7 = s+r-1.0;
850 t9 = t7*t7;
851 t13 = t9*t9;
852 f = (t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0)*(t5+20.0/9.0*t4*t*t7+5.0/
853 3.0*t4*t9+10.0/21.0*t*t9*t7+5.0/126.0*t13);
854 /* p=11 3 1 3 */
855 } else if( i==3 && j==1 && k==3 ) {
856 t1 = r*r;
857 t6 = t*t;
858 t8 = s+r-1.0;
859 t10 = t8*t8;
860 f = (t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0)*(s+2.0/11.0*r-2.0/11.0)*(
861 t6*t+12.0/7.0*t8*t6+6.0/7.0*t*t10+4.0/35.0*t10*t8);
862 /* p=11 3 2 2 */
863 } else if( i==3 && j==2 && k==2 ) {
864 t1 = r*r;
865 t4 = s*s;
866 t5 = r-1.0;
867 t7 = t5*t5;
868 t10 = t*t;
869 t11 = s+r-1.0;
870 t13 = t11*t11;
871 f = (t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0)*(t4+6.0/11.0*s*t5+3.0/55.0
872 *t7)*(t10+6.0/5.0*t*t11+3.0/10.0*t13);
873 /* p=11 3 3 1 */
874 } else if( i==3 && j==3 && k==1 ) {
875 t1 = r*r;
876 t4 = s*s;
877 t6 = r-1.0;
878 t8 = t6*t6;
879 f = (t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0)*(s*t4+12.0/11.0*t4*t6+18.0
880 /55.0*s*t8+4.0/165.0*t8*t6)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
881 /* p=11 3 4 0 */
882 } else if( i==3 && j==4 && k==0 ) {
883 t1 = r*r;
884 t4 = s*s;
885 t5 = t4*t4;
886 t7 = r-1.0;
887 t9 = t7*t7;
888 t13 = t9*t9;
889 f = (t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0)*(t5+20.0/11.0*t4*s*t7+12.0
890 /11.0*t4*t9+8.0/33.0*s*t9*t7+t13/66);
891 /* p=11 4 0 3 */
892 } else if( i==4 && j==0 && k==3 ) {
893 t1 = r*r;
894 t2 = t1*t1;
895 t5 = t*t;
896 t7 = s+r-1.0;
897 t9 = t7*t7;
898 f = (t2-20.0/19.0*t1*r+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0)*(t5*t+12.0
899 /7.0*t7*t5+6.0/7.0*t9*t+4.0/35.0*t9*t7);
900 /* p=11 4 1 2 */
901 } else if( i==4 && j==1 && k==2 ) {
902 t1 = r*r;
903 t2 = t1*t1;
904 t7 = t*t;
905 t8 = s+r-1.0;
906 t10 = t8*t8;
907 f = (t2-20.0/19.0*t1*r+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0)*(s+2.0/9.0
908 *r-2.0/9.0)*(t7+6.0/5.0*t*t8+3.0/10.0*t10);
909 /* p=11 4 2 1 */
910 } else if( i==4 && j==2 && k==1 ) {
911 t1 = r*r;
912 t2 = t1*t1;
913 t5 = s*s;
914 t6 = r-1.0;
915 t8 = t6*t6;
916 f = (t2-20.0/19.0*t1*r+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0)*(t5+2.0/
917 3.0*s*t6+t8/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
918 /* p=11 4 3 0 */
919 } else if( i==4 && j==3 && k==0 ) {
920 t1 = r*r;
921 t2 = t1*t1;
922 t5 = s*s;
923 t7 = r-1.0;
924 t9 = t7*t7;
925 f = (t2-20.0/19.0*t1*r+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0)*(t5*s+4.0/
926 3.0*t7*t5+s*t9/2+t9*t7/21);
927 /* p=11 5 0 2 */
928 } else if( i==5 && j==0 && k==2 ) {
929 t1 = r*r;
930 t2 = t1*t1;
931 t6 = t*t;
932 t7 = s+r-1.0;
933 t9 = t7*t7;
934 f = (t2*r-30.0/19.0*t2+50.0/57.0*t1*r-200.0/969.0*t1+25.0/1292.0*r-1.0/
935 1938.0)*(t6+6.0/5.0*t*t7+3.0/10.0*t9);
936 /* p=11 5 1 1 */
937 } else if( i==5 && j==1 && k==1 ) {
938 t1 = r*r;
939 t2 = t1*t1;
940 f = (t2*r-30.0/19.0*t2+50.0/57.0*t1*r-200.0/969.0*t1+25.0/1292.0*r-1.0/
941 1938.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
942 /* p=11 5 2 0 */
943 } else if( i==5 && j==2 && k==0 ) {
944 t1 = r*r;
945 t2 = t1*t1;
946 t6 = s*s;
947 t7 = r-1.0;
948 t9 = t7*t7;
949 f = (t2*r-30.0/19.0*t2+50.0/57.0*t1*r-200.0/969.0*t1+25.0/1292.0*r-1.0/
950 1938.0)*(t6+6.0/7.0*s*t7+t9/7);
951 /* p=11 6 0 1 */
952 } else if( i==6 && j==0 && k==1 ) {
953 t1 = r*r;
954 t2 = t1*t1;
955 f = (t2*t1-42.0/19.0*t2*r+35.0/19.0*t2-700.0/969.0*t1*r+175.0/1292.0*t1
956 -7.0/646.0*r+1.0/3876.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
957 /* p=11 6 1 0 */
958 } else if( i==6 && j==1 && k==0 ) {
959 t1 = r*r;
960 t2 = t1*t1;
961 f = (t2*t1-42.0/19.0*t2*r+35.0/19.0*t2-700.0/969.0*t1*r+175.0/1292.0*t1
962 -7.0/646.0*r+1.0/3876.0)*(s+2.0/5.0*r-2.0/5.0);
963 /* p=11 7 0 0 */
964 } else if( i==7 && j==0 && k==0 ) {
965 t1 = r*r;
966 t2 = t1*r;
967 t3 = t1*t1;
968 f = t3*t2-56.0/19.0*t3*t1+196.0/57.0*t3*r-1960.0/969.0*t3+1225.0/1938.0*
969 t2-98.0/969.0*t1+7.0/969.0*r-2.0/12597.0;
970 /* p=12 0 0 8 */
971 } else if( i==0 && j==0 && k==8 ) {
972 t1 = t*t;
973 t2 = t1*t1;
974 t3 = t2*t2;
975 t4 = t1*t;
976 t6 = s+r-1.0;
977 t9 = t6*t6;
978 t12 = t9*t6;
979 t14 = t9*t9;
980 t22 = t14*t14;
981 f = t3+72.0/17.0*t2*t4*t6+126.0/17.0*t2*t1*t9+588.0/85.0*t2*t*t12+63.0/
982 17.0*t2*t14+252.0/221.0*t4*t14*t6+42.0/221.0*t1*t14*t9+36.0/2431.0*t*t14*t12+
983 9.0/24310.0*t22;
984 /* p=12 0 1 7 */
985 } else if( i==0 && j==1 && k==7 ) {
986 t2 = t*t;
987 t3 = t2*t;
988 t4 = t2*t2;
989 t7 = s+r-1.0;
990 t10 = t7*t7;
991 t12 = t10*t7;
992 t14 = t10*t10;
993 f = (s+2.0/19.0*r-2.0/19.0)*(t4*t3+56.0/15.0*t4*t2*t7+28.0/5.0*t4*t*t10
994 +56.0/13.0*t4*t12+70.0/39.0*t3*t14+56.0/143.0*t2*t14*t7+28.0/715.0*t*t14*t10+
995 8.0/6435.0*t14*t12);
996 /* p=12 0 2 6 */
997 } else if( i==0 && j==2 && k==6 ) {
998 t1 = s*s;
999 t2 = r-1.0;
1000 t4 = t2*t2;
1001 t6 = t*t;
1002 t7 = t6*t6;
1003 t10 = s+r-1.0;
1004 t12 = t10*t10;
1005 t17 = t12*t12;
1006 f = (t1+6.0/19.0*s*t2+t4/57)*(t6*t7+42.0/13.0*t7*t*t10+105.0/26.0*t7*
1007 t12+350.0/143.0*t6*t*t12*t10+105.0/143.0*t6*t17+14.0/143.0*t*t17*t10+7.0/1716.0
1008 *t17*t12);
1009 /* p=12 0 3 5 */
1010 } else if( i==0 && j==3 && k==5 ) {
1011 t1 = s*s;
1012 t3 = r-1.0;
1013 t5 = t3*t3;
1014 t9 = t*t;
1015 t10 = t9*t9;
1016 t12 = s+r-1.0;
1017 t15 = t12*t12;
1018 t19 = t15*t15;
1019 f = (t1*s+12.0/19.0*t1*t3+2.0/19.0*s*t5+4.0/969.0*t5*t3)*(t10*t+30.0/
1020 11.0*t12*t10+30.0/11.0*t9*t*t15+40.0/33.0*t9*t15*t12+5.0/22.0*t*t19+t19*t12/77)
1021 ;
1022 /* p=12 0 4 4 */
1023 } else if( i==0 && j==4 && k==4 ) {
1024 t1 = s*s;
1025 t2 = t1*t1;
1026 t4 = r-1.0;
1027 t6 = t4*t4;
1028 t10 = t6*t6;
1029 t12 = t*t;
1030 t13 = t12*t12;
1031 t15 = s+r-1.0;
1032 t17 = t15*t15;
1033 t21 = t17*t17;
1034 f = (t2+20.0/19.0*t4*t1*s+20.0/57.0*t1*t6+40.0/969.0*s*t6*t4+5.0/3876.0
1035 *t10)*(t13+20.0/9.0*t12*t*t15+5.0/3.0*t17*t12+10.0/21.0*t*t17*t15+5.0/126.0*t21
1036 );
1037 /* p=12 0 5 3 */
1038 } else if( i==0 && j==5 && k==3 ) {
1039 t1 = s*s;
1040 t2 = t1*t1;
1041 t4 = r-1.0;
1042 t7 = t4*t4;
1043 t11 = t7*t7;
1044 t15 = t*t;
1045 t17 = s+r-1.0;
1046 t19 = t17*t17;
1047 f = (s*t2+30.0/19.0*t4*t2+50.0/57.0*t1*s*t7+200.0/969.0*t1*t7*t4+25.0/
1048 1292.0*s*t11+t11*t4/1938)*(t15*t+12.0/7.0*t17*t15+6.0/7.0*t*t19+4.0/35.0*t19*
1049 t17);
1050 /* p=12 0 6 2 */
1051 } else if( i==0 && j==6 && k==2 ) {
1052 t1 = s*s;
1053 t2 = t1*t1;
1054 t5 = r-1.0;
1055 t7 = t5*t5;
1056 t12 = t7*t7;
1057 t18 = t*t;
1058 t19 = s+r-1.0;
1059 t21 = t19*t19;
1060 f = (t2*t1+42.0/19.0*t2*s*t5+35.0/19.0*t2*t7+700.0/969.0*t1*s*t7*t5+
1061 175.0/1292.0*t1*t12+7.0/646.0*s*t12*t5+t7*t12/3876)*(t18+6.0/5.0*t*t19+3.0/10.0
1062 *t21);
1063 /* p=12 0 7 1 */
1064 } else if( i==0 && j==7 && k==1 ) {
1065 t1 = s*s;
1066 t2 = t1*s;
1067 t3 = t1*t1;
1068 t6 = r-1.0;
1069 t9 = t6*t6;
1070 t11 = t9*t6;
1071 t13 = t9*t9;
1072 f = (t3*t2+56.0/19.0*t1*t3*t6+196.0/57.0*t3*s*t9+1960.0/969.0*t3*t11+
1073 1225.0/1938.0*t2*t13+98.0/969.0*t1*t13*t6+7.0/969.0*s*t13*t9+2.0/12597.0*t13*
1074 t11)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
1075 /* p=12 0 8 0 */
1076 } else if( i==0 && j==8 && k==0 ) {
1077 t1 = s*s;
1078 t2 = t1*t1;
1079 t3 = t2*t2;
1080 t4 = t1*s;
1081 t6 = r-1.0;
1082 t9 = t6*t6;
1083 t12 = t9*t6;
1084 t14 = t9*t9;
1085 t22 = t14*t14;
1086 f = t3+72.0/19.0*t2*t4*t6+112.0/19.0*t2*t1*t9+1568.0/323.0*t2*s*t12+
1087 735.0/323.0*t2*t14+196.0/323.0*t4*t14*t6+28.0/323.0*t1*t14*t9+24.0/4199.0*s*t14
1088 *t12+t22/8398;
1089 /* p=12 1 0 7 */
1090 } else if( i==1 && j==0 && k==7 ) {
1091 t2 = t*t;
1092 t3 = t2*t;
1093 t4 = t2*t2;
1094 t7 = s+r-1.0;
1095 t10 = t7*t7;
1096 t12 = t10*t7;
1097 t14 = t10*t10;
1098 f = (r-2.0/21.0)*(t4*t3+56.0/15.0*t2*t4*t7+28.0/5.0*t4*t*t10+56.0/13.0*
1099 t4*t12+70.0/39.0*t3*t14+56.0/143.0*t2*t14*t7+28.0/715.0*t*t14*t10+8.0/6435.0*
1100 t14*t12);
1101 /* p=12 1 1 6 */
1102 } else if( i==1 && j==1 && k==6 ) {
1103 t4 = t*t;
1104 t5 = t4*t4;
1105 t8 = s+r-1.0;
1106 t10 = t8*t8;
1107 t15 = t10*t10;
1108 f = (s+2.0/17.0*r-2.0/17.0)*(r-2.0/21.0)*(t5*t4+42.0/13.0*t5*t*t8+105.0
1109 /26.0*t5*t10+350.0/143.0*t4*t*t10*t8+105.0/143.0*t4*t15+14.0/143.0*t*t15*t8+7.0
1110 /1716.0*t15*t10);
1111 /* p=12 1 2 5 */
1112 } else if( i==1 && j==2 && k==5 ) {
1113 t2 = s*s;
1114 t3 = r-1.0;
1115 t5 = t3*t3;
1116 t8 = t*t;
1117 t9 = t8*t8;
1118 t11 = s+r-1.0;
1119 t14 = t11*t11;
1120 t18 = t14*t14;
1121 f = (r-2.0/21.0)*(t2+6.0/17.0*s*t3+3.0/136.0*t5)*(t9*t+30.0/11.0*t9*t11
1122 +30.0/11.0*t8*t*t14+40.0/33.0*t8*t14*t11+5.0/22.0*t*t18+t18*t11/77);
1123 /* p=12 1 3 4 */
1124 } else if( i==1 && j==3 && k==4 ) {
1125 t2 = s*s;
1126 t4 = r-1.0;
1127 t6 = t4*t4;
1128 t11 = t*t;
1129 t12 = t11*t11;
1130 t14 = s+r-1.0;
1131 t16 = t14*t14;
1132 t20 = t16*t16;
1133 f = (r-2.0/21.0)*(t2*s+12.0/17.0*t2*t4+9.0/68.0*s*t6+t6*t4/170)*(t12+
1134 20.0/9.0*t11*t*t14+5.0/3.0*t11*t16+10.0/21.0*t*t16*t14+5.0/126.0*t20);
1135 /* p=12 1 4 3 */
1136 } else if( i==1 && j==4 && k==3 ) {
1137 t2 = s*s;
1138 t3 = t2*t2;
1139 t5 = r-1.0;
1140 t7 = t5*t5;
1141 t11 = t7*t7;
1142 t14 = t*t;
1143 t16 = s+r-1.0;
1144 t18 = t16*t16;
1145 f = (r-2.0/21.0)*(t3+20.0/17.0*t5*t2*s+15.0/34.0*t2*t7+s*t7*t5/17+t11
1146 /476)*(t14*t+12.0/7.0*t16*t14+6.0/7.0*t*t18+4.0/35.0*t18*t16);
1147 /* p=12 1 5 2 */
1148 } else if( i==1 && j==5 && k==2 ) {
1149 t2 = s*s;
1150 t3 = t2*t2;
1151 t5 = r-1.0;
1152 t8 = t5*t5;
1153 t12 = t8*t8;
1154 t17 = t*t;
1155 t18 = s+r-1.0;
1156 t20 = t18*t18;
1157 f = (r-2.0/21.0)*(s*t3+30.0/17.0*t3*t5+75.0/68.0*t2*s*t8+5.0/17.0*t2*t8
1158 *t5+15.0/476.0*s*t12+3.0/3094.0*t12*t5)*(t17+6.0/5.0*t*t18+3.0/10.0*t20);
1159 /* p=12 1 6 1 */
1160 } else if( i==1 && j==6 && k==1 ) {
1161 t2 = s*s;
1162 t3 = t2*t2;
1163 t6 = r-1.0;
1164 t8 = t6*t6;
1165 t13 = t8*t8;
1166 f = (r-2.0/21.0)*(t3*t2+42.0/17.0*t3*s*t6+315.0/136.0*t3*t8+35.0/34.0*
1167 t2*s*t8*t6+15.0/68.0*t2*t13+9.0/442.0*s*t13*t6+t13*t8/1768)*(t+2.0/3.0*s+2.0/
1168 3.0*r-2.0/3.0);
1169 /* p=12 1 7 0 */
1170 } else if( i==1 && j==7 && k==0 ) {
1171 t2 = s*s;
1172 t3 = t2*s;
1173 t4 = t2*t2;
1174 t7 = r-1.0;
1175 t10 = t7*t7;
1176 t12 = t10*t7;
1177 t14 = t10*t10;
1178 f = (r-2.0/21.0)*(t4*t3+56.0/17.0*t4*t2*t7+147.0/34.0*t4*s*t10+49.0/
1179 17.0*t4*t12+35.0/34.0*t3*t14+42.0/221.0*t2*t14*t7+7.0/442.0*s*t14*t10+t14*t12
1180 /2431);
1181 /* p=12 2 0 6 */
1182 } else if( i==2 && j==0 && k==6 ) {
1183 t1 = r*r;
1184 t3 = t*t;
1185 t4 = t3*t3;
1186 t7 = s+r-1.0;
1187 t9 = t7*t7;
1188 t14 = t9*t9;
1189 f = (t1-2.0/7.0*r+1.0/70.0)*(t4*t3+42.0/13.0*t4*t*t7+105.0/26.0*t4*t9+
1190 350.0/143.0*t3*t*t9*t7+105.0/143.0*t3*t14+14.0/143.0*t*t14*t7+7.0/1716.0*t14*t9
1191 );
1192 /* p=12 2 1 5 */
1193 } else if( i==2 && j==1 && k==5 ) {
1194 t1 = r*r;
1195 t5 = t*t;
1196 t6 = t5*t5;
1197 t8 = s+r-1.0;
1198 t11 = t8*t8;
1199 t15 = t11*t11;
1200 f = (s+2.0/15.0*r-2.0/15.0)*(t1-2.0/7.0*r+1.0/70.0)*(t6*t+30.0/11.0*t6*
1201 t8+30.0/11.0*t5*t*t11+40.0/33.0*t5*t11*t8+5.0/22.0*t*t15+t15*t8/77);
1202 /* p=12 2 2 4 */
1203 } else if( i==2 && j==2 && k==4 ) {
1204 t1 = r*r;
1205 t3 = s*s;
1206 t4 = r-1.0;
1207 t6 = t4*t4;
1208 t9 = t*t;
1209 t10 = t9*t9;
1210 t12 = s+r-1.0;
1211 t14 = t12*t12;
1212 t18 = t14*t14;
1213 f = (t1-2.0/7.0*r+1.0/70.0)*(t3+2.0/5.0*s*t4+t6/35)*(t10+20.0/9.0*t9*t*
1214 t12+5.0/3.0*t14*t9+10.0/21.0*t*t14*t12+5.0/126.0*t18);
1215 /* p=12 2 3 3 */
1216 } else if( i==2 && j==3 && k==3 ) {
1217 t1 = r*r;
1218 t3 = s*s;
1219 t5 = r-1.0;
1220 t7 = t5*t5;
1221 t12 = t*t;
1222 t14 = s+r-1.0;
1223 t16 = t14*t14;
1224 f = (t1-2.0/7.0*r+1.0/70.0)*(t3*s+4.0/5.0*t3*t5+6.0/35.0*s*t7+4.0/455.0
1225 *t7*t5)*(t12*t+12.0/7.0*t14*t12+6.0/7.0*t*t16+4.0/35.0*t16*t14);
1226 /* p=12 2 4 2 */
1227 } else if( i==2 && j==4 && k==2 ) {
1228 t1 = r*r;
1229 t3 = s*s;
1230 t4 = t3*t3;
1231 t6 = r-1.0;
1232 t8 = t6*t6;
1233 t12 = t8*t8;
1234 t15 = t*t;
1235 t16 = s+r-1.0;
1236 t18 = t16*t16;
1237 f = (t1-2.0/7.0*r+1.0/70.0)*(t4+4.0/3.0*t3*s*t6+4.0/7.0*t3*t8+8.0/91.0*
1238 s*t6*t8+t12/273)*(t15+6.0/5.0*t*t16+3.0/10.0*t18);
1239 /* p=12 2 5 1 */
1240 } else if( i==2 && j==5 && k==1 ) {
1241 t1 = r*r;
1242 t3 = s*s;
1243 t4 = t3*t3;
1244 t6 = r-1.0;
1245 t9 = t6*t6;
1246 t13 = t9*t9;
1247 f = (t1-2.0/7.0*r+1.0/70.0)*(s*t4+2.0*t4*t6+10.0/7.0*t3*s*t9+40.0/91.0*
1248 t3*t9*t6+5.0/91.0*s*t13+2.0/1001.0*t13*t6)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
1249 /* p=12 2 6 0 */
1250 } else if( i==2 && j==6 && k==0 ) {
1251 t1 = r*r;
1252 t3 = s*s;
1253 t4 = t3*t3;
1254 t7 = r-1.0;
1255 t9 = t7*t7;
1256 t14 = t9*t9;
1257 f = (t1-2.0/7.0*r+1.0/70.0)*(t4*t3+14.0/5.0*t4*s*t7+3.0*t4*t9+20.0/13.0
1258 *t3*s*t9*t7+5.0/13.0*t3*t14+6.0/143.0*s*t14*t7+t14*t9/715);
1259 /* p=12 3 0 5 */
1260 } else if( i==3 && j==0 && k==5 ) {
1261 t1 = r*r;
1262 t4 = t*t;
1263 t5 = t4*t4;
1264 t7 = s+r-1.0;
1265 t10 = t7*t7;
1266 t14 = t10*t10;
1267 f = (t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(t5*t+30.0/11.0*t5*t7+30.0/
1268 11.0*t4*t*t10+40.0/33.0*t4*t10*t7+5.0/22.0*t*t14+t14*t7/77);
1269 /* p=12 3 1 4 */
1270 } else if( i==3 && j==1 && k==4 ) {
1271 t1 = r*r;
1272 t6 = t*t;
1273 t7 = t6*t6;
1274 t9 = s+r-1.0;
1275 t11 = t9*t9;
1276 t15 = t11*t11;
1277 f = (s+2.0/13.0*r-2.0/13.0)*(t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(t7+
1278 20.0/9.0*t6*t*t9+5.0/3.0*t6*t11+10.0/21.0*t*t11*t9+5.0/126.0*t15);
1279 /* p=12 3 2 3 */
1280 } else if( i==3 && j==2 && k==3 ) {
1281 t1 = r*r;
1282 t4 = s*s;
1283 t5 = r-1.0;
1284 t7 = t5*t5;
1285 t10 = t*t;
1286 t12 = s+r-1.0;
1287 t14 = t12*t12;
1288 f = (t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(t4+6.0/13.0*s*t5+t7/26)*(
1289 t10*t+12.0/7.0*t10*t12+6.0/7.0*t*t14+4.0/35.0*t14*t12);
1290 /* p=12 3 3 2 */
1291 } else if( i==3 && j==3 && k==2 ) {
1292 t1 = r*r;
1293 t4 = s*s;
1294 t6 = r-1.0;
1295 t8 = t6*t6;
1296 t13 = t*t;
1297 t14 = s+r-1.0;
1298 t16 = t14*t14;
1299 f = (t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(t4*s+12.0/13.0*t4*t6+3.0/
1300 13.0*s*t8+2.0/143.0*t8*t6)*(t13+6.0/5.0*t*t14+3.0/10.0*t16);
1301 /* p=12 3 4 1 */
1302 } else if( i==3 && j==4 && k==1 ) {
1303 t1 = r*r;
1304 t4 = s*s;
1305 t5 = t4*t4;
1306 t7 = r-1.0;
1307 t9 = t7*t7;
1308 t13 = t9*t9;
1309 f = (t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(t5+20.0/13.0*t4*s*t7+10.0/
1310 13.0*t4*t9+20.0/143.0*s*t9*t7+t13/143)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
1311 /* p=12 3 5 0 */
1312 } else if( i==3 && j==5 && k==0 ) {
1313 t1 = r*r;
1314 t4 = s*s;
1315 t5 = t4*t4;
1316 t7 = r-1.0;
1317 t10 = t7*t7;
1318 t14 = t10*t10;
1319 f = (t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(s*t5+30.0/13.0*t5*t7+25.0/
1320 13.0*t4*s*t10+100.0/143.0*t4*t10*t7+15.0/143.0*s*t14+2.0/429.0*t14*t7);
1321 /* p=12 4 0 4 */
1322 } else if( i==4 && j==0 && k==4 ) {
1323 t1 = r*r;
1324 t2 = t1*t1;
1325 t5 = t*t;
1326 t6 = t5*t5;
1327 t8 = s+r-1.0;
1328 t10 = t8*t8;
1329 t14 = t10*t10;
1330 f = (t2-20.0/21.0*t1*r+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0)*(t6+20.0/9.0*
1331 t5*t*t8+5.0/3.0*t5*t10+10.0/21.0*t*t10*t8+5.0/126.0*t14);
1332 /* p=12 4 1 3 */
1333 } else if( i==4 && j==1 && k==3 ) {
1334 t1 = r*r;
1335 t2 = t1*t1;
1336 t7 = t*t;
1337 t9 = s+r-1.0;
1338 t11 = t9*t9;
1339 f = (t2-20.0/21.0*t1*r+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0)*(s+2.0/11.0*r
1340 -2.0/11.0)*(t7*t+12.0/7.0*t9*t7+6.0/7.0*t*t11+4.0/35.0*t11*t9);
1341 /* p=12 4 2 2 */
1342 } else if( i==4 && j==2 && k==2 ) {
1343 t1 = r*r;
1344 t2 = t1*t1;
1345 t5 = s*s;
1346 t6 = r-1.0;
1347 t8 = t6*t6;
1348 t11 = t*t;
1349 t12 = s+r-1.0;
1350 t14 = t12*t12;
1351 f = (t2-20.0/21.0*t1*r+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0)*(t5+6.0/11.0*
1352 s*t6+3.0/55.0*t8)*(t11+6.0/5.0*t*t12+3.0/10.0*t14);
1353 /* p=12 4 3 1 */
1354 } else if( i==4 && j==3 && k==1 ) {
1355 t1 = r*r;
1356 t2 = t1*t1;
1357 t5 = s*s;
1358 t7 = r-1.0;
1359 t9 = t7*t7;
1360 f = (t2-20.0/21.0*t1*r+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0)*(t5*s+12.0/
1361 11.0*t5*t7+18.0/55.0*s*t9+4.0/165.0*t9*t7)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
1362 /* p=12 4 4 0 */
1363 } else if( i==4 && j==4 && k==0 ) {
1364 t1 = r*r;
1365 t2 = t1*t1;
1366 t5 = s*s;
1367 t6 = t5*t5;
1368 t8 = r-1.0;
1369 t10 = t8*t8;
1370 t14 = t10*t10;
1371 f = (t2-20.0/21.0*t1*r+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0)*(t6+20.0/11.0
1372 *t5*s*t8+12.0/11.0*t5*t10+8.0/33.0*s*t10*t8+t14/66);
1373 /* p=12 5 0 3 */
1374 } else if( i==5 && j==0 && k==3 ) {
1375 t1 = r*r;
1376 t2 = t1*t1;
1377 t6 = t*t;
1378 t8 = s+r-1.0;
1379 t10 = t8*t8;
1380 f = (t2*r-10.0/7.0*t2+5.0/7.0*t1*r-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0
1381 )*(t6*t+12.0/7.0*t6*t8+6.0/7.0*t*t10+4.0/35.0*t10*t8);
1382 /* p=12 5 1 2 */
1383 } else if( i==5 && j==1 && k==2 ) {
1384 t1 = r*r;
1385 t2 = t1*t1;
1386 t8 = t*t;
1387 t9 = s+r-1.0;
1388 t11 = t9*t9;
1389 f = (t2*r-10.0/7.0*t2+5.0/7.0*t1*r-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0
1390 )*(s+2.0/9.0*r-2.0/9.0)*(t8+6.0/5.0*t*t9+3.0/10.0*t11);
1391 /* p=12 5 2 1 */
1392 } else if( i==5 && j==2 && k==1 ) {
1393 t1 = r*r;
1394 t2 = t1*t1;
1395 t6 = s*s;
1396 t7 = r-1.0;
1397 t9 = t7*t7;
1398 f = (t2*r-10.0/7.0*t2+5.0/7.0*t1*r-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0
1399 )*(t6+2.0/3.0*s*t7+t9/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
1400 /* p=12 5 3 0 */
1401 } else if( i==5 && j==3 && k==0 ) {
1402 t1 = r*r;
1403 t2 = t1*t1;
1404 t6 = s*s;
1405 t8 = r-1.0;
1406 t10 = t8*t8;
1407 f = (t2*r-10.0/7.0*t2+5.0/7.0*t1*r-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0
1408 )*(s*t6+4.0/3.0*t6*t8+s*t10/2+t10*t8/21);
1409 /* p=12 6 0 2 */
1410 } else if( i==6 && j==0 && k==2 ) {
1411 t1 = r*r;
1412 t2 = t1*t1;
1413 t7 = t*t;
1414 t8 = s+r-1.0;
1415 t10 = t8*t8;
1416 f = (t2*t1-2.0*t2*r+3.0/2.0*t2-10.0/19.0*t1*r+5.0/57.0*t1-2.0/323.0*r+
1417 1.0/7752.0)*(t7+6.0/5.0*t*t8+3.0/10.0*t10);
1418 /* p=12 6 1 1 */
1419 } else if( i==6 && j==1 && k==1 ) {
1420 t1 = r*r;
1421 t2 = t1*t1;
1422 f = (t2*t1-2.0*t2*r+3.0/2.0*t2-10.0/19.0*t1*r+5.0/57.0*t1-2.0/323.0*r+
1423 1.0/7752.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
1424 /* p=12 6 2 0 */
1425 } else if( i==6 && j==2 && k==0 ) {
1426 t1 = r*r;
1427 t2 = t1*t1;
1428 t7 = s*s;
1429 t8 = r-1.0;
1430 t10 = t8*t8;
1431 f = (t2*t1-2.0*t2*r+3.0/2.0*t2-10.0/19.0*t1*r+5.0/57.0*t1-2.0/323.0*r+
1432 1.0/7752.0)*(t7+6.0/7.0*s*t8+t10/7);
1433 /* p=12 7 0 1 */
1434 } else if( i==7 && j==0 && k==1 ) {
1435 t1 = r*r;
1436 t2 = t1*r;
1437 t3 = t1*t1;
1438 f = (t3*t2-8.0/3.0*t3*t1+14.0/5.0*t3*r-28.0/19.0*t3+70.0/171.0*t2-56.0/
1439 969.0*t1+7.0/1938.0*r-1.0/14535.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
1440 /* p=12 7 1 0 */
1441 } else if( i==7 && j==1 && k==0 ) {
1442 t1 = r*r;
1443 t2 = t1*r;
1444 t3 = t1*t1;
1445 f = (t3*t2-8.0/3.0*t3*t1+14.0/5.0*t3*r-28.0/19.0*t3+70.0/171.0*t2-56.0/
1446 969.0*t1+7.0/1938.0*r-1.0/14535.0)*(s+2.0/5.0*r-2.0/5.0);
1447 /* p=12 8 0 0 */
1448 } else if( i==8 && j==0 && k==0 ) {
1449 t1 = r*r;
1450 t2 = t1*t1;
1451 t3 = t2*t2;
1452 t4 = t1*r;
1453 f = t3-24.0/7.0*t2*t4+24.0/5.0*t2*t1-336.0/95.0*t2*r+28.0/19.0*t2-112.0/
1454 323.0*t4+14.0/323.0*t1-4.0/1615.0*r+1.0/22610.0;
1455 /* p=13 0 0 9 */
1456 } else if( i==0 && j==0 && k==9 ) {
1457 t1 = t*t;
1458 t2 = t1*t1;
1459 t3 = t2*t2;
1460 t5 = s+r-1.0;
1461 t7 = t1*t;
1462 t9 = t5*t5;
1463 t12 = t9*t5;
1464 t15 = t9*t9;
1465 t23 = t15*t15;
1466 f = t3*t+90.0/19.0*t3*t5+180.0/19.0*t2*t7*t9+3360.0/323.0*t2*t1*t12+
1467 2205.0/323.0*t2*t*t15+882.0/323.0*t2*t15*t5+210.0/323.0*t7*t15*t9+360.0/4199.0*
1468 t1*t15*t12+45.0/8398.0*t*t23+5.0/46189.0*t23*t5;
1469 /* p=13 0 1 8 */
1470 } else if( i==0 && j==1 && k==8 ) {
1471 t2 = t*t;
1472 t3 = t2*t2;
1473 t4 = t3*t3;
1474 t5 = t2*t;
1475 t7 = s+r-1.0;
1476 t10 = t7*t7;
1477 t13 = t10*t7;
1478 t15 = t10*t10;
1479 t23 = t15*t15;
1480 f = (s+2.0/21.0*r-2.0/21.0)*(t4+72.0/17.0*t3*t5*t7+126.0/17.0*t3*t2*t10
1481 +588.0/85.0*t3*t*t13+63.0/17.0*t3*t15+252.0/221.0*t5*t15*t7+42.0/221.0*t2*t15*
1482 t10+36.0/2431.0*t*t15*t13+9.0/24310.0*t23);
1483 /* p=13 0 2 7 */
1484 } else if( i==0 && j==2 && k==7 ) {
1485 t1 = s*s;
1486 t2 = r-1.0;
1487 t4 = t2*t2;
1488 t6 = t*t;
1489 t7 = t6*t;
1490 t8 = t6*t6;
1491 t11 = s+r-1.0;
1492 t14 = t11*t11;
1493 t16 = t14*t11;
1494 t18 = t14*t14;
1495 f = (t1+2.0/7.0*s*t2+t4/70)*(t7*t8+56.0/15.0*t8*t6*t11+28.0/5.0*t8*t*
1496 t14+56.0/13.0*t8*t16+70.0/39.0*t7*t18+56.0/143.0*t6*t18*t11+28.0/715.0*t*t18*
1497 t14+8.0/6435.0*t18*t16);
1498 /* p=13 0 3 6 */
1499 } else if( i==0 && j==3 && k==6 ) {
1500 t1 = s*s;
1501 t3 = r-1.0;
1502 t5 = t3*t3;
1503 t9 = t*t;
1504 t10 = t9*t9;
1505 t13 = s+r-1.0;
1506 t15 = t13*t13;
1507 t20 = t15*t15;
1508 f = (t1*s+4.0/7.0*t3*t1+3.0/35.0*s*t5+2.0/665.0*t3*t5)*(t9*t10+42.0/
1509 13.0*t10*t*t13+105.0/26.0*t15*t10+350.0/143.0*t9*t*t15*t13+105.0/143.0*t9*t20+
1510 14.0/143.0*t*t20*t13+7.0/1716.0*t20*t15);
1511 /* p=13 0 4 5 */
1512 } else if( i==0 && j==4 && k==5 ) {
1513 t1 = s*s;
1514 t2 = t1*t1;
1515 t4 = r-1.0;
1516 t6 = t4*t4;
1517 t10 = t6*t6;
1518 t12 = t*t;
1519 t13 = t12*t12;
1520 t15 = s+r-1.0;
1521 t18 = t15*t15;
1522 t22 = t18*t18;
1523 f = (t2+20.0/21.0*t1*s*t4+2.0/7.0*t1*t6+4.0/133.0*s*t6*t4+t10/1197)*(
1524 t13*t+30.0/11.0*t13*t15+30.0/11.0*t12*t*t18+40.0/33.0*t12*t18*t15+5.0/22.0*t*
1525 t22+t22*t15/77);
1526 /* p=13 0 5 4 */
1527 } else if( i==0 && j==5 && k==4 ) {
1528 t1 = s*s;
1529 t2 = t1*t1;
1530 t4 = r-1.0;
1531 t7 = t4*t4;
1532 t11 = t7*t7;
1533 t15 = t*t;
1534 t16 = t15*t15;
1535 t18 = s+r-1.0;
1536 t20 = t18*t18;
1537 t24 = t20*t20;
1538 f = (t2*s+10.0/7.0*t2*t4+5.0/7.0*t1*s*t7+20.0/133.0*t1*t7*t4+5.0/399.0*
1539 s*t11+2.0/6783.0*t11*t4)*(t16+20.0/9.0*t15*t*t18+5.0/3.0*t20*t15+10.0/21.0*t*
1540 t20*t18+5.0/126.0*t24);
1541 /* p=13 0 6 3 */
1542 } else if( i==0 && j==6 && k==3 ) {
1543 t1 = s*s;
1544 t2 = t1*t1;
1545 t5 = r-1.0;
1546 t7 = t5*t5;
1547 t12 = t7*t7;
1548 t18 = t*t;
1549 t20 = s+r-1.0;
1550 t22 = t20*t20;
1551 f = (t2*t1+2.0*t2*s*t5+3.0/2.0*t2*t7+10.0/19.0*t1*s*t7*t5+5.0/57.0*t1*
1552 t12+2.0/323.0*s*t12*t5+t12*t7/7752)*(t18*t+12.0/7.0*t20*t18+6.0/7.0*t*t22+4.0/
1553 35.0*t22*t20);
1554 /* p=13 0 7 2 */
1555 } else if( i==0 && j==7 && k==2 ) {
1556 t1 = s*s;
1557 t2 = t1*s;
1558 t3 = t1*t1;
1559 t6 = r-1.0;
1560 t9 = t6*t6;
1561 t11 = t9*t6;
1562 t13 = t9*t9;
1563 t21 = t*t;
1564 t22 = s+r-1.0;
1565 t24 = t22*t22;
1566 f = (t3*t2+8.0/3.0*t3*t1*t6+14.0/5.0*t3*s*t9+28.0/19.0*t3*t11+70.0/
1567 171.0*t2*t13+56.0/969.0*t1*t13*t6+7.0/1938.0*s*t13*t9+t13*t11/14535)*(t21+6.0/
1568 5.0*t*t22+3.0/10.0*t24);
1569 /* p=13 0 8 1 */
1570 } else if( i==0 && j==8 && k==1 ) {
1571 t1 = s*s;
1572 t2 = t1*t1;
1573 t3 = t2*t2;
1574 t4 = t1*s;
1575 t6 = r-1.0;
1576 t9 = t6*t6;
1577 t12 = t9*t6;
1578 t14 = t9*t9;
1579 t22 = t14*t14;
1580 f = (t3+24.0/7.0*t2*t4*t6+24.0/5.0*t2*t1*t9+336.0/95.0*t2*s*t12+28.0/
1581 19.0*t2*t14+112.0/323.0*t4*t14*t6+14.0/323.0*t1*t14*t9+4.0/1615.0*s*t14*t12+t22
1582 /22610)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
1583 /* p=13 0 9 0 */
1584 } else if( i==0 && j==9 && k==0 ) {
1585 t1 = s*s;
1586 t2 = t1*t1;
1587 t3 = t2*t2;
1588 t5 = r-1.0;
1589 t7 = t1*s;
1590 t9 = t5*t5;
1591 t12 = t9*t5;
1592 t15 = t9*t9;
1593 t23 = t15*t15;
1594 f = t3*s+30.0/7.0*t3*t5+54.0/7.0*t2*t7*t9+144.0/19.0*t2*t1*t12+84.0/
1595 19.0*t2*s*t15+504.0/323.0*t2*t15*t5+105.0/323.0*t7*t15*t9+12.0/323.0*t1*t15*t12
1596 +9.0/4522.0*s*t23+t23*t5/29393;
1597 /* p=13 1 0 8 */
1598 } else if( i==1 && j==0 && k==8 ) {
1599 t2 = t*t;
1600 t3 = t2*t2;
1601 t4 = t3*t3;
1602 t5 = t2*t;
1603 t7 = s+r-1.0;
1604 t10 = t7*t7;
1605 t13 = t10*t7;
1606 t15 = t10*t10;
1607 t23 = t15*t15;
1608 f = (r-2.0/23.0)*(t4+72.0/17.0*t3*t5*t7+126.0/17.0*t3*t2*t10+588.0/85.0
1609 *t3*t*t13+63.0/17.0*t3*t15+252.0/221.0*t5*t15*t7+42.0/221.0*t2*t15*t10+36.0/
1610 2431.0*t*t13*t15+9.0/24310.0*t23);
1611 /* p=13 1 1 7 */
1612 } else if( i==1 && j==1 && k==7 ) {
1613 t4 = t*t;
1614 t5 = t4*t;
1615 t6 = t4*t4;
1616 t9 = s+r-1.0;
1617 t12 = t9*t9;
1618 t14 = t12*t9;
1619 t16 = t12*t12;
1620 f = (r-2.0/23.0)*(s+2.0/19.0*r-2.0/19.0)*(t6*t5+56.0/15.0*t6*t4*t9+28.0
1621 /5.0*t6*t*t12+56.0/13.0*t6*t14+70.0/39.0*t5*t16+56.0/143.0*t4*t16*t9+28.0/715.0
1622 *t*t16*t12+8.0/6435.0*t16*t14);
1623 /* p=13 1 2 6 */
1624 } else if( i==1 && j==2 && k==6 ) {
1625 t2 = s*s;
1626 t3 = r-1.0;
1627 t5 = t3*t3;
1628 t8 = t*t;
1629 t9 = t8*t8;
1630 t12 = s+r-1.0;
1631 t14 = t12*t12;
1632 t19 = t14*t14;
1633 f = (r-2.0/23.0)*(t2+6.0/19.0*s*t3+t5/57)*(t8*t9+42.0/13.0*t9*t*t12+
1634 105.0/26.0*t9*t14+350.0/143.0*t8*t*t14*t12+105.0/143.0*t8*t19+14.0/143.0*t*t19*
1635 t12+7.0/1716.0*t19*t14);
1636 /* p=13 1 3 5 */
1637 } else if( i==1 && j==3 && k==5 ) {
1638 t2 = s*s;
1639 t4 = r-1.0;
1640 t6 = t4*t4;
1641 t11 = t*t;
1642 t12 = t11*t11;
1643 t14 = s+r-1.0;
1644 t17 = t14*t14;
1645 t21 = t17*t17;
1646 f = (r-2.0/23.0)*(t2*s+12.0/19.0*t2*t4+2.0/19.0*s*t6+4.0/969.0*t6*t4)*(
1647 t12*t+30.0/11.0*t14*t12+30.0/11.0*t11*t*t17+40.0/33.0*t11*t17*t14+5.0/22.0*t*
1648 t21+t21*t14/77);
1649 /* p=13 1 4 4 */
1650 } else if( i==1 && j==4 && k==4 ) {
1651 t2 = s*s;
1652 t3 = t2*t2;
1653 t5 = r-1.0;
1654 t7 = t5*t5;
1655 t11 = t7*t7;
1656 t14 = t*t;
1657 t15 = t14*t14;
1658 t17 = s+r-1.0;
1659 t19 = t17*t17;
1660 t23 = t19*t19;
1661 f = (r-2.0/23.0)*(t3+20.0/19.0*t2*s*t5+20.0/57.0*t2*t7+40.0/969.0*s*t7*
1662 t5+5.0/3876.0*t11)*(t15+20.0/9.0*t14*t*t17+5.0/3.0*t19*t14+10.0/21.0*t*t19*t17+
1663 5.0/126.0*t23);
1664 /* p=13 1 5 3 */
1665 } else if( i==1 && j==5 && k==3 ) {
1666 t2 = s*s;
1667 t3 = t2*t2;
1668 t5 = r-1.0;
1669 t8 = t5*t5;
1670 t12 = t8*t8;
1671 t17 = t*t;
1672 t19 = s+r-1.0;
1673 t21 = t19*t19;
1674 f = (r-2.0/23.0)*(s*t3+30.0/19.0*t3*t5+50.0/57.0*t2*s*t8+200.0/969.0*t2
1675 *t8*t5+25.0/1292.0*s*t12+t12*t5/1938)*(t17*t+12.0/7.0*t19*t17+6.0/7.0*t*t21+4.0
1676 /35.0*t21*t19);
1677 /* p=13 1 6 2 */
1678 } else if( i==1 && j==6 && k==2 ) {
1679 t2 = s*s;
1680 t3 = t2*t2;
1681 t6 = r-1.0;
1682 t8 = t6*t6;
1683 t13 = t8*t8;
1684 t20 = t*t;
1685 t21 = s+r-1.0;
1686 t23 = t21*t21;
1687 f = (r-2.0/23.0)*(t3*t2+42.0/19.0*t6*t3*s+35.0/19.0*t3*t8+700.0/969.0*
1688 t2*s*t8*t6+175.0/1292.0*t2*t13+7.0/646.0*s*t13*t6+t13*t8/3876)*(t20+6.0/5.0*t*
1689 t21+3.0/10.0*t23);
1690 /* p=13 1 7 1 */
1691 } else if( i==1 && j==7 && k==1 ) {
1692 t2 = s*s;
1693 t3 = t2*s;
1694 t4 = t2*t2;
1695 t7 = r-1.0;
1696 t10 = t7*t7;
1697 t12 = t10*t7;
1698 t14 = t10*t10;
1699 f = (r-2.0/23.0)*(t4*t3+56.0/19.0*t4*t2*t7+196.0/57.0*t4*s*t10+1960.0/
1700 969.0*t4*t12+1225.0/1938.0*t3*t14+98.0/969.0*t2*t14*t7+7.0/969.0*s*t14*t10+2.0/
1701 12597.0*t14*t12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
1702 /* p=13 1 8 0 */
1703 } else if( i==1 && j==8 && k==0 ) {
1704 t2 = s*s;
1705 t3 = t2*t2;
1706 t4 = t3*t3;
1707 t5 = t2*s;
1708 t7 = r-1.0;
1709 t10 = t7*t7;
1710 t13 = t10*t7;
1711 t15 = t10*t10;
1712 t23 = t15*t15;
1713 f = (r-2.0/23.0)*(t4+72.0/19.0*t3*t5*t7+112.0/19.0*t3*t2*t10+1568.0/
1714 323.0*t3*s*t13+735.0/323.0*t3*t15+196.0/323.0*t5*t15*t7+28.0/323.0*t2*t15*t10+
1715 24.0/4199.0*s*t15*t13+t23/8398);
1716 /* p=13 2 0 7 */
1717 } else if( i==2 && j==0 && k==7 ) {
1718 t1 = r*r;
1719 t3 = t*t;
1720 t4 = t3*t;
1721 t5 = t3*t3;
1722 t8 = s+r-1.0;
1723 t11 = t8*t8;
1724 t13 = t11*t8;
1725 t15 = t11*t11;
1726 f = (t1-6.0/23.0*r+3.0/253.0)*(t5*t4+56.0/15.0*t3*t5*t8+28.0/5.0*t5*t*
1727 t11+56.0/13.0*t5*t13+70.0/39.0*t4*t15+56.0/143.0*t3*t15*t8+28.0/715.0*t*t15*t11
1728 +8.0/6435.0*t15*t13);
1729 /* p=13 2 1 6 */
1730 } else if( i==2 && j==1 && k==6 ) {
1731 t1 = r*r;
1732 t5 = t*t;
1733 t6 = t5*t5;
1734 t9 = s+r-1.0;
1735 t11 = t9*t9;
1736 t16 = t11*t11;
1737 f = (s+2.0/17.0*r-2.0/17.0)*(t1-6.0/23.0*r+3.0/253.0)*(t6*t5+42.0/13.0*
1738 t6*t*t9+105.0/26.0*t6*t11+350.0/143.0*t5*t*t11*t9+105.0/143.0*t5*t16+14.0/143.0
1739 *t*t16*t9+7.0/1716.0*t16*t11);
1740 /* p=13 2 2 5 */
1741 } else if( i==2 && j==2 && k==5 ) {
1742 t1 = r*r;
1743 t3 = s*s;
1744 t4 = r-1.0;
1745 t6 = t4*t4;
1746 t9 = t*t;
1747 t10 = t9*t9;
1748 t12 = s+r-1.0;
1749 t15 = t12*t12;
1750 t19 = t15*t15;
1751 f = (t1-6.0/23.0*r+3.0/253.0)*(t3+6.0/17.0*s*t4+3.0/136.0*t6)*(t10*t+
1752 30.0/11.0*t10*t12+30.0/11.0*t9*t*t15+40.0/33.0*t9*t15*t12+5.0/22.0*t*t19+t19*
1753 t12/77);
1754 /* p=13 2 3 4 */
1755 } else if( i==2 && j==3 && k==4 ) {
1756 t1 = r*r;
1757 t3 = s*s;
1758 t5 = r-1.0;
1759 t7 = t5*t5;
1760 t12 = t*t;
1761 t13 = t12*t12;
1762 t15 = s+r-1.0;
1763 t17 = t15*t15;
1764 t21 = t17*t17;
1765 f = (t1-6.0/23.0*r+3.0/253.0)*(t3*s+12.0/17.0*t3*t5+9.0/68.0*s*t7+t7*t5
1766 /170)*(t13+20.0/9.0*t12*t*t15+5.0/3.0*t12*t17+10.0/21.0*t*t17*t15+5.0/126.0*t21
1767 );
1768 /* p=13 2 4 3 */
1769 } else if( i==2 && j==4 && k==3 ) {
1770 t1 = r*r;
1771 t3 = s*s;
1772 t4 = t3*t3;
1773 t6 = r-1.0;
1774 t8 = t6*t6;
1775 t12 = t8*t8;
1776 t15 = t*t;
1777 t17 = s+r-1.0;
1778 t19 = t17*t17;
1779 f = (t1-6.0/23.0*r+3.0/253.0)*(t4+20.0/17.0*t6*t3*s+15.0/34.0*t3*t8+s*
1780 t8*t6/17+t12/476)*(t15*t+12.0/7.0*t17*t15+6.0/7.0*t*t19+4.0/35.0*t19*t17);
1781 /* p=13 2 5 2 */
1782 } else if( i==2 && j==5 && k==2 ) {
1783 t1 = r*r;
1784 t3 = s*s;
1785 t4 = t3*t3;
1786 t6 = r-1.0;
1787 t9 = t6*t6;
1788 t13 = t9*t9;
1789 t18 = t*t;
1790 t19 = s+r-1.0;
1791 t21 = t19*t19;
1792 f = (t1-6.0/23.0*r+3.0/253.0)*(t4*s+30.0/17.0*t4*t6+75.0/68.0*t3*s*t9+
1793 5.0/17.0*t3*t9*t6+15.0/476.0*s*t13+3.0/3094.0*t13*t6)*(t18+6.0/5.0*t*t19+3.0/
1794 10.0*t21);
1795 /* p=13 2 6 1 */
1796 } else if( i==2 && j==6 && k==1 ) {
1797 t1 = r*r;
1798 t3 = s*s;
1799 t4 = t3*t3;
1800 t7 = r-1.0;
1801 t9 = t7*t7;
1802 t14 = t9*t9;
1803 f = (t1-6.0/23.0*r+3.0/253.0)*(t4*t3+42.0/17.0*t4*s*t7+315.0/136.0*t4*
1804 t9+35.0/34.0*t3*s*t9*t7+15.0/68.0*t3*t14+9.0/442.0*s*t14*t7+t14*t9/1768)*(t+2.0
1805 /3.0*s+2.0/3.0*r-2.0/3.0);
1806 /* p=13 2 7 0 */
1807 } else if( i==2 && j==7 && k==0 ) {
1808 t1 = r*r;
1809 t3 = s*s;
1810 t4 = t3*s;
1811 t5 = t3*t3;
1812 t8 = r-1.0;
1813 t11 = t8*t8;
1814 t13 = t11*t8;
1815 t15 = t11*t11;
1816 f = (t1-6.0/23.0*r+3.0/253.0)*(t5*t4+56.0/17.0*t5*t3*t8+147.0/34.0*t5*s
1817 *t11+49.0/17.0*t5*t13+35.0/34.0*t4*t15+42.0/221.0*t3*t15*t8+7.0/442.0*s*t15*t11
1818 +t15*t13/2431);
1819 /* p=13 3 0 6 */
1820 } else if( i==3 && j==0 && k==6 ) {
1821 t1 = r*r;
1822 t4 = t*t;
1823 t5 = t4*t4;
1824 t8 = s+r-1.0;
1825 t10 = t8*t8;
1826 t15 = t10*t10;
1827 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t5*t4+42.0/13.0*t5*t*
1828 t8+105.0/26.0*t5*t10+350.0/143.0*t4*t*t10*t8+105.0/143.0*t4*t15+14.0/143.0*t*
1829 t15*t8+7.0/1716.0*t15*t10);
1830 /* p=13 3 1 5 */
1831 } else if( i==3 && j==1 && k==5 ) {
1832 t1 = r*r;
1833 t6 = t*t;
1834 t7 = t6*t6;
1835 t9 = s+r-1.0;
1836 t12 = t9*t9;
1837 t16 = t12*t12;
1838 f = (s+2.0/15.0*r-2.0/15.0)*(t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)
1839 *(t7*t+30.0/11.0*t9*t7+30.0/11.0*t6*t*t12+40.0/33.0*t6*t12*t9+5.0/22.0*t*t16+
1840 t16*t9/77);
1841 /* p=13 3 2 4 */
1842 } else if( i==3 && j==2 && k==4 ) {
1843 t1 = r*r;
1844 t4 = s*s;
1845 t5 = r-1.0;
1846 t7 = t5*t5;
1847 t10 = t*t;
1848 t11 = t10*t10;
1849 t13 = s+r-1.0;
1850 t15 = t13*t13;
1851 t19 = t15*t15;
1852 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t4+2.0/5.0*s*t5+t7/35)
1853 *(t11+20.0/9.0*t10*t*t13+5.0/3.0*t15*t10+10.0/21.0*t*t15*t13+5.0/126.0*t19);
1854 /* p=13 3 3 3 */
1855 } else if( i==3 && j==3 && k==3 ) {
1856 t1 = r*r;
1857 t4 = s*s;
1858 t6 = r-1.0;
1859 t8 = t6*t6;
1860 t13 = t*t;
1861 t15 = s+r-1.0;
1862 t17 = t15*t15;
1863 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t4*s+4.0/5.0*t4*t6+6.0
1864 /35.0*s*t8+4.0/455.0*t8*t6)*(t13*t+12.0/7.0*t15*t13+6.0/7.0*t*t17+4.0/35.0*t17*
1865 t15);
1866 /* p=13 3 4 2 */
1867 } else if( i==3 && j==4 && k==2 ) {
1868 t1 = r*r;
1869 t4 = s*s;
1870 t5 = t4*t4;
1871 t7 = r-1.0;
1872 t9 = t7*t7;
1873 t13 = t9*t9;
1874 t16 = t*t;
1875 t17 = s+r-1.0;
1876 t19 = t17*t17;
1877 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t5+4.0/3.0*t4*s*t7+4.0
1878 /7.0*t4*t9+8.0/91.0*s*t9*t7+t13/273)*(t16+6.0/5.0*t*t17+3.0/10.0*t19);
1879 /* p=13 3 5 1 */
1880 } else if( i==3 && j==5 && k==1 ) {
1881 t1 = r*r;
1882 t4 = s*s;
1883 t5 = t4*t4;
1884 t7 = r-1.0;
1885 t10 = t7*t7;
1886 t14 = t10*t10;
1887 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t5*s+2.0*t5*t7+10.0/
1888 7.0*t4*s*t10+40.0/91.0*t4*t10*t7+5.0/91.0*s*t14+2.0/1001.0*t14*t7)*(t+2.0/3.0*s
1889 +2.0/3.0*r-2.0/3.0);
1890 /* p=13 3 6 0 */
1891 } else if( i==3 && j==6 && k==0 ) {
1892 t1 = r*r;
1893 t4 = s*s;
1894 t5 = t4*t4;
1895 t8 = r-1.0;
1896 t10 = t8*t8;
1897 t15 = t10*t10;
1898 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t5*t4+14.0/5.0*t5*s*t8
1899 +3.0*t5*t10+20.0/13.0*t4*s*t10*t8+5.0/13.0*t4*t15+6.0/143.0*s*t15*t8+t15*t10
1900 /715);
1901 /* p=13 4 0 5 */
1902 } else if( i==4 && j==0 && k==5 ) {
1903 t1 = r*r;
1904 t2 = t1*t1;
1905 t5 = t*t;
1906 t6 = t5*t5;
1907 t8 = s+r-1.0;
1908 t11 = t8*t8;
1909 t15 = t11*t11;
1910 f = (t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0)*(t6*t+
1911 30.0/11.0*t6*t8+30.0/11.0*t5*t*t11+40.0/33.0*t5*t11*t8+5.0/22.0*t*t15+t15*t8/77
1912 );
1913 /* p=13 4 1 4 */
1914 } else if( i==4 && j==1 && k==4 ) {
1915 t1 = r*r;
1916 t2 = t1*t1;
1917 t7 = t*t;
1918 t8 = t7*t7;
1919 t10 = s+r-1.0;
1920 t12 = t10*t10;
1921 t16 = t12*t12;
1922 f = (s+2.0/13.0*r-2.0/13.0)*(t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/
1923 1771.0*r+1.0/1771.0)*(t8+20.0/9.0*t7*t*t10+5.0/3.0*t7*t12+10.0/21.0*t*t12*t10+
1924 5.0/126.0*t16);
1925 /* p=13 4 2 3 */
1926 } else if( i==4 && j==2 && k==3 ) {
1927 t1 = r*r;
1928 t2 = t1*t1;
1929 t5 = s*s;
1930 t6 = r-1.0;
1931 t8 = t6*t6;
1932 t11 = t*t;
1933 t13 = s+r-1.0;
1934 t15 = t13*t13;
1935 f = (t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0)*(t5+6.0/
1936 13.0*s*t6+t8/26)*(t11*t+12.0/7.0*t11*t13+6.0/7.0*t*t15+4.0/35.0*t15*t13);
1937 /* p=13 4 3 2 */
1938 } else if( i==4 && j==3 && k==2 ) {
1939 t1 = r*r;
1940 t2 = t1*t1;
1941 t5 = s*s;
1942 t7 = r-1.0;
1943 t9 = t7*t7;
1944 t14 = t*t;
1945 t15 = s+r-1.0;
1946 t17 = t15*t15;
1947 f = (t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0)*(t5*s+
1948 12.0/13.0*t5*t7+3.0/13.0*s*t9+2.0/143.0*t9*t7)*(t14+6.0/5.0*t*t15+3.0/10.0*t17)
1949 ;
1950 /* p=13 4 4 1 */
1951 } else if( i==4 && j==4 && k==1 ) {
1952 t1 = r*r;
1953 t2 = t1*t1;
1954 t5 = s*s;
1955 t6 = t5*t5;
1956 t8 = r-1.0;
1957 t10 = t8*t8;
1958 t14 = t10*t10;
1959 f = (t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0)*(t6+20.0
1960 /13.0*t5*s*t8+10.0/13.0*t5*t10+20.0/143.0*s*t10*t8+t14/143)*(t+2.0/3.0*s+2.0/
1961 3.0*r-2.0/3.0);
1962 /* p=13 4 5 0 */
1963 } else if( i==4 && j==5 && k==0 ) {
1964 t1 = r*r;
1965 t2 = t1*t1;
1966 t5 = s*s;
1967 t6 = t5*t5;
1968 t8 = r-1.0;
1969 t11 = t8*t8;
1970 t15 = t11*t11;
1971 f = (t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0)*(t6*s+
1972 30.0/13.0*t6*t8+25.0/13.0*t5*s*t11+100.0/143.0*t5*t11*t8+15.0/143.0*s*t15+2.0/
1973 429.0*t15*t8);
1974 /* p=13 5 0 4 */
1975 } else if( i==5 && j==0 && k==4 ) {
1976 t1 = r*r;
1977 t2 = t1*t1;
1978 t6 = t*t;
1979 t7 = t6*t6;
1980 t9 = s+r-1.0;
1981 t11 = t9*t9;
1982 t15 = t11*t11;
1983 f = (t2*r-30.0/23.0*t2+150.0/253.0*t1*r-200.0/1771.0*t1+15.0/1771.0*r
1984 -6.0/33649.0)*(t7+20.0/9.0*t6*t*t9+5.0/3.0*t6*t11+10.0/21.0*t*t11*t9+5.0/126.0*
1985 t15);
1986 /* p=13 5 1 3 */
1987 } else if( i==5 && j==1 && k==3 ) {
1988 t1 = r*r;
1989 t2 = t1*t1;
1990 t8 = t*t;
1991 t10 = s+r-1.0;
1992 t12 = t10*t10;
1993 f = (t2*r-30.0/23.0*t2+150.0/253.0*t1*r-200.0/1771.0*t1+15.0/1771.0*r
1994 -6.0/33649.0)*(s+2.0/11.0*r-2.0/11.0)*(t8*t+12.0/7.0*t8*t10+6.0/7.0*t*t12+4.0/
1995 35.0*t12*t10);
1996 /* p=13 5 2 2 */
1997 } else if( i==5 && j==2 && k==2 ) {
1998 t1 = r*r;
1999 t2 = t1*t1;
2000 t6 = s*s;
2001 t7 = r-1.0;
2002 t9 = t7*t7;
2003 t12 = t*t;
2004 t13 = s+r-1.0;
2005 t15 = t13*t13;
2006 f = (t2*r-30.0/23.0*t2+150.0/253.0*t1*r-200.0/1771.0*t1+15.0/1771.0*r
2007 -6.0/33649.0)*(t6+6.0/11.0*s*t7+3.0/55.0*t9)*(t12+6.0/5.0*t*t13+3.0/10.0*t15);
2008 /* p=13 5 3 1 */
2009 } else if( i==5 && j==3 && k==1 ) {
2010 t1 = r*r;
2011 t2 = t1*t1;
2012 t6 = s*s;
2013 t8 = r-1.0;
2014 t10 = t8*t8;
2015 f = (t2*r-30.0/23.0*t2+150.0/253.0*t1*r-200.0/1771.0*t1+15.0/1771.0*r
2016 -6.0/33649.0)*(t6*s+12.0/11.0*t6*t8+18.0/55.0*s*t10+4.0/165.0*t8*t10)*(t+2.0/
2017 3.0*s+2.0/3.0*r-2.0/3.0);
2018 /* p=13 5 4 0 */
2019 } else if( i==5 && j==4 && k==0 ) {
2020 t1 = r*r;
2021 t2 = t1*t1;
2022 t6 = s*s;
2023 t7 = t6*t6;
2024 t9 = r-1.0;
2025 t11 = t9*t9;
2026 t15 = t11*t11;
2027 f = (t2*r-30.0/23.0*t2+150.0/253.0*t1*r-200.0/1771.0*t1+15.0/1771.0*r
2028 -6.0/33649.0)*(t7+20.0/11.0*t6*s*t9+12.0/11.0*t6*t11+8.0/33.0*s*t11*t9+t15/66);
2029 /* p=13 6 0 3 */
2030 } else if( i==6 && j==0 && k==3 ) {
2031 t1 = r*r;
2032 t2 = t1*t1;
2033 t7 = t*t;
2034 t9 = s+r-1.0;
2035 t11 = t9*t9;
2036 f = (t2*t1-42.0/23.0*t2*r+315.0/253.0*t2-100.0/253.0*t1*r+15.0/253.0*t1
2037 -18.0/4807.0*r+1.0/14421.0)*(t7*t+12.0/7.0*t7*t9+6.0/7.0*t*t11+4.0/35.0*t11*t9)
2038 ;
2039 /* p=13 6 1 2 */
2040 } else if( i==6 && j==1 && k==2 ) {
2041 t1 = r*r;
2042 t2 = t1*t1;
2043 t9 = t*t;
2044 t10 = s+r-1.0;
2045 t12 = t10*t10;
2046 f = (t2*t1-42.0/23.0*t2*r+315.0/253.0*t2-100.0/253.0*t1*r+15.0/253.0*t1
2047 -18.0/4807.0*r+1.0/14421.0)*(s+2.0/9.0*r-2.0/9.0)*(t9+6.0/5.0*t*t10+3.0/10.0*
2048 t12);
2049 /* p=13 6 2 1 */
2050 } else if( i==6 && j==2 && k==1 ) {
2051 t1 = r*r;
2052 t2 = t1*t1;
2053 t7 = s*s;
2054 t8 = r-1.0;
2055 t10 = t8*t8;
2056 f = (t2*t1-42.0/23.0*t2*r+315.0/253.0*t2-100.0/253.0*t1*r+15.0/253.0*t1
2057 -18.0/4807.0*r+1.0/14421.0)*(t7+2.0/3.0*s*t8+t10/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0
2058 /3.0);
2059 /* p=13 6 3 0 */
2060 } else if( i==6 && j==3 && k==0 ) {
2061 t1 = r*r;
2062 t2 = t1*t1;
2063 t7 = s*s;
2064 t9 = r-1.0;
2065 t11 = t9*t9;
2066 f = (t2*t1-42.0/23.0*t2*r+315.0/253.0*t2-100.0/253.0*t1*r+15.0/253.0*t1
2067 -18.0/4807.0*r+1.0/14421.0)*(t7*s+4.0/3.0*t7*t9+s*t11/2+t11*t9/21);
2068 /* p=13 7 0 2 */
2069 } else if( i==7 && j==0 && k==2 ) {
2070 t1 = r*r;
2071 t2 = t1*r;
2072 t3 = t1*t1;
2073 t8 = t*t;
2074 t9 = s+r-1.0;
2075 t11 = t9*t9;
2076 f = (t3*t2-56.0/23.0*t3*t1+588.0/253.0*t3*r-280.0/253.0*t3+70.0/253.0*
2077 t2-168.0/4807.0*t1+28.0/14421.0*r-8.0/245157.0)*(t8+6.0/5.0*t*t9+3.0/10.0*t11);
2078 /* p=13 7 1 1 */
2079 } else if( i==7 && j==1 && k==1 ) {
2080 t1 = r*r;
2081 t2 = t1*r;
2082 t3 = t1*t1;
2083 f = (t3*t2-56.0/23.0*t3*t1+588.0/253.0*t3*r-280.0/253.0*t3+70.0/253.0*
2084 t2-168.0/4807.0*t1+28.0/14421.0*r-8.0/245157.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/
2085 3.0*s+2.0/3.0*r-2.0/3.0);
2086 /* p=13 7 2 0 */
2087 } else if( i==7 && j==2 && k==0 ) {
2088 t1 = r*r;
2089 t2 = t1*r;
2090 t3 = t1*t1;
2091 t8 = s*s;
2092 t9 = r-1.0;
2093 t11 = t9*t9;
2094 f = (t3*t2-56.0/23.0*t3*t1+588.0/253.0*t3*r-280.0/253.0*t3+70.0/253.0*
2095 t2-168.0/4807.0*t1+28.0/14421.0*r-8.0/245157.0)*(t8+6.0/7.0*s*t9+t11/7);
2096 /* p=13 8 0 1 */
2097 } else if( i==8 && j==0 && k==1 ) {
2098 t1 = r*r;
2099 t2 = t1*t1;
2100 t3 = t2*t2;
2101 t4 = t1*r;
2102 f = (t3-72.0/23.0*t2*t4+1008.0/253.0*t2*t1-672.0/253.0*t2*r+252.0/253.0
2103 *t2-1008.0/4807.0*t4+112.0/4807.0*t1-96.0/81719.0*r+3.0/163438.0)*(t+2.0/3.0*s+
2104 2.0/3.0*r-2.0/3.0);
2105 /* p=13 8 1 0 */
2106 } else if( i==8 && j==1 && k==0 ) {
2107 t1 = r*r;
2108 t2 = t1*t1;
2109 t3 = t2*t2;
2110 t4 = t1*r;
2111 f = (t3-72.0/23.0*t2*t4+1008.0/253.0*t2*t1-672.0/253.0*t2*r+252.0/253.0
2112 *t2-1008.0/4807.0*t4+112.0/4807.0*t1-96.0/81719.0*r+3.0/163438.0)*(s+2.0/5.0*r
2113 -2.0/5.0);
2114 /* p=13 9 0 0 */
2115 } else if( i==9 && j==0 && k==0 ) {
2116 t1 = r*r;
2117 t2 = t1*t1;
2118 t3 = t2*t2;
2119 t5 = t1*r;
2120 f = t3*r-90.0/23.0*t3+1620.0/253.0*t2*t5-1440.0/253.0*t2*t1+756.0/253.0*
2121 t2*r-4536.0/4807.0*t2+840.0/4807.0*t5-1440.0/81719.0*t1+135.0/163438.0*r-1.0/
2122 81719.0;
2123 /* p=14 0 0 10 */
2124 } else if( i==0 && j==0 && k==10 ) {
2125 t1 = t*t;
2126 t2 = t1*t1;
2127 t3 = t2*t2;
2128 t6 = s+r-1.0;
2129 t8 = t6*t6;
2130 t10 = t1*t;
2131 t12 = t8*t6;
2132 t15 = t8*t8;
2133 t24 = t15*t15;
2134 f = t3*t1+110.0/21.0*t3*t*t6+165.0/14.0*t3*t8+1980.0/133.0*t2*t10*t12+
2135 220.0/19.0*t2*t1*t15+1848.0/323.0*t2*t*t15*t6+1155.0/646.0*t2*t15*t8+110.0/
2136 323.0*t10*t15*t12+165.0/4522.0*t1*t24+55.0/29393.0*t*t24*t6+11.0/352716.0*t24*
2137 t8;
2138 /* p=14 0 1 9 */
2139 } else if( i==0 && j==1 && k==9 ) {
2140 t2 = t*t;
2141 t3 = t2*t2;
2142 t4 = t3*t3;
2143 t6 = s+r-1.0;
2144 t8 = t2*t;
2145 t10 = t6*t6;
2146 t13 = t10*t6;
2147 t16 = t10*t10;
2148 t24 = t16*t16;
2149 f = (s+2.0/23.0*r-2.0/23.0)*(t4*t+90.0/19.0*t4*t6+180.0/19.0*t3*t8*t10+
2150 3360.0/323.0*t3*t2*t13+2205.0/323.0*t3*t*t16+882.0/323.0*t3*t16*t6+210.0/323.0*
2151 t8*t16*t10+360.0/4199.0*t2*t16*t13+45.0/8398.0*t*t24+5.0/46189.0*t24*t6);
2152 /* p=14 0 2 8 */
2153 } else if( i==0 && j==2 && k==8 ) {
2154 t1 = s*s;
2155 t2 = r-1.0;
2156 t4 = t2*t2;
2157 t6 = t*t;
2158 t7 = t6*t6;
2159 t8 = t7*t7;
2160 t9 = t6*t;
2161 t11 = s+r-1.0;
2162 t14 = t11*t11;
2163 t17 = t14*t11;
2164 t19 = t14*t14;
2165 t27 = t19*t19;
2166 f = (t1+6.0/23.0*s*t2+3.0/253.0*t4)*(t8+72.0/17.0*t7*t9*t11+126.0/17.0*
2167 t7*t6*t14+588.0/85.0*t7*t*t17+63.0/17.0*t7*t19+252.0/221.0*t9*t19*t11+42.0/
2168 221.0*t6*t19*t14+36.0/2431.0*t*t19*t17+9.0/24310.0*t27);
2169 /* p=14 0 3 7 */
2170 } else if( i==0 && j==3 && k==7 ) {
2171 t1 = s*s;
2172 t3 = r-1.0;
2173 t5 = t3*t3;
2174 t9 = t*t;
2175 t10 = t9*t;
2176 t11 = t9*t9;
2177 t14 = s+r-1.0;
2178 t17 = t14*t14;
2179 t19 = t17*t14;
2180 t21 = t17*t17;
2181 f = (t1*s+12.0/23.0*t3*t1+18.0/253.0*s*t5+4.0/1771.0*t5*t3)*(t10*t11+
2182 56.0/15.0*t11*t9*t14+28.0/5.0*t11*t*t17+56.0/13.0*t19*t11+70.0/39.0*t10*t21+
2183 56.0/143.0*t9*t21*t14+28.0/715.0*t*t21*t17+8.0/6435.0*t21*t19);
2184 /* p=14 0 4 6 */
2185 } else if( i==0 && j==4 && k==6 ) {
2186 t1 = s*s;
2187 t2 = t1*t1;
2188 t4 = r-1.0;
2189 t6 = t4*t4;
2190 t10 = t6*t6;
2191 t12 = t*t;
2192 t13 = t12*t12;
2193 t16 = s+r-1.0;
2194 t18 = t16*t16;
2195 t23 = t18*t18;
2196 f = (t2+20.0/23.0*t1*s*t4+60.0/253.0*t1*t6+40.0/1771.0*s*t4*t6+t10/1771
2197 )*(t12*t13+42.0/13.0*t13*t*t16+105.0/26.0*t13*t18+350.0/143.0*t12*t*t18*t16+
2198 105.0/143.0*t12*t23+14.0/143.0*t*t23*t16+7.0/1716.0*t23*t18);
2199 /* p=14 0 5 5 */
2200 } else if( i==0 && j==5 && k==5 ) {
2201 t1 = s*s;
2202 t2 = t1*t1;
2203 t4 = r-1.0;
2204 t7 = t4*t4;
2205 t11 = t7*t7;
2206 t15 = t*t;
2207 t16 = t15*t15;
2208 t18 = s+r-1.0;
2209 t21 = t18*t18;
2210 t25 = t21*t21;
2211 f = (s*t2+30.0/23.0*t2*t4+150.0/253.0*t7*t1*s+200.0/1771.0*t1*t7*t4+
2212 15.0/1771.0*s*t11+6.0/33649.0*t11*t4)*(t16*t+30.0/11.0*t18*t16+30.0/11.0*t15*t*
2213 t21+40.0/33.0*t15*t21*t18+5.0/22.0*t*t25+t25*t18/77);
2214 /* p=14 0 6 4 */
2215 } else if( i==0 && j==6 && k==4 ) {
2216 t1 = s*s;
2217 t2 = t1*t1;
2218 t5 = r-1.0;
2219 t7 = t5*t5;
2220 t12 = t7*t7;
2221 t18 = t*t;
2222 t19 = t18*t18;
2223 t21 = s+r-1.0;
2224 t23 = t21*t21;
2225 t27 = t23*t23;
2226 f = (t2*t1+42.0/23.0*t2*s*t5+315.0/253.0*t2*t7+100.0/253.0*t1*s*t7*t5+
2227 15.0/253.0*t1*t12+18.0/4807.0*s*t12*t5+t12*t7/14421)*(t19+20.0/9.0*t18*t*t21+
2228 5.0/3.0*t23*t18+10.0/21.0*t*t23*t21+5.0/126.0*t27);
2229 /* p=14 0 7 3 */
2230 } else if( i==0 && j==7 && k==3 ) {
2231 t1 = s*s;
2232 t2 = t1*s;
2233 t3 = t1*t1;
2234 t6 = r-1.0;
2235 t9 = t6*t6;
2236 t11 = t9*t6;
2237 t13 = t9*t9;
2238 t21 = t*t;
2239 t23 = s+r-1.0;
2240 t25 = t23*t23;
2241 f = (t3*t2+56.0/23.0*t3*t1*t6+588.0/253.0*t3*s*t9+280.0/253.0*t3*t11+
2242 70.0/253.0*t2*t13+168.0/4807.0*t1*t13*t6+28.0/14421.0*s*t13*t9+8.0/245157.0*t13
2243 *t11)*(t21*t+12.0/7.0*t21*t23+6.0/7.0*t*t25+4.0/35.0*t25*t23);
2244 /* p=14 0 8 2 */
2245 } else if( i==0 && j==8 && k==2 ) {
2246 t1 = s*s;
2247 t2 = t1*t1;
2248 t3 = t2*t2;
2249 t4 = t1*s;
2250 t6 = r-1.0;
2251 t9 = t6*t6;
2252 t12 = t9*t6;
2253 t14 = t9*t9;
2254 t22 = t14*t14;
2255 t24 = t*t;
2256 t25 = s+r-1.0;
2257 t27 = t25*t25;
2258 f = (t3+72.0/23.0*t2*t4*t6+1008.0/253.0*t2*t1*t9+672.0/253.0*t2*s*t12+
2259 252.0/253.0*t2*t14+1008.0/4807.0*t4*t14*t6+112.0/4807.0*t1*t14*t9+96.0/81719.0*
2260 s*t14*t12+3.0/163438.0*t22)*(t24+6.0/5.0*t*t25+3.0/10.0*t27);
2261 /* p=14 0 9 1 */
2262 } else if( i==0 && j==9 && k==1 ) {
2263 t1 = s*s;
2264 t2 = t1*t1;
2265 t3 = t2*t2;
2266 t5 = r-1.0;
2267 t7 = t1*s;
2268 t9 = t5*t5;
2269 t12 = t9*t5;
2270 t15 = t9*t9;
2271 t23 = t15*t15;
2272 f = (t3*s+90.0/23.0*t3*t5+1620.0/253.0*t2*t7*t9+1440.0/253.0*t2*t1*t12+
2273 756.0/253.0*t2*s*t15+4536.0/4807.0*t2*t15*t5+840.0/4807.0*t7*t15*t9+1440.0/
2274 81719.0*t1*t15*t12+135.0/163438.0*s*t23+t23*t5/81719)*(t+2.0/3.0*s+2.0/3.0*r
2275 -2.0/3.0);
2276 /* p=14 0 10 0 */
2277 } else if( i==0 && j==10 && k==0 ) {
2278 t1 = s*s;
2279 t2 = t1*t1;
2280 t3 = t2*t2;
2281 t6 = r-1.0;
2282 t8 = t6*t6;
2283 t10 = t1*s;
2284 t12 = t8*t6;
2285 t15 = t8*t8;
2286 t24 = t15*t15;
2287 f = t3*t1+110.0/23.0*t3*s*t6+225.0/23.0*t3*t8+1800.0/161.0*t2*t10*t12+
2288 180.0/23.0*t2*t1*t15+1512.0/437.0*t2*s*t15*t6+420.0/437.0*t2*t15*t8+1200.0/
2289 7429.0*t10*t15*t12+225.0/14858.0*t1*t24+5.0/7429.0*s*t24*t6+t24*t8/104006;
2290 /* p=14 1 0 9 */
2291 } else if( i==1 && j==0 && k==9 ) {
2292 t2 = t*t;
2293 t3 = t2*t2;
2294 t4 = t3*t3;
2295 t6 = s+r-1.0;
2296 t8 = t2*t;
2297 t10 = t6*t6;
2298 t13 = t10*t6;
2299 t16 = t10*t10;
2300 t24 = t16*t16;
2301 f = (r-2.0/25.0)*(t4*t+90.0/19.0*t4*t6+180.0/19.0*t3*t8*t10+3360.0/
2302 323.0*t3*t2*t13+2205.0/323.0*t3*t*t16+882.0/323.0*t3*t16*t6+210.0/323.0*t8*t16*
2303 t10+360.0/4199.0*t2*t16*t13+45.0/8398.0*t*t24+5.0/46189.0*t24*t6);
2304 /* p=14 1 1 8 */
2305 } else if( i==1 && j==1 && k==8 ) {
2306 t4 = t*t;
2307 t5 = t4*t4;
2308 t6 = t5*t5;
2309 t7 = t4*t;
2310 t9 = s+r-1.0;
2311 t12 = t9*t9;
2312 t15 = t12*t9;
2313 t17 = t12*t12;
2314 t25 = t17*t17;
2315 f = (s+2.0/21.0*r-2.0/21.0)*(r-2.0/25.0)*(t6+72.0/17.0*t5*t7*t9+126.0/
2316 17.0*t5*t4*t12+588.0/85.0*t5*t*t15+63.0/17.0*t5*t17+252.0/221.0*t7*t17*t9+42.0/
2317 221.0*t4*t17*t12+36.0/2431.0*t*t17*t15+9.0/24310.0*t25);
2318 /* p=14 1 2 7 */
2319 } else if( i==1 && j==2 && k==7 ) {
2320 t2 = s*s;
2321 t3 = r-1.0;
2322 t5 = t3*t3;
2323 t8 = t*t;
2324 t9 = t8*t;
2325 t10 = t8*t8;
2326 t13 = s+r-1.0;
2327 t16 = t13*t13;
2328 t18 = t16*t13;
2329 t20 = t16*t16;
2330 f = (r-2.0/25.0)*(t2+2.0/7.0*s*t3+t5/70)*(t10*t9+56.0/15.0*t10*t8*t13+
2331 28.0/5.0*t10*t*t16+56.0/13.0*t10*t18+70.0/39.0*t9*t20+56.0/143.0*t8*t20*t13+
2332 28.0/715.0*t*t20*t16+8.0/6435.0*t20*t18);
2333 /* p=14 1 3 6 */
2334 } else if( i==1 && j==3 && k==6 ) {
2335 t2 = s*s;
2336 t4 = r-1.0;
2337 t6 = t4*t4;
2338 t11 = t*t;
2339 t12 = t11*t11;
2340 t15 = s+r-1.0;
2341 t17 = t15*t15;
2342 t22 = t17*t17;
2343 f = (r-2.0/25.0)*(t2*s+4.0/7.0*t2*t4+3.0/35.0*s*t6+2.0/665.0*t6*t4)*(
2344 t11*t12+42.0/13.0*t12*t*t15+105.0/26.0*t17*t12+350.0/143.0*t11*t*t17*t15+105.0/
2345 143.0*t11*t22+14.0/143.0*t*t22*t15+7.0/1716.0*t22*t17);
2346 /* p=14 1 4 5 */
2347 } else if( i==1 && j==4 && k==5 ) {
2348 t2 = s*s;
2349 t3 = t2*t2;
2350 t5 = r-1.0;
2351 t7 = t5*t5;
2352 t11 = t7*t7;
2353 t14 = t*t;
2354 t15 = t14*t14;
2355 t17 = s+r-1.0;
2356 t20 = t17*t17;
2357 t24 = t20*t20;
2358 f = (r-2.0/25.0)*(t3+20.0/21.0*t5*t2*s+2.0/7.0*t2*t7+4.0/133.0*s*t5*t7+
2359 t11/1197)*(t15*t+30.0/11.0*t17*t15+30.0/11.0*t14*t*t20+40.0/33.0*t14*t20*t17+
2360 5.0/22.0*t*t24+t24*t17/77);
2361 /* p=14 1 5 4 */
2362 } else if( i==1 && j==5 && k==4 ) {
2363 t2 = s*s;
2364 t3 = t2*t2;
2365 t5 = r-1.0;
2366 t8 = t5*t5;
2367 t12 = t8*t8;
2368 t17 = t*t;
2369 t18 = t17*t17;
2370 t20 = s+r-1.0;
2371 t22 = t20*t20;
2372 t26 = t22*t22;
2373 f = (r-2.0/25.0)*(s*t3+10.0/7.0*t3*t5+5.0/7.0*t2*s*t8+20.0/133.0*t2*t8*
2374 t5+5.0/399.0*s*t12+2.0/6783.0*t12*t5)*(t18+20.0/9.0*t17*t*t20+5.0/3.0*t22*t17+
2375 10.0/21.0*t*t22*t20+5.0/126.0*t26);
2376 /* p=14 1 6 3 */
2377 } else if( i==1 && j==6 && k==3 ) {
2378 t2 = s*s;
2379 t3 = t2*t2;
2380 t6 = r-1.0;
2381 t8 = t6*t6;
2382 t13 = t8*t8;
2383 t20 = t*t;
2384 t22 = s+r-1.0;
2385 t24 = t22*t22;
2386 f = (r-2.0/25.0)*(t3*t2+2.0*t3*s*t6+3.0/2.0*t3*t8+10.0/19.0*t2*s*t8*t6+
2387 5.0/57.0*t2*t13+2.0/323.0*s*t13*t6+t13*t8/7752)*(t20*t+12.0/7.0*t22*t20+6.0/7.0
2388 *t*t24+4.0/35.0*t24*t22);
2389 /* p=14 1 7 2 */
2390 } else if( i==1 && j==7 && k==2 ) {
2391 t2 = s*s;
2392 t3 = t2*s;
2393 t4 = t2*t2;
2394 t7 = r-1.0;
2395 t10 = t7*t7;
2396 t12 = t10*t7;
2397 t14 = t10*t10;
2398 t23 = t*t;
2399 t24 = s+r-1.0;
2400 t26 = t24*t24;
2401 f = (r-2.0/25.0)*(t4*t3+8.0/3.0*t2*t4*t7+14.0/5.0*t10*t4*s+28.0/19.0*t4
2402 *t12+70.0/171.0*t3*t14+56.0/969.0*t2*t14*t7+7.0/1938.0*s*t14*t10+t14*t12/14535)
2403 *(t23+6.0/5.0*t*t24+3.0/10.0*t26);
2404 /* p=14 1 8 1 */
2405 } else if( i==1 && j==8 && k==1 ) {
2406 t2 = s*s;
2407 t3 = t2*t2;
2408 t4 = t3*t3;
2409 t5 = t2*s;
2410 t7 = r-1.0;
2411 t10 = t7*t7;
2412 t13 = t10*t7;
2413 t15 = t10*t10;
2414 t23 = t15*t15;
2415 f = (r-2.0/25.0)*(t4+24.0/7.0*t3*t5*t7+24.0/5.0*t3*t2*t10+336.0/95.0*t3
2416 *s*t13+28.0/19.0*t3*t15+112.0/323.0*t5*t15*t7+14.0/323.0*t2*t15*t10+4.0/1615.0*
2417 s*t15*t13+t23/22610)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
2418 /* p=14 1 9 0 */
2419 } else if( i==1 && j==9 && k==0 ) {
2420 t2 = s*s;
2421 t3 = t2*t2;
2422 t4 = t3*t3;
2423 t6 = r-1.0;
2424 t8 = t2*s;
2425 t10 = t6*t6;
2426 t13 = t10*t6;
2427 t16 = t10*t10;
2428 t24 = t16*t16;
2429 f = (r-2.0/25.0)*(t4*s+30.0/7.0*t4*t6+54.0/7.0*t3*t8*t10+144.0/19.0*t3*
2430 t2*t13+84.0/19.0*t3*s*t16+504.0/323.0*t3*t16*t6+105.0/323.0*t8*t16*t10+12.0/
2431 323.0*t2*t16*t13+9.0/4522.0*s*t24+t24*t6/29393);
2432 /* p=14 2 0 8 */
2433 } else if( i==2 && j==0 && k==8 ) {
2434 t1 = r*r;
2435 t3 = t*t;
2436 t4 = t3*t3;
2437 t5 = t4*t4;
2438 t6 = t3*t;
2439 t8 = s+r-1.0;
2440 t11 = t8*t8;
2441 t14 = t11*t8;
2442 t16 = t11*t11;
2443 t24 = t16*t16;
2444 f = (t1-6.0/25.0*r+1.0/100.0)*(t5+72.0/17.0*t4*t6*t8+126.0/17.0*t4*t3*
2445 t11+588.0/85.0*t4*t*t14+63.0/17.0*t4*t16+252.0/221.0*t6*t16*t8+42.0/221.0*t3*
2446 t16*t11+36.0/2431.0*t*t16*t14+9.0/24310.0*t24);
2447 /* p=14 2 1 7 */
2448 } else if( i==2 && j==1 && k==7 ) {
2449 t1 = r*r;
2450 t5 = t*t;
2451 t6 = t5*t;
2452 t7 = t5*t5;
2453 t10 = s+r-1.0;
2454 t13 = t10*t10;
2455 t15 = t13*t10;
2456 t17 = t13*t13;
2457 f = (s+2.0/19.0*r-2.0/19.0)*(t1-6.0/25.0*r+1.0/100.0)*(t6*t7+56.0/15.0*
2458 t7*t5*t10+28.0/5.0*t7*t*t13+56.0/13.0*t15*t7+70.0/39.0*t6*t17+56.0/143.0*t5*t17
2459 *t10+28.0/715.0*t*t17*t13+8.0/6435.0*t17*t15);
2460 /* p=14 2 2 6 */
2461 } else if( i==2 && j==2 && k==6 ) {
2462 t1 = r*r;
2463 t3 = s*s;
2464 t4 = r-1.0;
2465 t6 = t4*t4;
2466 t9 = t*t;
2467 t10 = t9*t9;
2468 t13 = s+r-1.0;
2469 t15 = t13*t13;
2470 t20 = t15*t15;
2471 f = (t1-6.0/25.0*r+1.0/100.0)*(t3+6.0/19.0*t4*s+t6/57)*(t9*t10+42.0/
2472 13.0*t10*t*t13+105.0/26.0*t15*t10+350.0/143.0*t9*t*t15*t13+105.0/143.0*t9*t20+
2473 14.0/143.0*t*t20*t13+7.0/1716.0*t20*t15);
2474 /* p=14 2 3 5 */
2475 } else if( i==2 && j==3 && k==5 ) {
2476 t1 = r*r;
2477 t3 = s*s;
2478 t5 = r-1.0;
2479 t7 = t5*t5;
2480 t12 = t*t;
2481 t13 = t12*t12;
2482 t15 = s+r-1.0;
2483 t18 = t15*t15;
2484 t22 = t18*t18;
2485 f = (t1-6.0/25.0*r+1.0/100.0)*(t3*s+12.0/19.0*t3*t5+2.0/19.0*s*t7+4.0/
2486 969.0*t7*t5)*(t13*t+30.0/11.0*t15*t13+30.0/11.0*t12*t*t18+40.0/33.0*t12*t18*t15
2487 +5.0/22.0*t*t22+t22*t15/77);
2488 /* p=14 2 4 4 */
2489 } else if( i==2 && j==4 && k==4 ) {
2490 t1 = r*r;
2491 t3 = s*s;
2492 t4 = t3*t3;
2493 t6 = r-1.0;
2494 t8 = t6*t6;
2495 t12 = t8*t8;
2496 t15 = t*t;
2497 t16 = t15*t15;
2498 t18 = s+r-1.0;
2499 t20 = t18*t18;
2500 t24 = t20*t20;
2501 f = (t1-6.0/25.0*r+1.0/100.0)*(t4+20.0/19.0*t3*s*t6+20.0/57.0*t3*t8+
2502 40.0/969.0*s*t8*t6+5.0/3876.0*t12)*(t16+20.0/9.0*t15*t*t18+5.0/3.0*t15*t20+10.0
2503 /21.0*t*t20*t18+5.0/126.0*t24);
2504 /* p=14 2 5 3 */
2505 } else if( i==2 && j==5 && k==3 ) {
2506 t1 = r*r;
2507 t3 = s*s;
2508 t4 = t3*t3;
2509 t6 = r-1.0;
2510 t9 = t6*t6;
2511 t13 = t9*t9;
2512 t18 = t*t;
2513 t20 = s+r-1.0;
2514 t22 = t20*t20;
2515 f = (t1-6.0/25.0*r+1.0/100.0)*(t4*s+30.0/19.0*t4*t6+50.0/57.0*t3*s*t9+
2516 200.0/969.0*t3*t9*t6+25.0/1292.0*s*t13+t13*t6/1938)*(t18*t+12.0/7.0*t20*t18+6.0
2517 /7.0*t*t22+4.0/35.0*t22*t20);
2518 /* p=14 2 6 2 */
2519 } else if( i==2 && j==6 && k==2 ) {
2520 t1 = r*r;
2521 t3 = s*s;
2522 t4 = t3*t3;
2523 t7 = r-1.0;
2524 t9 = t7*t7;
2525 t14 = t9*t9;
2526 t21 = t*t;
2527 t22 = s+r-1.0;
2528 t24 = t22*t22;
2529 f = (t1-6.0/25.0*r+1.0/100.0)*(t4*t3+42.0/19.0*t4*s*t7+35.0/19.0*t4*t9+
2530 700.0/969.0*t3*s*t9*t7+175.0/1292.0*t3*t14+7.0/646.0*s*t14*t7+t14*t9/3876)*(t21
2531 +6.0/5.0*t*t22+3.0/10.0*t24);
2532 /* p=14 2 7 1 */
2533 } else if( i==2 && j==7 && k==1 ) {
2534 t1 = r*r;
2535 t3 = s*s;
2536 t4 = t3*s;
2537 t5 = t3*t3;
2538 t8 = r-1.0;
2539 t11 = t8*t8;
2540 t13 = t11*t8;
2541 t15 = t11*t11;
2542 f = (t1-6.0/25.0*r+1.0/100.0)*(t5*t4+56.0/19.0*t5*t3*t8+196.0/57.0*t5*s
2543 *t11+1960.0/969.0*t5*t13+1225.0/1938.0*t4*t15+98.0/969.0*t3*t15*t8+7.0/969.0*s*
2544 t15*t11+2.0/12597.0*t15*t13)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
2545 /* p=14 2 8 0 */
2546 } else if( i==2 && j==8 && k==0 ) {
2547 t1 = r*r;
2548 t3 = s*s;
2549 t4 = t3*t3;
2550 t5 = t4*t4;
2551 t6 = t3*s;
2552 t8 = r-1.0;
2553 t11 = t8*t8;
2554 t14 = t11*t8;
2555 t16 = t11*t11;
2556 t24 = t16*t16;
2557 f = (t1-6.0/25.0*r+1.0/100.0)*(t5+72.0/19.0*t4*t6*t8+112.0/19.0*t4*t3*
2558 t11+1568.0/323.0*t4*s*t14+735.0/323.0*t4*t16+196.0/323.0*t6*t16*t8+28.0/323.0*
2559 t3*t16*t11+24.0/4199.0*s*t16*t14+t24/8398);
2560 /* p=14 3 0 7 */
2561 } else if( i==3 && j==0 && k==7 ) {
2562 t1 = r*r;
2563 t4 = t*t;
2564 t5 = t4*t;
2565 t6 = t4*t4;
2566 t9 = s+r-1.0;
2567 t12 = t9*t9;
2568 t14 = t12*t9;
2569 t16 = t12*t12;
2570 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t5*t6+56.0/15.0*t4*t6*t9+
2571 28.0/5.0*t6*t*t12+56.0/13.0*t6*t14+70.0/39.0*t5*t16+56.0/143.0*t4*t16*t9+28.0/
2572 715.0*t*t16*t12+8.0/6435.0*t16*t14);
2573 /* p=14 3 1 6 */
2574 } else if( i==3 && j==1 && k==6 ) {
2575 t1 = r*r;
2576 t6 = t*t;
2577 t7 = t6*t6;
2578 t10 = s+r-1.0;
2579 t12 = t10*t10;
2580 t17 = t12*t12;
2581 f = (s+2.0/17.0*r-2.0/17.0)*(t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(
2582 t6*t7+42.0/13.0*t7*t*t10+105.0/26.0*t7*t12+350.0/143.0*t6*t*t12*t10+105.0/143.0
2583 *t6*t17+14.0/143.0*t*t17*t10+7.0/1716.0*t17*t12);
2584 /* p=14 3 2 5 */
2585 } else if( i==3 && j==2 && k==5 ) {
2586 t1 = r*r;
2587 t4 = s*s;
2588 t5 = r-1.0;
2589 t7 = t5*t5;
2590 t10 = t*t;
2591 t11 = t10*t10;
2592 t13 = s+r-1.0;
2593 t16 = t13*t13;
2594 t20 = t16*t16;
2595 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t4+6.0/17.0*s*t5+3.0/
2596 136.0*t7)*(t11*t+30.0/11.0*t11*t13+30.0/11.0*t10*t*t16+40.0/33.0*t10*t16*t13+
2597 5.0/22.0*t*t20+t20*t13/77);
2598 /* p=14 3 3 4 */
2599 } else if( i==3 && j==3 && k==4 ) {
2600 t1 = r*r;
2601 t4 = s*s;
2602 t6 = r-1.0;
2603 t8 = t6*t6;
2604 t13 = t*t;
2605 t14 = t13*t13;
2606 t16 = s+r-1.0;
2607 t18 = t16*t16;
2608 t22 = t18*t18;
2609 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t4*s+12.0/17.0*t4*t6+9.0/
2610 68.0*s*t8+t8*t6/170)*(t14+20.0/9.0*t13*t*t16+5.0/3.0*t13*t18+10.0/21.0*t*t18*
2611 t16+5.0/126.0*t22);
2612 /* p=14 3 4 3 */
2613 } else if( i==3 && j==4 && k==3 ) {
2614 t1 = r*r;
2615 t4 = s*s;
2616 t5 = t4*t4;
2617 t7 = r-1.0;
2618 t9 = t7*t7;
2619 t13 = t9*t9;
2620 t16 = t*t;
2621 t18 = s+r-1.0;
2622 t20 = t18*t18;
2623 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t5+20.0/17.0*t4*s*t7+15.0
2624 /34.0*t4*t9+s*t9*t7/17+t13/476)*(t16*t+12.0/7.0*t18*t16+6.0/7.0*t*t20+4.0/35.0*
2625 t20*t18);
2626 /* p=14 3 5 2 */
2627 } else if( i==3 && j==5 && k==2 ) {
2628 t1 = r*r;
2629 t4 = s*s;
2630 t5 = t4*t4;
2631 t7 = r-1.0;
2632 t10 = t7*t7;
2633 t14 = t10*t10;
2634 t19 = t*t;
2635 t20 = s+r-1.0;
2636 t22 = t20*t20;
2637 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(s*t5+30.0/17.0*t5*t7+75.0
2638 /68.0*t4*s*t10+5.0/17.0*t4*t10*t7+15.0/476.0*s*t14+3.0/3094.0*t14*t7)*(t19+6.0/
2639 5.0*t*t20+3.0/10.0*t22);
2640 /* p=14 3 6 1 */
2641 } else if( i==3 && j==6 && k==1 ) {
2642 t1 = r*r;
2643 t4 = s*s;
2644 t5 = t4*t4;
2645 t8 = r-1.0;
2646 t10 = t8*t8;
2647 t15 = t10*t10;
2648 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t5*t4+42.0/17.0*t5*s*t8+
2649 315.0/136.0*t5*t10+35.0/34.0*t4*s*t10*t8+15.0/68.0*t4*t15+9.0/442.0*s*t15*t8+
2650 t15*t10/1768)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
2651 /* p=14 3 7 0 */
2652 } else if( i==3 && j==7 && k==0 ) {
2653 t1 = r*r;
2654 t4 = s*s;
2655 t5 = t4*s;
2656 t6 = t4*t4;
2657 t9 = r-1.0;
2658 t12 = t9*t9;
2659 t14 = t12*t9;
2660 t16 = t12*t12;
2661 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t6*t5+56.0/17.0*t4*t6*t9+
2662 147.0/34.0*t6*s*t12+49.0/17.0*t6*t14+35.0/34.0*t5*t16+42.0/221.0*t4*t16*t9+7.0/
2663 442.0*s*t16*t12+t16*t14/2431);
2664 /* p=14 4 0 6 */
2665 } else if( i==4 && j==0 && k==6 ) {
2666 t1 = r*r;
2667 t2 = t1*t1;
2668 t5 = t*t;
2669 t6 = t5*t5;
2670 t9 = s+r-1.0;
2671 t11 = t9*t9;
2672 t16 = t11*t11;
2673 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(t6*t5+42.0/13.0*t6*t
2674 *t9+105.0/26.0*t6*t11+350.0/143.0*t5*t*t11*t9+105.0/143.0*t5*t16+14.0/143.0*t*
2675 t16*t9+7.0/1716.0*t16*t11);
2676 /* p=14 4 1 5 */
2677 } else if( i==4 && j==1 && k==5 ) {
2678 t1 = r*r;
2679 t2 = t1*t1;
2680 t7 = t*t;
2681 t8 = t7*t7;
2682 t10 = s+r-1.0;
2683 t13 = t10*t10;
2684 t17 = t13*t13;
2685 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(s+2.0/15.0*r-2.0/
2686 15.0)*(t8*t+30.0/11.0*t8*t10+30.0/11.0*t7*t*t13+40.0/33.0*t7*t13*t10+5.0/22.0*t
2687 *t17+t17*t10/77);
2688 /* p=14 4 2 4 */
2689 } else if( i==4 && j==2 && k==4 ) {
2690 t1 = r*r;
2691 t2 = t1*t1;
2692 t5 = s*s;
2693 t6 = r-1.0;
2694 t8 = t6*t6;
2695 t11 = t*t;
2696 t12 = t11*t11;
2697 t14 = s+r-1.0;
2698 t16 = t14*t14;
2699 t20 = t16*t16;
2700 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(t5+2.0/5.0*s*t6+t8
2701 /35)*(t12+20.0/9.0*t11*t*t14+5.0/3.0*t16*t11+10.0/21.0*t*t16*t14+5.0/126.0*t20)
2702 ;
2703 /* p=14 4 3 3 */
2704 } else if( i==4 && j==3 && k==3 ) {
2705 t1 = r*r;
2706 t2 = t1*t1;
2707 t5 = s*s;
2708 t7 = r-1.0;
2709 t9 = t7*t7;
2710 t14 = t*t;
2711 t16 = s+r-1.0;
2712 t18 = t16*t16;
2713 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(t5*s+4.0/5.0*t5*t7+
2714 6.0/35.0*s*t9+4.0/455.0*t9*t7)*(t14*t+12.0/7.0*t16*t14+6.0/7.0*t*t18+4.0/35.0*
2715 t18*t16);
2716 /* p=14 4 4 2 */
2717 } else if( i==4 && j==4 && k==2 ) {
2718 t1 = r*r;
2719 t2 = t1*t1;
2720 t5 = s*s;
2721 t6 = t5*t5;
2722 t8 = r-1.0;
2723 t10 = t8*t8;
2724 t14 = t10*t10;
2725 t17 = t*t;
2726 t18 = s+r-1.0;
2727 t20 = t18*t18;
2728 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(t6+4.0/3.0*t5*s*t8+
2729 4.0/7.0*t5*t10+8.0/91.0*s*t8*t10+t14/273)*(t17+6.0/5.0*t*t18+3.0/10.0*t20);
2730 /* p=14 4 5 1 */
2731 } else if( i==4 && j==5 && k==1 ) {
2732 t1 = r*r;
2733 t2 = t1*t1;
2734 t5 = s*s;
2735 t6 = t5*t5;
2736 t8 = r-1.0;
2737 t11 = t8*t8;
2738 t15 = t11*t11;
2739 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(s*t6+2.0*t6*t8+10.0/
2740 7.0*t5*s*t11+40.0/91.0*t5*t11*t8+5.0/91.0*s*t15+2.0/1001.0*t15*t8)*(t+2.0/3.0*s
2741 +2.0/3.0*r-2.0/3.0);
2742 /* p=14 4 6 0 */
2743 } else if( i==4 && j==6 && k==0 ) {
2744 t1 = r*r;
2745 t2 = t1*t1;
2746 t5 = s*s;
2747 t6 = t5*t5;
2748 t9 = r-1.0;
2749 t11 = t9*t9;
2750 t16 = t11*t11;
2751 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(t6*t5+14.0/5.0*t6*s*
2752 t9+3.0*t6*t11+20.0/13.0*t5*s*t11*t9+5.0/13.0*t5*t16+6.0/143.0*s*t16*t9+t16*t11
2753 /715);
2754 /* p=14 5 0 5 */
2755 } else if( i==5 && j==0 && k==5 ) {
2756 t1 = r*r;
2757 t2 = t1*t1;
2758 t6 = t*t;
2759 t7 = t6*t6;
2760 t9 = s+r-1.0;
2761 t12 = t9*t9;
2762 t16 = t12*t12;
2763 f = (t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0)*(t7*t+
2764 30.0/11.0*t7*t9+30.0/11.0*t6*t*t12+40.0/33.0*t6*t12*t9+5.0/22.0*t*t16+t16*t9/77
2765 );
2766 /* p=14 5 1 4 */
2767 } else if( i==5 && j==1 && k==4 ) {
2768 t1 = r*r;
2769 t2 = t1*t1;
2770 t8 = t*t;
2771 t9 = t8*t8;
2772 t11 = s+r-1.0;
2773 t13 = t11*t11;
2774 t17 = t13*t13;
2775 f = (s+2.0/13.0*r-2.0/13.0)*(t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/
2776 506.0*r-1.0/8855.0)*(t9+20.0/9.0*t8*t*t11+5.0/3.0*t8*t13+10.0/21.0*t*t13*t11+
2777 5.0/126.0*t17);
2778 /* p=14 5 2 3 */
2779 } else if( i==5 && j==2 && k==3 ) {
2780 t1 = r*r;
2781 t2 = t1*t1;
2782 t6 = s*s;
2783 t7 = r-1.0;
2784 t9 = t7*t7;
2785 t12 = t*t;
2786 t14 = s+r-1.0;
2787 t16 = t14*t14;
2788 f = (t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0)*(t6+6.0
2789 /13.0*s*t7+t9/26)*(t12*t+12.0/7.0*t12*t14+6.0/7.0*t*t16+4.0/35.0*t16*t14);
2790 /* p=14 5 3 2 */
2791 } else if( i==5 && j==3 && k==2 ) {
2792 t1 = r*r;
2793 t2 = t1*t1;
2794 t6 = s*s;
2795 t8 = r-1.0;
2796 t10 = t8*t8;
2797 t15 = t*t;
2798 t16 = s+r-1.0;
2799 t18 = t16*t16;
2800 f = (t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0)*(t6*s+
2801 12.0/13.0*t6*t8+3.0/13.0*s*t10+2.0/143.0*t10*t8)*(t15+6.0/5.0*t*t16+3.0/10.0*
2802 t18);
2803 /* p=14 5 4 1 */
2804 } else if( i==5 && j==4 && k==1 ) {
2805 t1 = r*r;
2806 t2 = t1*t1;
2807 t6 = s*s;
2808 t7 = t6*t6;
2809 t9 = r-1.0;
2810 t11 = t9*t9;
2811 t15 = t11*t11;
2812 f = (t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0)*(t7+
2813 20.0/13.0*t6*s*t9+10.0/13.0*t6*t11+20.0/143.0*s*t11*t9+t15/143)*(t+2.0/3.0*s+
2814 2.0/3.0*r-2.0/3.0);
2815 /* p=14 5 5 0 */
2816 } else if( i==5 && j==5 && k==0 ) {
2817 t1 = r*r;
2818 t2 = t1*t1;
2819 t6 = s*s;
2820 t7 = t6*t6;
2821 t9 = r-1.0;
2822 t12 = t9*t9;
2823 t16 = t12*t12;
2824 f = (t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0)*(s*t7+
2825 30.0/13.0*t7*t9+25.0/13.0*t6*s*t12+100.0/143.0*t6*t12*t9+15.0/143.0*s*t16+2.0/
2826 429.0*t16*t9);
2827 /* p=14 6 0 4 */
2828 } else if( i==6 && j==0 && k==4 ) {
2829 t1 = r*r;
2830 t2 = t1*t1;
2831 t7 = t*t;
2832 t8 = t7*t7;
2833 t10 = s+r-1.0;
2834 t12 = t10*t10;
2835 t16 = t12*t12;
2836 f = (t2*t1-42.0/25.0*t2*r+21.0/20.0*t2-7.0/23.0*t1*r+21.0/506.0*t1-3.0/
2837 1265.0*r+1.0/25300.0)*(t8+20.0/9.0*t7*t*t10+5.0/3.0*t7*t12+10.0/21.0*t*t12*t10+
2838 5.0/126.0*t16);
2839 /* p=14 6 1 3 */
2840 } else if( i==6 && j==1 && k==3 ) {
2841 t1 = r*r;
2842 t2 = t1*t1;
2843 t9 = t*t;
2844 t11 = s+r-1.0;
2845 t13 = t11*t11;
2846 f = (t2*t1-42.0/25.0*t2*r+21.0/20.0*t2-7.0/23.0*t1*r+21.0/506.0*t1-3.0/
2847 1265.0*r+1.0/25300.0)*(s+2.0/11.0*r-2.0/11.0)*(t9*t+12.0/7.0*t11*t9+6.0/7.0*t*
2848 t13+4.0/35.0*t13*t11);
2849 /* p=14 6 2 2 */
2850 } else if( i==6 && j==2 && k==2 ) {
2851 t1 = r*r;
2852 t2 = t1*t1;
2853 t7 = s*s;
2854 t8 = r-1.0;
2855 t10 = t8*t8;
2856 t13 = t*t;
2857 t14 = s+r-1.0;
2858 t16 = t14*t14;
2859 f = (t2*t1-42.0/25.0*t2*r+21.0/20.0*t2-7.0/23.0*t1*r+21.0/506.0*t1-3.0/
2860 1265.0*r+1.0/25300.0)*(t7+6.0/11.0*s*t8+3.0/55.0*t10)*(t13+6.0/5.0*t*t14+3.0/
2861 10.0*t16);
2862 /* p=14 6 3 1 */
2863 } else if( i==6 && j==3 && k==1 ) {
2864 t1 = r*r;
2865 t2 = t1*t1;
2866 t7 = s*s;
2867 t9 = r-1.0;
2868 t11 = t9*t9;
2869 f = (t2*t1-42.0/25.0*t2*r+21.0/20.0*t2-7.0/23.0*t1*r+21.0/506.0*t1-3.0/
2870 1265.0*r+1.0/25300.0)*(t7*s+12.0/11.0*t7*t9+18.0/55.0*s*t11+4.0/165.0*t11*t9)*(
2871 t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
2872 /* p=14 6 4 0 */
2873 } else if( i==6 && j==4 && k==0 ) {
2874 t1 = r*r;
2875 t2 = t1*t1;
2876 t7 = s*s;
2877 t8 = t7*t7;
2878 t10 = r-1.0;
2879 t12 = t10*t10;
2880 t16 = t12*t12;
2881 f = (t2*t1-42.0/25.0*t2*r+21.0/20.0*t2-7.0/23.0*t1*r+21.0/506.0*t1-3.0/
2882 1265.0*r+1.0/25300.0)*(t8+20.0/11.0*t7*s*t10+12.0/11.0*t7*t12+8.0/33.0*s*t12*
2883 t10+t16/66);
2884 /* p=14 7 0 3 */
2885 } else if( i==7 && j==0 && k==3 ) {
2886 t1 = r*r;
2887 t2 = t1*r;
2888 t3 = t1*t1;
2889 t8 = t*t;
2890 t10 = s+r-1.0;
2891 t12 = t10*t10;
2892 f = (t3*t2-56.0/25.0*t3*t1+49.0/25.0*t3*r-98.0/115.0*t3+49.0/253.0*t2
2893 -28.0/1265.0*t1+7.0/6325.0*r-2.0/120175.0)*(t8*t+12.0/7.0*t8*t10+6.0/7.0*t*t12+
2894 4.0/35.0*t12*t10);
2895 /* p=14 7 1 2 */
2896 } else if( i==7 && j==1 && k==2 ) {
2897 t1 = r*r;
2898 t2 = t1*r;
2899 t3 = t1*t1;
2900 t10 = t*t;
2901 t11 = s+r-1.0;
2902 t13 = t11*t11;
2903 f = (t3*t2-56.0/25.0*t3*t1+49.0/25.0*t3*r-98.0/115.0*t3+49.0/253.0*t2
2904 -28.0/1265.0*t1+7.0/6325.0*r-2.0/120175.0)*(s+2.0/9.0*r-2.0/9.0)*(t10+6.0/5.0*t
2905 *t11+3.0/10.0*t13);
2906 /* p=14 7 2 1 */
2907 } else if( i==7 && j==2 && k==1 ) {
2908 t1 = r*r;
2909 t2 = t1*r;
2910 t3 = t1*t1;
2911 t8 = s*s;
2912 t9 = r-1.0;
2913 t11 = t9*t9;
2914 f = (t3*t2-56.0/25.0*t3*t1+49.0/25.0*t3*r-98.0/115.0*t3+49.0/253.0*t2
2915 -28.0/1265.0*t1+7.0/6325.0*r-2.0/120175.0)*(t8+2.0/3.0*s*t9+t11/12)*(t+2.0/3.0*
2916 s+2.0/3.0*r-2.0/3.0);
2917 /* p=14 7 3 0 */
2918 } else if( i==7 && j==3 && k==0 ) {
2919 t1 = r*r;
2920 t2 = t1*r;
2921 t3 = t1*t1;
2922 t8 = s*s;
2923 t10 = r-1.0;
2924 t12 = t10*t10;
2925 f = (t3*t2-56.0/25.0*t3*t1+49.0/25.0*t3*r-98.0/115.0*t3+49.0/253.0*t2
2926 -28.0/1265.0*t1+7.0/6325.0*r-2.0/120175.0)*(s*t8+4.0/3.0*t8*t10+s*t12/2+t12*t10
2927 /21);
2928 /* p=14 8 0 2 */
2929 } else if( i==8 && j==0 && k==2 ) {
2930 t1 = r*r;
2931 t2 = t1*t1;
2932 t3 = t2*t2;
2933 t4 = t1*r;
2934 t9 = t*t;
2935 t10 = s+r-1.0;
2936 t12 = t10*t10;
2937 f = (t3-72.0/25.0*t2*t4+84.0/25.0*t2*t1-1176.0/575.0*t2*r+882.0/1265.0*
2938 t2-168.0/1265.0*t4+84.0/6325.0*t1-72.0/120175.0*r+1.0/120175.0)*(t9+6.0/5.0*t*
2939 t10+3.0/10.0*t12);
2940 /* p=14 8 1 1 */
2941 } else if( i==8 && j==1 && k==1 ) {
2942 t1 = r*r;
2943 t2 = t1*t1;
2944 t3 = t2*t2;
2945 t4 = t1*r;
2946 f = (t3-72.0/25.0*t2*t4+84.0/25.0*t2*t1-1176.0/575.0*t2*r+882.0/1265.0*
2947 t2-168.0/1265.0*t4+84.0/6325.0*t1-72.0/120175.0*r+1.0/120175.0)*(s+2.0/7.0*r
2948 -2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
2949 /* p=14 8 2 0 */
2950 } else if( i==8 && j==2 && k==0 ) {
2951 t1 = r*r;
2952 t2 = t1*t1;
2953 t3 = t2*t2;
2954 t4 = t1*r;
2955 t9 = s*s;
2956 t10 = r-1.0;
2957 t12 = t10*t10;
2958 f = (t3-72.0/25.0*t2*t4+84.0/25.0*t2*t1-1176.0/575.0*t2*r+882.0/1265.0*
2959 t2-168.0/1265.0*t4+84.0/6325.0*t1-72.0/120175.0*r+1.0/120175.0)*(t9+6.0/7.0*s*
2960 t10+t12/7);
2961 /* p=14 9 0 1 */
2962 } else if( i==9 && j==0 && k==1 ) {
2963 t1 = r*r;
2964 t2 = t1*t1;
2965 t3 = t2*t2;
2966 t5 = t1*r;
2967 f = (t3*r-18.0/5.0*t3+27.0/5.0*t2*t5-504.0/115.0*t2*t1+2646.0/1265.0*t2
2968 *r-756.0/1265.0*t2+126.0/1265.0*t5-216.0/24035.0*t1+9.0/24035.0*r-2.0/408595.0)
2969 *(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
2970 /* p=14 9 1 0 */
2971 } else if( i==9 && j==1 && k==0 ) {
2972 t1 = r*r;
2973 t2 = t1*t1;
2974 t3 = t2*t2;
2975 t5 = t1*r;
2976 f = (t3*r-18.0/5.0*t3+27.0/5.0*t2*t5-504.0/115.0*t2*t1+2646.0/1265.0*t2
2977 *r-756.0/1265.0*t2+126.0/1265.0*t5-216.0/24035.0*t1+9.0/24035.0*r-2.0/408595.0)
2978 *(s+2.0/5.0*r-2.0/5.0);
2979 /* p=14 10 0 0 */
2980 } else if( i==10 && j==0 && k==0 ) {
2981 t1 = r*r;
2982 t2 = t1*t1;
2983 t3 = t2*t2;
2984 t6 = t1*r;
2985 f = t3*t1-22.0/5.0*t3*r+33.0/4.0*t3-198.0/23.0*t2*t6+126.0/23.0*t2*t1
2986 -252.0/115.0*t2*r+63.0/115.0*t2-36.0/437.0*t6+3.0/437.0*t1-2.0/7429.0*r+1.0/
2987 297160.0;
2988 /* p=15 0 0 11 */
2989 } else if( i==0 && j==0 && k==11 ) {
2990 t1 = t*t;
2991 t2 = t1*t;
2992 t3 = t1*t1;
2993 t4 = t3*t3;
2994 t7 = s+r-1.0;
2995 t10 = t7*t7;
2996 t12 = t10*t7;
2997 t15 = t10*t10;
2998 t25 = t15*t15;
2999 f = t2*t4+132.0/23.0*t4*t1*t7+330.0/23.0*t4*t*t10+3300.0/161.0*t4*t12+
3000 2970.0/161.0*t3*t2*t15+4752.0/437.0*t3*t1*t15*t7+1848.0/437.0*t3*t*t15*t10+
3001 7920.0/7429.0*t3*t15*t12+2475.0/14858.0*t2*t25+110.0/7429.0*t1*t25*t7+33.0/
3002 52003.0*t*t25*t10+6.0/676039.0*t25*t12;
3003 /* p=15 0 1 10 */
3004 } else if( i==0 && j==1 && k==10 ) {
3005 t2 = t*t;
3006 t3 = t2*t2;
3007 t4 = t3*t3;
3008 t7 = s+r-1.0;
3009 t9 = t7*t7;
3010 t11 = t2*t;
3011 t13 = t9*t7;
3012 t16 = t9*t9;
3013 t25 = t16*t16;
3014 t30 = t2*t4+110.0/21.0*t4*t*t7+165.0/14.0*t4*t9+1980.0/133.0*t3*t11*t13+
3015 220.0/19.0*t3*t2*t16+1848.0/323.0*t3*t*t16*t7+1155.0/646.0*t3*t16*t9+110.0/
3016 323.0*t11*t16*t13+165.0/4522.0*t2*t25+55.0/29393.0*t*t25*t7+11.0/352716.0*t25*
3017 t9;
3018 f = (s+2.0/25.0*r-2.0/25.0)*t30;
3019 /* p=15 0 2 9 */
3020 } else if( i==0 && j==2 && k==9 ) {
3021 t1 = s*s;
3022 t2 = r-1.0;
3023 t4 = t2*t2;
3024 t6 = t*t;
3025 t7 = t6*t6;
3026 t8 = t7*t7;
3027 t10 = s+r-1.0;
3028 t12 = t6*t;
3029 t14 = t10*t10;
3030 t17 = t14*t10;
3031 t20 = t14*t14;
3032 t28 = t20*t20;
3033 f = (t1+6.0/25.0*s*t2+t4/100)*(t8*t+90.0/19.0*t8*t10+180.0/19.0*t7*t12*
3034 t14+3360.0/323.0*t7*t6*t17+2205.0/323.0*t7*t*t20+882.0/323.0*t7*t20*t10+210.0/
3035 323.0*t12*t20*t14+360.0/4199.0*t6*t20*t17+45.0/8398.0*t*t28+5.0/46189.0*t28*t10
3036 );
3037 /* p=15 0 3 8 */
3038 } else if( i==0 && j==3 && k==8 ) {
3039 t1 = s*s;
3040 t3 = r-1.0;
3041 t5 = t3*t3;
3042 t9 = t*t;
3043 t10 = t9*t9;
3044 t11 = t10*t10;
3045 t12 = t9*t;
3046 t14 = s+r-1.0;
3047 t17 = t14*t14;
3048 t20 = t17*t14;
3049 t22 = t17*t17;
3050 t30 = t22*t22;
3051 f = (t1*s+12.0/25.0*t3*t1+3.0/50.0*s*t5+t5*t3/575)*(t11+72.0/17.0*t10*
3052 t12*t14+126.0/17.0*t10*t9*t17+588.0/85.0*t10*t*t20+63.0/17.0*t10*t22+252.0/
3053 221.0*t12*t22*t14+42.0/221.0*t9*t22*t17+36.0/2431.0*t*t22*t20+9.0/24310.0*t30);
3054 /* p=15 0 4 7 */
3055 } else if( i==0 && j==4 && k==7 ) {
3056 t1 = s*s;
3057 t2 = t1*t1;
3058 t4 = r-1.0;
3059 t6 = t4*t4;
3060 t10 = t6*t6;
3061 t12 = t*t;
3062 t13 = t12*t;
3063 t14 = t12*t12;
3064 t17 = s+r-1.0;
3065 t20 = t17*t17;
3066 t22 = t20*t17;
3067 t24 = t20*t20;
3068 f = (t2+4.0/5.0*t1*s*t4+t1*t6/5+2.0/115.0*s*t6*t4+t10/2530)*(t13*t14+
3069 56.0/15.0*t14*t12*t17+28.0/5.0*t14*t*t20+56.0/13.0*t22*t14+70.0/39.0*t13*t24+
3070 56.0/143.0*t12*t24*t17+28.0/715.0*t*t24*t20+8.0/6435.0*t24*t22);
3071 /* p=15 0 5 6 */
3072 } else if( i==0 && j==5 && k==6 ) {
3073 t1 = s*s;
3074 t2 = t1*t1;
3075 t4 = r-1.0;
3076 t7 = t4*t4;
3077 t11 = t7*t7;
3078 t15 = t*t;
3079 t16 = t15*t15;
3080 t19 = s+r-1.0;
3081 t21 = t19*t19;
3082 t26 = t21*t21;
3083 f = (s*t2+6.0/5.0*t2*t4+t7*t1*s/2+2.0/23.0*t1*t7*t4+3.0/506.0*s*t11+t11
3084 *t4/8855)*(t16*t15+42.0/13.0*t16*t*t19+105.0/26.0*t16*t21+350.0/143.0*t15*t*t21
3085 *t19+105.0/143.0*t15*t26+14.0/143.0*t*t26*t19+7.0/1716.0*t26*t21);
3086 /* p=15 0 6 5 */
3087 } else if( i==0 && j==6 && k==5 ) {
3088 t1 = s*s;
3089 t2 = t1*t1;
3090 t5 = r-1.0;
3091 t7 = t5*t5;
3092 t12 = t7*t7;
3093 t18 = t*t;
3094 t19 = t18*t18;
3095 t21 = s+r-1.0;
3096 t24 = t21*t21;
3097 t28 = t24*t24;
3098 f = (t2*t1+42.0/25.0*t2*s*t5+21.0/20.0*t2*t7+7.0/23.0*t7*t5*t1*s+21.0/
3099 506.0*t1*t12+3.0/1265.0*s*t12*t5+t7*t12/25300)*(t19*t+30.0/11.0*t21*t19+30.0/
3100 11.0*t18*t*t24+40.0/33.0*t18*t24*t21+5.0/22.0*t*t28+t28*t21/77);
3101 /* p=15 0 7 4 */
3102 } else if( i==0 && j==7 && k==4 ) {
3103 t1 = s*s;
3104 t2 = t1*s;
3105 t3 = t1*t1;
3106 t6 = r-1.0;
3107 t9 = t6*t6;
3108 t11 = t9*t6;
3109 t13 = t9*t9;
3110 t21 = t*t;
3111 t22 = t21*t21;
3112 t24 = s+r-1.0;
3113 t26 = t24*t24;
3114 t30 = t26*t26;
3115 f = (t3*t2+56.0/25.0*t3*t1*t6+49.0/25.0*t3*s*t9+98.0/115.0*t3*t11+49.0/
3116 253.0*t2*t13+28.0/1265.0*t1*t13*t6+7.0/6325.0*s*t13*t9+2.0/120175.0*t13*t11)*(
3117 t22+20.0/9.0*t21*t*t24+5.0/3.0*t21*t26+10.0/21.0*t*t26*t24+5.0/126.0*t30);
3118 /* p=15 0 8 3 */
3119 } else if( i==0 && j==8 && k==3 ) {
3120 t1 = s*s;
3121 t2 = t1*t1;
3122 t3 = t2*t2;
3123 t4 = t1*s;
3124 t6 = r-1.0;
3125 t9 = t6*t6;
3126 t12 = t9*t6;
3127 t14 = t9*t9;
3128 t22 = t14*t14;
3129 t24 = t*t;
3130 t26 = s+r-1.0;
3131 t28 = t26*t26;
3132 f = (t3+72.0/25.0*t2*t4*t6+84.0/25.0*t2*t1*t9+1176.0/575.0*t2*s*t12+
3133 882.0/1265.0*t2*t14+168.0/1265.0*t4*t14*t6+84.0/6325.0*t1*t14*t9+72.0/120175.0*
3134 s*t14*t12+t22/120175)*(t24*t+12.0/7.0*t26*t24+6.0/7.0*t*t28+4.0/35.0*t28*t26);
3135 /* p=15 0 9 2 */
3136 } else if( i==0 && j==9 && k==2 ) {
3137 t1 = s*s;
3138 t2 = t1*t1;
3139 t3 = t2*t2;
3140 t5 = r-1.0;
3141 t7 = t1*s;
3142 t9 = t5*t5;
3143 t12 = t9*t5;
3144 t15 = t9*t9;
3145 t23 = t15*t15;
3146 t27 = t*t;
3147 t28 = s+r-1.0;
3148 t30 = t28*t28;
3149 f = (t3*s+18.0/5.0*t3*t5+27.0/5.0*t2*t7*t9+504.0/115.0*t2*t1*t12+2646.0
3150 /1265.0*t2*s*t15+756.0/1265.0*t2*t15*t5+126.0/1265.0*t7*t15*t9+216.0/24035.0*t1
3151 *t15*t12+9.0/24035.0*s*t23+2.0/408595.0*t23*t5)*(t27+6.0/5.0*t*t28+3.0/10.0*t30
3152 );
3153 /* p=15 0 10 1 */
3154 } else if( i==0 && j==10 && k==1 ) {
3155 t1 = s*s;
3156 t2 = t1*t1;
3157 t3 = t2*t2;
3158 t6 = r-1.0;
3159 t8 = t6*t6;
3160 t10 = t1*s;
3161 t12 = t8*t6;
3162 t15 = t8*t8;
3163 t24 = t15*t15;
3164 t29 = t3*t1+22.0/5.0*t3*s*t6+33.0/4.0*t3*t8+198.0/23.0*t2*t10*t12+126.0/
3165 23.0*t2*t1*t15+252.0/115.0*t2*s*t15*t6+63.0/115.0*t2*t15*t8+36.0/437.0*t10*t15*
3166 t12+3.0/437.0*t1*t24+2.0/7429.0*s*t24*t6+t24*t8/297160;
3167 f = t29*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
3168 /* p=15 0 11 0 */
3169 } else if( i==0 && j==11 && k==0 ) {
3170 t1 = s*s;
3171 t2 = t1*s;
3172 t3 = t1*t1;
3173 t4 = t3*t3;
3174 t7 = r-1.0;
3175 t10 = t7*t7;
3176 t12 = t10*t7;
3177 t15 = t10*t10;
3178 t25 = t15*t15;
3179 f = t2*t4+132.0/25.0*t4*t1*t7+121.0/10.0*t4*s*t10+363.0/23.0*t4*t12+
3180 297.0/23.0*t3*t2*t15+792.0/115.0*t3*t1*t15*t7+1386.0/575.0*t3*s*t15*t10+1188.0/
3181 2185.0*t3*t15*t12+33.0/437.0*t2*t25+44.0/7429.0*t1*t25*t7+33.0/148580.0*s*t25*
3182 t10+t25*t12/371450;
3183 /* p=15 1 0 10 */
3184 } else if( i==1 && j==0 && k==10 ) {
3185 t2 = t*t;
3186 t3 = t2*t2;
3187 t4 = t3*t3;
3188 t7 = s+r-1.0;
3189 t9 = t7*t7;
3190 t11 = t2*t;
3191 t13 = t9*t7;
3192 t16 = t9*t9;
3193 t25 = t16*t16;
3194 t30 = t2*t4+110.0/21.0*t4*t*t7+165.0/14.0*t4*t9+1980.0/133.0*t3*t11*t13+
3195 220.0/19.0*t3*t2*t16+1848.0/323.0*t3*t*t16*t7+1155.0/646.0*t3*t16*t9+110.0/
3196 323.0*t11*t16*t13+165.0/4522.0*t2*t25+55.0/29393.0*t*t25*t7+11.0/352716.0*t25*
3197 t9;
3198 f = (r-2.0/27.0)*t30;
3199 /* p=15 1 1 9 */
3200 } else if( i==1 && j==1 && k==9 ) {
3201 t4 = t*t;
3202 t5 = t4*t4;
3203 t6 = t5*t5;
3204 t8 = s+r-1.0;
3205 t10 = t4*t;
3206 t12 = t8*t8;
3207 t15 = t12*t8;
3208 t18 = t12*t12;
3209 t26 = t18*t18;
3210 f = (s+2.0/23.0*r-2.0/23.0)*(r-2.0/27.0)*(t6*t+90.0/19.0*t6*t8+180.0/
3211 19.0*t5*t10*t12+3360.0/323.0*t5*t4*t15+2205.0/323.0*t5*t*t18+882.0/323.0*t5*t18
3212 *t8+210.0/323.0*t10*t18*t12+360.0/4199.0*t4*t18*t15+45.0/8398.0*t*t26+5.0/
3213 46189.0*t26*t8);
3214 /* p=15 1 2 8 */
3215 } else if( i==1 && j==2 && k==8 ) {
3216 t2 = s*s;
3217 t3 = r-1.0;
3218 t5 = t3*t3;
3219 t8 = t*t;
3220 t9 = t8*t8;
3221 t10 = t9*t9;
3222 t11 = t8*t;
3223 t13 = s+r-1.0;
3224 t16 = t13*t13;
3225 t19 = t16*t13;
3226 t21 = t16*t16;
3227 t29 = t21*t21;
3228 f = (r-2.0/27.0)*(t2+6.0/23.0*s*t3+3.0/253.0*t5)*(t10+72.0/17.0*t9*t11*
3229 t13+126.0/17.0*t9*t8*t16+588.0/85.0*t9*t*t19+63.0/17.0*t9*t21+252.0/221.0*t11*
3230 t21*t13+42.0/221.0*t8*t21*t16+36.0/2431.0*t*t21*t19+9.0/24310.0*t29);
3231 /* p=15 1 3 7 */
3232 } else if( i==1 && j==3 && k==7 ) {
3233 t2 = s*s;
3234 t4 = r-1.0;
3235 t6 = t4*t4;
3236 t11 = t*t;
3237 t12 = t11*t;
3238 t13 = t11*t11;
3239 t16 = s+r-1.0;
3240 t19 = t16*t16;
3241 t21 = t19*t16;
3242 t23 = t19*t19;
3243 f = (r-2.0/27.0)*(t2*s+12.0/23.0*t2*t4+18.0/253.0*s*t6+4.0/1771.0*t6*t4
3244 )*(t12*t13+56.0/15.0*t11*t16*t13+28.0/5.0*t13*t*t19+56.0/13.0*t21*t13+70.0/39.0
3245 *t12*t23+56.0/143.0*t11*t23*t16+28.0/715.0*t*t23*t19+8.0/6435.0*t23*t21);
3246 /* p=15 1 4 6 */
3247 } else if( i==1 && j==4 && k==6 ) {
3248 t2 = s*s;
3249 t3 = t2*t2;
3250 t5 = r-1.0;
3251 t7 = t5*t5;
3252 t11 = t7*t7;
3253 t14 = t*t;
3254 t15 = t14*t14;
3255 t18 = s+r-1.0;
3256 t20 = t18*t18;
3257 t25 = t20*t20;
3258 f = (r-2.0/27.0)*(t3+20.0/23.0*t5*t2*s+60.0/253.0*t2*t7+40.0/1771.0*s*
3259 t7*t5+t11/1771)*(t14*t15+42.0/13.0*t15*t*t18+105.0/26.0*t15*t20+350.0/143.0*t14
3260 *t*t20*t18+105.0/143.0*t14*t25+14.0/143.0*t*t25*t18+7.0/1716.0*t25*t20);
3261 /* p=15 1 5 5 */
3262 } else if( i==1 && j==5 && k==5 ) {
3263 t2 = s*s;
3264 t3 = t2*t2;
3265 t5 = r-1.0;
3266 t8 = t5*t5;
3267 t12 = t8*t8;
3268 t17 = t*t;
3269 t18 = t17*t17;
3270 t20 = s+r-1.0;
3271 t23 = t20*t20;
3272 t27 = t23*t23;
3273 f = (r-2.0/27.0)*(s*t3+30.0/23.0*t3*t5+150.0/253.0*t2*s*t8+200.0/1771.0
3274 *t2*t8*t5+15.0/1771.0*s*t12+6.0/33649.0*t12*t5)*(t18*t+30.0/11.0*t20*t18+30.0/
3275 11.0*t17*t*t23+40.0/33.0*t17*t23*t20+5.0/22.0*t*t27+t27*t20/77);
3276 /* p=15 1 6 4 */
3277 } else if( i==1 && j==6 && k==4 ) {
3278 t2 = s*s;
3279 t3 = t2*t2;
3280 t6 = r-1.0;
3281 t8 = t6*t6;
3282 t13 = t8*t8;
3283 t20 = t*t;
3284 t21 = t20*t20;
3285 t23 = s+r-1.0;
3286 t25 = t23*t23;
3287 t29 = t25*t25;
3288 f = (r-2.0/27.0)*(t3*t2+42.0/23.0*t3*s*t6+315.0/253.0*t3*t8+100.0/253.0
3289 *t2*s*t8*t6+15.0/253.0*t2*t13+18.0/4807.0*s*t13*t6+t13*t8/14421)*(t21+20.0/9.0*
3290 t20*t*t23+5.0/3.0*t20*t25+10.0/21.0*t*t25*t23+5.0/126.0*t29);
3291 /* p=15 1 7 3 */
3292 } else if( i==1 && j==7 && k==3 ) {
3293 t2 = s*s;
3294 t3 = t2*s;
3295 t4 = t2*t2;
3296 t7 = r-1.0;
3297 t10 = t7*t7;
3298 t12 = t10*t7;
3299 t14 = t10*t10;
3300 t23 = t*t;
3301 t25 = s+r-1.0;
3302 t27 = t25*t25;
3303 f = (r-2.0/27.0)*(t4*t3+56.0/23.0*t4*t2*t7+588.0/253.0*t4*s*t10+280.0/
3304 253.0*t4*t12+70.0/253.0*t3*t14+168.0/4807.0*t2*t14*t7+28.0/14421.0*s*t14*t10+
3305 8.0/245157.0*t14*t12)*(t23*t+12.0/7.0*t25*t23+6.0/7.0*t*t27+4.0/35.0*t27*t25);
3306 /* p=15 1 8 2 */
3307 } else if( i==1 && j==8 && k==2 ) {
3308 t2 = s*s;
3309 t3 = t2*t2;
3310 t4 = t3*t3;
3311 t5 = t2*s;
3312 t7 = r-1.0;
3313 t10 = t7*t7;
3314 t13 = t10*t7;
3315 t15 = t10*t10;
3316 t23 = t15*t15;
3317 t26 = t*t;
3318 t27 = s+r-1.0;
3319 t29 = t27*t27;
3320 f = (r-2.0/27.0)*(t4+72.0/23.0*t3*t5*t7+1008.0/253.0*t3*t2*t10+672.0/
3321 253.0*t3*s*t13+252.0/253.0*t3*t15+1008.0/4807.0*t5*t15*t7+112.0/4807.0*t2*t15*
3322 t10+96.0/81719.0*s*t15*t13+3.0/163438.0*t23)*(t26+6.0/5.0*t*t27+3.0/10.0*t29);
3323 /* p=15 1 9 1 */
3324 } else if( i==1 && j==9 && k==1 ) {
3325 t2 = s*s;
3326 t3 = t2*t2;
3327 t4 = t3*t3;
3328 t6 = r-1.0;
3329 t8 = t2*s;
3330 t10 = t6*t6;
3331 t13 = t10*t6;
3332 t16 = t10*t10;
3333 t24 = t16*t16;
3334 f = (r-2.0/27.0)*(t4*s+90.0/23.0*t4*t6+1620.0/253.0*t3*t8*t10+1440.0/
3335 253.0*t3*t2*t13+756.0/253.0*t3*s*t16+4536.0/4807.0*t3*t16*t6+840.0/4807.0*t8*
3336 t16*t10+1440.0/81719.0*t2*t16*t13+135.0/163438.0*s*t24+t24*t6/81719)*(t+2.0/3.0
3337 *s+2.0/3.0*r-2.0/3.0);
3338 /* p=15 1 10 0 */
3339 } else if( i==1 && j==10 && k==0 ) {
3340 t2 = s*s;
3341 t3 = t2*t2;
3342 t4 = t3*t3;
3343 t7 = r-1.0;
3344 t9 = t7*t7;
3345 t11 = t2*s;
3346 t13 = t9*t7;
3347 t16 = t9*t9;
3348 t25 = t16*t16;
3349 t30 = t4*t2+110.0/23.0*t4*s*t7+225.0/23.0*t4*t9+1800.0/161.0*t3*t11*t13+
3350 180.0/23.0*t3*t2*t16+1512.0/437.0*t3*s*t16*t7+420.0/437.0*t3*t16*t9+1200.0/
3351 7429.0*t11*t16*t13+225.0/14858.0*t2*t25+5.0/7429.0*s*t25*t7+t25*t9/104006;
3352 f = (r-2.0/27.0)*t30;
3353 /* p=15 2 0 9 */
3354 } else if( i==2 && j==0 && k==9 ) {
3355 t1 = r*r;
3356 t3 = t*t;
3357 t4 = t3*t3;
3358 t5 = t4*t4;
3359 t7 = s+r-1.0;
3360 t9 = t3*t;
3361 t11 = t7*t7;
3362 t14 = t11*t7;
3363 t17 = t11*t11;
3364 t25 = t17*t17;
3365 f = (t1-2.0/9.0*r+1.0/117.0)*(t5*t+90.0/19.0*t5*t7+180.0/19.0*t4*t9*t11
3366 +3360.0/323.0*t4*t3*t14+2205.0/323.0*t4*t*t17+882.0/323.0*t4*t17*t7+210.0/323.0
3367 *t9*t17*t11+360.0/4199.0*t3*t17*t14+45.0/8398.0*t*t25+5.0/46189.0*t25*t7);
3368 /* p=15 2 1 8 */
3369 } else if( i==2 && j==1 && k==8 ) {
3370 t1 = r*r;
3371 t5 = t*t;
3372 t6 = t5*t5;
3373 t7 = t6*t6;
3374 t8 = t5*t;
3375 t10 = s+r-1.0;
3376 t13 = t10*t10;
3377 t16 = t13*t10;
3378 t18 = t13*t13;
3379 t26 = t18*t18;
3380 f = (t1-2.0/9.0*r+1.0/117.0)*(s+2.0/21.0*r-2.0/21.0)*(t7+72.0/17.0*t6*
3381 t8*t10+126.0/17.0*t6*t5*t13+588.0/85.0*t6*t*t16+63.0/17.0*t6*t18+252.0/221.0*t8
3382 *t18*t10+42.0/221.0*t5*t18*t13+36.0/2431.0*t*t18*t16+9.0/24310.0*t26);
3383 /* p=15 2 2 7 */
3384 } else if( i==2 && j==2 && k==7 ) {
3385 t1 = r*r;
3386 t3 = s*s;
3387 t4 = r-1.0;
3388 t6 = t4*t4;
3389 t9 = t*t;
3390 t10 = t9*t;
3391 t11 = t9*t9;
3392 t14 = s+r-1.0;
3393 t17 = t14*t14;
3394 t19 = t17*t14;
3395 t21 = t17*t17;
3396 f = (t1-2.0/9.0*r+1.0/117.0)*(t3+2.0/7.0*s*t4+t6/70)*(t11*t10+56.0/15.0
3397 *t11*t9*t14+28.0/5.0*t11*t*t17+56.0/13.0*t11*t19+70.0/39.0*t10*t21+56.0/143.0*
3398 t9*t21*t14+28.0/715.0*t*t21*t17+8.0/6435.0*t21*t19);
3399 /* p=15 2 3 6 */
3400 } else if( i==2 && j==3 && k==6 ) {
3401 t1 = r*r;
3402 t3 = s*s;
3403 t5 = r-1.0;
3404 t7 = t5*t5;
3405 t12 = t*t;
3406 t13 = t12*t12;
3407 t16 = s+r-1.0;
3408 t18 = t16*t16;
3409 t23 = t18*t18;
3410 f = (t1-2.0/9.0*r+1.0/117.0)*(t3*s+4.0/7.0*t3*t5+3.0/35.0*s*t7+2.0/
3411 665.0*t5*t7)*(t12*t13+42.0/13.0*t13*t*t16+105.0/26.0*t18*t13+350.0/143.0*t12*t*
3412 t18*t16+105.0/143.0*t12*t23+14.0/143.0*t*t23*t16+7.0/1716.0*t23*t18);
3413 /* p=15 2 4 5 */
3414 } else if( i==2 && j==4 && k==5 ) {
3415 t1 = r*r;
3416 t3 = s*s;
3417 t4 = t3*t3;
3418 t6 = r-1.0;
3419 t8 = t6*t6;
3420 t12 = t8*t8;
3421 t15 = t*t;
3422 t16 = t15*t15;
3423 t18 = s+r-1.0;
3424 t21 = t18*t18;
3425 t25 = t21*t21;
3426 f = (t1-2.0/9.0*r+1.0/117.0)*(t4+20.0/21.0*t6*t3*s+2.0/7.0*t3*t8+4.0/
3427 133.0*s*t6*t8+t12/1197)*(t16*t+30.0/11.0*t18*t16+30.0/11.0*t15*t*t21+40.0/33.0*
3428 t15*t21*t18+5.0/22.0*t*t25+t25*t18/77);
3429 /* p=15 2 5 4 */
3430 } else if( i==2 && j==5 && k==4 ) {
3431 t1 = r*r;
3432 t3 = s*s;
3433 t4 = t3*t3;
3434 t6 = r-1.0;
3435 t9 = t6*t6;
3436 t13 = t9*t9;
3437 t18 = t*t;
3438 t19 = t18*t18;
3439 t21 = s+r-1.0;
3440 t23 = t21*t21;
3441 t27 = t23*t23;
3442 f = (t1-2.0/9.0*r+1.0/117.0)*(s*t4+10.0/7.0*t4*t6+5.0/7.0*t3*s*t9+20.0/
3443 133.0*t3*t9*t6+5.0/399.0*s*t13+2.0/6783.0*t13*t6)*(t19+20.0/9.0*t18*t*t21+5.0/
3444 3.0*t23*t18+10.0/21.0*t*t23*t21+5.0/126.0*t27);
3445 /* p=15 2 6 3 */
3446 } else if( i==2 && j==6 && k==3 ) {
3447 t1 = r*r;
3448 t3 = s*s;
3449 t4 = t3*t3;
3450 t7 = r-1.0;
3451 t9 = t7*t7;
3452 t14 = t9*t9;
3453 t21 = t*t;
3454 t23 = s+r-1.0;
3455 t25 = t23*t23;
3456 f = (t1-2.0/9.0*r+1.0/117.0)*(t4*t3+2.0*t4*s*t7+3.0/2.0*t4*t9+10.0/19.0
3457 *t3*s*t9*t7+5.0/57.0*t3*t14+2.0/323.0*s*t14*t7+t14*t9/7752)*(t21*t+12.0/7.0*t23
3458 *t21+6.0/7.0*t*t25+4.0/35.0*t25*t23);
3459 /* p=15 2 7 2 */
3460 } else if( i==2 && j==7 && k==2 ) {
3461 t1 = r*r;
3462 t3 = s*s;
3463 t4 = t3*s;
3464 t5 = t3*t3;
3465 t8 = r-1.0;
3466 t11 = t8*t8;
3467 t13 = t11*t8;
3468 t15 = t11*t11;
3469 t24 = t*t;
3470 t25 = s+r-1.0;
3471 t27 = t25*t25;
3472 f = (t1-2.0/9.0*r+1.0/117.0)*(t5*t4+8.0/3.0*t5*t3*t8+14.0/5.0*t5*s*t11+
3473 28.0/19.0*t5*t13+70.0/171.0*t4*t15+56.0/969.0*t3*t15*t8+7.0/1938.0*s*t15*t11+
3474 t15*t13/14535)*(t24+6.0/5.0*t*t25+3.0/10.0*t27);
3475 /* p=15 2 8 1 */
3476 } else if( i==2 && j==8 && k==1 ) {
3477 t1 = r*r;
3478 t3 = s*s;
3479 t4 = t3*t3;
3480 t5 = t4*t4;
3481 t6 = t3*s;
3482 t8 = r-1.0;
3483 t11 = t8*t8;
3484 t14 = t11*t8;
3485 t16 = t11*t11;
3486 t24 = t16*t16;
3487 f = (t1-2.0/9.0*r+1.0/117.0)*(t5+24.0/7.0*t4*t6*t8+24.0/5.0*t4*t3*t11+
3488 336.0/95.0*t4*s*t14+28.0/19.0*t4*t16+112.0/323.0*t6*t16*t8+14.0/323.0*t3*t16*
3489 t11+4.0/1615.0*s*t16*t14+t24/22610)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
3490 /* p=15 2 9 0 */
3491 } else if( i==2 && j==9 && k==0 ) {
3492 t1 = r*r;
3493 t3 = s*s;
3494 t4 = t3*t3;
3495 t5 = t4*t4;
3496 t7 = r-1.0;
3497 t9 = t3*s;
3498 t11 = t7*t7;
3499 t14 = t11*t7;
3500 t17 = t11*t11;
3501 t25 = t17*t17;
3502 f = (t1-2.0/9.0*r+1.0/117.0)*(t5*s+30.0/7.0*t5*t7+54.0/7.0*t4*t9*t11+
3503 144.0/19.0*t4*t3*t14+84.0/19.0*t4*s*t17+504.0/323.0*t4*t17*t7+105.0/323.0*t9*
3504 t17*t11+12.0/323.0*t3*t17*t14+9.0/4522.0*s*t25+t25*t7/29393);
3505 /* p=15 3 0 8 */
3506 } else if( i==3 && j==0 && k==8 ) {
3507 t1 = r*r;
3508 t4 = t*t;
3509 t5 = t4*t4;
3510 t6 = t5*t5;
3511 t7 = t4*t;
3512 t9 = s+r-1.0;
3513 t12 = t9*t9;
3514 t15 = t12*t9;
3515 t17 = t12*t12;
3516 t25 = t17*t17;
3517 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t6+72.0/17.0*t5*t7*t9+
3518 126.0/17.0*t5*t4*t12+588.0/85.0*t5*t*t15+63.0/17.0*t5*t17+252.0/221.0*t7*t17*t9
3519 +42.0/221.0*t4*t17*t12+36.0/2431.0*t*t17*t15+9.0/24310.0*t25);
3520 /* p=15 3 1 7 */
3521 } else if( i==3 && j==1 && k==7 ) {
3522 t1 = r*r;
3523 t6 = t*t;
3524 t7 = t6*t;
3525 t8 = t6*t6;
3526 t11 = s+r-1.0;
3527 t14 = t11*t11;
3528 t16 = t14*t11;
3529 t18 = t14*t14;
3530 f = (s+2.0/19.0*r-2.0/19.0)*(t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t7
3531 *t8+56.0/15.0*t8*t6*t11+28.0/5.0*t8*t*t14+56.0/13.0*t16*t8+70.0/39.0*t7*t18+
3532 56.0/143.0*t6*t18*t11+28.0/715.0*t*t18*t14+8.0/6435.0*t18*t16);
3533 /* p=15 3 2 6 */
3534 } else if( i==3 && j==2 && k==6 ) {
3535 t1 = r*r;
3536 t4 = s*s;
3537 t5 = r-1.0;
3538 t7 = t5*t5;
3539 t10 = t*t;
3540 t11 = t10*t10;
3541 t14 = s+r-1.0;
3542 t16 = t14*t14;
3543 t21 = t16*t16;
3544 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t4+6.0/19.0*s*t5+t7/57)*(
3545 t10*t11+42.0/13.0*t11*t*t14+105.0/26.0*t16*t11+350.0/143.0*t10*t*t16*t14+105.0/
3546 143.0*t10*t21+14.0/143.0*t*t21*t14+7.0/1716.0*t21*t16);
3547 /* p=15 3 3 5 */
3548 } else if( i==3 && j==3 && k==5 ) {
3549 t1 = r*r;
3550 t4 = s*s;
3551 t6 = r-1.0;
3552 t8 = t6*t6;
3553 t13 = t*t;
3554 t14 = t13*t13;
3555 t16 = s+r-1.0;
3556 t19 = t16*t16;
3557 t23 = t19*t19;
3558 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t4*s+12.0/19.0*t4*t6+2.0/
3559 19.0*s*t8+4.0/969.0*t8*t6)*(t14*t+30.0/11.0*t16*t14+30.0/11.0*t13*t*t19+40.0/
3560 33.0*t13*t19*t16+5.0/22.0*t*t23+t23*t16/77);
3561 /* p=15 3 4 4 */
3562 } else if( i==3 && j==4 && k==4 ) {
3563 t1 = r*r;
3564 t4 = s*s;
3565 t5 = t4*t4;
3566 t7 = r-1.0;
3567 t9 = t7*t7;
3568 t13 = t9*t9;
3569 t16 = t*t;
3570 t17 = t16*t16;
3571 t19 = s+r-1.0;
3572 t21 = t19*t19;
3573 t25 = t21*t21;
3574 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t5+20.0/19.0*t4*s*t7+20.0/
3575 57.0*t4*t9+40.0/969.0*s*t9*t7+5.0/3876.0*t13)*(t17+20.0/9.0*t16*t*t19+5.0/3.0*
3576 t21*t16+10.0/21.0*t*t21*t19+5.0/126.0*t25);
3577 /* p=15 3 5 3 */
3578 } else if( i==3 && j==5 && k==3 ) {
3579 t1 = r*r;
3580 t4 = s*s;
3581 t5 = t4*t4;
3582 t7 = r-1.0;
3583 t10 = t7*t7;
3584 t14 = t10*t10;
3585 t19 = t*t;
3586 t21 = s+r-1.0;
3587 t23 = t21*t21;
3588 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(s*t5+30.0/19.0*t5*t7+50.0/
3589 57.0*t4*s*t10+200.0/969.0*t4*t10*t7+25.0/1292.0*s*t14+t14*t7/1938)*(t19*t+12.0/
3590 7.0*t21*t19+6.0/7.0*t*t23+4.0/35.0*t23*t21);
3591 /* p=15 3 6 2 */
3592 } else if( i==3 && j==6 && k==2 ) {
3593 t1 = r*r;
3594 t4 = s*s;
3595 t5 = t4*t4;
3596 t8 = r-1.0;
3597 t10 = t8*t8;
3598 t15 = t10*t10;
3599 t22 = t*t;
3600 t23 = s+r-1.0;
3601 t25 = t23*t23;
3602 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t5*t4+42.0/19.0*t5*s*t8+
3603 35.0/19.0*t5*t10+700.0/969.0*t4*s*t10*t8+175.0/1292.0*t4*t15+7.0/646.0*s*t15*t8
3604 +t15*t10/3876)*(t22+6.0/5.0*t*t23+3.0/10.0*t25);
3605 /* p=15 3 7 1 */
3606 } else if( i==3 && j==7 && k==1 ) {
3607 t1 = r*r;
3608 t4 = s*s;
3609 t5 = t4*s;
3610 t6 = t4*t4;
3611 t9 = r-1.0;
3612 t12 = t9*t9;
3613 t14 = t12*t9;
3614 t16 = t12*t12;
3615 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t6*t5+56.0/19.0*t4*t6*t9+
3616 196.0/57.0*t6*s*t12+1960.0/969.0*t6*t14+1225.0/1938.0*t5*t16+98.0/969.0*t4*t16*
3617 t9+7.0/969.0*s*t16*t12+2.0/12597.0*t16*t14)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
3618 /* p=15 3 8 0 */
3619 } else if( i==3 && j==8 && k==0 ) {
3620 t1 = r*r;
3621 t4 = s*s;
3622 t5 = t4*t4;
3623 t6 = t5*t5;
3624 t7 = t4*s;
3625 t9 = r-1.0;
3626 t12 = t9*t9;
3627 t15 = t12*t9;
3628 t17 = t12*t12;
3629 t25 = t17*t17;
3630 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t6+72.0/19.0*t5*t7*t9+
3631 112.0/19.0*t5*t4*t12+1568.0/323.0*t5*s*t15+735.0/323.0*t5*t17+196.0/323.0*t7*
3632 t17*t9+28.0/323.0*t4*t17*t12+24.0/4199.0*s*t17*t15+t25/8398);
3633 /* p=15 4 0 7 */
3634 } else if( i==4 && j==0 && k==7 ) {
3635 t1 = r*r;
3636 t2 = t1*t1;
3637 t5 = t*t;
3638 t6 = t5*t;
3639 t7 = t5*t5;
3640 t10 = s+r-1.0;
3641 t13 = t10*t10;
3642 t15 = t13*t10;
3643 t17 = t13*t13;
3644 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t7*t6+
3645 56.0/15.0*t5*t7*t10+28.0/5.0*t7*t*t13+56.0/13.0*t7*t15+70.0/39.0*t6*t17+56.0/
3646 143.0*t5*t17*t10+28.0/715.0*t*t17*t13+8.0/6435.0*t17*t15);
3647 /* p=15 4 1 6 */
3648 } else if( i==4 && j==1 && k==6 ) {
3649 t1 = r*r;
3650 t2 = t1*t1;
3651 t7 = t*t;
3652 t8 = t7*t7;
3653 t11 = s+r-1.0;
3654 t13 = t11*t11;
3655 t18 = t13*t13;
3656 f = (s+2.0/17.0*r-2.0/17.0)*(t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*
3657 r+1.0/3510.0)*(t8*t7+42.0/13.0*t8*t*t11+105.0/26.0*t8*t13+350.0/143.0*t7*t*t13*
3658 t11+105.0/143.0*t7*t18+14.0/143.0*t*t18*t11+7.0/1716.0*t18*t13);
3659 /* p=15 4 2 5 */
3660 } else if( i==4 && j==2 && k==5 ) {
3661 t1 = r*r;
3662 t2 = t1*t1;
3663 t5 = s*s;
3664 t6 = r-1.0;
3665 t8 = t6*t6;
3666 t11 = t*t;
3667 t12 = t11*t11;
3668 t14 = s+r-1.0;
3669 t17 = t14*t14;
3670 t21 = t17*t17;
3671 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t5+6.0/
3672 17.0*s*t6+3.0/136.0*t8)*(t12*t+30.0/11.0*t12*t14+30.0/11.0*t11*t*t17+40.0/33.0*
3673 t11*t17*t14+5.0/22.0*t*t21+t21*t14/77);
3674 /* p=15 4 3 4 */
3675 } else if( i==4 && j==3 && k==4 ) {
3676 t1 = r*r;
3677 t2 = t1*t1;
3678 t5 = s*s;
3679 t7 = r-1.0;
3680 t9 = t7*t7;
3681 t14 = t*t;
3682 t15 = t14*t14;
3683 t17 = s+r-1.0;
3684 t19 = t17*t17;
3685 t23 = t19*t19;
3686 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t5*s+12.0
3687 /17.0*t5*t7+9.0/68.0*s*t9+t9*t7/170)*(t15+20.0/9.0*t14*t*t17+5.0/3.0*t14*t19+
3688 10.0/21.0*t*t19*t17+5.0/126.0*t23);
3689 /* p=15 4 4 3 */
3690 } else if( i==4 && j==4 && k==3 ) {
3691 t1 = r*r;
3692 t2 = t1*t1;
3693 t5 = s*s;
3694 t6 = t5*t5;
3695 t8 = r-1.0;
3696 t10 = t8*t8;
3697 t14 = t10*t10;
3698 t17 = t*t;
3699 t19 = s+r-1.0;
3700 t21 = t19*t19;
3701 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t6+20.0/
3702 17.0*t8*t5*s+15.0/34.0*t5*t10+s*t10*t8/17+t14/476)*(t17*t+12.0/7.0*t19*t17+6.0/
3703 7.0*t*t21+4.0/35.0*t21*t19);
3704 /* p=15 4 5 2 */
3705 } else if( i==4 && j==5 && k==2 ) {
3706 t1 = r*r;
3707 t2 = t1*t1;
3708 t5 = s*s;
3709 t6 = t5*t5;
3710 t8 = r-1.0;
3711 t11 = t8*t8;
3712 t15 = t11*t11;
3713 t20 = t*t;
3714 t21 = s+r-1.0;
3715 t23 = t21*t21;
3716 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(s*t6+30.0
3717 /17.0*t6*t8+75.0/68.0*t5*s*t11+5.0/17.0*t5*t11*t8+15.0/476.0*s*t15+3.0/3094.0*
3718 t15*t8)*(t20+6.0/5.0*t*t21+3.0/10.0*t23);
3719 /* p=15 4 6 1 */
3720 } else if( i==4 && j==6 && k==1 ) {
3721 t1 = r*r;
3722 t2 = t1*t1;
3723 t5 = s*s;
3724 t6 = t5*t5;
3725 t9 = r-1.0;
3726 t11 = t9*t9;
3727 t16 = t11*t11;
3728 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t6*t5+
3729 42.0/17.0*t6*s*t9+315.0/136.0*t6*t11+35.0/34.0*t5*s*t11*t9+15.0/68.0*t5*t16+9.0
3730 /442.0*s*t16*t9+t16*t11/1768)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
3731 /* p=15 4 7 0 */
3732 } else if( i==4 && j==7 && k==0 ) {
3733 t1 = r*r;
3734 t2 = t1*t1;
3735 t5 = s*s;
3736 t6 = t5*s;
3737 t7 = t5*t5;
3738 t10 = r-1.0;
3739 t13 = t10*t10;
3740 t15 = t13*t10;
3741 t17 = t13*t13;
3742 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t7*t6+
3743 56.0/17.0*t7*t5*t10+147.0/34.0*t7*s*t13+49.0/17.0*t7*t15+35.0/34.0*t6*t17+42.0/
3744 221.0*t5*t17*t10+7.0/442.0*s*t17*t13+t17*t15/2431);
3745 /* p=15 5 0 6 */
3746 } else if( i==5 && j==0 && k==6 ) {
3747 t1 = r*r;
3748 t2 = t1*t1;
3749 t6 = t*t;
3750 t7 = t6*t6;
3751 t10 = s+r-1.0;
3752 t12 = t10*t10;
3753 t17 = t12*t12;
3754 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*(
3755 t7*t6+42.0/13.0*t7*t*t10+105.0/26.0*t7*t12+350.0/143.0*t6*t*t12*t10+105.0/143.0
3756 *t6*t17+14.0/143.0*t*t17*t10+7.0/1716.0*t17*t12);
3757 /* p=15 5 1 5 */
3758 } else if( i==5 && j==1 && k==5 ) {
3759 t1 = r*r;
3760 t2 = t1*t1;
3761 t8 = t*t;
3762 t9 = t8*t8;
3763 t11 = s+r-1.0;
3764 t14 = t11*t11;
3765 t18 = t14*t14;
3766 f = (s+2.0/15.0*r-2.0/15.0)*(t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0
3767 *t1+r/234-1.0/13455.0)*(t9*t+30.0/11.0*t11*t9+30.0/11.0*t8*t*t14+40.0/33.0*t8*
3768 t14*t11+5.0/22.0*t*t18+t18*t11/77);
3769 /* p=15 5 2 4 */
3770 } else if( i==5 && j==2 && k==4 ) {
3771 t1 = r*r;
3772 t2 = t1*t1;
3773 t6 = s*s;
3774 t7 = r-1.0;
3775 t9 = t7*t7;
3776 t12 = t*t;
3777 t13 = t12*t12;
3778 t15 = s+r-1.0;
3779 t17 = t15*t15;
3780 t21 = t17*t17;
3781 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*(
3782 t6+2.0/5.0*s*t7+t9/35)*(t13+20.0/9.0*t12*t*t15+5.0/3.0*t17*t12+10.0/21.0*t*t17*
3783 t15+5.0/126.0*t21);
3784 /* p=15 5 3 3 */
3785 } else if( i==5 && j==3 && k==3 ) {
3786 t1 = r*r;
3787 t2 = t1*t1;
3788 t6 = s*s;
3789 t8 = r-1.0;
3790 t10 = t8*t8;
3791 t15 = t*t;
3792 t17 = s+r-1.0;
3793 t19 = t17*t17;
3794 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*(
3795 s*t6+4.0/5.0*t6*t8+6.0/35.0*s*t10+4.0/455.0*t10*t8)*(t15*t+12.0/7.0*t17*t15+6.0
3796 /7.0*t*t19+4.0/35.0*t19*t17);
3797 /* p=15 5 4 2 */
3798 } else if( i==5 && j==4 && k==2 ) {
3799 t1 = r*r;
3800 t2 = t1*t1;
3801 t6 = s*s;
3802 t7 = t6*t6;
3803 t9 = r-1.0;
3804 t11 = t9*t9;
3805 t15 = t11*t11;
3806 t18 = t*t;
3807 t19 = s+r-1.0;
3808 t21 = t19*t19;
3809 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*(
3810 t7+4.0/3.0*t6*s*t9+4.0/7.0*t6*t11+8.0/91.0*s*t11*t9+t15/273)*(t18+6.0/5.0*t*t19
3811 +3.0/10.0*t21);
3812 /* p=15 5 5 1 */
3813 } else if( i==5 && j==5 && k==1 ) {
3814 t1 = r*r;
3815 t2 = t1*t1;
3816 t6 = s*s;
3817 t7 = t6*t6;
3818 t9 = r-1.0;
3819 t12 = t9*t9;
3820 t16 = t12*t12;
3821 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*(
3822 t7*s+2.0*t7*t9+10.0/7.0*t6*s*t12+40.0/91.0*t6*t12*t9+5.0/91.0*s*t16+2.0/1001.0*
3823 t16*t9)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
3824 /* p=15 5 6 0 */
3825 } else if( i==5 && j==6 && k==0 ) {
3826 t1 = r*r;
3827 t2 = t1*t1;
3828 t6 = s*s;
3829 t7 = t6*t6;
3830 t10 = r-1.0;
3831 t12 = t10*t10;
3832 t17 = t12*t12;
3833 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*(
3834 t7*t6+14.0/5.0*t7*s*t10+3.0*t7*t12+20.0/13.0*t6*s*t12*t10+5.0/13.0*t6*t17+6.0/
3835 143.0*s*t17*t10+t17*t12/715);
3836 /* p=15 6 0 5 */
3837 } else if( i==6 && j==0 && k==5 ) {
3838 t1 = r*r;
3839 t2 = t1*t1;
3840 t7 = t*t;
3841 t8 = t7*t7;
3842 t10 = s+r-1.0;
3843 t13 = t10*t10;
3844 t17 = t13*t13;
3845 f = (t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/117.0*t1*r+7.0/234.0*t1-7.0/
3846 4485.0*r+7.0/296010.0)*(t8*t+30.0/11.0*t8*t10+30.0/11.0*t7*t*t13+40.0/33.0*t7*
3847 t13*t10+5.0/22.0*t*t17+t17*t10/77);
3848 /* p=15 6 1 4 */
3849 } else if( i==6 && j==1 && k==4 ) {
3850 t1 = r*r;
3851 t2 = t1*t1;
3852 t9 = t*t;
3853 t10 = t9*t9;
3854 t12 = s+r-1.0;
3855 t14 = t12*t12;
3856 t18 = t14*t14;
3857 f = (s+2.0/13.0*r-2.0/13.0)*(t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/
3858 117.0*t1*r+7.0/234.0*t1-7.0/4485.0*r+7.0/296010.0)*(t10+20.0/9.0*t9*t*t12+5.0/
3859 3.0*t9*t14+10.0/21.0*t*t14*t12+5.0/126.0*t18);
3860 /* p=15 6 2 3 */
3861 } else if( i==6 && j==2 && k==3 ) {
3862 t1 = r*r;
3863 t2 = t1*t1;
3864 t7 = s*s;
3865 t8 = r-1.0;
3866 t10 = t8*t8;
3867 t13 = t*t;
3868 t15 = s+r-1.0;
3869 t17 = t15*t15;
3870 f = (t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/117.0*t1*r+7.0/234.0*t1-7.0/
3871 4485.0*r+7.0/296010.0)*(t7+6.0/13.0*s*t8+t10/26)*(t13*t+12.0/7.0*t13*t15+6.0/
3872 7.0*t*t17+4.0/35.0*t17*t15);
3873 /* p=15 6 3 2 */
3874 } else if( i==6 && j==3 && k==2 ) {
3875 t1 = r*r;
3876 t2 = t1*t1;
3877 t7 = s*s;
3878 t9 = r-1.0;
3879 t11 = t9*t9;
3880 t16 = t*t;
3881 t17 = s+r-1.0;
3882 t19 = t17*t17;
3883 f = (t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/117.0*t1*r+7.0/234.0*t1-7.0/
3884 4485.0*r+7.0/296010.0)*(t7*s+12.0/13.0*t7*t9+3.0/13.0*s*t11+2.0/143.0*t11*t9)*(
3885 t16+6.0/5.0*t*t17+3.0/10.0*t19);
3886 /* p=15 6 4 1 */
3887 } else if( i==6 && j==4 && k==1 ) {
3888 t1 = r*r;
3889 t2 = t1*t1;
3890 t7 = s*s;
3891 t8 = t7*t7;
3892 t10 = r-1.0;
3893 t12 = t10*t10;
3894 t16 = t12*t12;
3895 f = (t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/117.0*t1*r+7.0/234.0*t1-7.0/
3896 4485.0*r+7.0/296010.0)*(t8+20.0/13.0*t7*s*t10+10.0/13.0*t7*t12+20.0/143.0*s*t12
3897 *t10+t16/143)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
3898 /* p=15 6 5 0 */
3899 } else if( i==6 && j==5 && k==0 ) {
3900 t1 = r*r;
3901 t2 = t1*t1;
3902 t7 = s*s;
3903 t8 = t7*t7;
3904 t10 = r-1.0;
3905 t13 = t10*t10;
3906 t17 = t13*t13;
3907 f = (t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/117.0*t1*r+7.0/234.0*t1-7.0/
3908 4485.0*r+7.0/296010.0)*(t8*s+30.0/13.0*t8*t10+25.0/13.0*t7*s*t13+100.0/143.0*t7
3909 *t13*t10+15.0/143.0*s*t17+2.0/429.0*t17*t10);
3910 /* p=15 7 0 4 */
3911 } else if( i==7 && j==0 && k==4 ) {
3912 t1 = r*r;
3913 t2 = t1*r;
3914 t3 = t1*t1;
3915 t8 = t*t;
3916 t9 = t8*t8;
3917 t11 = s+r-1.0;
3918 t13 = t11*t11;
3919 t17 = t13*t13;
3920 f = (t3*t2-56.0/27.0*t3*t1+196.0/117.0*t3*r-392.0/585.0*t3+49.0/351.0*
3921 t2-196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0)*(t9+20.0/9.0*t8*t*t11+5.0/3.0
3922 *t8*t13+10.0/21.0*t*t13*t11+5.0/126.0*t17);
3923 /* p=15 7 1 3 */
3924 } else if( i==7 && j==1 && k==3 ) {
3925 t1 = r*r;
3926 t2 = t1*r;
3927 t3 = t1*t1;
3928 t10 = t*t;
3929 t12 = s+r-1.0;
3930 t14 = t12*t12;
3931 f = (t3*t2-56.0/27.0*t3*t1+196.0/117.0*t3*r-392.0/585.0*t3+49.0/351.0*
3932 t2-196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0)*(s+2.0/11.0*r-2.0/11.0)*(t10*
3933 t+12.0/7.0*t12*t10+6.0/7.0*t*t14+4.0/35.0*t14*t12);
3934 /* p=15 7 2 2 */
3935 } else if( i==7 && j==2 && k==2 ) {
3936 t1 = r*r;
3937 t2 = t1*r;
3938 t3 = t1*t1;
3939 t8 = s*s;
3940 t9 = r-1.0;
3941 t11 = t9*t9;
3942 t14 = t*t;
3943 t15 = s+r-1.0;
3944 t17 = t15*t15;
3945 f = (t3*t2-56.0/27.0*t3*t1+196.0/117.0*t3*r-392.0/585.0*t3+49.0/351.0*
3946 t2-196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0)*(t8+6.0/11.0*s*t9+3.0/55.0*
3947 t11)*(t14+6.0/5.0*t*t15+3.0/10.0*t17);
3948 /* p=15 7 3 1 */
3949 } else if( i==7 && j==3 && k==1 ) {
3950 t1 = r*r;
3951 t2 = t1*r;
3952 t3 = t1*t1;
3953 t8 = s*s;
3954 t10 = r-1.0;
3955 t12 = t10*t10;
3956 f = (t3*t2-56.0/27.0*t3*t1+196.0/117.0*t3*r-392.0/585.0*t3+49.0/351.0*
3957 t2-196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0)*(t8*s+12.0/11.0*t8*t10+18.0/
3958 55.0*s*t12+4.0/165.0*t12*t10)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
3959 /* p=15 7 4 0 */
3960 } else if( i==7 && j==4 && k==0 ) {
3961 t1 = r*r;
3962 t2 = t1*r;
3963 t3 = t1*t1;
3964 t8 = s*s;
3965 t9 = t8*t8;
3966 t11 = r-1.0;
3967 t13 = t11*t11;
3968 t17 = t13*t13;
3969 f = (t3*t2-56.0/27.0*t3*t1+196.0/117.0*t3*r-392.0/585.0*t3+49.0/351.0*
3970 t2-196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0)*(t9+20.0/11.0*t8*s*t11+12.0/
3971 11.0*t8*t13+8.0/33.0*s*t13*t11+t17/66);
3972 /* p=15 8 0 3 */
3973 } else if( i==8 && j==0 && k==3 ) {
3974 t1 = r*r;
3975 t2 = t1*t1;
3976 t3 = t2*t2;
3977 t4 = t1*r;
3978 t9 = t*t;
3979 t11 = s+r-1.0;
3980 t13 = t11*t11;
3981 f = (t3-8.0/3.0*t2*t4+112.0/39.0*t2*t1-1568.0/975.0*t2*r+98.0/195.0*t2
3982 -392.0/4485.0*t4+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0)*(t9*t+12.0/7.0*
3983 t9*t11+6.0/7.0*t*t13+4.0/35.0*t13*t11);
3984 /* p=15 8 1 2 */
3985 } else if( i==8 && j==1 && k==2 ) {
3986 t1 = r*r;
3987 t2 = t1*t1;
3988 t3 = t2*t2;
3989 t4 = t1*r;
3990 t11 = t*t;
3991 t12 = s+r-1.0;
3992 t14 = t12*t12;
3993 f = (t3-8.0/3.0*t2*t4+112.0/39.0*t2*t1-1568.0/975.0*t2*r+98.0/195.0*t2
3994 -392.0/4485.0*t4+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0)*(s+2.0/9.0*r-2.0
3995 /9.0)*(t11+6.0/5.0*t*t12+3.0/10.0*t14);
3996 /* p=15 8 2 1 */
3997 } else if( i==8 && j==2 && k==1 ) {
3998 t1 = r*r;
3999 t2 = t1*t1;
4000 t3 = t2*t2;
4001 t4 = t1*r;
4002 t9 = s*s;
4003 t10 = r-1.0;
4004 t12 = t10*t10;
4005 f = (t3-8.0/3.0*t2*t4+112.0/39.0*t2*t1-1568.0/975.0*t2*r+98.0/195.0*t2
4006 -392.0/4485.0*t4+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0)*(t9+2.0/3.0*s*
4007 t10+t12/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
4008 /* p=15 8 3 0 */
4009 } else if( i==8 && j==3 && k==0 ) {
4010 t1 = r*r;
4011 t2 = t1*t1;
4012 t3 = t2*t2;
4013 t4 = t1*r;
4014 t9 = s*s;
4015 t11 = r-1.0;
4016 t13 = t11*t11;
4017 f = (t3-8.0/3.0*t2*t4+112.0/39.0*t2*t1-1568.0/975.0*t2*r+98.0/195.0*t2
4018 -392.0/4485.0*t4+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0)*(t9*s+4.0/3.0*t9
4019 *t11+s*t13/2+t13*t11/21);
4020 /* p=15 9 0 2 */
4021 } else if( i==9 && j==0 && k==2 ) {
4022 t1 = r*r;
4023 t2 = t1*t1;
4024 t3 = t2*t2;
4025 t5 = t1*r;
4026 t10 = t*t;
4027 t11 = s+r-1.0;
4028 t13 = t11*t11;
4029 f = (t3*r-10.0/3.0*t3+60.0/13.0*t2*t5-224.0/65.0*t2*t1+98.0/65.0*t2*r
4030 -588.0/1495.0*t2+196.0/3289.0*t5-16.0/3289.0*t1+3.0/16445.0*r-2.0/937365.0)*(
4031 t10+6.0/5.0*t*t11+3.0/10.0*t13);
4032 /* p=15 9 1 1 */
4033 } else if( i==9 && j==1 && k==1 ) {
4034 t1 = r*r;
4035 t2 = t1*t1;
4036 t3 = t2*t2;
4037 t5 = t1*r;
4038 f = (t3*r-10.0/3.0*t3+60.0/13.0*t2*t5-224.0/65.0*t2*t1+98.0/65.0*t2*r
4039 -588.0/1495.0*t2+196.0/3289.0*t5-16.0/3289.0*t1+3.0/16445.0*r-2.0/937365.0)*(s+
4040 2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
4041 /* p=15 9 2 0 */
4042 } else if( i==9 && j==2 && k==0 ) {
4043 t1 = r*r;
4044 t2 = t1*t1;
4045 t3 = t2*t2;
4046 t5 = t1*r;
4047 t10 = s*s;
4048 t11 = r-1.0;
4049 t13 = t11*t11;
4050 f = (t3*r-10.0/3.0*t3+60.0/13.0*t2*t5-224.0/65.0*t2*t1+98.0/65.0*t2*r
4051 -588.0/1495.0*t2+196.0/3289.0*t5-16.0/3289.0*t1+3.0/16445.0*r-2.0/937365.0)*(
4052 t10+6.0/7.0*s*t11+t13/7);
4053 /* p=15 10 0 1 */
4054 } else if( i==10 && j==0 && k==1 ) {
4055 t1 = r*r;
4056 t2 = t1*t1;
4057 t3 = t2*t2;
4058 t6 = t1*r;
4059 t10 = t3*t1-110.0/27.0*t3*r+275.0/39.0*t3-88.0/13.0*t2*t6+154.0/39.0*t2*
4060 t1-2156.0/1495.0*t2*r+98.0/299.0*t2-40.0/897.0*t6+t1/299-2.0/17043.0*r+1.0/
4061 766935.0;
4062 f = t10*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);
4063 /* p=15 10 1 0 */
4064 } else if( i==10 && j==1 && k==0 ) {
4065 t1 = r*r;
4066 t2 = t1*t1;
4067 t3 = t2*t2;
4068 t6 = t1*r;
4069 t10 = t3*t1-110.0/27.0*t3*r+275.0/39.0*t3-88.0/13.0*t2*t6+154.0/39.0*t2*
4070 t1-2156.0/1495.0*t2*r+98.0/299.0*t2-40.0/897.0*t6+t1/299-2.0/17043.0*r+1.0/
4071 766935.0;
4072 f = t10*(s+2.0/5.0*r-2.0/5.0);
4073 /* p=15 11 0 0 */
4074 } else if( i==11 && j==0 && k==0 ) {
4075 t1 = r*r;
4076 t2 = t1*r;
4077 t3 = t1*t1;
4078 t4 = t3*t3;
4079 f = t4*t2-44.0/9.0*t4*t1+1210.0/117.0*t4*r-484.0/39.0*t4+121.0/13.0*t3*
4080 t2-6776.0/1495.0*t3*t1+2156.0/1495.0*t3*r-88.0/299.0*t3+11.0/299.0*t2-44.0/
4081 17043.0*t1+22.0/255645.0*r-4.0/4345965.0;
4082 }
4083 return f ;
4084 }
4085
4086 #ifdef __cplusplus
4087 }
4088 #endif
4089