xref: /phasta/shapeFunction/src/FnDrv.c (revision 1e99f302ca5103688ae35115c2fefb7cfa6714f1) 
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              derivative of a tri. face mode shape. i,j are the highest poly.
11              order of r,s in the polynomial. maple generated code.
12 -------------------------------------------------------------------------*/
13 
14 #ifdef __cplusplus
15 extern "C" {
16 #endif
17 
18 int FnDrv(int i,int j, double r, double s, double drv[2]) {
19    double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18;
20    double t19,t20,t21,t22,t23,t24,t25,t26,t27,t28,t29,t30,t31,t32,t33,t34;
21    double t35,t36,t37,t38,t39,t40,t41,t42,t43,t44,t45,t46,t47,t48,t50,t52;
22    double t53,t57;
23 
24    /* p=3 0 0 */
25    if( i==0 && j==0 ) {
26       drv[0] = 0.0;
27       drv[1] = 0.0;
28    /* p=4 0 1 */
29    } else if( i==0 && j==1 ) {
30       drv[0] = 0.0;
31       drv[1] = 1.0;
32    /* p=4 1 0 */
33    } else if( i==1 && j==0 ) {
34       drv[0] = 1.0;
35       drv[1] = 0.0;
36    /* p=5 0 2 */
37    } else if( i==0 && j==2 ) {
38       drv[0] = 0.0;
39       drv[1] = 2.0*s-3.0/4.0;
40    /* p=5 1 1 */
41    } else if( i==1 && j==1 ) {
42       drv[0] = s-1.0/4.0;
43       drv[1] = r-1.0/4.0;
44    /* p=5 2 0 */
45    } else if( i==2 && j==0 ) {
46       drv[0] = 2.0*r-3.0/4.0;
47       drv[1] = 0.0;
48    /* p=6 0 3 */
49    } else if( i==0 && j==3 ) {
50       t1 = s*s;
51       drv[0] = 0.0;
52       drv[1] = 3.0*t1-12.0/5.0*s+2.0/5.0;
53    /* p=6 1 2 */
54    } else if( i==1 && j==2 ) {
55       t1 = s*s;
56       drv[0] = t1-3.0/5.0*s+1.0/15.0;
57       drv[1] = 2.0*r*s-3.0/5.0*r-2.0/5.0*s+2.0/15.0;
58    /* p=6 2 1 */
59    } else if( i==2 && j==1 ) {
60       t3 = r*r;
61       drv[0] = 2.0*r*s-2.0/5.0*r-3.0/5.0*s+2.0/15.0;
62       drv[1] = t3-3.0/5.0*r+1.0/15.0;
63    /* p=6 3 0 */
64    } else if( i==3 && j==0 ) {
65       t1 = r*r;
66       drv[0] = 3.0*t1-12.0/5.0*r+2.0/5.0;
67       drv[1] = 0.0;
68    /* p=7 0 4 */
69    } else if( i==0 && j==4 ) {
70       t1 = s*s;
71       drv[0] = 0.0;
72       drv[1] = 4.0*t1*s-5.0*t1+20.0/11.0*s-2.0/11.0;
73    /* p=7 1 3 */
74    } else if( i==1 && j==3 ) {
75       t1 = s*s;
76       drv[0] = t1*s-t1+3.0/11.0*s-1.0/55.0;
77       drv[1] = 3.0*r*t1-2.0*r*s+3.0/11.0*r-t1/2+4.0/11.0*s-3.0/55.0;
78    /* p=7 2 2 */
79    } else if( i==2 && j==2 ) {
80       t1 = s*s;
81       t3 = r*s;
82       t5 = r*r;
83       drv[0] = 2.0*r*t1-t3+r/11-t1/2+3.0/11.0*s-3.0/110.0;
84       drv[1] = 2.0*t5*s-t5/2-t3+3.0/11.0*r+s/11-3.0/110.0;
85    /* p=7 3 1 */
86    } else if( i==3 && j==1 ) {
87       t1 = r*r;
88       drv[0] = 3.0*t1*s-t1/2-2.0*r*s+4.0/11.0*r+3.0/11.0*s-3.0/55.0;
89       drv[1] = r*t1-t1+3.0/11.0*r-1.0/55.0;
90    /* p=7 4 0 */
91    } else if( i==4 && j==0 ) {
92       t1 = r*r;
93       drv[0] = 4.0*r*t1-5.0*t1+20.0/11.0*r-2.0/11.0;
94       drv[1] = 0.0;
95    /* p=8 0 5 */
96    } else if( i==0 && j==5 ) {
97       t1 = s*s;
98       t2 = t1*t1;
99       drv[0] = 0.0;
100       drv[1] = 5.0*t2-60.0/7.0*t1*s+450.0/91.0*t1-100.0/91.0*s+75.0/1001.0;
101    /* p=8 1 4 */
102    } else if( i==1 && j==4 ) {
103       t1 = s*s;
104       t2 = t1*t1;
105       t3 = t1*s;
106       drv[0] = t2-10.0/7.0*t3+60.0/91.0*t1-10.0/91.0*s+5.0/1001.0;
107       drv[1] = 4.0*r*t3-30.0/7.0*r*t1+120.0/91.0*r*s-10.0/91.0*r-4.0/7.0*t3+
108 60.0/91.0*t1-20.0/91.0*s+20.0/1001.0;
109    /* p=8 2 3 */
110    } else if( i==2 && j==3 ) {
111       t1 = s*s;
112       t2 = t1*s;
113       t4 = r*t1;
114       t5 = r*s;
115       t7 = r*r;
116       drv[0] = 2.0*r*t2-12.0/7.0*t4+36.0/91.0*t5-2.0/91.0*r-3.0/7.0*t2+36.0/
117 91.0*t1-9.0/91.0*s+6.0/1001.0;
118       drv[1] = 3.0*t7*t1-12.0/7.0*t7*s+18.0/91.0*t7-9.0/7.0*t4+72.0/91.0*t5-9.0
119 /91.0*r+9.0/91.0*t1-6.0/91.0*s+9.0/1001.0;
120    /* p=8 3 2 */
121    } else if( i==3 && j==2 ) {
122       t1 = r*r;
123       t2 = s*s;
124       t4 = t1*s;
125       t6 = r*s;
126       t8 = r*t1;
127       drv[0] = 3.0*t1*t2-9.0/7.0*t4+9.0/91.0*t1-12.0/7.0*r*t2+72.0/91.0*t6-6.0/
128 91.0*r+18.0/91.0*t2-9.0/91.0*s+9.0/1001.0;
129       drv[1] = 2.0*t8*s-3.0/7.0*t8-12.0/7.0*t4+36.0/91.0*t1+36.0/91.0*t6-9.0/
130 91.0*r-2.0/91.0*s+6.0/1001.0;
131    /* p=8 4 1 */
132    } else if( i==4 && j==1 ) {
133       t1 = r*r;
134       t2 = r*t1;
135       t7 = t1*t1;
136       drv[0] = 4.0*t2*s-4.0/7.0*t2-30.0/7.0*t1*s+60.0/91.0*t1+120.0/91.0*r*s
137 -20.0/91.0*r-10.0/91.0*s+20.0/1001.0;
138       drv[1] = t7-10.0/7.0*t2+60.0/91.0*t1-10.0/91.0*r+5.0/1001.0;
139    /* p=8 5 0 */
140    } else if( i==5 && j==0 ) {
141       t1 = r*r;
142       t2 = t1*t1;
143       drv[0] = 5.0*t2-60.0/7.0*r*t1+450.0/91.0*t1-100.0/91.0*r+75.0/1001.0;
144       drv[1] = 0.0;
145    /* p=9 0 6 */
146    } else if( i==0 && j==6 ) {
147       t1 = s*s;
148       t2 = t1*t1;
149       drv[0] = 0.0;
150       drv[1] = 6.0*t2*s-105.0/8.0*t2+21.0/2.0*t1*s-15.0/4.0*t1+15.0/26.0*s-3.0/
151 104.0;
152    /* p=9 1 5 */
153    } else if( i==1 && j==5 ) {
154       t1 = s*s;
155       t2 = t1*t1;
156       t4 = t1*s;
157       drv[0] = t2*s-15.0/8.0*t2+5.0/4.0*t4-5.0/14.0*t1+15.0/364.0*s-1.0/728.0;
158       drv[1] = 5.0*r*t2-15.0/2.0*r*t4+15.0/4.0*r*t1-5.0/7.0*r*s+15.0/364.0*r
159 -5.0/8.0*t2+t4-15.0/28.0*t1+10.0/91.0*s-5.0/728.0;
160    /* p=9 2 4 */
161    } else if( i==2 && j==4 ) {
162       t1 = s*s;
163       t2 = t1*t1;
164       t4 = t1*s;
165       t5 = r*t4;
166       t6 = r*t1;
167       t7 = r*s;
168       t9 = r*r;
169       t13 = 4.0*t9*t4-15.0/4.0*t9*t1+t9*s-t9/14-3.0/2.0*t5+3.0/2.0*t6-3.0/7.0*
170 t7+3.0/91.0*r+t4/10-3.0/28.0*t1+3.0/91.0*s-1.0/364.0;
171       drv[0] = 2.0*r*t2-5.0/2.0*t5+t6-t7/7+r/182-3.0/8.0*t2+t4/2-3.0/14.0*t1+
172 3.0/91.0*s-1.0/728.0;
173       drv[1] = t13;
174    /* p=9 3 3 */
175    } else if( i==3 && j==3 ) {
176       t1 = r*r;
177       t2 = s*s;
178       t3 = t2*s;
179       t5 = t1*t2;
180       t6 = t1*s;
181       t8 = r*t2;
182       t9 = r*s;
183       t10 = 3.0*t1*t3-9.0/4.0*t5+9.0/20.0*t6-3.0/140.0*t1-3.0/2.0*r*t3+6.0/5.0*
184 t8-9.0/35.0*t9+6.0/455.0*r+3.0/20.0*t3-9.0/70.0*t2+27.0/910.0*s-3.0/1820.0;
185       t11 = r*t1;
186       t14 = 3.0*t11*t2-3.0/2.0*t11*s+3.0/20.0*t11-9.0/4.0*t5+6.0/5.0*t6-9.0/
187 70.0*t1+9.0/20.0*t8-9.0/35.0*t9+27.0/910.0*r-3.0/140.0*t2+6.0/455.0*s-3.0/
188 1820.0;
189       drv[0] = t10;
190       drv[1] = t14;
191    /* p=9 4 2 */
192    } else if( i==4 && j==2 ) {
193       t1 = r*r;
194       t2 = r*t1;
195       t3 = s*s;
196       t5 = t2*s;
197       t7 = t1*s;
198       t9 = r*s;
199       t10 = 4.0*t2*t3-3.0/2.0*t5+t2/10-15.0/4.0*t1*t3+3.0/2.0*t7-3.0/28.0*t1+r*
200 t3-3.0/7.0*t9+3.0/91.0*r-t3/14+3.0/91.0*s-1.0/364.0;
201       t11 = t1*t1;
202       drv[0] = t10;
203       drv[1] = 2.0*t11*s-3.0/8.0*t11-5.0/2.0*t5+t2/2+t7-3.0/14.0*t1-t9/7+3.0/
204 91.0*r+s/182-1.0/728.0;
205    /* p=9 5 1 */
206    } else if( i==5 && j==1 ) {
207       t1 = r*r;
208       t2 = t1*t1;
209       t4 = r*t1;
210       drv[0] = 5.0*t2*s-5.0/8.0*t2-15.0/2.0*t4*s+t4+15.0/4.0*t1*s-15.0/28.0*t1
211 -5.0/7.0*r*s+10.0/91.0*r+15.0/364.0*s-5.0/728.0;
212       drv[1] = r*t2-15.0/8.0*t2+5.0/4.0*t4-5.0/14.0*t1+15.0/364.0*r-1.0/728.0;
213    /* p=9 6 0 */
214    } else if( i==6 && j==0 ) {
215       t1 = r*r;
216       t2 = t1*t1;
217       drv[0] = 6.0*r*t2-105.0/8.0*t2+21.0/2.0*r*t1-15.0/4.0*t1+15.0/26.0*r-3.0/
218 104.0;
219       drv[1] = 0.0;
220    /* p=10 0 7 */
221    } else if( i==0 && j==7 ) {
222       t1 = s*s;
223       t2 = t1*t1;
224       drv[0] = 0.0;
225       drv[1] = 7.0*t1*t2-56.0/3.0*t2*s+980.0/51.0*t2-490.0/51.0*t1*s+245.0/
226 102.0*t1-14.0/51.0*s+7.0/663.0;
227    /* p=10 1 6 */
228    } else if( i==1 && j==6 ) {
229       t1 = s*s;
230       t2 = t1*t1;
231       t4 = t2*s;
232       t5 = t1*s;
233       t12 = 6.0*r*t4-35.0/3.0*r*t2+140.0/17.0*r*t5-175.0/68.0*r*t1+35.0/102.0*r
234 *s-r/68-2.0/3.0*t4+70.0/51.0*t2-35.0/34.0*t5+35.0/102.0*t1-5.0/102.0*s+1.0/
235 442.0;
236       drv[0] = t1*t2-7.0/3.0*t4+35.0/17.0*t2-175.0/204.0*t5+35.0/204.0*t1-s/68+
237 1.0/2652.0;
238       drv[1] = t12;
239    /* p=10 2 5 */
240    } else if( i==2 && j==5 ) {
241       t1 = s*s;
242       t2 = t1*t1;
243       t3 = t2*s;
244       t5 = r*t2;
245       t6 = t1*s;
246       t7 = r*t6;
247       t8 = r*t1;
248       t9 = r*s;
249       t10 = 2.0*r*t3-10.0/3.0*t5+100.0/51.0*t7-25.0/51.0*t8+5.0/102.0*t9-r/714-
250 t3/3+10.0/17.0*t2-25.0/68.0*t6+5.0/51.0*t1-5.0/476.0*s+1.0/3094.0;
251       t11 = r*r;
252       t16 = 5.0*t11*t2-20.0/3.0*t11*t6+50.0/17.0*t11*t1-25.0/51.0*t11*s+5.0/
253 204.0*t11-5.0/3.0*t5+40.0/17.0*t7-75.0/68.0*t8+10.0/51.0*t9-5.0/476.0*r+5.0/
254 51.0*t2-5.0/34.0*t6+5.0/68.0*t1-5.0/357.0*s+5.0/6188.0;
255       drv[0] = t10;
256       drv[1] = t16;
257    /* p=10 3 4 */
258    } else if( i==3 && j==4 ) {
259       t1 = r*r;
260       t2 = s*s;
261       t3 = t2*t2;
262       t5 = t2*s;
263       t6 = t1*t5;
264       t7 = t1*t2;
265       t8 = t1*s;
266       t10 = r*t5;
267       t11 = r*t2;
268       t12 = r*s;
269       t13 = 3.0*t1*t3-10.0/3.0*t6+20.0/17.0*t7-5.0/34.0*t8+t1/204-4.0/3.0*r*t3+
270 80.0/51.0*t10-10.0/17.0*t11+4.0/51.0*t12-r/357+2.0/17.0*t3-5.0/34.0*t5+t2/17-s
271 /119+1.0/3094.0;
272       t14 = r*t1;
273       t18 = 4.0*t14*t5-10.0/3.0*t14*t2+40.0/51.0*t14*s-5.0/102.0*t14-8.0/3.0*t6
274 +40.0/17.0*t7-10.0/17.0*t8+2.0/51.0*t1+8.0/17.0*t10-15.0/34.0*t11+2.0/17.0*t12-
275 r/119-t5/51+t2/51-2.0/357.0*s+2.0/4641.0;
276       drv[0] = t13;
277       drv[1] = t18;
278    /* p=10 4 3 */
279    } else if( i==4 && j==3 ) {
280       t1 = r*r;
281       t2 = r*t1;
282       t3 = s*s;
283       t4 = t3*s;
284       t6 = t2*t3;
285       t7 = t2*s;
286       t9 = t1*t3;
287       t10 = t1*s;
288       t12 = r*t3;
289       t13 = r*s;
290       t14 = 4.0*t2*t4-8.0/3.0*t6+8.0/17.0*t7-t2/51-10.0/3.0*t1*t4+40.0/17.0*t9
291 -15.0/34.0*t10+t1/51+40.0/51.0*r*t4-10.0/17.0*t12+2.0/17.0*t13-2.0/357.0*r-5.0/
292 102.0*t4+2.0/51.0*t3-s/119+2.0/4641.0;
293       t15 = t1*t1;
294       t18 = 3.0*t15*t3-4.0/3.0*t15*s+2.0/17.0*t15-10.0/3.0*t6+80.0/51.0*t7-5.0/
295 34.0*t2+20.0/17.0*t9-10.0/17.0*t10+t1/17-5.0/34.0*t12+4.0/51.0*t13-r/119+t3/204
296 -s/357+1.0/3094.0;
297       drv[0] = t14;
298       drv[1] = t18;
299    /* p=10 5 2 */
300    } else if( i==5 && j==2 ) {
301       t1 = r*r;
302       t2 = t1*t1;
303       t3 = s*s;
304       t5 = t2*s;
305       t6 = t1*r;
306       t8 = t6*s;
307       t10 = t1*s;
308       t12 = r*s;
309       t13 = 5.0*t2*t3-5.0/3.0*t5+5.0/51.0*t2-20.0/3.0*t6*t3+40.0/17.0*t8-5.0/
310 34.0*t6+50.0/17.0*t1*t3-75.0/68.0*t10+5.0/68.0*t1-25.0/51.0*r*t3+10.0/51.0*t12
311 -5.0/357.0*r+5.0/204.0*t3-5.0/476.0*s+5.0/6188.0;
312       t14 = t2*r;
313       t16 = 2.0*t14*s-t14/3-10.0/3.0*t5+10.0/17.0*t2+100.0/51.0*t8-25.0/68.0*t6
314 -25.0/51.0*t10+5.0/51.0*t1+5.0/102.0*t12-5.0/476.0*r-s/714+1.0/3094.0;
315       drv[0] = t13;
316       drv[1] = t16;
317    /* p=10 6 1 */
318    } else if( i==6 && j==1 ) {
319       t1 = r*r;
320       t2 = t1*t1;
321       t3 = t2*r;
322       t6 = t1*r;
323       t10 = 6.0*t3*s-2.0/3.0*t3-35.0/3.0*t2*s+70.0/51.0*t2+140.0/17.0*t6*s-35.0
324 /34.0*t6-175.0/68.0*t1*s+35.0/102.0*t1+35.0/102.0*r*s-5.0/102.0*r-s/68+1.0/
325 442.0;
326       drv[0] = t10;
327       drv[1] = t2*t1-7.0/3.0*t3+35.0/17.0*t2-175.0/204.0*t6+35.0/204.0*t1-r/68+
328 1.0/2652.0;
329    /* p=10 7 0 */
330    } else if( i==7 && j==0 ) {
331       t1 = r*r;
332       t2 = t1*t1;
333       drv[0] = 7.0*t2*t1-56.0/3.0*t2*r+980.0/51.0*t2-490.0/51.0*t1*r+245.0/
334 102.0*t1-14.0/51.0*r+7.0/663.0;
335       drv[1] = 0.0;
336    /* p=11 0 8 */
337    } else if( i==0 && j==8 ) {
338       t1 = s*s;
339       t2 = t1*s;
340       t3 = t1*t1;
341       drv[0] = 0.0;
342       drv[1] = 8.0*t2*t3-126.0/5.0*t1*t3+3024.0/95.0*t3*s-392.0/19.0*t3+2352.0/
343 323.0*t2-441.0/323.0*t1+196.0/1615.0*s-6.0/1615.0;
344    /* p=11 1 7 */
345    } else if( i==1 && j==7 ) {
346       t1 = s*s;
347       t2 = t1*s;
348       t3 = t1*t1;
349       t5 = t1*t3;
350       t6 = t3*s;
351       t14 = 7.0*r*t5-84.0/5.0*r*t6+294.0/19.0*r*t3-392.0/57.0*t2*r+490.0/323.0*
352 t1*r-49.0/323.0*r*s+49.0/9690.0*r-7.0/10.0*t5+168.0/95.0*t6-98.0/57.0*t3+784.0/
353 969.0*t2-245.0/1292.0*t1+98.0/4845.0*s-7.0/9690.0;
354       drv[0] = t2*t3-14.0/5.0*t5+294.0/95.0*t6-98.0/57.0*t3+490.0/969.0*t2-49.0
355 /646.0*t1+49.0/9690.0*s-1.0/9690.0;
356       drv[1] = t14;
357    /* p=11 2 6 */
358    } else if( i==2 && j==6 ) {
359       t1 = s*s;
360       t2 = t1*t1;
361       t3 = t2*t1;
362       t5 = t2*s;
363       t6 = r*t5;
364       t7 = r*t2;
365       t8 = t1*s;
366       t9 = r*t8;
367       t10 = r*t1;
368       t11 = r*s;
369       t12 = 2.0*r*t3-21.0/5.0*t6+63.0/19.0*t7-70.0/57.0*t9+70.0/323.0*t10-21.0/
370 1292.0*t11+7.0/19380.0*r-3.0/10.0*t3+63.0/95.0*t5-21.0/38.0*t2+70.0/323.0*t8
371 -105.0/2584.0*t1+21.0/6460.0*s-1.0/12920.0;
372       t13 = r*r;
373       t19 = 6.0*t13*t5-21.0/2.0*t13*t2+126.0/19.0*t13*t8-35.0/19.0*t13*t1+70.0/
374 323.0*t13*s-21.0/2584.0*t13-9.0/5.0*t6+63.0/19.0*t7-42.0/19.0*t9+210.0/323.0*
375 t10-105.0/1292.0*t11+21.0/6460.0*r+9.0/95.0*t5-7.0/38.0*t2+42.0/323.0*t8-105.0/
376 2584.0*t1+7.0/1292.0*s-3.0/12920.0;
377       drv[0] = t12;
378       drv[1] = t19;
379    /* p=11 3 5 */
380    } else if( i==3 && j==5 ) {
381       t1 = r*r;
382       t2 = s*s;
383       t3 = t2*t2;
384       t4 = t3*s;
385       t6 = t1*t3;
386       t7 = t2*s;
387       t8 = t1*t7;
388       t9 = t2*t1;
389       t10 = t1*s;
390       t12 = r*t3;
391       t13 = r*t7;
392       t14 = r*t2;
393       t15 = r*s;
394       t16 = 3.0*t1*t4-9.0/2.0*t6+45.0/19.0*t8-10.0/19.0*t9+15.0/323.0*t10-3.0/
395 2584.0*t1-6.0/5.0*r*t4+36.0/19.0*t12-20.0/19.0*t13+80.0/323.0*t14-15.0/646.0*
396 t15+r/1615+9.0/95.0*t4-3.0/19.0*t3+30.0/323.0*t7-15.0/646.0*t2+3.0/1292.0*s-3.0
397 /45220.0;
398       t17 = r*t1;
399       t22 = 5.0*t17*t3-6.0*t17*t7+45.0/19.0*t17*t2-20.0/57.0*t17*s+5.0/323.0*
400 t17-3.0*t6+72.0/19.0*t8-30.0/19.0*t9+80.0/323.0*t10-15.0/1292.0*t1+9.0/19.0*t12
401 -12.0/19.0*t13+90.0/323.0*t14-15.0/323.0*t15+3.0/1292.0*r-t3/57+8.0/323.0*t7
402 -15.0/1292.0*t2+2.0/969.0*s-1.0/9044.0;
403       drv[0] = t16;
404       drv[1] = t22;
405    /* p=11 4 4 */
406    } else if( i==4 && j==4 ) {
407       t1 = r*r;
408       t2 = t1*s;
409       t3 = r*s;
410       t4 = s*s;
411       t5 = r*t4;
412       t6 = t1*t4;
413       t7 = r*t1;
414       t8 = t7*s;
415       t9 = t4*s;
416       t10 = t1*t9;
417       t11 = r*t9;
418       t12 = t7*t4;
419       t13 = t4*t4;
420       t16 = t7*t9;
421       t18 = 40.0/323.0*t2-1.0/13566.0-10.0/323.0*t3+80.0/323.0*t5+2.0/969.0*s
422 -20.0/19.0*t6+r/969-8.0/57.0*t8+40.0/969.0*t9+60.0/19.0*t10-5.0/323.0*t4-5.0/
423 1292.0*t1-40.0/57.0*t11+24.0/19.0*t12+4.0/969.0*t7+4.0*t7*t13-2.0/57.0*t13-3.0*
424 t13*t1-4.0*t16+12.0/19.0*r*t13;
425       t19 = t1*t1;
426       t23 = 40.0/969.0*t7-1.0/13566.0-5.0/323.0*t1+80.0/323.0*t2+s/969+2.0/
427 969.0*r-40.0/57.0*t8-2.0/57.0*t19-5.0/1292.0*t4+40.0/323.0*t5+4.0/969.0*t9+60.0
428 /19.0*t12-10.0/323.0*t3-20.0/19.0*t6+12.0/19.0*t19*s+4.0*t19*t9-8.0/57.0*t11
429 -4.0*t16-3.0*t19*t4+24.0/19.0*t10;
430       drv[0] = t18;
431       drv[1] = t23;
432    /* p=11 5 3 */
433    } else if( i==5 && j==3 ) {
434       t1 = r*r;
435       t2 = t1*t1;
436       t3 = s*s;
437       t4 = t3*s;
438       t6 = t2*t3;
439       t7 = t2*s;
440       t8 = r*t1;
441       t10 = t8*t3;
442       t11 = t8*s;
443       t13 = t1*t3;
444       t14 = t1*s;
445       t16 = r*t3;
446       t17 = r*s;
447       t18 = 5.0*t2*t4-3.0*t6+9.0/19.0*t7-t2/57-6.0*t8*t4+72.0/19.0*t10-12.0/
448 19.0*t11+8.0/323.0*t8+45.0/19.0*t1*t4-30.0/19.0*t13+90.0/323.0*t14-15.0/1292.0*
449 t1-20.0/57.0*r*t4+80.0/323.0*t16-15.0/323.0*t17+2.0/969.0*r+5.0/323.0*t4-15.0/
450 1292.0*t3+3.0/1292.0*s-1.0/9044.0;
451       t19 = t2*r;
452       t22 = 3.0*t19*t3-6.0/5.0*t19*s+9.0/95.0*t19-9.0/2.0*t6+36.0/19.0*t7-3.0/
453 19.0*t2+45.0/19.0*t10-20.0/19.0*t11+30.0/323.0*t8-10.0/19.0*t13+80.0/323.0*t14
454 -15.0/646.0*t1+15.0/323.0*t16-15.0/646.0*t17+3.0/1292.0*r-3.0/2584.0*t3+s/1615
455 -3.0/45220.0;
456       drv[0] = t18;
457       drv[1] = t22;
458    /* p=11 6 2 */
459    } else if( i==6 && j==2 ) {
460       t1 = r*r;
461       t2 = t1*t1;
462       t3 = t2*r;
463       t4 = s*s;
464       t6 = t3*s;
465       t8 = t2*s;
466       t9 = r*t1;
467       t11 = t9*s;
468       t13 = t1*s;
469       t15 = r*s;
470       t16 = 6.0*t3*t4-9.0/5.0*t6+9.0/95.0*t3-21.0/2.0*t2*t4+63.0/19.0*t8-7.0/
471 38.0*t2+126.0/19.0*t9*t4-42.0/19.0*t11+42.0/323.0*t9-35.0/19.0*t1*t4+210.0/
472 323.0*t13-105.0/2584.0*t1+70.0/323.0*r*t4-105.0/1292.0*t15+7.0/1292.0*r-21.0/
473 2584.0*t4+21.0/6460.0*s-3.0/12920.0;
474       t17 = t2*t1;
475       t19 = 2.0*t17*s-3.0/10.0*t17-21.0/5.0*t6+63.0/95.0*t3+63.0/19.0*t8-21.0/
476 38.0*t2-70.0/57.0*t11+70.0/323.0*t9+70.0/323.0*t13-105.0/2584.0*t1-21.0/1292.0*
477 t15+21.0/6460.0*r+7.0/19380.0*s-1.0/12920.0;
478       drv[0] = t16;
479       drv[1] = t19;
480    /* p=11 7 1 */
481    } else if( i==7 && j==1 ) {
482       t1 = r*r;
483       t2 = t1*t1;
484       t3 = t2*t1;
485       t5 = t2*r;
486       t8 = r*t1;
487       t12 = 7.0*t3*s-7.0/10.0*t3-84.0/5.0*t5*s+168.0/95.0*t5+294.0/19.0*t2*s
488 -98.0/57.0*t2-392.0/57.0*t8*s+784.0/969.0*t8+490.0/323.0*t1*s-245.0/1292.0*t1
489 -49.0/323.0*r*s+98.0/4845.0*r+49.0/9690.0*s-7.0/9690.0;
490       drv[0] = t12;
491       drv[1] = t2*t8-14.0/5.0*t3+294.0/95.0*t5-98.0/57.0*t2+490.0/969.0*t8-49.0
492 /646.0*t1+49.0/9690.0*r-1.0/9690.0;
493    /* p=11 8 0 */
494    } else if( i==8 && j==0 ) {
495       t1 = r*r;
496       t2 = r*t1;
497       t3 = t1*t1;
498       drv[0] = 8.0*t3*t2-126.0/5.0*t3*t1+3024.0/95.0*t3*r-392.0/19.0*t3+2352.0/
499 323.0*t2-441.0/323.0*t1+196.0/1615.0*r-6.0/1615.0;
500       drv[1] = 0.0;
501    /* p=12 0 9 */
502    } else if( i==0 && j==9 ) {
503       t1 = s*s;
504       t2 = t1*t1;
505       t3 = t2*t2;
506       t4 = t1*s;
507       drv[0] = 0.0;
508       drv[1] = 9.0*t3-360.0/11.0*t2*t4+540.0/11.0*t2*t1-432.0/11.0*t2*s+3780.0/
509 209.0*t2-1008.0/209.0*t4+2520.0/3553.0*t1-180.0/3553.0*s+9.0/7106.0;
510    /* p=12 1 8 */
511    } else if( i==1 && j==8 ) {
512       t1 = s*s;
513       t2 = t1*t1;
514       t3 = t2*t2;
515       t4 = t1*s;
516       t5 = t2*t4;
517       t6 = t2*t1;
518       t7 = t2*s;
519       t16 = 8.0*r*t5-252.0/11.0*r*t6+288.0/11.0*r*t7-168.0/11.0*t2*r+1008.0/
520 209.0*r*t4-168.0/209.0*r*t1+224.0/3553.0*r*s-6.0/3553.0*r-8.0/11.0*t5+24.0/11.0
521 *t6-144.0/55.0*t7+336.0/209.0*t2-112.0/209.0*t4+336.0/3553.0*t1-28.0/3553.0*s+
522 4.0/17765.0;
523       drv[0] = t3-36.0/11.0*t5+48.0/11.0*t6-168.0/55.0*t7+252.0/209.0*t2-56.0/
524 209.0*t4+112.0/3553.0*t1-6.0/3553.0*s+1.0/35530.0;
525       drv[1] = t16;
526    /* p=12 2 7 */
527    } else if( i==2 && j==7 ) {
528       t1 = s*s;
529       t2 = t1*s;
530       t3 = t1*t1;
531       t4 = t3*t2;
532       t6 = t3*t1;
533       t7 = r*t6;
534       t8 = t3*s;
535       t9 = r*t8;
536       t10 = t3*r;
537       t11 = t2*r;
538       t12 = r*t1;
539       t13 = r*s;
540       t14 = 2.0*r*t4-56.0/11.0*t7+56.0/11.0*t9-28.0/11.0*t10+140.0/209.0*t11
541 -56.0/627.0*t12+56.0/10659.0*t13-r/10659-3.0/11.0*t4+8.0/11.0*t6-42.0/55.0*t8+
542 84.0/209.0*t3-70.0/627.0*t2+56.0/3553.0*t1-7.0/7106.0*s+1.0/53295.0;
543       t15 = r*r;
544       t23 = 336.0/209.0*t11-70.0/209.0*t12+112.0/3553.0*t13-7.0/7106.0*r+t6/11
545 -12.0/55.0*t8+42.0/209.0*t3-56.0/627.0*t2+70.0/3553.0*t1-7.0/3553.0*s+7.0/
546 106590.0;
547       drv[0] = t14;
548       drv[1] = 7.0*t15*t6-168.0/11.0*t15*t8+140.0/11.0*t15*t3-56.0/11.0*t15*t2+
549 210.0/209.0*t15*t1-56.0/627.0*t15*s+28.0/10659.0*t15-21.0/11.0*t7+48.0/11.0*t9
550 -42.0/11.0*t10+t23;
551    /* p=12 3 6 */
552    } else if( i==3 && j==6 ) {
553       t1 = r*s;
554       t2 = s*s;
555       t3 = t2*r;
556       t4 = t2*t2;
557       t5 = r*t4;
558       t6 = r*r;
559       t7 = t4*s;
560       t8 = r*t7;
561       t9 = t2*s;
562       t10 = t6*t9;
563       t12 = t4*t2;
564       t13 = t6*s;
565       t14 = t6*t7;
566       t15 = r*t9;
567       t17 = t6*t2;
568       t19 = t6*t4;
569       t20 = 6.0/77.0*t12-3.0/209.0*t13-63.0/11.0*t14+120.0/209.0*t15+3.0*t6*t12
570 -10.0/209.0*t9+27.0/209.0*t4-9.0/55.0*t7+45.0/209.0*t17-12.0/11.0*r*t12+45.0/
571 11.0*t19;
572       t22 = 1.0/35530.0+12.0/3553.0*t6-20.0/209.0*t13-5.0/7106.0*s-6.0/385.0*t7
573 -36.0/11.0*t10-9.0/14212.0*r+20.0/3553.0*t2-4.0/209.0*t9+6.0/209.0*t4+60.0/
574 3553.0*t1+60.0/11.0*t19;
575       t23 = r*t6;
576       t29 = -15.0/11.0*t23*t2+36.0/77.0*t8-t23/209-105.0/11.0*t23*t4+180.0/
577 209.0*t17+6.0*t23*t7-30.0/209.0*t3+108.0/209.0*t15-9.0/11.0*t5+30.0/209.0*t23*s
578 -36.0/11.0*t14+60.0/11.0*t23*t9;
579       drv[0] = 1.0/71060.0+24.0/3553.0*t1-20.0/209.0*t3-18.0/11.0*t5-9.0/
580 14212.0*s-r/7106+30.0/3553.0*t2+t6/3553+24.0/11.0*t8-15.0/11.0*t10+t20;
581       drv[1] = t22+t29;
582    /* p=12 4 5 */
583    } else if( i==4 && j==5 ) {
584       t1 = r*s;
585       t2 = s*s;
586       t3 = t2*r;
587       t4 = r*r;
588       t5 = r*t4;
589       t6 = t2*t2;
590       t7 = t6*s;
591       t9 = t5*t6;
592       t10 = r*t6;
593       t12 = t2*s;
594       t13 = t4*t12;
595       t14 = 1.0/74613.0+200.0/24871.0*t1-400.0/4389.0*t3+4.0*t5*t7-60.0/11.0*t9
596 -60.0/77.0*t10-25.0/49742.0*s-5.0/24871.0*r+400.0/74613.0*t2+20.0/24871.0*t4+
597 40.0/77.0*r*t7-150.0/77.0*t13;
598       t15 = t5*s;
599       t16 = t4*s;
600       t18 = r*t12;
601       t19 = t5*t12;
602       t20 = t4*t2;
603       t21 = t6*t4;
604       t22 = t5*t2;
605       t23 = 60.0/1463.0*t15-50.0/1463.0*t16-30.0/11.0*t4*t7+600.0/1463.0*t18
606 -4.0/4389.0*t5-100.0/4389.0*t12+60.0/1463.0*t6-2.0/77.0*t7+200.0/77.0*t19+600.0
607 /1463.0*t20+300.0/77.0*t21-40.0/77.0*t22;
608       t25 = t4*t4;
609       t27 = 5.0/298452.0+800.0/74613.0*t1-50.0/11.0*t9-10.0/77.0*t10+5.0/1463.0
610 *t6-25.0/74613.0*s+50.0/24871.0*t2-50.0/4389.0*t5-100.0/1463.0*t3+15.0/1463.0*
611 t25+5.0*t25*t6-20.0/4389.0*t12;
612       t31 = -120.0/77.0*t13+400.0/1463.0*t15-400.0/4389.0*t16+240.0/1463.0*t18+
613 400.0/77.0*t19+900.0/1463.0*t20+100.0/77.0*t21-150.0/77.0*t22-20.0/77.0*t25*s
614 -60.0/11.0*t25*t12+150.0/77.0*t25*t2+100.0/24871.0*t4-25.0/49742.0*r;
615       drv[0] = t14+t23;
616       drv[1] = t27+t31;
617    /* p=12 5 4 */
618    } else if( i==5 && j==4 ) {
619       t1 = r*r;
620       t2 = t1*r;
621       t3 = t2*s;
622       t4 = s*s;
623       t5 = t4*s;
624       t6 = t5*t2;
625       t7 = t1*t1;
626       t8 = t4*t4;
627       t9 = t1*t4;
628       t10 = t1*s;
629       t11 = 240.0/1463.0*t3+400.0/77.0*t6-20.0/4389.0*t2-50.0/4389.0*t5+5.0/
630 1463.0*t7+15.0/1463.0*t8+100.0/24871.0*t4+900.0/1463.0*t9+50.0/24871.0*t1-100.0
631 /1463.0*t10-25.0/74613.0*r-25.0/49742.0*s;
632       t14 = r*s;
633       t15 = r*t4;
634       t16 = t7*s;
635       t17 = t1*t5;
636       t18 = r*t5;
637       t20 = t5*t7;
638       t21 = t4*t2;
639       t23 = t4*t7;
640       t24 = 5.0*t7*t8+150.0/77.0*t1*t8+800.0/74613.0*t14-400.0/4389.0*t15-10.0/
641 77.0*t16-150.0/77.0*t17+400.0/1463.0*t18-20.0/77.0*r*t8+5.0/298452.0-50.0/11.0*
642 t20-120.0/77.0*t21-60.0/11.0*t2*t8+100.0/77.0*t23;
643       t26 = r*t7;
644       t29 = 1.0/74613.0+400.0/74613.0*t1-400.0/4389.0*t10-60.0/11.0*t20-5.0/
645 24871.0*s-40.0/77.0*t17-25.0/49742.0*r+20.0/24871.0*t4-4.0/4389.0*t5+40.0/77.0*
646 t26*s+200.0/24871.0*t14+4.0*t26*t5;
647       t31 = -30.0/11.0*t26*t4-150.0/77.0*t21+60.0/1463.0*t7-100.0/4389.0*t2-2.0
648 /77.0*t26+600.0/1463.0*t9-50.0/1463.0*t15+60.0/1463.0*t18+300.0/77.0*t23+600.0/
649 1463.0*t3+200.0/77.0*t6-60.0/77.0*t16;
650       drv[0] = t11+t24;
651       drv[1] = t29+t31;
652    /* p=12 6 3 */
653    } else if( i==6 && j==3 ) {
654       t1 = r*s;
655       t2 = s*s;
656       t3 = t2*r;
657       t4 = r*r;
658       t5 = t4*t4;
659       t6 = r*t5;
660       t7 = t5*t2;
661       t8 = t6*t2;
662       t9 = t2*s;
663       t12 = 1.0/35530.0+60.0/3553.0*t1-20.0/209.0*t3-6.0/385.0*t6-9.0/14212.0*s
664 -5.0/7106.0*r+60.0/11.0*t7+12.0/3553.0*t2-36.0/11.0*t8-105.0/11.0*t5*t9+20.0/
665 3553.0*t4-15.0/11.0*t4*t9;
666       t13 = r*t4;
667       t14 = t13*s;
668       t15 = t4*s;
669       t16 = t5*s;
670       t20 = t6*s;
671       t21 = t4*t2;
672       t22 = t13*t2;
673       t23 = 108.0/209.0*t14-30.0/209.0*t15-9.0/11.0*t16+6.0*t6*t9+30.0/209.0*r*
674 t9-4.0/209.0*t13-t9/209+60.0/11.0*t13*t9+36.0/77.0*t20+180.0/209.0*t21-36.0/
675 11.0*t22+6.0/209.0*t5;
676       t25 = t5*t4;
677       t29 = -15.0/11.0*t22+27.0/209.0*t5-10.0/209.0*t13-9.0/55.0*t6+3.0*t25*t2+
678 45.0/209.0*t21-3.0/209.0*t3+45.0/11.0*t7+6.0/77.0*t25+120.0/209.0*t14-18.0/11.0
679 *t16;
680       drv[0] = t12+t23;
681       drv[1] = 1.0/71060.0+30.0/3553.0*t4-20.0/209.0*t15-s/7106-9.0/14212.0*r+
682 t2/3553+24.0/11.0*t20+24.0/3553.0*t1-12.0/11.0*t25*s-63.0/11.0*t8+t29;
683    /* p=12 7 2 */
684    } else if( i==7 && j==2 ) {
685       t1 = r*r;
686       t2 = t1*t1;
687       t3 = t2*t1;
688       t4 = s*s;
689       t6 = t3*s;
690       t7 = t2*r;
691       t9 = t7*s;
692       t11 = t2*s;
693       t12 = t1*r;
694       t15 = t12*s;
695       t17 = t1*s;
696       t19 = r*s;
697       t20 = 336.0/209.0*t15-56.0/627.0*t12+210.0/209.0*t1*t4-70.0/209.0*t17+
698 70.0/3553.0*t1-56.0/627.0*r*t4+112.0/3553.0*t19-7.0/3553.0*r+28.0/10659.0*t4
699 -7.0/7106.0*s+7.0/106590.0;
700       t22 = t2*t12;
701       t24 = 2.0*t22*s-3.0/11.0*t22-56.0/11.0*t6+8.0/11.0*t3+56.0/11.0*t9-42.0/
702 55.0*t7-28.0/11.0*t11+84.0/209.0*t2+140.0/209.0*t15-70.0/627.0*t12-56.0/627.0*
703 t17+56.0/3553.0*t1+56.0/10659.0*t19-7.0/7106.0*r-s/10659+1.0/53295.0;
704       drv[0] = 7.0*t3*t4-21.0/11.0*t6+t3/11-168.0/11.0*t7*t4+48.0/11.0*t9-12.0/
705 55.0*t7+140.0/11.0*t4*t2-42.0/11.0*t11+42.0/209.0*t2-56.0/11.0*t12*t4+t20;
706       drv[1] = t24;
707    /* p=12 8 1 */
708    } else if( i==8 && j==1 ) {
709       t1 = r*r;
710       t2 = t1*r;
711       t3 = t1*t1;
712       t4 = t3*t2;
713       t6 = t3*t1;
714       t8 = t3*r;
715       t14 = 8.0*t4*s-8.0/11.0*t4-252.0/11.0*t6*s+24.0/11.0*t6+288.0/11.0*t8*s
716 -144.0/55.0*t8-168.0/11.0*t3*s+336.0/209.0*t3+1008.0/209.0*t2*s-112.0/209.0*t2
717 -168.0/209.0*t1*s+336.0/3553.0*t1+224.0/3553.0*r*s-28.0/3553.0*r-6.0/3553.0*s+
718 4.0/17765.0;
719       t15 = t3*t3;
720       drv[0] = t14;
721       drv[1] = t15-36.0/11.0*t4+48.0/11.0*t6-168.0/55.0*t8+252.0/209.0*t3-56.0/
722 209.0*t2+112.0/3553.0*t1-6.0/3553.0*r+1.0/35530.0;
723    /* p=12 9 0 */
724    } else if( i==9 && j==0 ) {
725       t1 = r*r;
726       t2 = t1*t1;
727       t3 = t2*t2;
728       t4 = t1*r;
729       drv[0] = 9.0*t3-360.0/11.0*t4*t2+540.0/11.0*t2*t1-432.0/11.0*t2*r+3780.0/
730 209.0*t2-1008.0/209.0*t4+2520.0/3553.0*t1-180.0/3553.0*r+9.0/7106.0;
731       drv[1] = 0.0;
732    /* p=13 0 10 */
733    } else if( i==0 && j==10 ) {
734       t1 = s*s;
735       t2 = t1*t1;
736       t3 = t2*t2;
737       t5 = t1*s;
738       drv[0] = 0.0;
739       drv[1] = 10.0*t3*s-165.0/4.0*t3+1650.0/23.0*t2*t5-1575.0/23.0*t2*t1+900.0
740 /23.0*t2*s-315.0/23.0*t2+1260.0/437.0*t5-150.0/437.0*t1+150.0/7429.0*s-25.0/
741 59432.0;
742    /* p=13 1 9 */
743    } else if( i==1 && j==9 ) {
744       t1 = s*s;
745       t2 = t1*t1;
746       t3 = t2*t2;
747       t5 = t1*s;
748       t6 = t2*t5;
749       t7 = t2*t1;
750       t8 = t2*s;
751       t18 = 9.0*r*t3-30.0*r*t6+945.0/23.0*r*t7-7560.0/253.0*r*t8+3150.0/253.0*r
752 *t2-756.0/253.0*r*t5+1890.0/4807.0*r*t1-120.0/4807.0*r*s+45.0/81719.0*r-3.0/4.0
753 *t3+60.0/23.0*t6-945.0/253.0*t7+720.0/253.0*t8-315.0/253.0*t2+1512.0/4807.0*t5
754 -210.0/4807.0*t1+240.0/81719.0*s-45.0/653752.0;
755       drv[0] = t3*s-15.0/4.0*t3+135.0/23.0*t6-1260.0/253.0*t7+630.0/253.0*t8
756 -189.0/253.0*t2+630.0/4807.0*t5-60.0/4807.0*t1+45.0/81719.0*s-5.0/653752.0;
757       drv[1] = t18;
758    /* p=13 2 8 */
759    } else if( i==2 && j==8 ) {
760       t1 = s*s;
761       t2 = t1*t1;
762       t3 = t2*s;
763       t4 = r*t2;
764       t5 = r*t3;
765       t6 = t2*t1;
766       t7 = r*t6;
767       t8 = t2*t2;
768       t10 = t1*s;
769       t11 = t2*t10;
770       t12 = r*t11;
771       t13 = r*s;
772       t14 = r*t1;
773       t15 = r*t10;
774       t16 = -3.0/653752.0-28.0/4807.0*t1+24.0/81719.0*s+2.0/81719.0*r+168.0/
775 253.0*t3+420.0/253.0*t4-1176.0/253.0*t5+168.0/23.0*t7+2.0*r*t8-6.0*t12+252.0/
776 4807.0*t10-8.0/4807.0*t13-252.0/253.0*t6+168.0/4807.0*t14-84.0/253.0*t15-63.0/
777 253.0*t2+18.0/23.0*t11-t8/4;
778       t17 = r*r;
779       t23 = 56.0/81719.0*s-3.0/163438.0-56.0/4807.0*t13+24.0/81719.0*r+252.0/
780 4807.0*t10-42.0/4807.0*t1+840.0/253.0*t4+840.0/253.0*t17*t10-2940.0/253.0*t17*
781 t2+504.0/23.0*t17*t3+8.0*t17*t11-21.0*t17*t6;
782       t26 = 756.0/4807.0*t14-4.0/4807.0*t17-1512.0/253.0*t5+168.0/4807.0*t17*s
783 -126.0/253.0*t17*t1-252.0/253.0*t15+126.0/23.0*t7+2.0/23.0*t11-63.0/253.0*t6
784 -2.0*t12+72.0/253.0*t3-42.0/253.0*t2;
785       drv[0] = t16;
786       drv[1] = t23+t26;
787    /* p=13 3 7 */
788    } else if( i==3 && j==7 ) {
789       t1 = s*s;
790       t2 = r*r;
791       t3 = t1*t1;
792       t4 = t3*t1;
793       t5 = t2*t4;
794       t6 = t3*s;
795       t7 = r*t3;
796       t8 = t2*t3;
797       t9 = r*t6;
798       t10 = r*t4;
799       t11 = t1*s;
800       t12 = t3*t11;
801       t14 = -1.0/326876.0-14.0/4807.0*t1-7.0*t5+14.0/81719.0*s+8.0/245157.0*r+
802 42.0/253.0*t6-t2/14421+280.0/253.0*t7-735.0/253.0*t8-588.0/253.0*t9+56.0/23.0*
803 t10-r*t12;
804       t15 = r*s;
805       t16 = r*t1;
806       t17 = r*t11;
807       t18 = t2*t11;
808       t19 = t2*t1;
809       t20 = t2*s;
810       t22 = t2*t6;
811       t23 = 105.0/4807.0*t11-28.0/14421.0*t15-42.0/253.0*t4+168.0/4807.0*t16
812 -70.0/253.0*t17-21.0/253.0*t3+175.0/253.0*t18-21.0/253.0*t19+3.0/46.0*t12+21.0/
813 4807.0*t20+3.0*t2*t12+147.0/23.0*t22;
814       t25 = t2*r;
815       t30 = -7.0/506.0*t4-14.0/14421.0*t2+14.0/81719.0*r+56.0/245157.0*s-7.0/
816 2.0*t5-35.0/14421.0*t1+210.0/253.0*t7-1470.0/253.0*t8-252.0/253.0*t9+21.0/46.0*
817 t10-14.0*t25*t6-980.0/253.0*t25*t11+175.0/253.0*t25*t1-14.0/253.0*t25*s;
818       t33 = 7.0/4807.0*t25+168.0/4807.0*t20+168.0/23.0*t22+7.0*t25*t4+245.0/
819 23.0*t25*t3+56.0/4807.0*t11-28.0/4807.0*t15+315.0/4807.0*t16-84.0/253.0*t17+
820 560.0/253.0*t18-105.0/253.0*t19-7.0/253.0*t3+8.0/253.0*t6-7.0/980628.0;
821       drv[0] = t14+t23;
822       drv[1] = t30+t33;
823    /* p=13 4 6 */
824    } else if( i==4 && j==6 ) {
825       t1 = r*r;
826       t2 = s*s;
827       t3 = t2*r;
828       t4 = r*t1;
829       t5 = t2*t2;
830       t6 = t5*t2;
831       t8 = r*s;
832       t9 = t1*t5;
833       t10 = t4*t5;
834       t11 = t5*s;
835       t12 = t1*t11;
836       t14 = t2*s;
837       t15 = r*t14;
838       t16 = t1*s;
839       t18 = t2*t1;
840       t19 = t1*t14;
841       t20 = r*t5;
842       t21 = -5.0/28842.0*t1+150.0/4807.0*t3+4.0*t4*t6-10.0/4807.0*t8-1575.0/
843 506.0*t9+105.0/23.0*t10+105.0/23.0*t12-5.0/2.0*t1*t6-50.0/253.0*t15+45.0/4807.0
844 *t16+10.0/23.0*r*t6-75.0/506.0*t18+250.0/253.0*t19+150.0/253.0*t20;
845       t22 = t4*t14;
846       t23 = t2*t4;
847       t24 = t4*s;
848       t25 = t4*t11;
849       t26 = r*t11;
850       t27 = -350.0/253.0*t22+50.0/253.0*t23-3.0/253.0*t24-7.0*t25+10.0/245157.0
851 *r+10.0/253.0*t11-210.0/253.0*t26-5.0/253.0*t6-15.0/506.0*t5+t4/4807-25.0/
852 14421.0*t2+10.0/81719.0*s+50.0/4807.0*t14-5.0/1961256.0;
853       t29 = t1*t1;
854       t31 = -5.0/4807.0*t1+25.0/253.0*t29*s+150.0/4807.0*t3-50.0/14421.0*t8
855 -525.0/253.0*t9+175.0/23.0*t10+30.0/23.0*t12-30.0/253.0*t15+150.0/4807.0*t16
856 -75.0/253.0*t18+300.0/253.0*t19+50.0/253.0*t20-1050.0/253.0*t22+250.0/253.0*t23
857 -25.0/253.0*t24;
858       t36 = -5.0*t25-525.0/506.0*t29*t2+10.0/81719.0*r+105.0/23.0*t29*t14+6.0*
859 t29*t11-35.0/4.0*t29*t5-3.0/1012.0*t29+5.0/1771.0*t11-30.0/253.0*t26-5.0/1012.0
860 *t5+15.0/4807.0*t4-25.0/28842.0*t2+25.0/245157.0*s+15.0/4807.0*t14-5.0/
861 1307504.0;
862       drv[0] = t21+t27;
863       drv[1] = t31+t36;
864    /* p=13 5 5 */
865    } else if( i==5 && j==5 ) {
866       t1 = r*r;
867       t2 = r*t1;
868       t3 = s*s;
869       t4 = t2*t3;
870       t5 = t1*t1;
871       t6 = t3*t3;
872       t7 = t6*s;
873       t8 = t3*s;
874       t9 = r*t3;
875       t10 = r*s;
876       t12 = t1*t6;
877       t13 = t5*t3;
878       t14 = t2*t6;
879       t15 = t5*s;
880       t16 = -25.0/9152528.0+1000.0/1771.0*t4-5.0/7084.0*t5+30.0/33649.0*t2+25.0
881 /3542.0*t7+375.0/67298.0*t8-25.0/67298.0*t1-75.0/7084.0*t6+1000.0/33649.0*t9
882 -250.0/100947.0*t10-50.0/253.0*r*t7-1125.0/506.0*t12-125.0/253.0*t13+150.0/23.0
883 *t14+125.0/3542.0*t15;
884       t18 = r*t8;
885       t19 = t5*t8;
886       t20 = t1*s;
887       t21 = t1*t3;
888       t22 = t1*t8;
889       t24 = t5*t6;
890       t25 = r*t6;
891       t26 = t2*t8;
892       t27 = t2*s;
893       t29 = 75.0/46.0*t1*t7-250.0/1771.0*t18+125.0/46.0*t19+1125.0/67298.0*t20
894 -375.0/1771.0*t21+1875.0/1771.0*t22+5.0*t5*t7-25.0/4.0*t24+500.0/1771.0*t25
895 -750.0/253.0*t26-75.0/1771.0*t27-5.0*t2*t7+100.0/1716099.0*r-125.0/100947.0*t3+
896 125.0/1144066.0*s;
897       t31 = -25.0/9152528.0+1875.0/1771.0*t4-75.0/7084.0*t5+375.0/67298.0*t2+
898 30.0/33649.0*t8-125.0/100947.0*t1-5.0/7084.0*t6+1125.0/67298.0*t9-250.0/
899 100947.0*t10-125.0/253.0*t12-1125.0/506.0*t13+125.0/46.0*t14+500.0/1771.0*t15
900 -75.0/1771.0*t18+150.0/23.0*t19;
901       t32 = r*t5;
902       t37 = 1000.0/33649.0*t20-375.0/1771.0*t21+1000.0/1771.0*t22-25.0/4.0*t24+
903 125.0/3542.0*t25-750.0/253.0*t26-250.0/1771.0*t27+125.0/1144066.0*r-25.0/
904 67298.0*t3+100.0/1716099.0*s-50.0/253.0*t32*s+75.0/46.0*t32*t3-5.0*t32*t8+25.0/
905 3542.0*t32+5.0*t32*t6;
906       drv[0] = t16+t29;
907       drv[1] = t31+t37;
908    /* p=13 6 4 */
909    } else if( i==6 && j==4 ) {
910       t1 = s*s;
911       t2 = r*r;
912       t3 = t2*t2;
913       t4 = t3*r;
914       t5 = t1*t1;
915       t7 = r*t1;
916       t8 = r*s;
917       t9 = t2*r;
918       t10 = t1*s;
919       t11 = t3*t10;
920       t12 = r*t10;
921       t13 = t3*t1;
922       t14 = 25.0/245157.0*r+10.0/81719.0*s-5.0/4807.0*t1+6.0*t4*t5+150.0/4807.0
923 *t7-25.0/28842.0*t2-50.0/14421.0*t8-3.0/1012.0*t5+5.0/1771.0*t4+15.0/4807.0*t9+
924 15.0/4807.0*t10+175.0/23.0*t11-25.0/253.0*t12-5.0/1012.0*t3-525.0/253.0*t13;
925       t16 = t9*s;
926       t17 = t9*t1;
927       t18 = t9*t10;
928       t20 = t3*s;
929       t21 = t4*s;
930       t22 = t2*s;
931       t23 = t2*t1;
932       t24 = t2*t10;
933       t27 = t4*t1;
934       t28 = t4*t10;
935       t29 = -525.0/506.0*t2*t5-30.0/253.0*t16+300.0/253.0*t17-1050.0/253.0*t18+
936 105.0/23.0*t9*t5+50.0/253.0*t20-30.0/253.0*t21+150.0/4807.0*t22-75.0/253.0*t23+
937 250.0/253.0*t24-35.0/4.0*t3*t5+25.0/253.0*r*t5+30.0/23.0*t27-5.0*t28-5.0/
938 1307504.0;
939       t31 = 10.0/81719.0*r+10.0/245157.0*s-5.0/28842.0*t1+45.0/4807.0*t7-25.0/
940 14421.0*t2-10.0/4807.0*t8+10.0/253.0*t4+50.0/4807.0*t9+t10/4807+105.0/23.0*t11
941 -3.0/253.0*t12-15.0/506.0*t3-1575.0/506.0*t13-50.0/253.0*t16;
942       t32 = t3*t2;
943       t36 = 250.0/253.0*t17-350.0/253.0*t18-5.0/253.0*t32+4.0*t32*t10-5.0/2.0*
944 t32*t1+150.0/253.0*t20-210.0/253.0*t21+150.0/4807.0*t22-75.0/506.0*t23+50.0/
945 253.0*t24+105.0/23.0*t27-7.0*t28+10.0/23.0*t32*s-5.0/1961256.0;
946       drv[0] = t14+t29;
947       drv[1] = t31+t36;
948    /* p=13 7 3 */
949    } else if( i==7 && j==3 ) {
950       t1 = s*s;
951       t2 = r*t1;
952       t3 = r*s;
953       t4 = r*r;
954       t5 = t4*t4;
955       t6 = t4*t5;
956       t7 = t6*s;
957       t8 = r*t5;
958       t9 = t8*s;
959       t10 = t1*s;
960       t13 = t4*s;
961       t14 = t8*t1;
962       t17 = t4*r;
963       t19 = t17*t1;
964       t20 = t17*s;
965       t21 = 168.0/4807.0*t2-28.0/4807.0*t3+21.0/46.0*t7-252.0/253.0*t9-14.0/
966 253.0*r*t10+245.0/23.0*t5*t10+315.0/4807.0*t13+168.0/23.0*t14-14.0*t8*t10+7.0*
967 t6*t10-980.0/253.0*t17*t10+560.0/253.0*t19-84.0/253.0*t20+14.0/81719.0*s;
968       t23 = t4*t1;
969       t24 = t5*t1;
970       t25 = t5*s;
971       t26 = t6*t1;
972       t27 = 56.0/245157.0*r-14.0/14421.0*t1+175.0/253.0*t4*t10-105.0/253.0*t23
973 -1470.0/253.0*t24+8.0/253.0*t8+210.0/253.0*t25-7.0/253.0*t5+7.0/4807.0*t10+56.0
974 /4807.0*t17-7.0/506.0*t6-35.0/14421.0*t4-7.0/980628.0-7.0/2.0*t26;
975       t29 = t5*t17;
976       t31 = 8.0/245157.0*s-1.0/326876.0-28.0/14421.0*t3+14.0/81719.0*r-t1/14421
977 -42.0/253.0*t6-t29*s+280.0/253.0*t25-21.0/253.0*t5+42.0/253.0*t8+3.0/46.0*t29+
978 21.0/4807.0*t2;
979       t33 = -735.0/253.0*t24-14.0/4807.0*t4+168.0/4807.0*t13+3.0*t29*t1-21.0/
980 253.0*t23-588.0/253.0*t9+147.0/23.0*t14-7.0*t26+175.0/253.0*t19-70.0/253.0*t20+
981 56.0/23.0*t7+105.0/4807.0*t17;
982       drv[0] = t21+t27;
983       drv[1] = t31+t33;
984    /* p=13 8 2 */
985    } else if( i==8 && j==2 ) {
986       t1 = s*s;
987       t2 = r*r;
988       t3 = t2*t2;
989       t4 = t3*r;
990       t5 = t2*r;
991       t6 = t3*t5;
992       t8 = t6*s;
993       t9 = t4*s;
994       t10 = t3*t2;
995       t13 = t5*s;
996       t14 = -3.0/163438.0-4.0/4807.0*t1+72.0/253.0*t4+24.0/81719.0*s+8.0*t6*t1+
997 56.0/81719.0*r-2.0*t8-1512.0/253.0*t9-42.0/4807.0*t2-21.0*t10*t1+840.0/253.0*t5
998 *t1-252.0/253.0*t13;
999       t15 = t3*s;
1000       t16 = t10*s;
1001       t17 = r*s;
1002       t22 = t2*s;
1003       t23 = 840.0/253.0*t15+252.0/4807.0*t5+126.0/23.0*t16-56.0/4807.0*t17+
1004 168.0/4807.0*r*t1-2940.0/253.0*t3*t1-63.0/253.0*t10-42.0/253.0*t3-126.0/253.0*
1005 t2*t1+2.0/23.0*t6+504.0/23.0*t4*t1+756.0/4807.0*t22;
1006       t25 = t3*t3;
1007       t27 = 2.0/81719.0*s-3.0/653752.0-8.0/4807.0*t17+24.0/81719.0*r+2.0*t25*s-
1008 t25/4-252.0/253.0*t10-6.0*t8+420.0/253.0*t15-63.0/253.0*t3+168.0/253.0*t4+18.0/
1009 23.0*t6-28.0/4807.0*t2+168.0/4807.0*t22-1176.0/253.0*t9-84.0/253.0*t13+168.0/
1010 23.0*t16+252.0/4807.0*t5;
1011       drv[0] = t14+t23;
1012       drv[1] = t27;
1013    /* p=13 9 1 */
1014    } else if( i==9 && j==1 ) {
1015       t1 = r*r;
1016       t2 = t1*t1;
1017       t3 = t2*t2;
1018       t5 = r*t1;
1019       t6 = t2*t5;
1020       t8 = t2*t1;
1021       t10 = t2*r;
1022       t16 = 9.0*t3*s-3.0/4.0*t3-30.0*t6*s+60.0/23.0*t6+945.0/23.0*t8*s-945.0/
1023 253.0*t8-7560.0/253.0*t10*s+720.0/253.0*t10+3150.0/253.0*t2*s-315.0/253.0*t2
1024 -756.0/253.0*t5*s+1512.0/4807.0*t5+1890.0/4807.0*t1*s-210.0/4807.0*t1-120.0/
1025 4807.0*r*s+240.0/81719.0*r+45.0/81719.0*s-45.0/653752.0;
1026       drv[0] = t16;
1027       drv[1] = t3*r-15.0/4.0*t3+135.0/23.0*t6-1260.0/253.0*t8+630.0/253.0*t10
1028 -189.0/253.0*t2+630.0/4807.0*t5-60.0/4807.0*t1+45.0/81719.0*r-5.0/653752.0;
1029    /* p=13 10 0 */
1030    } else if( i==10 && j==0 ) {
1031       t1 = r*r;
1032       t2 = t1*t1;
1033       t3 = t2*t2;
1034       t5 = t1*r;
1035       drv[0] = 10.0*t3*r-165.0/4.0*t3+1650.0/23.0*t2*t5-1575.0/23.0*t2*t1+900.0
1036 /23.0*t2*r-315.0/23.0*t2+1260.0/437.0*t5-150.0/437.0*t1+150.0/7429.0*r-25.0/
1037 59432.0;
1038       drv[1] = 0.0;
1039    /* p=14 0 11 */
1040    } else if( i==0 && j==11 ) {
1041       t1 = s*s;
1042       t2 = t1*t1;
1043       t3 = t2*t2;
1044       t6 = t1*s;
1045       t10 = 11.0*t3*t1-660.0/13.0*t3*s+6534.0/65.0*t3-1452.0/13.0*t2*t6+22869.0
1046 /299.0*t2*t1-49896.0/1495.0*t2*s+2772.0/299.0*t2-2376.0/1495.0*t6+891.0/5681.0*
1047 t1-44.0/5681.0*s+66.0/482885.0;
1048       drv[0] = 0.0;
1049       drv[1] = t10;
1050    /* p=14 1 10 */
1051    } else if( i==1 && j==10 ) {
1052       t1 = s*s;
1053       t2 = t1*t1;
1054       t3 = t2*t2;
1055       t5 = t3*s;
1056       t6 = t1*s;
1057       t7 = t2*t6;
1058       t8 = t2*t1;
1059       t9 = t2*s;
1060       t10 = t3*t1-55.0/13.0*t5+99.0/13.0*t3-99.0/13.0*t7+1386.0/299.0*t8-2646.0
1061 /1495.0*t9+126.0/299.0*t2-18.0/299.0*t6+27.0/5681.0*t1-s/5681+1.0/482885.0;
1062       t20 = 10.0*r*t5-495.0/13.0*t3*r+792.0/13.0*r*t7-693.0/13.0*r*t8+8316.0/
1063 299.0*r*t9-2646.0/299.0*t2*r+504.0/299.0*r*t6-54.0/299.0*t1*r+54.0/5681.0*r*s-r
1064 /5681-10.0/13.0*t5+198.0/65.0*t3-66.0/13.0*t7+1386.0/299.0*t8-756.0/299.0*t9+
1065 252.0/299.0*t2-252.0/1495.0*t6+108.0/5681.0*t1-6.0/5681.0*s+2.0/96577.0;
1066       drv[0] = t10;
1067       drv[1] = t20;
1068    /* p=14 2 9 */
1069    } else if( i==2 && j==9 ) {
1070       t1 = s*s;
1071       t2 = t1*t1;
1072       t3 = t2*t2;
1073       t4 = t2*s;
1074       t5 = t2*t1;
1075       t6 = t3*s;
1076       t7 = t1*s;
1077       t8 = r*t5;
1078       t9 = t2*t7;
1079       t10 = r*t9;
1080       t11 = r*s;
1081       t12 = t1*r;
1082       t13 = r*t7;
1083       t14 = t2*r;
1084       t15 = r*t4;
1085       t17 = t3*r;
1086       t18 = 6.0/5311735.0+54.0/65.0*t3-7938.0/16445.0*t4+1512.0/1495.0*t5-3.0/
1087 13.0*t6-27.0/312455.0*s+648.0/312455.0*t1-378.0/16445.0*t7+2268.0/16445.0*t2
1088 -2.0/312455.0*r-504.0/65.0*t8+648.0/65.0*t10+162.0/312455.0*t11-216.0/16445.0*
1089 t12+504.0/3289.0*t13-15876.0/16445.0*t14+5292.0/1495.0*t15+2.0*r*t6-81.0/65.0*
1090 t9-90.0/13.0*t17;
1091       t19 = r*r;
1092       t25 = 2646.0/299.0*t19*t2-27.0/13.0*t17+27.0/5311735.0+432.0/65.0*t10+
1093 1296.0/312455.0*t11-1134.0/16445.0*t12+9072.0/16445.0*t13-7938.0/3289.0*t14+
1094 9072.0/1495.0*t15-31752.0/16445.0*t19*t7-360.0/13.0*t19*t9-216.0/16445.0*t19*s+
1095 756.0/3289.0*t19*t1;
1096       t29 = -1512.0/65.0*t19*t4+9.0*t19*t3+2268.0/65.0*t19*t5+81.0/312455.0*t19
1097 -18.0/65.0*t9-567.0/65.0*t8+1134.0/312455.0*t1-2268.0/82225.0*t7+378.0/3289.0*
1098 t2-27.0/312455.0*r+27.0/325.0*t3-72.0/312455.0*s-4536.0/16445.0*t4+567.0/1495.0
1099 *t5;
1100       drv[0] = t18;
1101       drv[1] = t25+t29;
1102    /* p=14 3 8 */
1103    } else if( i==3 && j==8 ) {
1104       t1 = r*r;
1105       t2 = s*s;
1106       t3 = t2*t2;
1107       t4 = t3*t3;
1108       t6 = t1*t3;
1109       t7 = t1*s;
1110       t8 = t2*s;
1111       t9 = t3*t8;
1112       t10 = t3*t2;
1113       t11 = r*t10;
1114       t12 = t3*s;
1115       t13 = r*t12;
1116       t14 = 18.0/26558675.0+3.0*t1*t4+2646.0/1495.0*t6-108.0/82225.0*t7+756.0/
1117 16445.0*t3-54.0/325.0*t9+27.0/1562275.0*t1-756.0/82225.0*t8+1512.0/1562275.0*t2
1118 +1512.0/7475.0*t10-1008.0/325.0*t11-10584.0/82225.0*t12+14112.0/7475.0*t13;
1119       t16 = r*s;
1120       t17 = t2*r;
1121       t18 = r*t8;
1122       t19 = t3*r;
1123       t20 = t1*t10;
1124       t21 = t1*t9;
1125       t22 = t2*t1;
1126       t23 = t1*t8;
1127       t24 = r*t9;
1128       t25 = t1*t12;
1129       t26 = 18.0/325.0*t4-72.0/1562275.0*s-12.0/1562275.0*r-12.0/13.0*r*t4+
1130 864.0/1562275.0*t16-1008.0/82225.0*t17+2016.0/16445.0*t18-10584.0/16445.0*t19+
1131 3024.0/325.0*t20-108.0/13.0*t21+504.0/16445.0*t22-5292.0/16445.0*t23+864.0/
1132 325.0*t24-1764.0/325.0*t25;
1133       t28 = t1*r;
1134       t31 = 48.0/26558675.0+6048.0/325.0*t28*t12-252.0/13.0*t28*t10+7056.0/
1135 1495.0*t6-1008.0/82225.0*t7+336.0/16445.0*t3-4.0/325.0*t9+432.0/1562275.0*t1
1136 -504.0/82225.0*t8+1512.0/1562275.0*t2+252.0/7475.0*t10-378.0/325.0*t11-3024.0/
1137 82225.0*t12+9072.0/7475.0*t13-112.0/1562275.0*s-72.0/1562275.0*r;
1138       t37 = 3024.0/1562275.0*t16-2268.0/82225.0*t17+3024.0/16445.0*t18-10584.0/
1139 16445.0*t19+3024.0/325.0*t20-48.0/13.0*t21+3024.0/16445.0*t22-21168.0/16445.0*
1140 t23+144.0/325.0*t24-3024.0/325.0*t25+8.0*t28*t9+3528.0/1495.0*t28*t8+336.0/
1141 16445.0*t28*s-5292.0/16445.0*t28*t2-36.0/82225.0*t28-588.0/65.0*t28*t3;
1142       drv[0] = t14+t26;
1143       drv[1] = t31+t37;
1144    /* p=14 4 7 */
1145    } else if( i==4 && j==7 ) {
1146       t1 = r*s;
1147       t2 = s*s;
1148       t3 = t2*t2;
1149       t4 = t3*t2;
1150       t5 = r*t4;
1151       t6 = r*r;
1152       t7 = t6*r;
1153       t8 = t2*s;
1154       t9 = t7*t8;
1155       t10 = t3*s;
1156       t11 = r*t10;
1157       t12 = t3*t8;
1158       t13 = t7*t3;
1159       t14 = t7*t10;
1160       t15 = t3*r;
1161       t16 = r*t8;
1162       t17 = 168.0/312455.0*t1-56.0/65.0*t5+196.0/299.0*t9+56.0/1495.0*t4+1176.0
1163 /1495.0*t11-588.0/16445.0*t10-28.0/937365.0*s-8.0/937365.0*r-t12/65-196.0/65.0*
1164 t13-4.0/82225.0*t7+2352.0/325.0*t14+12.0/312455.0*t6-1176.0/3289.0*t15+56.0/
1165 3289.0*t3+280.0/3289.0*t16;
1166       t18 = t2*r;
1167       t19 = t6*t4;
1168       t21 = t7*t4;
1169       t22 = t6*t10;
1170       t23 = t6*s;
1171       t25 = t6*t2;
1172       t26 = t7*t2;
1173       t28 = t6*t3;
1174       t29 = t6*t8;
1175       t30 = t7*s;
1176       t31 = -14.0/3289.0*t8-168.0/16445.0*t18+336.0/65.0*t19+4.0*t7*t12-112.0/
1177 13.0*t21+168.0/312455.0*t2-294.0/65.0*t22-42.0/16445.0*t23-30.0/13.0*t6*t12+
1178 168.0/3289.0*t25-1176.0/16445.0*t26+24.0/65.0*r*t12+588.0/299.0*t28-1470.0/
1179 3289.0*t29+56.0/16445.0*t30+8.0/15935205.0;
1180       t33 = 336.0/312455.0*t1-7.0/65.0*t5+784.0/299.0*t9+7.0/2990.0*t4+336.0/
1181 1495.0*t11-84.0/16445.0*t10-28.0/937365.0*s-28.0/937365.0*r-98.0/13.0*t13-14.0/
1182 16445.0*t7+672.0/65.0*t14+84.0/312455.0*t6-588.0/3289.0*t15+14.0/3289.0*t3+
1183 224.0/3289.0*t16-28.0/16445.0*t8-42.0/3289.0*t18;
1184       t34 = t6*t6;
1185       t41 = 84.0/65.0*t19-70.0/13.0*t21+21.0/62491.0*t2-168.0/65.0*t22-168.0/
1186 16445.0*t23+420.0/3289.0*t25-1470.0/3289.0*t26+588.0/299.0*t28+14.0/16445.0*t34
1187 +7.0*t34*t4-196.0/65.0*t34*t8-588.0/16445.0*t34*s+147.0/299.0*t34*t2+588.0/65.0
1188 *t34*t3-168.0/13.0*t34*t10-2352.0/3289.0*t29+112.0/3289.0*t30+14.0/15935205.0;
1189       drv[0] = t17+t31;
1190       drv[1] = t33+t41;
1191    /* p=14 5 6 */
1192    } else if( i==5 && j==6 ) {
1193       t1 = r*r;
1194       t2 = t1*t1;
1195       t3 = s*s;
1196       t4 = t3*t3;
1197       t5 = t4*t3;
1198       t7 = t1*r;
1199       t9 = t3*s;
1200       t10 = t7*t9;
1201       t11 = t3*t2;
1202       t12 = t4*s;
1203       t13 = r*t12;
1204       t15 = r*s;
1205       t16 = t3*r;
1206       t17 = r*t9;
1207       t18 = r*t4;
1208       t19 = t7*t4;
1209       t20 = t2*t12;
1210       t21 = t7*t12;
1211       t23 = t7*s;
1212       t24 = t7*t3;
1213       t25 = 5.0*t2*t5-60.0/13.0*t7*t5+420.0/299.0*t10+105.0/598.0*t11+84.0/
1214 299.0*t13+1.0/2124694.0-2.0/13.0*r*t5+36.0/62491.0*t15-30.0/3289.0*t16+200.0/
1215 3289.0*t17-630.0/3289.0*t18-63.0/13.0*t19-105.0/13.0*t20+504.0/65.0*t21+18.0/
1216 13.0*t1*t5+36.0/3289.0*t23-630.0/3289.0*t24;
1217       t26 = t1*s;
1218       t27 = t1*t3;
1219       t28 = t1*t9;
1220       t29 = t1*t4;
1221       t30 = t1*t12;
1222       t31 = t2*t9;
1223       t32 = t2*s;
1224       t33 = t2*t4;
1225       t34 = -27.0/6578.0*t26+225.0/3289.0*t27-1575.0/3289.0*t28+945.0/598.0*t29
1226 -63.0/26.0*t30-35.0/26.0*t31-63.0/6578.0*t32+63.0/13.0*t33+t2/6578-3.0/16445.0*
1227 t7-63.0/6578.0*t12+3.0/598.0*t5+45.0/124982.0*t3-15.0/6578.0*t9+45.0/6578.0*t4+
1228 9.0/124982.0*t1-3.0/124982.0*s-2.0/187473.0*r;
1229       t36 = 630.0/299.0*t10+315.0/299.0*t11+9.0/299.0*t13+3.0/5311735.0+45.0/
1230 62491.0*t15-45.0/6578.0*t16+90.0/3289.0*t17-315.0/6578.0*t18-105.0/26.0*t19
1231 -90.0/13.0*t20+36.0/13.0*t21+150.0/3289.0*t23-1575.0/3289.0*t24-30.0/3289.0*t26
1232 +300.0/3289.0*t27-1260.0/3289.0*t28+210.0/299.0*t29-6.0/13.0*t30;
1233       t37 = t2*r;
1234       t43 = -63.0/13.0*t31-315.0/3289.0*t32+126.0/13.0*t33+21.0/299.0*t37*s
1235 -105.0/13.0*t37*t4-21.0/26.0*t37*t3+6.0*t37*t12-63.0/32890.0*t37+252.0/65.0*t37
1236 *t9+9.0/3289.0*t2-9.0/6578.0*t7-9.0/16445.0*t12+9.0/62491.0*t3-9.0/16445.0*t9+
1237 3.0/3289.0*t4+18.0/62491.0*t1-s/62491-3.0/124982.0*r;
1238       drv[0] = t25+t34;
1239       drv[1] = t36+t43;
1240    /* p=14 6 5 */
1241    } else if( i==6 && j==5 ) {
1242       t1 = r*r;
1243       t2 = t1*t1;
1244       t3 = s*s;
1245       t4 = t3*t2;
1246       t5 = t2*r;
1247       t6 = t3*t3;
1248       t7 = t6*s;
1249       t9 = t1*r;
1250       t10 = t9*s;
1251       t11 = t5*t6;
1252       t12 = t5*t3;
1253       t13 = t3*s;
1254       t14 = t9*t13;
1255       t15 = t9*t3;
1256       t16 = t1*s;
1257       t17 = t1*t3;
1258       t18 = r*t13;
1259       t19 = 210.0/299.0*t4+6.0*t7*t5+3.0/5311735.0+90.0/3289.0*t10-90.0/13.0*
1260 t11-6.0/13.0*t12+630.0/299.0*t14-1260.0/3289.0*t15-45.0/6578.0*t16+300.0/3289.0
1261 *t17+3.0/3289.0*t2-9.0/16445.0*t9+9.0/62491.0*t1-9.0/16445.0*t5+9.0/3289.0*t6
1262 -9.0/6578.0*t13+150.0/3289.0*t18+18.0/62491.0*t3;
1263       t21 = r*s;
1264       t22 = t3*r;
1265       t23 = r*t6;
1266       t24 = t5*t13;
1267       t25 = t9*t6;
1268       t27 = t5*s;
1269       t29 = t1*t13;
1270       t30 = t1*t6;
1271       t32 = t2*t13;
1272       t33 = t2*s;
1273       t34 = t2*t6;
1274       t35 = -63.0/32890.0*t7-3.0/124982.0*s-r/62491+21.0/299.0*r*t7+45.0/
1275 62491.0*t21-30.0/3289.0*t22-315.0/3289.0*t23+36.0/13.0*t24-63.0/13.0*t25-105.0/
1276 13.0*t2*t7+9.0/299.0*t27+252.0/65.0*t7*t9-1575.0/3289.0*t29+315.0/299.0*t30
1277 -21.0/26.0*t1*t7-105.0/26.0*t32-315.0/6578.0*t33+126.0/13.0*t34;
1278       t37 = t2*t1;
1279       t39 = 5.0*t37*t6+945.0/598.0*t4+1.0/2124694.0+200.0/3289.0*t10-105.0/13.0
1280 *t11-63.0/26.0*t12+420.0/299.0*t14-1575.0/3289.0*t15-30.0/3289.0*t16+225.0/
1281 3289.0*t17+45.0/6578.0*t2-15.0/6578.0*t9+45.0/124982.0*t1-63.0/6578.0*t5+t6
1282 /6578-3.0/16445.0*t13+36.0/3289.0*t18;
1283       t43 = 9.0/124982.0*t3-2.0/187473.0*s-3.0/124982.0*r+36.0/62491.0*t21-27.0
1284 /6578.0*t22-63.0/6578.0*t23+504.0/65.0*t24-35.0/26.0*t25+84.0/299.0*t27+18.0/
1285 13.0*t37*t3-60.0/13.0*t37*t13-2.0/13.0*t37*s+3.0/598.0*t37-630.0/3289.0*t29+
1286 105.0/598.0*t30-63.0/13.0*t32-630.0/3289.0*t33+63.0/13.0*t34;
1287       drv[0] = t19+t35;
1288       drv[1] = t39+t43;
1289    /* p=14 7 4 */
1290    } else if( i==7 && j==4 ) {
1291       t1 = s*s;
1292       t2 = t1*r;
1293       t3 = r*s;
1294       t4 = t1*t1;
1295       t6 = r*r;
1296       t7 = t6*t6;
1297       t8 = r*t7;
1298       t9 = t1*s;
1299       t10 = t8*t9;
1300       t11 = t7*t6;
1301       t12 = t11*s;
1302       t13 = r*t6;
1303       t15 = t8*s;
1304       t17 = r*t9;
1305       t18 = 14.0/15935205.0+84.0/312455.0*t1-168.0/16445.0*t2+336.0/312455.0*t3
1306 -588.0/16445.0*r*t4+672.0/65.0*t10-7.0/65.0*t12-196.0/65.0*t13*t4+336.0/1495.0*
1307 t15+7.0*t11*t4-84.0/16445.0*t8+112.0/3289.0*t17-14.0/16445.0*t9-28.0/937365.0*r
1308 -28.0/937365.0*s+7.0/2990.0*t11+14.0/3289.0*t7;
1309       t20 = t8*t1;
1310       t21 = t9*t13;
1311       t22 = t13*s;
1312       t23 = t1*t13;
1313       t24 = t6*s;
1314       t25 = t9*t6;
1315       t27 = t1*t7;
1316       t28 = t7*s;
1317       t29 = t11*t1;
1318       t30 = t11*t9;
1319       t32 = t7*t9;
1320       t33 = t1*t6;
1321       t34 = -168.0/13.0*t8*t4-168.0/65.0*t20+784.0/299.0*t21+224.0/3289.0*t22
1322 -2352.0/3289.0*t23-42.0/3289.0*t24-1470.0/3289.0*t25+147.0/299.0*t6*t4+588.0/
1323 299.0*t27-588.0/3289.0*t28+84.0/65.0*t29-28.0/16445.0*t13+14.0/16445.0*t4-70.0/
1324 13.0*t30+21.0/62491.0*t6+588.0/65.0*t7*t4-98.0/13.0*t32+420.0/3289.0*t33;
1325       t36 = 8.0/15935205.0+12.0/312455.0*t1-42.0/16445.0*t2+168.0/312455.0*t3+
1326 2352.0/325.0*t10-56.0/65.0*t12+1176.0/1495.0*t15-588.0/16445.0*t8+56.0/16445.0*
1327 t17-4.0/82225.0*t9-28.0/937365.0*r-8.0/937365.0*s+56.0/1495.0*t11+56.0/3289.0*
1328 t7-294.0/65.0*t20+196.0/299.0*t21;
1329       t37 = t7*t13;
1330       t41 = 280.0/3289.0*t22-1470.0/3289.0*t23-168.0/16445.0*t24-1176.0/16445.0
1331 *t25+588.0/299.0*t27-1176.0/3289.0*t28+336.0/65.0*t29-14.0/3289.0*t13-112.0/
1332 13.0*t30+168.0/312455.0*t6+24.0/65.0*t37*s-196.0/65.0*t32+168.0/3289.0*t33-30.0
1333 /13.0*t37*t1+4.0*t37*t9-t37/65;
1334       drv[0] = t18+t34;
1335       drv[1] = t36+t41;
1336    /* p=14 8 3 */
1337    } else if( i==8 && j==3 ) {
1338       t1 = r*s;
1339       t2 = s*s;
1340       t3 = t2*r;
1341       t4 = r*r;
1342       t5 = t4*t4;
1343       t6 = t5*r;
1344       t7 = t2*s;
1345       t9 = t5*t4;
1346       t11 = t9*t2;
1347       t12 = t4*r;
1348       t13 = t5*t12;
1349       t15 = t12*t2;
1350       t16 = t12*s;
1351       t17 = t4*s;
1352       t18 = t4*t2;
1353       t20 = t13*t2;
1354       t22 = t2*t5;
1355       t23 = t5*s;
1356       t24 = 3024.0/1562275.0*t1-1008.0/82225.0*t3+6048.0/325.0*t6*t7-252.0/13.0
1357 *t9*t7+3024.0/325.0*t11+8.0*t13*t7-21168.0/16445.0*t15+3024.0/16445.0*t16
1358 -2268.0/82225.0*t17+3024.0/16445.0*t18-5292.0/16445.0*t4*t7-48.0/13.0*t20-588.0
1359 /65.0*t5*t7+7056.0/1495.0*t22-10584.0/16445.0*t23-4.0/325.0*t13;
1360       t25 = t6*s;
1361       t26 = t6*t2;
1362       t28 = t13*s;
1363       t29 = t9*s;
1364       t31 = 9072.0/7475.0*t25-3024.0/325.0*t26+3528.0/1495.0*t12*t7+144.0/325.0
1365 *t28-504.0/82225.0*t12+336.0/16445.0*t5-378.0/325.0*t29+1512.0/1562275.0*t4
1366 -3024.0/82225.0*t6-36.0/82225.0*t7+336.0/16445.0*r*t7+432.0/1562275.0*t2+252.0/
1367 7475.0*t9-112.0/1562275.0*r-72.0/1562275.0*s+48.0/26558675.0;
1368       t33 = t5*t5;
1369       t36 = 864.0/1562275.0*t1-108.0/82225.0*t3+3.0*t33*t2+18.0/325.0*t33+
1370 3024.0/325.0*t11-5292.0/16445.0*t15+2016.0/16445.0*t16-1008.0/82225.0*t17+504.0
1371 /16445.0*t18-108.0/13.0*t20+2646.0/1495.0*t22-10584.0/16445.0*t23-12.0/13.0*t33
1372 *s;
1373       t37 = -54.0/325.0*t13+14112.0/7475.0*t25-1764.0/325.0*t26+864.0/325.0*t28
1374 -756.0/82225.0*t12+756.0/16445.0*t5-1008.0/325.0*t29+1512.0/1562275.0*t4
1375 -10584.0/82225.0*t6+27.0/1562275.0*t2+1512.0/7475.0*t9-72.0/1562275.0*r-12.0/
1376 1562275.0*s+18.0/26558675.0;
1377       drv[0] = t24+t31;
1378       drv[1] = t36+t37;
1379    /* p=14 9 2 */
1380    } else if( i==9 && j==2 ) {
1381       t1 = r*r;
1382       t2 = t1*t1;
1383       t3 = s*s;
1384       t5 = t1*r;
1385       t6 = t5*s;
1386       t7 = r*s;
1387       t9 = t2*t2;
1388       t10 = t9*s;
1389       t11 = t2*r;
1390       t14 = t1*s;
1391       t16 = t2*t5;
1392       t18 = t2*s;
1393       t19 = 1134.0/312455.0*t1+2646.0/299.0*t2*t3+9072.0/16445.0*t6+1296.0/
1394 312455.0*t7-216.0/16445.0*t3*r-27.0/13.0*t10-1512.0/65.0*t11*t3+27.0/5311735.0
1395 -31752.0/16445.0*t5*t3-1134.0/16445.0*t14+756.0/3289.0*t1*t3-360.0/13.0*t16*t3
1396 -7938.0/3289.0*t18;
1397       t20 = t2*t1;
1398       t23 = t16*s;
1399       t24 = t20*s;
1400       t25 = t11*s;
1401       t26 = 2268.0/65.0*t20*t3+9.0*t9*t3+432.0/65.0*t23+378.0/3289.0*t2+27.0/
1402 325.0*t9-18.0/65.0*t16-2268.0/82225.0*t5-4536.0/16445.0*t11-567.0/65.0*t24+
1403 9072.0/1495.0*t25-27.0/312455.0*s+567.0/1495.0*t20+81.0/312455.0*t3-72.0/
1404 312455.0*r;
1405       t28 = t9*r;
1406       t30 = 6.0/5311735.0-27.0/312455.0*r+162.0/312455.0*t7+1512.0/1495.0*t20
1407 -7938.0/16445.0*t11-2.0/312455.0*s+54.0/65.0*t9+5292.0/1495.0*t25-81.0/65.0*t16
1408 -378.0/16445.0*t5-3.0/13.0*t28+504.0/3289.0*t6+648.0/312455.0*t1-216.0/16445.0*
1409 t14+648.0/65.0*t23+2.0*t28*s-15876.0/16445.0*t18+2268.0/16445.0*t2-504.0/65.0*
1410 t24-90.0/13.0*t10;
1411       drv[0] = t19+t26;
1412       drv[1] = t30;
1413    /* p=14 10 1 */
1414    } else if( i==10 && j==1 ) {
1415       t1 = r*r;
1416       t2 = t1*t1;
1417       t3 = t2*t2;
1418       t4 = t3*r;
1419       t7 = t1*r;
1420       t8 = t2*t7;
1421       t10 = t2*t1;
1422       t12 = t2*r;
1423       t18 = 10.0*t4*s-10.0/13.0*t4-495.0/13.0*t3*s+198.0/65.0*t3+792.0/13.0*t8*
1424 s-66.0/13.0*t8-693.0/13.0*t10*s+1386.0/299.0*t10+8316.0/299.0*t12*s-756.0/299.0
1425 *t12-2646.0/299.0*t2*s+252.0/299.0*t2+504.0/299.0*t7*s-252.0/1495.0*t7-54.0/
1426 299.0*t1*s+108.0/5681.0*t1+54.0/5681.0*r*s-6.0/5681.0*r-s/5681+2.0/96577.0;
1427       t20 = t1*t3-55.0/13.0*t4+99.0/13.0*t3-99.0/13.0*t8+1386.0/299.0*t10
1428 -2646.0/1495.0*t12+126.0/299.0*t2-18.0/299.0*t7+27.0/5681.0*t1-r/5681+1.0/
1429 482885.0;
1430       drv[0] = t18;
1431       drv[1] = t20;
1432    /* p=14 11 0 */
1433    } else if( i==11 && j==0 ) {
1434       t1 = r*r;
1435       t2 = t1*t1;
1436       t3 = t2*t2;
1437       t6 = t1*r;
1438       t10 = 11.0*t1*t3-660.0/13.0*t3*r+6534.0/65.0*t3-1452.0/13.0*t2*t6+22869.0
1439 /299.0*t2*t1-49896.0/1495.0*t2*r+2772.0/299.0*t2-2376.0/1495.0*t6+891.0/5681.0*
1440 t1-44.0/5681.0*r+66.0/482885.0;
1441       drv[0] = t10;
1442       drv[1] = 0.0;
1443    /* p=15 0 12 */
1444    } else if( i==0 && j==12 ) {
1445       t1 = s*s;
1446       t2 = t1*s;
1447       t3 = t1*t1;
1448       t4 = t3*t3;
1449       t11 = 12.0*t4*t2-429.0/7.0*t4*t1+2860.0/21.0*t4*s-1210.0/7.0*t4+968.0/7.0
1450 *t2*t3-363.0/5.0*t1*t3+2904.0/115.0*t3*s-132.0/23.0*t3+132.0/161.0*t2-11.0/
1451 161.0*t1+44.0/15295.0*s-2.0/45885.0;
1452       drv[0] = 0.0;
1453       drv[1] = t11;
1454    /* p=15 1 11 */
1455    } else if( i==1 && j==11 ) {
1456       t1 = s*s;
1457       t2 = t1*s;
1458       t3 = t1*t1;
1459       t4 = t3*t3;
1460       t6 = t4*t1;
1461       t7 = t4*s;
1462       t8 = t3*t2;
1463       t9 = t3*t1;
1464       t10 = t3*s;
1465       t11 = t4*t2-33.0/7.0*t6+605.0/63.0*t7-3025.0/273.0*t4+726.0/91.0*t8-242.0
1466 /65.0*t9+1694.0/1495.0*t10-66.0/299.0*t3+55.0/2093.0*t2-11.0/6279.0*t1+11.0/
1467 198835.0*s-1.0/1789515.0;
1468       t22 = 11.0*r*t6-330.0/7.0*r*t7+605.0/7.0*r*t4-24200.0/273.0*r*t8+726.0/
1469 13.0*r*t9-1452.0/65.0*r*t10+1694.0/299.0*r*t3-264.0/299.0*r*t2+165.0/2093.0*r*
1470 t1-22.0/6279.0*r*s+11.0/198835.0*r;
1471       t23 = -11.0/14.0*t6+220.0/63.0*t7-605.0/91.0*t4+1936.0/273.0*t8-121.0/
1472 26.0*t9+2904.0/1495.0*t10-154.0/299.0*t3+176.0/2093.0*t2-33.0/4186.0*t1+44.0/
1473 119301.0*s-11.0/1789515.0;
1474       drv[0] = t11;
1475       drv[1] = t22+t23;
1476    /* p=15 2 10 */
1477    } else if( i==2 && j==10 ) {
1478       t1 = s*s;
1479       t2 = t1*t1;
1480       t3 = t2*t2;
1481       t4 = t3*s;
1482       t5 = r*t4;
1483       t6 = r*t3;
1484       t7 = t1*s;
1485       t8 = t2*t1;
1486       t9 = r*t8;
1487       t10 = t2*t7;
1488       t11 = r*t10;
1489       t12 = t3*t1;
1490       t14 = -1.0/3579030.0-55.0/7.0*t5+275.0/21.0*t6+s/39767-3.0/4186.0*t1+20.0
1491 /2093.0*t7+88.0/13.0*t9-1100.0/91.0*t11+2.0*r*t12+r/596505+55.0/63.0*t4;
1492       t15 = r*t2;
1493       t16 = r*s;
1494       t17 = r*t1;
1495       t18 = r*t7;
1496       t19 = t2*s;
1497       t20 = r*t19;
1498       t21 = 154.0/299.0*t15-t16/6279-3.0/14.0*t12-275.0/182.0*t3+10.0/2093.0*
1499 t17-20.0/299.0*t18-154.0/65.0*t20-21.0/299.0*t2+462.0/1495.0*t19+132.0/91.0*t10
1500 -11.0/13.0*t8;
1501       t23 = r*r;
1502       t31 = -77.0/13.0*t23*t2-1.0/715806.0+3.0/39767.0*s+55.0/7.0*t6-t23/12558
1503 -66.0/13.0*t20+10.0*t23*t4-495.0/14.0*t23*t3-550.0/13.0*t23*t8+1100.0/21.0*t23*
1504 t10+10.0/2093.0*t23*s-30.0/299.0*t23*t1+66.0/299.0*t19-21.0/299.0*t2+462.0/
1505 299.0*t15;
1506       t34 = -1100.0/91.0*t11+44.0/91.0*t10-11.0/26.0*t8+132.0/13.0*t9+4.0/299.0
1507 *t7-3.0/2093.0*t1-3.0/2093.0*t16+60.0/2093.0*t17-84.0/299.0*t18-55.0/182.0*t3
1508 -15.0/7.0*t5+r/39767+5.0/63.0*t4+308.0/299.0*t23*t7+264.0/13.0*t23*t19;
1509       drv[0] = t14+t21;
1510       drv[1] = t31+t34;
1511    /* p=15 3 9 */
1512    } else if( i==3 && j==9 ) {
1513       t1 = r*r;
1514       t2 = t1*s;
1515       t3 = s*s;
1516       t4 = t3*t3;
1517       t5 = r*t4;
1518       t6 = t3*s;
1519       t7 = t4*t6;
1520       t8 = r*t7;
1521       t9 = t4*t3;
1522       t10 = r*t9;
1523       t11 = t4*s;
1524       t12 = t1*t7;
1525       t13 = r*t11;
1526       t14 = t4*t4;
1527       t15 = t14*s;
1528       t17 = r*t14;
1529       t18 = t1*t3;
1530       t19 = t1*t6;
1531       t20 = t1*t11;
1532       t21 = t1*t4;
1533       t22 = 9.0/23023.0*t2-378.0/16445.0*t4+504.0/1495.0*t5-360.0/91.0*t8+192.0
1534 /65.0*t10+126.0/1495.0*t11+108.0/455.0*t7+90.0/7.0*t12-84.0/65.0*t13-6.0/7.0*r*
1535 t15+20.0/7.0*t17-36.0/3289.0*t18+42.0/299.0*t19+252.0/65.0*t20-63.0/65.0*t21;
1536       t23 = r*t3;
1537       t24 = r*s;
1538       t26 = t1*t14;
1539       t27 = r*t6;
1540       t28 = t1*t9;
1541       t29 = 96.0/23023.0*t23-18.0/115115.0*t24+4.0/2187185.0*r-15.0/91.0*t14+
1542 t15/21+3.0*t1*t15-135.0/14.0*t26-t1/230230+27.0/2187185.0*s-36.0/115115.0*t3
1543 -168.0/3289.0*t27-120.0/13.0*t28+12.0/3289.0*t6-1.0/6561555.0-12.0/65.0*t9;
1544       t31 = t1*r;
1545       t34 = 96.0/23023.0*t2-42.0/3289.0*t4+126.0/299.0*t5-120.0/91.0*t8+108.0/
1546 65.0*t10+48.0/1495.0*t11+30.0*t31*t9-180.0/7.0*t31*t7-1.0/2187185.0+3.0/23023.0
1547 *t31+16.0/455.0*t7+80.0/7.0*t12-72.0/65.0*t13+3.0/7.0*t17-252.0/3289.0*t18+
1548 1008.0/1495.0*t19+576.0/65.0*t20-42.0/13.0*t21;
1549       t41 = -84.0/65.0*t31*t6+84.0/13.0*t31*t4+9.0*t31*t14+42.0/299.0*t31*t3
1550 -24.0/3289.0*t31*s-240.0/13.0*t31*t11+36.0/3289.0*t23-72.0/115115.0*t24+27.0/
1551 2187185.0*r-t14/91-27.0/7.0*t26-9.0/115115.0*t1+48.0/2187185.0*s-6.0/16445.0*t3
1552 -1512.0/16445.0*t27-180.0/13.0*t28+48.0/16445.0*t6-3.0/65.0*t9;
1553       drv[0] = t22+t29;
1554       drv[1] = t34+t41;
1555    /* p=15 4 8 */
1556    } else if( i==4 && j==8 ) {
1557       t1 = s*s;
1558       t2 = t1*t1;
1559       t3 = t2*t1;
1560       t4 = r*r;
1561       t5 = t4*t3;
1562       t6 = t1*s;
1563       t7 = r*t6;
1564       t8 = r*t1;
1565       t9 = r*s;
1566       t10 = t2*t2;
1567       t13 = t4*s;
1568       t14 = t4*t2;
1569       t15 = t4*t1;
1570       t16 = t4*r;
1571       t17 = -8.0/195.0*t3-80.0/13.0*t5+16.0/9867.0*t6-112.0/3289.0*t7-8.0/
1572 49335.0*t1+32.0/9867.0*t8-16.0/115115.0*t9-15.0/7.0*t4*t10-10.0/819.0*t10-t4
1573 /115115+16.0/2187185.0*s+4.0/2187185.0*r+20.0/63.0*r*t10+16.0/23023.0*t13-14.0/
1574 13.0*t14-28.0/3289.0*t2-56.0/3289.0*t15+4.0/345345.0*t16;
1575       t18 = t2*t6;
1576       t19 = t4*t18;
1577       t20 = r*t18;
1578       t21 = t16*t18;
1579       t22 = t2*s;
1580       t23 = t4*t22;
1581       t24 = t16*s;
1582       t25 = t4*t6;
1583       t26 = t16*t3;
1584       t27 = t16*t2;
1585       t28 = t16*t22;
1586       t29 = r*t3;
1587       t30 = t16*t1;
1588       t32 = r*t2;
1589       t33 = r*t22;
1590       t34 = t16*t6;
1591       t35 = 40.0/7.0*t19+16.0/455.0*t18-80.0/91.0*t20-72.0/7.0*t21+224.0/65.0*
1592 t23+112.0/4485.0*t22-16.0/16445.0*t24+56.0/299.0*t25+32.0/3.0*t26+112.0/65.0*
1593 t27-224.0/39.0*t28+64.0/65.0*t29+112.0/4485.0*t30+4.0*t16*t10+56.0/299.0*t32
1594 -112.0/195.0*t33-56.0/195.0*t34-2.0/19684665.0;
1595       t37 = -t3/195-40.0/13.0*t5+8.0/9867.0*t6-112.0/3289.0*t7-2.0/16445.0*t1+
1596 16.0/3289.0*t8-16.0/49335.0*t9-8.0/115115.0*t4+8.0/937365.0*s+16.0/2187185.0*r+
1597 32.0/9867.0*t13-56.0/39.0*t14-28.0/9867.0*t2-168.0/3289.0*t15+16.0/69069.0*t16+
1598 80.0/63.0*t19+8.0/4095.0*t18-80.0/819.0*t20-40.0/7.0*t21+192.0/65.0*t23;
1599       t38 = t4*t4;
1600       t46 = 8.0/1495.0*t22-112.0/9867.0*t24+112.0/299.0*t25+40.0/3.0*t26+224.0/
1601 39.0*t27-160.0/13.0*t28+16.0/65.0*t29+56.0/299.0*t30+112.0/897.0*t32-16.0/65.0*
1602 t33-4.0/16445.0*t38+8.0*t38*t18-56.0/39.0*t34-280.0/39.0*t38*t2+112.0/65.0*t38*
1603 t6-14.0/65.0*t38*t1+56.0/4485.0*t38*s+16.0*t38*t22-18.0*t38*t3-4.0/19684665.0;
1604       drv[0] = t17+t35;
1605       drv[1] = t37+t46;
1606    /* p=15 5 7 */
1607    } else if( i==5 && j==7 ) {
1608       t1 = s*s;
1609       t2 = t1*s;
1610       t3 = t1*t1;
1611       t4 = t3*t2;
1612       t6 = t3*t1;
1613       t7 = r*t6;
1614       t8 = r*r;
1615       t9 = t8*t8;
1616       t10 = t9*t1;
1617       t11 = r*t1;
1618       t12 = r*s;
1619       t13 = t8*t2;
1620       t14 = t3*s;
1621       t15 = t8*t14;
1622       t16 = t9*t6;
1623       t17 = t9*t14;
1624       t18 = t8*r;
1625       t20 = t9*s;
1626       t21 = r*t3;
1627       t22 = r*t14;
1628       t23 = t18*t1;
1629       t24 = t18*s;
1630       t25 = t18*t14;
1631       t26 = t18*t3;
1632       t27 = t18*t2;
1633       t28 = t8*s;
1634       t29 = -100.0/819.0*r*t4+32.0/117.0*t7-7.0/117.0*t10+80.0/29601.0*t11-4.0/
1635 29601.0*t12-1.0/11810799.0+175.0/897.0*t13+28.0/13.0*t15-10.0*t16+70.0/9.0*t17
1636 -30.0/7.0*t4*t18+7.0/2691.0*t20+280.0/2691.0*t21-28.0/117.0*t22+56.0/897.0*t23
1637 -28.0/9867.0*t24-280.0/39.0*t25+112.0/39.0*t26-70.0/117.0*t27+10.0/9867.0*t28;
1638       t30 = t9*t3;
1639       t31 = r*t2;
1640       t32 = t8*t6;
1641       t35 = t8*t3;
1642       t36 = t9*t2;
1643       t37 = t8*t1;
1644       t38 = t18*t6;
1645       t39 = -t1/9867+25.0/29601.0*t2+8.0/3936933.0*r+s/187473-t8/69069-t6/117+
1646 t4/273-t9/29601+8.0/207207.0*t18-35.0/9867.0*t3+7.0/897.0*t14-350.0/117.0*t30
1647 -700.0/29601.0*t31-100.0/39.0*t32+25.0/21.0*t8*t4+5.0*t9*t4-35.0/39.0*t35+70.0/
1648 117.0*t36-70.0/3289.0*t37+80.0/9.0*t38;
1649       t41 = t9*r;
1650       t46 = 28.0/897.0*t20+35.0/897.0*t21-2.0/39.0*t22+175.0/897.0*t23-140.0/
1651 9867.0*t24-200.0/39.0*t25+140.0/39.0*t26-140.0/117.0*t27-14.0/585.0*t41*s+14.0/
1652 39.0*t41*t1-12.0*t41*t14;
1653       t52 = -7.0/9867.0*t9+10.0/29601.0*t18-7.0/9867.0*t3+4.0/4485.0*t14-350.0/
1654 39.0*t30-140.0/9867.0*t31-50.0/117.0*t32-70.0/117.0*t35+112.0/39.0*t36-350.0/
1655 9867.0*t37+25.0/9.0*t38;
1656       drv[0] = t29+t39;
1657       drv[1] = t7/39-35.0/78.0*t10+25.0/9867.0*t11-2.0/9867.0*t12-1.0/8436285.0
1658 +7.0/13455.0*t41+560.0/2691.0*t13+32.0/39.0*t15-15.0/2.0*t16+40.0/3.0*t17+t46+
1659 7.0*t41*t6-280.0/117.0*t41*t2+70.0/9.0*t41*t3+80.0/29601.0*t28-t1/19734+8.0/
1660 29601.0*t2+r/187473+4.0/937365.0*s-2.0/29601.0*t8-t6/2340+t52;
1661    /* p=15 6 6 */
1662    } else if( i==6 && j==6 ) {
1663       t1 = r*r;
1664       t2 = t1*t1;
1665       t3 = t2*r;
1666       t4 = s*s;
1667       t5 = t4*t4;
1668       t6 = t3*t5;
1669       t7 = t4*s;
1670       t8 = t1*t7;
1671       t9 = r*t5;
1672       t10 = t3*s;
1673       t11 = r*s;
1674       t12 = r*t7;
1675       t13 = t5*t4;
1676       t15 = t1*t5;
1677       t16 = t3*t4;
1678       t17 = t1*s;
1679       t19 = t2*t4;
1680       t20 = t2*t7;
1681       t21 = t1*t4;
1682       t22 = t1*r;
1683       t24 = t5*s;
1684       t25 = t3*t24;
1685       t27 = t2*t5;
1686       t28 = t3*t7;
1687       t29 = t1*t24;
1688       t32 = -35.0/156.0*t19+70.0/39.0*t20-175.0/6578.0*t21+10.0/3.0*t22*t13-9.0
1689 *t25+6.0*t3*t13-175.0/26.0*t27-50.0/39.0*t28+14.0/13.0*t29-15.0/2.0*t2*t13+2.0/
1690 39.0*r*t13;
1691       t34 = t22*t4;
1692       t35 = t22*s;
1693       t36 = t22*t24;
1694       t37 = t22*t5;
1695       t38 = t22*t7;
1696       t39 = r*t4;
1697       t41 = r*t24;
1698       t42 = t2*s;
1699       t43 = t2*t24;
1700       t44 = 3.0/624910.0*s+t3/8970-t4/13156+r/374946+7.0/2990.0*t24-7.0/39468.0
1701 *t2+t22/9867-21.0/13156.0*t5-7.0/78.0*t41+7.0/598.0*t42+35.0/3.0*t43;
1702       t48 = 42.0/13.0*t20-175.0/6578.0*t21-9.0*t25-175.0/26.0*t27-70.0/13.0*t28
1703 +2.0/13.0*t29+175.0/897.0*t34-175.0/9867.0*t35-50.0/39.0*t36+70.0/39.0*t37-35.0
1704 /39.0*t38;
1705       t50 = t2*t1;
1706       t57 = t24/8970-21.0/13156.0*t2+5.0/9867.0*t22-7.0/39468.0*t5-t50/780-t41
1707 /130+35.0/598.0*t42+5.0*t43-15.0/2.0*t50*t5+10.0/3.0*t50*t7+6.0*t50*t24;
1708       drv[0] = 5.0*t6+175.0/897.0*t8+35.0/598.0*t9-t10/130-t11/6578-175.0/
1709 9867.0*t12-25.0/39.0*t1*t13-35.0/52.0*t15+2.0/13.0*t16+5.0/3289.0*t17+t32+35.0/
1710 299.0*t34-21.0/3289.0*t35-70.0/13.0*t36+42.0/13.0*t37-35.0/39.0*t38-1.0/
1711 11248380.0+5.0/9867.0*t7-t1/39468-t13/780+25.0/9867.0*t39+t44;
1712       drv[1] = 35.0/3.0*t6+35.0/299.0*t8+7.0/598.0*t9-7.0/78.0*t10-t11/6578
1713 -21.0/3289.0*t12-35.0/156.0*t15+14.0/13.0*t16+25.0/9867.0*t17-35.0/52.0*t19+t48
1714 +2.0/39.0*t50*s-25.0/39.0*t50*t4-1.0/11248380.0+t7/9867-t1/13156+5.0/3289.0*t39
1715 +s/374946+7.0/2990.0*t3-t4/39468+3.0/624910.0*r+t57;
1716    /* p=15 7 5 */
1717    } else if( i==7 && j==5 ) {
1718       t1 = r*r;
1719       t2 = t1*t1;
1720       t3 = s*s;
1721       t4 = t3*t3;
1722       t5 = t4*s;
1723       t8 = t3*s;
1724       t9 = t2*t8;
1725       t10 = t2*r;
1726       t11 = t10*s;
1727       t12 = r*s;
1728       t13 = r*t3;
1729       t14 = t2*t1;
1730       t15 = t14*s;
1731       t16 = t14*t3;
1732       t17 = t1*t4;
1733       t18 = t10*t3;
1734       t20 = t14*t4;
1735       t21 = t14*t8;
1736       t22 = t2*t3;
1737       t23 = t1*t3;
1738       t26 = t2*t4;
1739       t27 = t10*t8;
1740       t28 = t10*t4;
1741       t29 = t1*t8;
1742       t30 = t2*s;
1743       t31 = -15.0/2.0*t20+25.0/9.0*t21-70.0/117.0*t22-350.0/9867.0*t23+7.0*t14*
1744 t5-12.0*t10*t5-350.0/39.0*t26-200.0/39.0*t27+40.0/3.0*t28+175.0/897.0*t29+35.0/
1745 897.0*t30;
1746       t33 = r*t4;
1747       t35 = t1*r;
1748       t36 = t35*t3;
1749       t37 = t35*s;
1750       t39 = t35*t4;
1751       t40 = t35*t8;
1752       t42 = r*t8;
1753       t43 = t1*s;
1754       t44 = -t1/19734+8.0/29601.0*t35-140.0/9867.0*t42-t14/2340+4.0/4485.0*t10
1755 -2.0/29601.0*t3+7.0/13455.0*t5-7.0/9867.0*t2+25.0/9867.0*t43-7.0/9867.0*t4-1.0/
1756 8436285.0;
1757       t47 = t2*t35;
1758       t50 = 5.0*t47*t4+25.0/21.0*t47*t3+112.0/39.0*t9-1.0/11810799.0-28.0/117.0
1759 *t11-4.0/29601.0*t12+10.0/9867.0*t13+32.0/117.0*t15-100.0/39.0*t16-7.0/117.0*
1760 t17+28.0/13.0*t18-10.0*t20+80.0/9.0*t21-35.0/39.0*t22-70.0/3289.0*t23-350.0/
1761 117.0*t26-280.0/39.0*t27+70.0/9.0*t28+56.0/897.0*t29+280.0/2691.0*t30;
1762       t53 = 7.0/2691.0*t33+175.0/897.0*t36-700.0/29601.0*t37+70.0/117.0*t39
1763 -70.0/117.0*t40-100.0/819.0*t47*s-30.0/7.0*t47*t8+8.0/3936933.0*s+r/187473+8.0/
1764 207207.0*t8-t1/9867+25.0/29601.0*t35-28.0/9867.0*t42-t14/117+7.0/897.0*t10-t3
1765 /69069-35.0/9867.0*t2+80.0/29601.0*t43-t4/29601+t47/273;
1766       drv[0] = 70.0/9.0*t2*t5+14.0/39.0*t1*t5+140.0/39.0*t9-2.0/39.0*t11-2.0/
1767 9867.0*t12+80.0/29601.0*t13+t15/39-50.0/117.0*t16-35.0/78.0*t17+32.0/39.0*t18+
1768 t31+28.0/897.0*t33-14.0/585.0*r*t5+560.0/2691.0*t36-140.0/9867.0*t37-280.0/
1769 117.0*t35*t5+112.0/39.0*t39-140.0/117.0*t40+s/187473+4.0/937365.0*r+10.0/
1770 29601.0*t8+t44;
1771       drv[1] = t50+t53;
1772    /* p=15 8 4 */
1773    } else if( i==8 && j==4 ) {
1774       t1 = r*s;
1775       t2 = s*s;
1776       t3 = t2*r;
1777       t4 = r*r;
1778       t5 = r*t4;
1779       t6 = t4*t4;
1780       t7 = t6*t5;
1781       t8 = t7*s;
1782       t9 = t2*s;
1783       t10 = r*t9;
1784       t11 = t6*t4;
1785       t12 = t11*s;
1786       t13 = t7*t9;
1787       t14 = t11*t2;
1788       t15 = t6*r;
1789       t16 = t15*t2;
1790       t17 = t2*t2;
1791       t19 = t2*t4;
1792       t20 = t6*t9;
1793       t21 = t11*t9;
1794       t22 = t6*t2;
1795       t23 = t4*s;
1796       t26 = t15*t9;
1797       t27 = -16.0/49335.0*t1-4.0/19684665.0+8.0/937365.0*r+16.0/2187185.0*s+
1798 32.0/9867.0*t3-80.0/819.0*t8-112.0/9867.0*t10+16.0/65.0*t12-40.0/7.0*t13-40.0/
1799 13.0*t14+192.0/65.0*t16+8.0*t7*t17-168.0/3289.0*t19+224.0/39.0*t20+40.0/3.0*t21
1800 -56.0/39.0*t22+16.0/3289.0*t23-18.0*t11*t17-280.0/39.0*t6*t17-160.0/13.0*t26;
1801       t29 = t7*t2;
1802       t31 = t15*s;
1803       t32 = t5*t9;
1804       t33 = t5*s;
1805       t34 = t2*t5;
1806       t36 = t4*t9;
1807       t37 = t6*s;
1808       t39 = 16.0*t15*t17-28.0/9867.0*t6+80.0/63.0*t29-14.0/65.0*t4*t17-4.0/
1809 16445.0*t17+8.0/9867.0*t5-2.0/16445.0*t4+8.0/4095.0*t7+16.0/69069.0*t9-8.0/
1810 115115.0*t2+8.0/1495.0*t15-16.0/65.0*t31-t11/195-56.0/39.0*t32-112.0/3289.0*t33
1811 +112.0/299.0*t34+56.0/4485.0*r*t17+56.0/299.0*t36+112.0/897.0*t37+112.0/65.0*t5
1812 *t17;
1813       t41 = t6*t6;
1814       t43 = -16.0/115115.0*t1+4.0*t41*t9-2.0/19684665.0+16.0/2187185.0*r+4.0/
1815 2187185.0*s+16.0/23023.0*t3-80.0/91.0*t8-16.0/16445.0*t10+64.0/65.0*t12-72.0/
1816 7.0*t13-80.0/13.0*t14+224.0/65.0*t16-56.0/3289.0*t19+112.0/65.0*t20+32.0/3.0*
1817 t21-14.0/13.0*t22+32.0/9867.0*t23-224.0/39.0*t26;
1818       t46 = -28.0/3289.0*t6+40.0/7.0*t29+16.0/9867.0*t5-8.0/49335.0*t4+16.0/
1819 455.0*t7+4.0/345345.0*t9-t2/115115+112.0/4485.0*t15-112.0/195.0*t31-8.0/195.0*
1820 t11+20.0/63.0*t41*s-10.0/819.0*t41-15.0/7.0*t41*t2-56.0/195.0*t32-112.0/3289.0*
1821 t33+56.0/299.0*t34+112.0/4485.0*t36+56.0/299.0*t37;
1822       drv[0] = t27+t39;
1823       drv[1] = t43+t46;
1824    /* p=15 9 3 */
1825    } else if( i==9 && j==3 ) {
1826       t1 = r*r;
1827       t2 = t1*r;
1828       t3 = t1*t1;
1829       t4 = t3*t2;
1830       t5 = s*s;
1831       t6 = r*t5;
1832       t7 = t4*s;
1833       t8 = t3*t3;
1834       t9 = t8*t5;
1835       t10 = t3*r;
1836       t11 = t3*t1;
1837       t12 = t10*s;
1838       t13 = r*s;
1839       t14 = t3*s;
1840       t15 = t2*t5;
1841       t16 = t2*s;
1842       t17 = t5*s;
1843       t19 = t8*s;
1844       t20 = -1.0/2187185.0+16.0/455.0*t4+96.0/23023.0*t6-120.0/91.0*t7-t8/91
1845 -9.0/115115.0*t5-27.0/7.0*t9+48.0/1495.0*t10-3.0/65.0*t11+27.0/2187185.0*s+48.0
1846 /2187185.0*r-72.0/65.0*t12-72.0/115115.0*t13+126.0/299.0*t14+1008.0/1495.0*t15
1847 -1512.0/16445.0*t16-84.0/65.0*t2*t17+3.0/7.0*t19;
1848       t25 = t11*t5;
1849       t26 = t10*t5;
1850       t27 = t4*t5;
1851       t28 = t1*s;
1852       t30 = t3*t5;
1853       t32 = t1*t5;
1854       t34 = t11*s;
1855       t35 = -240.0/13.0*t10*t17-42.0/3289.0*t3+48.0/16445.0*t2-180.0/7.0*t4*t17
1856 -6.0/16445.0*t1-24.0/3289.0*r*t17+9.0*t8*t17-180.0/13.0*t25+576.0/65.0*t26+80.0
1857 /7.0*t27+36.0/3289.0*t28+30.0*t11*t17-42.0/13.0*t30+84.0/13.0*t3*t17-252.0/
1858 3289.0*t32+42.0/299.0*t1*t17+108.0/65.0*t34+3.0/23023.0*t17;
1859       t37 = t8*r;
1860       t39 = -1.0/6561555.0+108.0/455.0*t4+9.0/23023.0*t6-360.0/91.0*t7-15.0/
1861 91.0*t8-t5/230230-135.0/14.0*t9+126.0/1495.0*t10-12.0/65.0*t11+4.0/2187185.0*s+
1862 27.0/2187185.0*r-84.0/65.0*t12-18.0/115115.0*t13-6.0/7.0*t37*s+504.0/1495.0*t14
1863 ;
1864       t41 = 42.0/299.0*t15-168.0/3289.0*t16+20.0/7.0*t19-378.0/16445.0*t3+12.0/
1865 3289.0*t2-36.0/115115.0*t1-120.0/13.0*t25+252.0/65.0*t26+90.0/7.0*t27+96.0/
1866 23023.0*t28-63.0/65.0*t30-36.0/3289.0*t32+192.0/65.0*t34+t37/21+3.0*t37*t5;
1867       drv[0] = t20+t35;
1868       drv[1] = t39+t41;
1869    /* p=15 10 2 */
1870    } else if( i==10 && j==2 ) {
1871       t1 = r*r;
1872       t2 = s*s;
1873       t4 = t1*t1;
1874       t5 = t4*t4;
1875       t6 = t5*r;
1876       t7 = t6*s;
1877       t8 = t4*r;
1878       t9 = r*s;
1879       t10 = t4*t1;
1880       t12 = t1*r;
1881       t13 = t4*t12;
1882       t14 = t13*s;
1883       t17 = -1.0/715806.0-30.0/299.0*t1*t2-15.0/7.0*t7+66.0/299.0*t8-3.0/2093.0
1884 *t9-11.0/26.0*t10+s/39767+3.0/39767.0*r-550.0/13.0*t10*t2+44.0/91.0*t13-1100.0/
1885 91.0*t14+10.0/2093.0*r*t2-t2/12558-495.0/14.0*t5*t2-55.0/182.0*t5;
1886       t18 = t8*s;
1887       t21 = t10*s;
1888       t22 = t1*s;
1889       t25 = t5*s;
1890       t26 = t4*s;
1891       t28 = t12*s;
1892       t29 = -66.0/13.0*t18+264.0/13.0*t8*t2+1100.0/21.0*t13*t2+4.0/299.0*t12
1893 -3.0/2093.0*t1+132.0/13.0*t21+60.0/2093.0*t22+10.0*t6*t2+5.0/63.0*t6-77.0/13.0*
1894 t4*t2-21.0/299.0*t4+55.0/7.0*t25+462.0/299.0*t26+308.0/299.0*t12*t2-84.0/299.0*
1895 t28;
1896       t31 = -1.0/3579030.0+462.0/1495.0*t8+r/39767+275.0/21.0*t25-t9/6279+132.0
1897 /91.0*t13-11.0/13.0*t10+88.0/13.0*t21-154.0/65.0*t18-1100.0/91.0*t14-55.0/7.0*
1898 t7;
1899       t32 = t5*t1;
1900       t34 = -3.0/4186.0*t1+154.0/299.0*t26+s/596505+55.0/63.0*t6+10.0/2093.0*
1901 t22-3.0/14.0*t32-21.0/299.0*t4-20.0/299.0*t28+20.0/2093.0*t12-275.0/182.0*t5+
1902 2.0*t32*s;
1903       drv[0] = t17+t29;
1904       drv[1] = t31+t34;
1905    /* p=15 11 1 */
1906    } else if( i==11 && j==1 ) {
1907       t1 = r*r;
1908       t2 = t1*t1;
1909       t3 = t2*t2;
1910       t4 = t3*t1;
1911       t6 = t3*r;
1912       t9 = t1*r;
1913       t10 = t2*t9;
1914       t12 = t1*t2;
1915       t14 = r*t2;
1916       t16 = 11.0*t4*s-11.0/14.0*t4-330.0/7.0*t6*s+220.0/63.0*t6+605.0/7.0*t3*s
1917 -605.0/91.0*t3-24200.0/273.0*t10*s+1936.0/273.0*t10+726.0/13.0*t12*s-121.0/26.0
1918 *t12-1452.0/65.0*t14*s;
1919       t21 = 2904.0/1495.0*t14+1694.0/299.0*t2*s-154.0/299.0*t2-264.0/299.0*t9*s
1920 +176.0/2093.0*t9+165.0/2093.0*t1*s-33.0/4186.0*t1-22.0/6279.0*r*s+44.0/119301.0
1921 *r+11.0/198835.0*s-11.0/1789515.0;
1922       t24 = t3*t9-33.0/7.0*t4+605.0/63.0*t6-3025.0/273.0*t3+726.0/91.0*t10
1923 -242.0/65.0*t12+1694.0/1495.0*t14-66.0/299.0*t2+55.0/2093.0*t9-11.0/6279.0*t1+
1924 11.0/198835.0*r-1.0/1789515.0;
1925       drv[0] = t16+t21;
1926       drv[1] = t24;
1927    /* p=15 12 0 */
1928    } else if( i==12 && j==0 ) {
1929       t1 = r*r;
1930       t2 = t1*r;
1931       t3 = t1*t1;
1932       t4 = t3*t3;
1933       t11 = 12.0*t4*t2-429.0/7.0*t4*t1+2860.0/21.0*t4*r-1210.0/7.0*t4+968.0/7.0
1934 *t3*t2-363.0/5.0*t3*t1+2904.0/115.0*t3*r-132.0/23.0*t3+132.0/161.0*t2-11.0/
1935 161.0*t1+44.0/15295.0*r-2.0/45885.0;
1936       drv[0] = t11;
1937       drv[1] = 0.0;
1938     } else
1939       return 0;
1940     return 1;
1941 }
1942 #ifdef __cplusplus
1943 }
1944 #endif
1945