19 int BnDrv(int i,int j, int k, double r, double s, double t, double drv[3]) {  in BnDrv()  argument
46       t1 = 6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0;  in BnDrv()
49       drv[2] = 2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0;  in BnDrv()
52       drv[0] = 2.0/7.0*t+6.0/7.0*s+8.0/21.0*r-8.0/21.0;  in BnDrv()
53       drv[1] = t+4.0/3.0*s+6.0/7.0*r-6.0/7.0;  in BnDrv()
54       drv[2] = s+2.0/7.0*r-2.0/7.0;  in BnDrv()
57       drv[0] = 6.0/7.0*s+2.0/7.0*r-2.0/7.0;  in BnDrv()
58       drv[1] = 2.0*s+6.0/7.0*r-6.0/7.0;  in BnDrv()
62       drv[0] = t+2.0/3.0*s+4.0/3.0*r-22.0/27.0;  in BnDrv()
63       drv[1] = 2.0/3.0*r-4.0/27.0;  in BnDrv()
64       drv[2] = r-2.0/9.0;  in BnDrv()
67       drv[0] = s+4.0/5.0*r-22.0/45.0;  in BnDrv()
68       drv[1] = r-2.0/9.0;  in BnDrv()
72       drv[0] = 2.0*r-2.0/3.0;  in BnDrv()
78       t2 = s+r-1.0;  in BnDrv()
88       t2 = s+r-1.0;  in BnDrv()
91       t5 = s+2.0/9.0*r-2.0/9.0;  in BnDrv()
92       t7 = t5*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
95       drv[2] = t5*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
98       t2 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
100       t5 = r-1.0;  in BnDrv()
103       drv[0] = (2.0/3.0*s+r/6-1.0/6.0)*t2+2.0/3.0*t4+4.0/9.0*t6+t7/18;  in BnDrv()
104       drv[1] = (2.0*s+2.0/3.0*r-2.0/3.0)*t2+2.0/3.0*t4+4.0/9.0*t6+t7/18;  in BnDrv()
109       t2 = r-1.0;  in BnDrv()
118       t2 = s+r-1.0;  in BnDrv()
120       t5 = r-2.0/11.0;  in BnDrv()
121       t7 = t5*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
124       drv[2] = t5*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
127       t1 = s+2.0/7.0*r-2.0/7.0;  in BnDrv()
128       t2 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
129       t4 = r-2.0/11.0;  in BnDrv()
138       t2 = r-1.0;  in BnDrv()
140       t5 = r-2.0/11.0;  in BnDrv()
141       drv[0] = t1+6.0/7.0*s*t2+t4/7+t5*(6.0/7.0*s+2.0/7.0*r-2.0/7.0);  in BnDrv()
142       drv[1] = t5*(2.0*s+6.0/7.0*r-6.0/7.0);  in BnDrv()
146       t4 = r*r;  in BnDrv()
147       drv[0] = (2.0*r-6.0/11.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0)+2.0/3.0*t4-4.0/  in BnDrv()
148 11.0*r+2.0/55.0;  in BnDrv()
149       drv[1] = 2.0/3.0*t4-4.0/11.0*r+2.0/55.0;  in BnDrv()
150       drv[2] = t4-6.0/11.0*r+3.0/55.0;  in BnDrv()
153       t4 = r*r;  in BnDrv()
154       drv[0] = (2.0*r-6.0/11.0)*(s+2.0/5.0*r-2.0/5.0)+2.0/5.0*t4-12.0/55.0*r+  in BnDrv()
156       drv[1] = t4-6.0/11.0*r+3.0/55.0;  in BnDrv()
160       t1 = r*r;  in BnDrv()
161       drv[0] = 3.0*t1-24.0/11.0*r+18.0/55.0;  in BnDrv()
168       t3 = s+r-1.0;  in BnDrv()
181       t3 = s+r-1.0;  in BnDrv()
186       t8 = s+2.0/11.0*r-2.0/11.0;  in BnDrv()
195       t3 = s+r-1.0;  in BnDrv()
199       t9 = r-1.0;  in BnDrv()
202       t14 = t12*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
203       drv[0] = (6.0/11.0*s+6.0/55.0*r-6.0/55.0)*t6+t14;  in BnDrv()
204       drv[1] = (2.0*s+6.0/11.0*r-6.0/11.0)*t6+t14;  in BnDrv()
205       drv[2] = t12*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
209       t2 = r-1.0;  in BnDrv()
212       t6 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
226       t3 = r-1.0;  in BnDrv()
237       t3 = s+r-1.0;  in BnDrv()
239       t8 = r-2.0/13.0;  in BnDrv()
247       t1 = s+2.0/9.0*r-2.0/9.0;  in BnDrv()
249       t3 = s+r-1.0;  in BnDrv()
252       t8 = r-2.0/13.0;  in BnDrv()
255       t12 = t10*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
258       drv[2] = t10*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
262       t2 = r-1.0;  in BnDrv()
265       t6 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
266       t8 = r-2.0/13.0;  in BnDrv()
268       drv[0] = t5*t6+t8*(2.0/3.0*s+r/6-1.0/6.0)*t6+2.0/3.0*t12;  in BnDrv()
269       drv[1] = t8*(2.0*s+2.0/3.0*r-2.0/3.0)*t6+2.0/3.0*t12;  in BnDrv()
274       t3 = r-1.0;  in BnDrv()
276       t8 = r-2.0/13.0;  in BnDrv()
284       t3 = s+r-1.0;  in BnDrv()
286       t8 = r*r;  in BnDrv()
287       t9 = t8-6.0/13.0*r+1.0/26.0;  in BnDrv()
288       t11 = t9*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
289       drv[0] = (2.0*r-6.0/13.0)*(t2+6.0/5.0*t*t3+3.0/10.0*t5)+t11;  in BnDrv()
291       drv[2] = t9*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
294       t2 = s+2.0/7.0*r-2.0/7.0;  in BnDrv()
295       t4 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
296       t6 = r*r;  in BnDrv()
297       t7 = t6-6.0/13.0*r+1.0/26.0;  in BnDrv()
300       drv[0] = (2.0*r-6.0/13.0)*t2*t4+2.0/7.0*t8+2.0/3.0*t9;  in BnDrv()
306       t3 = r-1.0;  in BnDrv()
308       t8 = r*r;  in BnDrv()
309       t9 = t8-6.0/13.0*r+1.0/26.0;  in BnDrv()
310       drv[0] = (2.0*r-6.0/13.0)*(t2+6.0/7.0*s*t3+t5/7)+t9*(6.0/7.0*s+2.0/7.0*r  in BnDrv()
312       drv[1] = t9*(2.0*s+6.0/7.0*r-6.0/7.0);  in BnDrv()
316       t1 = r*r;  in BnDrv()
317       t5 = t1*r;  in BnDrv()
318       drv[0] = (t+2.0/3.0*s+2.0/3.0*r-2.0/3.0)*(3.0*t1-24.0/13.0*r+3.0/13.0)+  in BnDrv()
319 2.0/3.0*t5-8.0/13.0*t1+2.0/13.0*r-4.0/429.0;  in BnDrv()
320       drv[1] = 2.0/3.0*t5-8.0/13.0*t1+2.0/13.0*r-4.0/429.0;  in BnDrv()
321       drv[2] = t5-12.0/13.0*t1+3.0/13.0*r-2.0/143.0;  in BnDrv()
324       t1 = r*r;  in BnDrv()
325       t5 = t1*r;  in BnDrv()
326       drv[0] = (s+2.0/5.0*r-2.0/5.0)*(3.0*t1-24.0/13.0*r+3.0/13.0)+2.0/5.0*t5  in BnDrv()
327 -24.0/65.0*t1+6.0/65.0*r-4.0/715.0;  in BnDrv()
328       drv[1] = t5-12.0/13.0*t1+3.0/13.0*r-2.0/143.0;  in BnDrv()
332       t1 = r*r;  in BnDrv()
333       drv[0] = 4.0*t1*r-60.0/13.0*t1+20.0/13.0*r-20.0/143.0;  in BnDrv()
340       t4 = s+r-1.0;  in BnDrv()
355       t4 = s+r-1.0;  in BnDrv()
362       t11 = s+2.0/13.0*r-2.0/13.0;  in BnDrv()
373       t4 = s+r-1.0;  in BnDrv()
377       t12 = r-1.0;  in BnDrv()
382       drv[0] = (6.0/13.0*s+r/13-1.0/13.0)*t9+t18;  in BnDrv()
383       drv[1] = (2.0*s+6.0/13.0*r-6.0/13.0)*t9+t18;  in BnDrv()
388       t2 = r-1.0;  in BnDrv()
392       t7 = s+r-1.0;  in BnDrv()
396       t18 = t16*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
399       drv[2] = t16*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
404       t3 = r-1.0;  in BnDrv()
409       t9 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
424       t4 = r-1.0;  in BnDrv()
439       t4 = s+r-1.0;  in BnDrv()
443       t11 = r-2.0/15.0;  in BnDrv()
453       t1 = s+2.0/11.0*r-2.0/11.0;  in BnDrv()
455       t4 = s+r-1.0;  in BnDrv()
458       t11 = r-2.0/15.0;  in BnDrv()
469       t2 = r-1.0;  in BnDrv()
473       t7 = s+r-1.0;  in BnDrv()
476       t12 = r-2.0/15.0;  in BnDrv()
478       t18 = t16*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
479       drv[0] = t5*t10+t12*(6.0/11.0*s+6.0/55.0*r-6.0/55.0)*t10+t18;  in BnDrv()
480       drv[1] = t12*(2.0*s+6.0/11.0*r-6.0/11.0)*t10+t18;  in BnDrv()
481       drv[2] = t16*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
485       t3 = r-1.0;  in BnDrv()
488       t9 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
489       t11 = r-2.0/15.0;  in BnDrv()
501       t4 = r-1.0;  in BnDrv()
505       t11 = r-2.0/15.0;  in BnDrv()
515       t4 = s+r-1.0;  in BnDrv()
517       t11 = r*r;  in BnDrv()
518       t12 = t11-2.0/5.0*r+1.0/35.0;  in BnDrv()
521       drv[0] = (2.0*r-2.0/5.0)*(t2*t+12.0/7.0*t4*t2+6.0/7.0*t*t6+4.0/35.0*t6*t4  in BnDrv()
527       t2 = s+2.0/9.0*r-2.0/9.0;  in BnDrv()
529       t5 = s+r-1.0;  in BnDrv()
532       t10 = r*r;  in BnDrv()
533       t11 = t10-2.0/5.0*r+1.0/35.0;  in BnDrv()
536       t15 = t13*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
537       drv[0] = (2.0*r-2.0/5.0)*t2*t8+2.0/9.0*t12+t15;  in BnDrv()
539       drv[2] = t13*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
543       t3 = r-1.0;  in BnDrv()
546       t8 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
547       t10 = r*r;  in BnDrv()
548       t11 = t10-2.0/5.0*r+1.0/35.0;  in BnDrv()
550       drv[0] = (2.0*r-2.0/5.0)*t6*t8+t11*(2.0/3.0*s+r/6-1.0/6.0)*t8+2.0/3.0*t15  in BnDrv()
552       drv[1] = t11*(2.0*s+2.0/3.0*r-2.0/3.0)*t8+2.0/3.0*t15;  in BnDrv()
557       t4 = r-1.0;  in BnDrv()
559       t11 = r*r;  in BnDrv()
560       t12 = t11-2.0/5.0*r+1.0/35.0;  in BnDrv()
562       drv[0] = (2.0*r-2.0/5.0)*(t2*s+4.0/3.0*t4*t2+s*t6/2+t6*t4/21)+t12*(4.0/  in BnDrv()
568       t1 = r*r;  in BnDrv()
570       t4 = s+r-1.0;  in BnDrv()
572       t10 = t1*r-4.0/5.0*t1+6.0/35.0*r-4.0/455.0;  in BnDrv()
573       t12 = t10*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
574       drv[0] = (3.0*t1-8.0/5.0*r+6.0/35.0)*(t3+6.0/5.0*t*t4+3.0/10.0*t6)+t12;  in BnDrv()
576       drv[2] = t10*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
579       t1 = r*r;  in BnDrv()
580       t3 = s+2.0/7.0*r-2.0/7.0;  in BnDrv()
581       t5 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
582       t8 = t1*r-4.0/5.0*t1+6.0/35.0*r-4.0/455.0;  in BnDrv()
585       drv[0] = t3*(3.0*t1-8.0/5.0*r+6.0/35.0)*t5+2.0/7.0*t9+2.0/3.0*t10;  in BnDrv()
590       t1 = r*r;  in BnDrv()
592       t4 = r-1.0;  in BnDrv()
594       t10 = t1*r-4.0/5.0*t1+6.0/35.0*r-4.0/455.0;  in BnDrv()
595       drv[0] = (3.0*t1-8.0/5.0*r+6.0/35.0)*(t3+6.0/7.0*s*t4+t6/7)+t10*(6.0/7.0*  in BnDrv()
596 s+2.0/7.0*r-2.0/7.0);  in BnDrv()
597       drv[1] = t10*(2.0*s+6.0/7.0*r-6.0/7.0);  in BnDrv()
601       t1 = r*r;  in BnDrv()
602       t2 = t1*r;  in BnDrv()
604       drv[0] = (4.0*t2-4.0*t1+8.0/7.0*r-8.0/91.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/  in BnDrv()
605 3.0)+2.0/3.0*t6-8.0/9.0*t2+8.0/21.0*t1-16.0/273.0*r+2.0/819.0;  in BnDrv()
606       drv[1] = 2.0/3.0*t6-8.0/9.0*t2+8.0/21.0*t1-16.0/273.0*r+2.0/819.0;  in BnDrv()
607       drv[2] = t6-4.0/3.0*t2+4.0/7.0*t1-8.0/91.0*r+1.0/273.0;  in BnDrv()
610       t1 = r*r;  in BnDrv()
611       t2 = t1*r;  in BnDrv()
613       drv[0] = (4.0*t2-4.0*t1+8.0/7.0*r-8.0/91.0)*(s+2.0/5.0*r-2.0/5.0)+2.0/5.0  in BnDrv()
614 *t6-8.0/15.0*t2+8.0/35.0*t1-16.0/455.0*r+2.0/1365.0;  in BnDrv()
615       drv[1] = t6-4.0/3.0*t2+4.0/7.0*t1-8.0/91.0*r+1.0/273.0;  in BnDrv()
619       t1 = r*r;  in BnDrv()
621       drv[0] = 5.0*t2-8.0*t1*r+30.0/7.0*t1-80.0/91.0*r+5.0/91.0;  in BnDrv()
629       t4 = s+r-1.0;  in BnDrv()
648       t4 = s+r-1.0;  in BnDrv()
658       t14 = s+2.0/15.0*r-2.0/15.0;  in BnDrv()
675       t5 = s+r-1.0;  in BnDrv()
681       t15 = r-1.0;  in BnDrv()
687       drv[0] = (2.0/5.0*s+2.0/35.0*r-2.0/35.0)*t12+t22;  in BnDrv()
688       drv[1] = (2.0*s+2.0/5.0*r-2.0/5.0)*t12+t22;  in BnDrv()
693       t2 = r-1.0;  in BnDrv()
697       t8 = s+r-1.0;  in BnDrv()
710       t3 = r-1.0;  in BnDrv()
716       t10 = s+r-1.0;  in BnDrv()
722       t22 = t20*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
725       drv[2] = t20*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
731       t4 = r-1.0;  in BnDrv()
738       t12 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
757       t4 = r-1.0;  in BnDrv()
774       t4 = s+r-1.0;  in BnDrv()
779       t14 = r-2.0/17.0;  in BnDrv()
792       t1 = s+2.0/13.0*r-2.0/13.0;  in BnDrv()
796       t5 = s+r-1.0;  in BnDrv()
801       t14 = r-2.0/17.0;  in BnDrv()
813       t2 = r-1.0;  in BnDrv()
817       t8 = s+r-1.0;  in BnDrv()
820       t15 = r-2.0/17.0;  in BnDrv()
824       drv[0] = t5*t13+t15*(6.0/13.0*s+r/13-1.0/13.0)*t13+t22;  in BnDrv()
825       drv[1] = t15*(2.0*s+6.0/13.0*r-6.0/13.0)*t13+t22;  in BnDrv()
830       t3 = r-1.0;  in BnDrv()
834       t10 = s+r-1.0;  in BnDrv()
837       t15 = r-2.0/17.0;  in BnDrv()
840       t22 = t20*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
843       drv[2] = t20*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
849       t4 = r-1.0;  in BnDrv()
854       t12 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
855       t14 = r-2.0/17.0;  in BnDrv()
868       t4 = r-1.0;  in BnDrv()
873       t14 = r-2.0/17.0;  in BnDrv()
888       t5 = s+r-1.0;  in BnDrv()
892       t14 = r*r;  in BnDrv()
893       t15 = t14-6.0/17.0*r+3.0/136.0;  in BnDrv()
897       drv[0] = (2.0*r-6.0/17.0)*(t3+20.0/9.0*t4*t5+5.0/3.0*t2*t7+10.0/21.0*t*t9  in BnDrv()
903       t2 = s+2.0/11.0*r-2.0/11.0;  in BnDrv()
905       t6 = s+r-1.0;  in BnDrv()
908       t13 = r*r;  in BnDrv()
909       t14 = t13-6.0/17.0*r+3.0/136.0;  in BnDrv()
914       drv[0] = (2.0*r-6.0/17.0)*t2*t11+2.0/11.0*t15+t19;  in BnDrv()
920       t3 = r-1.0;  in BnDrv()
924       t9 = s+r-1.0;  in BnDrv()
927       t14 = r*r;  in BnDrv()
928       t15 = t14-6.0/17.0*r+3.0/136.0;  in BnDrv()
930       t21 = t19*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
931       drv[0] = (2.0*r-6.0/17.0)*t6*t12+t15*(6.0/11.0*s+6.0/55.0*r-6.0/55.0)*t12  in BnDrv()
933       drv[1] = t15*(2.0*s+6.0/11.0*r-6.0/11.0)*t12+t21;  in BnDrv()
934       drv[2] = t19*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
938       t4 = r-1.0;  in BnDrv()
941       t11 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
942       t13 = r*r;  in BnDrv()
943       t14 = t13-6.0/17.0*r+3.0/136.0;  in BnDrv()
946       drv[0] = (2.0*r-6.0/17.0)*t9*t11+t14*(12.0/11.0*t2+36.0/55.0*t15+4.0/55.0  in BnDrv()
955       t5 = r-1.0;  in BnDrv()
959       t14 = r*r;  in BnDrv()
960       t15 = t14-6.0/17.0*r+3.0/136.0;  in BnDrv()
963       drv[0] = (2.0*r-6.0/17.0)*(t3+20.0/11.0*t4*t5+12.0/11.0*t2*t7+8.0/33.0*s*  in BnDrv()
969       t1 = r*r;  in BnDrv()
971       t5 = s+r-1.0;  in BnDrv()
973       t13 = t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0;  in BnDrv()
976       drv[0] = (3.0*t1-24.0/17.0*r+9.0/68.0)*(t3*t+12.0/7.0*t3*t5+6.0/7.0*t*t7+  in BnDrv()
982       t1 = r*r;  in BnDrv()
983       t3 = s+2.0/9.0*r-2.0/9.0;  in BnDrv()
985       t6 = s+r-1.0;  in BnDrv()
988       t12 = t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0;  in BnDrv()
991       t16 = t14*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
992       drv[0] = t3*(3.0*t1-24.0/17.0*r+9.0/68.0)*t9+2.0/9.0*t13+t16;  in BnDrv()
994       drv[2] = t14*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
997       t1 = r*r;  in BnDrv()
999       t4 = r-1.0;  in BnDrv()
1002       t9 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
1003       t12 = t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0;  in BnDrv()
1005       drv[0] = (3.0*t1-24.0/17.0*r+9.0/68.0)*t7*t9+t12*(2.0/3.0*s+r/6-1.0/6.0)*  in BnDrv()
1007       drv[1] = t12*(2.0*s+2.0/3.0*r-2.0/3.0)*t9+2.0/3.0*t16;  in BnDrv()
1011       t1 = r*r;  in BnDrv()
1013       t5 = r-1.0;  in BnDrv()
1015       t13 = t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0;  in BnDrv()
1017       drv[0] = (3.0*t1-24.0/17.0*r+9.0/68.0)*(t3*s+4.0/3.0*t3*t5+s*t7/2+t7*t5  in BnDrv()
1023       t1 = r*r;  in BnDrv()
1024       t2 = t1*r;  in BnDrv()
1026       t5 = s+r-1.0;  in BnDrv()
1029       t11 = t10-20.0/17.0*t2+15.0/34.0*t1-r/17+1.0/476.0;  in BnDrv()
1030       t13 = t11*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
1031       drv[0] = (4.0*t2-60.0/17.0*t1+15.0/17.0*r-1.0/17.0)*(t4+6.0/5.0*t*t5+3.0/  in BnDrv()
1034       drv[2] = t11*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
1037       t1 = r*r;  in BnDrv()
1038       t2 = t1*r;  in BnDrv()
1039       t4 = s+2.0/7.0*r-2.0/7.0;  in BnDrv()
1040       t6 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
1042       t9 = t8-20.0/17.0*t2+15.0/34.0*t1-r/17+1.0/476.0;  in BnDrv()
1045       drv[0] = (4.0*t2-60.0/17.0*t1+15.0/17.0*r-1.0/17.0)*t4*t6+2.0/7.0*t10+2.0  in BnDrv()
1051       t1 = r*r;  in BnDrv()
1052       t2 = t1*r;  in BnDrv()
1054       t5 = r-1.0;  in BnDrv()
1057       t11 = t10-20.0/17.0*t2+15.0/34.0*t1-r/17+1.0/476.0;  in BnDrv()
1058       drv[0] = (4.0*t2-60.0/17.0*t1+15.0/17.0*r-1.0/17.0)*(t4+6.0/7.0*s*t5+t7/7  in BnDrv()
1059 )+t11*(6.0/7.0*s+2.0/7.0*r-2.0/7.0);  in BnDrv()
1060       drv[1] = t11*(2.0*s+6.0/7.0*r-6.0/7.0);  in BnDrv()
1064       t1 = r*r;  in BnDrv()
1066       t3 = t1*r;  in BnDrv()
1067       t7 = t2*r;  in BnDrv()
1068       drv[0] = (5.0*t2-120.0/17.0*t3+225.0/68.0*t1-10.0/17.0*r+15.0/476.0)*(t+  in BnDrv()
1069 2.0/3.0*s+2.0/3.0*r-2.0/3.0)+2.0/3.0*t7-20.0/17.0*t2+25.0/34.0*t3-10.0/51.0*t1+  in BnDrv()
1070 5.0/238.0*r-1.0/1547.0;  in BnDrv()
1071       drv[1] = 2.0/3.0*t7-20.0/17.0*t2+25.0/34.0*t3-10.0/51.0*t1+5.0/238.0*r  in BnDrv()
1073       drv[2] = t7-30.0/17.0*t2+75.0/68.0*t3-5.0/17.0*t1+15.0/476.0*r-3.0/3094.0  in BnDrv()
1077       t1 = r*r;  in BnDrv()
1079       t3 = t1*r;  in BnDrv()
1080       t7 = t2*r;  in BnDrv()
1081       drv[0] = (5.0*t2-120.0/17.0*t3+225.0/68.0*t1-10.0/17.0*r+15.0/476.0)*(s+  in BnDrv()
1082 2.0/5.0*r-2.0/5.0)+2.0/5.0*t7-12.0/17.0*t2+15.0/34.0*t3-2.0/17.0*t1+3.0/238.0*r  in BnDrv()
1084       drv[1] = t7-30.0/17.0*t2+75.0/68.0*t3-5.0/17.0*t1+15.0/476.0*r-3.0/3094.0  in BnDrv()
1089       t1 = r*r;  in BnDrv()
1091       drv[0] = 6.0*t2*r-210.0/17.0*t2+315.0/34.0*t1*r-105.0/34.0*t1+15.0/34.0*r  in BnDrv()
1100       t5 = s+r-1.0;  in BnDrv()
1121       t5 = s+r-1.0;  in BnDrv()
1133       t17 = s+2.0/17.0*r-2.0/17.0;  in BnDrv()
1150       t5 = s+r-1.0;  in BnDrv()
1158       t18 = r-1.0;  in BnDrv()
1166       drv[0] = (6.0/17.0*s+3.0/68.0*r-3.0/68.0)*t15+t26;  in BnDrv()
1167       drv[1] = (2.0*s+6.0/17.0*r-6.0/17.0)*t15+t26;  in BnDrv()
1173       t2 = r-1.0;  in BnDrv()
1179       t9 = s+r-1.0;  in BnDrv()
1195       t3 = r-1.0;  in BnDrv()
1201       t11 = s+r-1.0;  in BnDrv()
1217       t4 = r-1.0;  in BnDrv()
1225       t13 = s+r-1.0;  in BnDrv()
1230       t26 = t24*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
1235       drv[2] = t24*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
1241       t4 = r-1.0;  in BnDrv()
1251       t15 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
1272       t5 = r-1.0;  in BnDrv()
1291       t5 = s+r-1.0;  in BnDrv()
1297       t17 = r-2.0/19.0;  in BnDrv()
1311       t1 = s+2.0/15.0*r-2.0/15.0;  in BnDrv()
1314       t5 = s+r-1.0;  in BnDrv()
1321       t17 = r-2.0/19.0;  in BnDrv()
1336       t2 = r-1.0;  in BnDrv()
1342       t9 = s+r-1.0;  in BnDrv()
1347       t18 = r-2.0/19.0;  in BnDrv()
1352       drv[0] = t5*t16+t18*(2.0/5.0*s+2.0/35.0*r-2.0/35.0)*t16+t26;  in BnDrv()
1353       drv[1] = t18*(2.0*s+2.0/5.0*r-2.0/5.0)*t16+t26;  in BnDrv()
1358       t3 = r-1.0;  in BnDrv()
1362       t11 = s+r-1.0;  in BnDrv()
1365       t18 = r-2.0/19.0;  in BnDrv()
1378       t4 = r-1.0;  in BnDrv()
1384       t13 = s+r-1.0;  in BnDrv()
1387       t18 = r-2.0/19.0;  in BnDrv()
1391       t26 = t24*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
1395       drv[2] = t24*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
1400       t4 = r-1.0;  in BnDrv()
1407       t15 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
1408       t17 = r-2.0/19.0;  in BnDrv()
1423       t5 = r-1.0;  in BnDrv()
1429       t17 = r-2.0/19.0;  in BnDrv()
1444       t5 = s+r-1.0;  in BnDrv()
1449       t17 = r*r;  in BnDrv()
1450       t18 = t17-6.0/19.0*r+1.0/57.0;  in BnDrv()
1456       drv[0] = (2.0*r-6.0/19.0)*(t3*t+30.0/11.0*t3*t5+30.0/11.0*t7*t8+40.0/33.0  in BnDrv()
1463       t2 = s+2.0/13.0*r-2.0/13.0;  in BnDrv()
1467       t7 = s+r-1.0;  in BnDrv()
1472       t16 = r*r;  in BnDrv()
1473       t17 = t16-6.0/19.0*r+1.0/57.0;  in BnDrv()
1479       drv[0] = t2*(2.0*r-6.0/19.0)*t14+2.0/13.0*t18+t23;  in BnDrv()
1485       t3 = r-1.0;  in BnDrv()
1489       t10 = s+r-1.0;  in BnDrv()
1492       t17 = r*r;  in BnDrv()
1493       t18 = t17-6.0/19.0*r+1.0/57.0;  in BnDrv()
1497       drv[0] = (2.0*r-6.0/19.0)*t6*t15+t18*(6.0/13.0*s+r/13-1.0/13.0)*t15+t25;  in BnDrv()
1498       drv[1] = t18*(2.0*s+6.0/13.0*r-6.0/13.0)*t15+t25;  in BnDrv()
1503       t4 = r-1.0;  in BnDrv()
1507       t12 = s+r-1.0;  in BnDrv()
1510       t17 = r*r;  in BnDrv()
1511       t18 = t17-6.0/19.0*r+1.0/57.0;  in BnDrv()
1514       t25 = t23*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
1515       drv[0] = (2.0*r-6.0/19.0)*t9*t15+t18*(12.0/13.0*t2+6.0/13.0*t19+6.0/143.0  in BnDrv()
1518       drv[2] = t23*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
1524       t5 = r-1.0;  in BnDrv()
1529       t14 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
1530       t16 = r*r;  in BnDrv()
1531       t17 = t16-6.0/19.0*r+1.0/57.0;  in BnDrv()
1535       drv[0] = (2.0*r-6.0/19.0)*t12*t14+t17*(20.0/13.0*t4+20.0/13.0*t18+60.0/  in BnDrv()
1544       t5 = r-1.0;  in BnDrv()
1549       t17 = r*r;  in BnDrv()
1550       t18 = t17-6.0/19.0*r+1.0/57.0;  in BnDrv()
1554       drv[0] = (2.0*r-6.0/19.0)*(t3*s+30.0/13.0*t3*t5+25.0/13.0*t7*t8+100.0/  in BnDrv()
1562       t1 = r*r;  in BnDrv()
1566       t6 = s+r-1.0;  in BnDrv()
1570       t16 = t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0;  in BnDrv()
1574       drv[0] = (3.0*t1-24.0/19.0*r+2.0/19.0)*(t4+20.0/9.0*t5*t6+5.0/3.0*t3*t8+  in BnDrv()
1580       t1 = r*r;  in BnDrv()
1581       t3 = s+2.0/11.0*r-2.0/11.0;  in BnDrv()
1583       t7 = s+r-1.0;  in BnDrv()
1586       t15 = t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0;  in BnDrv()
1591       drv[0] = t3*(3.0*t1-24.0/19.0*r+2.0/19.0)*t12+2.0/11.0*t16+t20;  in BnDrv()
1596       t1 = r*r;  in BnDrv()
1598       t4 = r-1.0;  in BnDrv()
1602       t10 = s+r-1.0;  in BnDrv()
1605       t16 = t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0;  in BnDrv()
1607       t22 = t20*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
1608       drv[0] = (3.0*t1-24.0/19.0*r+2.0/19.0)*t7*t13+t16*(6.0/11.0*s+6.0/55.0*r  in BnDrv()
1610       drv[1] = t16*(2.0*s+6.0/11.0*r-6.0/11.0)*t13+t22;  in BnDrv()
1611       drv[2] = t20*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
1614       t1 = r*r;  in BnDrv()
1616       t5 = r-1.0;  in BnDrv()
1619       t12 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
1620       t15 = t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0;  in BnDrv()
1623       drv[0] = (3.0*t1-24.0/19.0*r+2.0/19.0)*t10*t12+t15*(12.0/11.0*t3+36.0/  in BnDrv()
1629       t1 = r*r;  in BnDrv()
1633       t6 = r-1.0;  in BnDrv()
1637       t16 = t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0;  in BnDrv()
1640       drv[0] = (3.0*t1-24.0/19.0*r+2.0/19.0)*(t4+20.0/11.0*t5*t6+12.0/11.0*t3*  in BnDrv()
1647       t1 = r*r;  in BnDrv()
1648       t2 = t1*r;  in BnDrv()
1650       t6 = s+r-1.0;  in BnDrv()
1653       t14 = t13-20.0/19.0*t2+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0;  in BnDrv()
1656       drv[0] = (4.0*t2-60.0/19.0*t1+40.0/57.0*r-40.0/969.0)*(t4*t+12.0/7.0*t4*  in BnDrv()
1662       t1 = r*r;  in BnDrv()
1663       t2 = t1*r;  in BnDrv()
1664       t4 = s+2.0/9.0*r-2.0/9.0;  in BnDrv()
1666       t7 = s+r-1.0;  in BnDrv()
1670       t13 = t12-20.0/19.0*t2+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0;  in BnDrv()
1673       t17 = t15*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
1674       drv[0] = (4.0*t2-60.0/19.0*t1+40.0/57.0*r-40.0/969.0)*t4*t10+2.0/9.0*t14+  in BnDrv()
1677       drv[2] = t15*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
1680       t1 = r*r;  in BnDrv()
1681       t2 = t1*r;  in BnDrv()
1683       t5 = r-1.0;  in BnDrv()
1686       t10 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
1688       t13 = t12-20.0/19.0*t2+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0;  in BnDrv()
1690       drv[0] = (4.0*t2-60.0/19.0*t1+40.0/57.0*r-40.0/969.0)*t8*t10+t13*(2.0/3.0  in BnDrv()
1691 *s+r/6-1.0/6.0)*t10+2.0/3.0*t17;  in BnDrv()
1692       drv[1] = t13*(2.0*s+2.0/3.0*r-2.0/3.0)*t10+2.0/3.0*t17;  in BnDrv()
1696       t1 = r*r;  in BnDrv()
1697       t2 = t1*r;  in BnDrv()
1699       t6 = r-1.0;  in BnDrv()
1702       t14 = t13-20.0/19.0*t2+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0;  in BnDrv()
1704       drv[0] = (4.0*t2-60.0/19.0*t1+40.0/57.0*r-40.0/969.0)*(t4*s+4.0/3.0*t4*t6  in BnDrv()
1710       t1 = r*r;  in BnDrv()
1712       t3 = t1*r;  in BnDrv()
1714       t6 = s+r-1.0;  in BnDrv()
1716       t12 = t2*r-30.0/19.0*t2+50.0/57.0*t3-200.0/969.0*t1+25.0/1292.0*r-1.0/  in BnDrv()
1718       t14 = t12*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
1719       drv[0] = (5.0*t2-120.0/19.0*t3+50.0/19.0*t1-400.0/969.0*r+25.0/1292.0)*(  in BnDrv()
1722       drv[2] = t12*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
1725       t1 = r*r;  in BnDrv()
1727       t3 = t1*r;  in BnDrv()
1728       t5 = s+2.0/7.0*r-2.0/7.0;  in BnDrv()
1729       t7 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
1730       t10 = t2*r-30.0/19.0*t2+50.0/57.0*t3-200.0/969.0*t1+25.0/1292.0*r-1.0/  in BnDrv()
1734       drv[0] = (5.0*t2-120.0/19.0*t3+50.0/19.0*t1-400.0/969.0*r+25.0/1292.0)*t5  in BnDrv()
1740       t1 = r*r;  in BnDrv()
1742       t3 = t1*r;  in BnDrv()
1744       t6 = r-1.0;  in BnDrv()
1746       t12 = t2*r-30.0/19.0*t2+50.0/57.0*t3-200.0/969.0*t1+25.0/1292.0*r-1.0/  in BnDrv()
1748       drv[0] = (5.0*t2-120.0/19.0*t3+50.0/19.0*t1-400.0/969.0*r+25.0/1292.0)*(  in BnDrv()
1749 t5+6.0/7.0*s*t6+t8/7)+t12*(6.0/7.0*s+2.0/7.0*r-2.0/7.0);  in BnDrv()
1750       drv[1] = t12*(2.0*s+6.0/7.0*r-6.0/7.0);  in BnDrv()
1754       t1 = r*r;  in BnDrv()
1756       t3 = t2*r;  in BnDrv()
1757       t4 = t1*r;  in BnDrv()
1759       drv[0] = (6.0*t3-210.0/19.0*t2+140.0/19.0*t4-700.0/323.0*t1+175.0/646.0*r  in BnDrv()
1760 -7.0/646.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0)+2.0/3.0*t8-28.0/19.0*t3+70.0/57.0*  in BnDrv()
1761 t2-1400.0/2907.0*t4+175.0/1938.0*t1-7.0/969.0*r+1.0/5814.0;  in BnDrv()
1763 1938.0*t1-7.0/969.0*r+1.0/5814.0;  in BnDrv()
1765 646.0*r+1.0/3876.0;  in BnDrv()
1768       t1 = r*r;  in BnDrv()
1770       t3 = t2*r;  in BnDrv()
1771       t4 = t1*r;  in BnDrv()
1773       drv[0] = (6.0*t3-210.0/19.0*t2+140.0/19.0*t4-700.0/323.0*t1+175.0/646.0*r  in BnDrv()
1774 -7.0/646.0)*(s+2.0/5.0*r-2.0/5.0)+2.0/5.0*t8-84.0/95.0*t3+14.0/19.0*t2-280.0/  in BnDrv()
1775 969.0*t4+35.0/646.0*t1-7.0/1615.0*r+1.0/9690.0;  in BnDrv()
1777 646.0*r+1.0/3876.0;  in BnDrv()
1781       t1 = r*r;  in BnDrv()
1783       drv[0] = 7.0*t2*t1-336.0/19.0*t2*r+980.0/57.0*t2-7840.0/969.0*t1*r+1225.0  in BnDrv()
1784 /646.0*t1-196.0/969.0*r+7.0/969.0;  in BnDrv()
1793       t6 = s+r-1.0;  in BnDrv()
1817       t6 = s+r-1.0;  in BnDrv()
1831       t20 = s+2.0/19.0*r-2.0/19.0;  in BnDrv()
1850       t6 = s+r-1.0;  in BnDrv()
1859       t21 = r-1.0;  in BnDrv()
1868       drv[0] = (6.0/19.0*s+2.0/57.0*r-2.0/57.0)*t18+t30;  in BnDrv()
1869       drv[1] = (2.0*s+6.0/19.0*r-6.0/19.0)*t18+t30;  in BnDrv()
1875       t2 = r-1.0;  in BnDrv()
1880       t9 = s+r-1.0;  in BnDrv()
1901       t3 = r-1.0;  in BnDrv()
1909       t12 = s+r-1.0;  in BnDrv()
1928       t4 = r-1.0;  in BnDrv()
1936       t14 = s+r-1.0;  in BnDrv()
1953       t4 = r-1.0;  in BnDrv()
1964       t16 = s+r-1.0;  in BnDrv()
1969       t30 = t28*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
1974       drv[2] = t28*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
1981       t5 = r-1.0;  in BnDrv()
1993       t18 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
2018       t6 = r-1.0;  in BnDrv()
2040       t6 = s+r-1.0;  in BnDrv()
2047       t20 = r-2.0/21.0;  in BnDrv()
2062       t1 = s+2.0/17.0*r-2.0/17.0;  in BnDrv()
2066       t6 = s+r-1.0;  in BnDrv()
2074       t20 = r-2.0/21.0;  in BnDrv()
2090       t2 = r-1.0;  in BnDrv()
2095       t9 = s+r-1.0;  in BnDrv()
2102       t21 = r-2.0/21.0;  in BnDrv()
2109       drv[0] = t5*t19+t21*(6.0/17.0*s+3.0/68.0*r-3.0/68.0)*t19+t30;  in BnDrv()
2110       drv[1] = t21*(2.0*s+6.0/17.0*r-6.0/17.0)*t19+t30;  in BnDrv()
2116       t3 = r-1.0;  in BnDrv()
2122       t12 = s+r-1.0;  in BnDrv()
2127       t21 = r-2.0/21.0;  in BnDrv()
2141       t4 = r-1.0;  in BnDrv()
2147       t14 = s+r-1.0;  in BnDrv()
2150       t21 = r-2.0/21.0;  in BnDrv()
2164       t4 = r-1.0;  in BnDrv()
2172       t16 = s+r-1.0;  in BnDrv()
2175       t21 = r-2.0/21.0;  in BnDrv()
2180       t30 = t28*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
2185       drv[2] = t28*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
2191       t5 = r-1.0;  in BnDrv()
2199       t18 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
2200       t20 = r-2.0/21.0;  in BnDrv()
2217       t6 = r-1.0;  in BnDrv()
2224       t20 = r-2.0/21.0;  in BnDrv()
2242       t6 = s+r-1.0;  in BnDrv()
2248       t20 = r*r;  in BnDrv()
2249       t21 = t20-2.0/7.0*r+1.0/70.0;  in BnDrv()
2256       drv[0] = (2.0*r-2.0/7.0)*(t3*t2+42.0/13.0*t5*t6+105.0/26.0*t3*t8+350.0/  in BnDrv()
2263       t2 = s+2.0/15.0*r-2.0/15.0;  in BnDrv()
2266       t7 = s+r-1.0;  in BnDrv()
2273       t19 = r*r;  in BnDrv()
2274       t20 = t19-2.0/7.0*r+1.0/70.0;  in BnDrv()
2282       drv[0] = (2.0*r-2.0/7.0)*t2*t17+2.0/15.0*t21+t27;  in BnDrv()
2289       t3 = r-1.0;  in BnDrv()
2295       t11 = s+r-1.0;  in BnDrv()
2300       t20 = r*r;  in BnDrv()
2301       t21 = t20-2.0/7.0*r+1.0/70.0;  in BnDrv()
2306       drv[0] = (2.0*r-2.0/7.0)*t6*t18+t21*(2.0/5.0*s+2.0/35.0*r-2.0/35.0)*t18+  in BnDrv()
2308       drv[1] = t21*(2.0*s+2.0/5.0*r-2.0/5.0)*t18+t29;  in BnDrv()
2313       t4 = r-1.0;  in BnDrv()
2317       t13 = s+r-1.0;  in BnDrv()
2320       t20 = r*r;  in BnDrv()
2321       t21 = t20-2.0/7.0*r+1.0/70.0;  in BnDrv()
2326       drv[0] = (2.0*r-2.0/7.0)*t9*t18+t21*(4.0/5.0*t2+12.0/35.0*t22+12.0/455.0*  in BnDrv()
2335       t5 = r-1.0;  in BnDrv()
2341       t15 = s+r-1.0;  in BnDrv()
2344       t20 = r*r;  in BnDrv()
2345       t21 = t20-2.0/7.0*r+1.0/70.0;  in BnDrv()
2349       t29 = t27*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
2350       drv[0] = (2.0*r-2.0/7.0)*t12*t18+t21*(4.0/3.0*t4+8.0/7.0*t22+24.0/91.0*  in BnDrv()
2353       drv[2] = t27*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
2358       t5 = r-1.0;  in BnDrv()
2365       t17 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
2366       t19 = r*r;  in BnDrv()
2367       t20 = t19-2.0/7.0*r+1.0/70.0;  in BnDrv()
2372       drv[0] = (2.0*r-2.0/7.0)*t15*t17+t20*(2.0*t3+20.0/7.0*t21+120.0/91.0*t22+  in BnDrv()
2382       t6 = r-1.0;  in BnDrv()
2388       t20 = r*r;  in BnDrv()
2389       t21 = t20-2.0/7.0*r+1.0/70.0;  in BnDrv()
2394       drv[0] = (2.0*r-2.0/7.0)*(t3*t2+14.0/5.0*t5*t6+3.0*t3*t8+20.0/13.0*t10*  in BnDrv()
2402       t1 = r*r;  in BnDrv()
2405       t6 = s+r-1.0;  in BnDrv()
2410       t19 = t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0;  in BnDrv()
2416       drv[0] = (3.0*t1-8.0/7.0*r+3.0/35.0)*(t4*t+30.0/11.0*t4*t6+30.0/11.0*t8*  in BnDrv()
2423       t1 = r*r;  in BnDrv()
2424       t3 = s+2.0/13.0*r-2.0/13.0;  in BnDrv()
2428       t8 = s+r-1.0;  in BnDrv()
2433       t18 = t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0;  in BnDrv()
2439       drv[0] = t3*(3.0*t1-8.0/7.0*r+3.0/35.0)*t15+2.0/13.0*t19+t24;  in BnDrv()
2444       t1 = r*r;  in BnDrv()
2446       t4 = r-1.0;  in BnDrv()
2450       t11 = s+r-1.0;  in BnDrv()
2453       t19 = t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0;  in BnDrv()
2457       drv[0] = (3.0*t1-8.0/7.0*r+3.0/35.0)*t7*t16+t19*(6.0/13.0*s+r/13-1.0/13.0  in BnDrv()
2459       drv[1] = t19*(2.0*s+6.0/13.0*r-6.0/13.0)*t16+t26;  in BnDrv()
2463       t1 = r*r;  in BnDrv()
2465       t5 = r-1.0;  in BnDrv()
2469       t13 = s+r-1.0;  in BnDrv()
2472       t19 = t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0;  in BnDrv()
2475       t26 = t24*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
2476       drv[0] = (3.0*t1-8.0/7.0*r+3.0/35.0)*t10*t16+t19*(12.0/13.0*t3+6.0/13.0*  in BnDrv()
2479       drv[2] = t24*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
2482       t1 = r*r;  in BnDrv()
2486       t6 = r-1.0;  in BnDrv()
2491       t15 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
2492       t18 = t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0;  in BnDrv()
2496       drv[0] = (3.0*t1-8.0/7.0*r+3.0/35.0)*t13*t15+t18*(20.0/13.0*t5+20.0/13.0*  in BnDrv()
2503       t1 = r*r;  in BnDrv()
2506       t6 = r-1.0;  in BnDrv()
2511       t19 = t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0;  in BnDrv()
2515       drv[0] = (3.0*t1-8.0/7.0*r+3.0/35.0)*(t4*s+30.0/13.0*t4*t6+25.0/13.0*t8*  in BnDrv()
2523       t1 = r*r;  in BnDrv()
2524       t2 = t1*r;  in BnDrv()
2528       t7 = s+r-1.0;  in BnDrv()
2533       t17 = t16-20.0/21.0*t2+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0;  in BnDrv()
2537       drv[0] = (4.0*t2-20.0/7.0*t1+4.0/7.0*r-4.0/133.0)*(t5+20.0/9.0*t6*t7+5.0/  in BnDrv()
2543       t1 = r*r;  in BnDrv()
2544       t2 = t1*r;  in BnDrv()
2545       t4 = s+2.0/11.0*r-2.0/11.0;  in BnDrv()
2547       t8 = s+r-1.0;  in BnDrv()
2551       t16 = t15-20.0/21.0*t2+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0;  in BnDrv()
2556       drv[0] = (4.0*t2-20.0/7.0*t1+4.0/7.0*r-4.0/133.0)*t4*t13+2.0/11.0*t17+t21  in BnDrv()
2562       t1 = r*r;  in BnDrv()
2563       t2 = t1*r;  in BnDrv()
2565       t5 = r-1.0;  in BnDrv()
2569       t11 = s+r-1.0;  in BnDrv()
2573       t17 = t16-20.0/21.0*t2+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0;  in BnDrv()
2575       t23 = t21*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
2576       drv[0] = (4.0*t2-20.0/7.0*t1+4.0/7.0*r-4.0/133.0)*t8*t14+t17*(6.0/11.0*s+  in BnDrv()
2577 6.0/55.0*r-6.0/55.0)*t14+t23;  in BnDrv()
2578       drv[1] = t17*(2.0*s+6.0/11.0*r-6.0/11.0)*t14+t23;  in BnDrv()
2579       drv[2] = t21*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
2582       t1 = r*r;  in BnDrv()
2583       t2 = t1*r;  in BnDrv()
2585       t6 = r-1.0;  in BnDrv()
2588       t13 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
2590       t16 = t15-20.0/21.0*t2+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0;  in BnDrv()
2593       drv[0] = (4.0*t2-20.0/7.0*t1+4.0/7.0*r-4.0/133.0)*t11*t13+t16*(12.0/11.0*  in BnDrv()
2599       t1 = r*r;  in BnDrv()
2600       t2 = t1*r;  in BnDrv()
2604       t7 = r-1.0;  in BnDrv()
2609       t17 = t16-20.0/21.0*t2+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0;  in BnDrv()
2612       drv[0] = (4.0*t2-20.0/7.0*t1+4.0/7.0*r-4.0/133.0)*(t5+20.0/11.0*t6*t7+  in BnDrv()
2619       t1 = r*r;  in BnDrv()
2621       t3 = t1*r;  in BnDrv()
2623       t7 = s+r-1.0;  in BnDrv()
2625       t15 = t2*r-10.0/7.0*t2+5.0/7.0*t3-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0;  in BnDrv()
2628       drv[0] = (5.0*t2-40.0/7.0*t3+15.0/7.0*t1-40.0/133.0*r+5.0/399.0)*(t5*t+  in BnDrv()
2634       t1 = r*r;  in BnDrv()
2636       t3 = t1*r;  in BnDrv()
2637       t5 = s+2.0/9.0*r-2.0/9.0;  in BnDrv()
2639       t8 = s+r-1.0;  in BnDrv()
2642       t14 = t2*r-10.0/7.0*t2+5.0/7.0*t3-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0;  in BnDrv()
2645       t18 = t16*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
2646       drv[0] = (5.0*t2-40.0/7.0*t3+15.0/7.0*t1-40.0/133.0*r+5.0/399.0)*t5*t11+  in BnDrv()
2649       drv[2] = t16*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
2652       t1 = r*r;  in BnDrv()
2654       t3 = t1*r;  in BnDrv()
2656       t6 = r-1.0;  in BnDrv()
2659       t11 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
2660       t14 = t2*r-10.0/7.0*t2+5.0/7.0*t3-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0;  in BnDrv()
2662       drv[0] = (5.0*t2-40.0/7.0*t3+15.0/7.0*t1-40.0/133.0*r+5.0/399.0)*t9*t11+  in BnDrv()
2663 t14*(2.0/3.0*s+r/6-1.0/6.0)*t11+2.0/3.0*t18;  in BnDrv()
2664       drv[1] = t14*(2.0*s+2.0/3.0*r-2.0/3.0)*t11+2.0/3.0*t18;  in BnDrv()
2668       t1 = r*r;  in BnDrv()
2670       t3 = t1*r;  in BnDrv()
2672       t7 = r-1.0;  in BnDrv()
2674       t15 = t2*r-10.0/7.0*t2+5.0/7.0*t3-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0;  in BnDrv()
2676       drv[0] = (5.0*t2-40.0/7.0*t3+15.0/7.0*t1-40.0/133.0*r+5.0/399.0)*(t5*s+  in BnDrv()
2682       t1 = r*r;  in BnDrv()
2684       t3 = t2*r;  in BnDrv()
2685       t4 = t1*r;  in BnDrv()
2687       t7 = s+r-1.0;  in BnDrv()
2689       t13 = t2*t1-2.0*t3+3.0/2.0*t2-10.0/19.0*t4+5.0/57.0*t1-2.0/323.0*r+1.0/  in BnDrv()
2691       t15 = t13*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
2692       drv[0] = (6.0*t3-10.0*t2+6.0*t4-30.0/19.0*t1+10.0/57.0*r-2.0/323.0)*(t6+  in BnDrv()
2695       drv[2] = t13*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
2698       t1 = r*r;  in BnDrv()
2700       t3 = t2*r;  in BnDrv()
2701       t4 = t1*r;  in BnDrv()
2702       t6 = s+2.0/7.0*r-2.0/7.0;  in BnDrv()
2703       t8 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
2704       t11 = t2*t1-2.0*t3+3.0/2.0*t2-10.0/19.0*t4+5.0/57.0*t1-2.0/323.0*r+1.0/  in BnDrv()
2708       drv[0] = (6.0*t3-10.0*t2+6.0*t4-30.0/19.0*t1+10.0/57.0*r-2.0/323.0)*t6*t8  in BnDrv()
2714       t1 = r*r;  in BnDrv()
2716       t3 = t2*r;  in BnDrv()
2717       t4 = t1*r;  in BnDrv()
2719       t7 = r-1.0;  in BnDrv()
2721       t13 = t2*t1-2.0*t3+3.0/2.0*t2-10.0/19.0*t4+5.0/57.0*t1-2.0/323.0*r+1.0/  in BnDrv()
2723       drv[0] = (6.0*t3-10.0*t2+6.0*t4-30.0/19.0*t1+10.0/57.0*r-2.0/323.0)*(t6+  in BnDrv()
2724 6.0/7.0*s*t7+t9/7)+t13*(6.0/7.0*s+2.0/7.0*r-2.0/7.0);  in BnDrv()
2725       drv[1] = t13*(2.0*s+6.0/7.0*r-6.0/7.0);  in BnDrv()
2729       t1 = r*r;  in BnDrv()
2732       t4 = t2*r;  in BnDrv()
2733       t5 = t1*r;  in BnDrv()
2735       drv[0] = (7.0*t3-16.0*t4+14.0*t2-112.0/19.0*t5+70.0/57.0*t1-112.0/969.0*r  in BnDrv()
2736 +7.0/1938.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0)+2.0/3.0*t9-16.0/9.0*t3+28.0/15.0*  in BnDrv()
2737 t4-56.0/57.0*t2+140.0/513.0*t5-112.0/2907.0*t1+7.0/2907.0*r-2.0/43605.0;  in BnDrv()
2739 -112.0/2907.0*t1+7.0/2907.0*r-2.0/43605.0;  in BnDrv()
2741 t1+7.0/1938.0*r-1.0/14535.0;  in BnDrv()
2744       t1 = r*r;  in BnDrv()
2747       t4 = t2*r;  in BnDrv()
2748       t5 = t1*r;  in BnDrv()
2750       drv[0] = (7.0*t3-16.0*t4+14.0*t2-112.0/19.0*t5+70.0/57.0*t1-112.0/969.0*r  in BnDrv()
2751 +7.0/1938.0)*(s+2.0/5.0*r-2.0/5.0)+2.0/5.0*t9-16.0/15.0*t3+28.0/25.0*t4-56.0/  in BnDrv()
2752 95.0*t2+28.0/171.0*t5-112.0/4845.0*t1+7.0/4845.0*r-2.0/72675.0;  in BnDrv()
2754 t1+7.0/1938.0*r-1.0/14535.0;  in BnDrv()
2758       t1 = r*r;  in BnDrv()
2759       t2 = t1*r;  in BnDrv()
2761       drv[0] = 8.0*t3*t2-24.0*t3*t1+144.0/5.0*t3*r-336.0/19.0*t3+112.0/19.0*t2  in BnDrv()
2762 -336.0/323.0*t1+28.0/323.0*r-4.0/1615.0;  in BnDrv()
2771       t6 = s+r-1.0;  in BnDrv()
2797       t6 = s+r-1.0;  in BnDrv()
2814       t23 = s+2.0/21.0*r-2.0/21.0;  in BnDrv()
2835       t7 = s+r-1.0;  in BnDrv()
2845       t24 = r-1.0;  in BnDrv()
2855       drv[0] = (2.0/7.0*s+r/35-1.0/35.0)*t21+t34;  in BnDrv()
2856       drv[1] = (2.0*s+2.0/7.0*r-2.0/7.0)*t21+t34;  in BnDrv()
2862       t2 = r-1.0;  in BnDrv()
2868       t10 = s+r-1.0;  in BnDrv()
2891       t3 = r-1.0;  in BnDrv()
2898       t12 = s+r-1.0;  in BnDrv()
2922       t4 = r-1.0;  in BnDrv()
2932       t15 = s+r-1.0;  in BnDrv()
2952       t4 = r-1.0;  in BnDrv()
2963       t17 = s+r-1.0;  in BnDrv()
2981       t5 = r-1.0;  in BnDrv()
2994       t19 = s+r-1.0;  in BnDrv()
2999       t34 = t32*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
3004       drv[2] = t32*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
3012       t6 = r-1.0;  in BnDrv()
3026       t21 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
3052       t6 = r-1.0;  in BnDrv()
3076       t6 = s+r-1.0;  in BnDrv()
3086       t23 = r-2.0/23.0;  in BnDrv()
3103       t1 = s+2.0/19.0*r-2.0/19.0;  in BnDrv()
3108       t7 = s+r-1.0;  in BnDrv()
3117       t23 = r-2.0/23.0;  in BnDrv()
3134       t2 = r-1.0;  in BnDrv()
3140       t10 = s+r-1.0;  in BnDrv()
3148       t24 = r-2.0/23.0;  in BnDrv()
3156       drv[0] = t5*t22+t24*(6.0/19.0*s+2.0/57.0*r-2.0/57.0)*t22+t34;  in BnDrv()
3157       drv[1] = t24*(2.0*s+6.0/19.0*r-6.0/19.0)*t22+t34;  in BnDrv()
3163       t3 = r-1.0;  in BnDrv()
3168       t12 = s+r-1.0;  in BnDrv()
3175       t24 = r-2.0/23.0;  in BnDrv()
3192       t4 = r-1.0;  in BnDrv()
3200       t15 = s+r-1.0;  in BnDrv()
3205       t24 = r-2.0/23.0;  in BnDrv()
3220       t4 = r-1.0;  in BnDrv()
3228       t17 = s+r-1.0;  in BnDrv()
3231       t24 = r-2.0/23.0;  in BnDrv()
3248       t5 = r-1.0;  in BnDrv()
3257       t19 = s+r-1.0;  in BnDrv()
3260       t24 = r-2.0/23.0;  in BnDrv()
3266       t34 = t32*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
3271       drv[2] = t32*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
3278       t6 = r-1.0;  in BnDrv()
3287       t21 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
3288       t23 = r-2.0/23.0;  in BnDrv()
3307       t6 = r-1.0;  in BnDrv()
3317       t23 = r-2.0/23.0;  in BnDrv()
3337       t7 = s+r-1.0;  in BnDrv()
3344       t23 = r*r;  in BnDrv()
3345       t24 = t23-6.0/23.0*r+3.0/253.0;  in BnDrv()
3353       drv[0] = (2.0*r-6.0/23.0)*(t4*t3+56.0/15.0*t6*t7+28.0/5.0*t9*t10+56.0/  in BnDrv()
3361       t2 = s+2.0/17.0*r-2.0/17.0;  in BnDrv()
3365       t8 = s+r-1.0;  in BnDrv()
3373       t22 = r*r;  in BnDrv()
3374       t23 = t22-6.0/23.0*r+3.0/253.0;  in BnDrv()
3383       drv[0] = t2*(2.0*r-6.0/23.0)*t20+2.0/17.0*t24+t31;  in BnDrv()
3390       t3 = r-1.0;  in BnDrv()
3395       t11 = s+r-1.0;  in BnDrv()
3402       t23 = r*r;  in BnDrv()
3403       t24 = t23-6.0/23.0*r+3.0/253.0;  in BnDrv()
3410       drv[0] = (2.0*r-6.0/23.0)*t6*t21+t24*(6.0/17.0*s+3.0/68.0*r-3.0/68.0)*t21  in BnDrv()
3412       drv[1] = t24*(2.0*s+6.0/17.0*r-6.0/17.0)*t21+t33;  in BnDrv()
3418       t4 = r-1.0;  in BnDrv()
3424       t14 = s+r-1.0;  in BnDrv()
3429       t23 = r*r;  in BnDrv()
3430       t24 = t23-6.0/23.0*r+3.0/253.0;  in BnDrv()
3436       drv[0] = (2.0*r-6.0/23.0)*t9*t21+t24*(12.0/17.0*t2+9.0/34.0*t25+3.0/170.0  in BnDrv()
3445       t5 = r-1.0;  in BnDrv()
3451       t16 = s+r-1.0;  in BnDrv()
3454       t23 = r*r;  in BnDrv()
3455       t24 = t23-6.0/23.0*r+3.0/253.0;  in BnDrv()
3461       drv[0] = (2.0*r-6.0/23.0)*t12*t21+t24*(20.0/17.0*t4+15.0/17.0*t25+3.0/  in BnDrv()
3469       t5 = r-1.0;  in BnDrv()
3477       t18 = s+r-1.0;  in BnDrv()
3480       t23 = r*r;  in BnDrv()
3481       t24 = t23-6.0/23.0*r+3.0/253.0;  in BnDrv()
3486       t33 = t31*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
3487       drv[0] = (2.0*r-6.0/23.0)*t15*t21+t24*(30.0/17.0*t3+75.0/34.0*t25+15.0/  in BnDrv()
3491       drv[2] = t31*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
3497       t6 = r-1.0;  in BnDrv()
3505       t20 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
3506       t22 = r*r;  in BnDrv()
3507       t23 = t22-6.0/23.0*r+3.0/253.0;  in BnDrv()
3513       drv[0] = (2.0*r-6.0/23.0)*t18*t20+t23*(42.0/17.0*t5+315.0/68.0*t24+105.0/  in BnDrv()
3524       t7 = r-1.0;  in BnDrv()
3531       t23 = r*r;  in BnDrv()
3532       t24 = t23-6.0/23.0*r+3.0/253.0;  in BnDrv()
3538       drv[0] = (2.0*r-6.0/23.0)*(t4*t3+56.0/17.0*t6*t7+147.0/34.0*t9*t10+49.0/  in BnDrv()
3547       t1 = r*r;  in BnDrv()
3551       t7 = s+r-1.0;  in BnDrv()
3557       t22 = t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0;  in BnDrv()
3564       drv[0] = (3.0*t1-24.0/23.0*r+18.0/253.0)*(t4*t3+42.0/13.0*t6*t7+105.0/  in BnDrv()
3572       t1 = r*r;  in BnDrv()
3573       t3 = s+2.0/15.0*r-2.0/15.0;  in BnDrv()
3576       t8 = s+r-1.0;  in BnDrv()
3583       t21 = t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0;  in BnDrv()
3591       drv[0] = t3*(3.0*t1-24.0/23.0*r+18.0/253.0)*t18+2.0/15.0*t22+t28;  in BnDrv()
3597       t1 = r*r;  in BnDrv()
3599       t4 = r-1.0;  in BnDrv()
3605       t12 = s+r-1.0;  in BnDrv()
3610       t22 = t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0;  in BnDrv()
3615       drv[0] = (3.0*t1-24.0/23.0*r+18.0/253.0)*t7*t19+t22*(2.0/5.0*s+2.0/35.0*r  in BnDrv()
3617       drv[1] = t22*(2.0*s+2.0/5.0*r-2.0/5.0)*t19+t30;  in BnDrv()
3621       t1 = r*r;  in BnDrv()
3623       t5 = r-1.0;  in BnDrv()
3627       t14 = s+r-1.0;  in BnDrv()
3630       t22 = t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0;  in BnDrv()
3635       drv[0] = (3.0*t1-24.0/23.0*r+18.0/253.0)*t10*t19+t22*(4.0/5.0*t3+12.0/  in BnDrv()
3641       t1 = r*r;  in BnDrv()
3645       t6 = r-1.0;  in BnDrv()
3651       t16 = s+r-1.0;  in BnDrv()
3654       t22 = t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0;  in BnDrv()
3658       t30 = t28*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
3659       drv[0] = (3.0*t1-24.0/23.0*r+18.0/253.0)*t13*t19+t22*(4.0/3.0*t5+8.0/7.0*  in BnDrv()
3662       drv[2] = t28*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
3665       t1 = r*r;  in BnDrv()
3668       t6 = r-1.0;  in BnDrv()
3675       t18 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
3676       t21 = t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0;  in BnDrv()
3681       drv[0] = (3.0*t1-24.0/23.0*r+18.0/253.0)*t16*t18+t21*(2.0*t4+20.0/7.0*t22  in BnDrv()
3688       t1 = r*r;  in BnDrv()
3692       t7 = r-1.0;  in BnDrv()
3698       t22 = t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0;  in BnDrv()
3703       drv[0] = (3.0*t1-24.0/23.0*r+18.0/253.0)*(t4*t3+14.0/5.0*t6*t7+3.0*t4*t9+  in BnDrv()
3711       t1 = r*r;  in BnDrv()
3712       t2 = t1*r;  in BnDrv()
3715       t7 = s+r-1.0;  in BnDrv()
3721       t20 = t19-20.0/23.0*t2+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0;  in BnDrv()
3727       drv[0] = (4.0*t2-60.0/23.0*t1+120.0/253.0*r-40.0/1771.0)*(t5*t+30.0/11.0*  in BnDrv()
3734       t1 = r*r;  in BnDrv()
3735       t2 = t1*r;  in BnDrv()
3736       t4 = s+2.0/13.0*r-2.0/13.0;  in BnDrv()
3740       t9 = s+r-1.0;  in BnDrv()
3746       t19 = t18-20.0/23.0*t2+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0;  in BnDrv()
3752       drv[0] = t4*(4.0*t2-60.0/23.0*t1+120.0/253.0*r-40.0/1771.0)*t16+2.0/13.0*  in BnDrv()
3758       t1 = r*r;  in BnDrv()
3759       t2 = t1*r;  in BnDrv()
3761       t5 = r-1.0;  in BnDrv()
3765       t12 = s+r-1.0;  in BnDrv()
3769       t20 = t19-20.0/23.0*t2+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0;  in BnDrv()
3773       drv[0] = (4.0*t2-60.0/23.0*t1+120.0/253.0*r-40.0/1771.0)*t8*t17+t20*(6.0/  in BnDrv()
3774 13.0*s+r/13-1.0/13.0)*t17+t27;  in BnDrv()
3775       drv[1] = t20*(2.0*s+6.0/13.0*r-6.0/13.0)*t17+t27;  in BnDrv()
3779       t1 = r*r;  in BnDrv()
3780       t2 = t1*r;  in BnDrv()
3782       t6 = r-1.0;  in BnDrv()
3786       t14 = s+r-1.0;  in BnDrv()
3790       t20 = t19-20.0/23.0*t2+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0;  in BnDrv()
3793       t27 = t25*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
3794       drv[0] = (4.0*t2-60.0/23.0*t1+120.0/253.0*r-40.0/1771.0)*t11*t17+t20*(  in BnDrv()
3797       drv[2] = t25*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
3800       t1 = r*r;  in BnDrv()
3801       t2 = t1*r;  in BnDrv()
3805       t7 = r-1.0;  in BnDrv()
3810       t16 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
3812       t19 = t18-20.0/23.0*t2+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0;  in BnDrv()
3816       drv[0] = (4.0*t2-60.0/23.0*t1+120.0/253.0*r-40.0/1771.0)*t14*t16+t19*(  in BnDrv()
3823       t1 = r*r;  in BnDrv()
3824       t2 = t1*r;  in BnDrv()
3827       t7 = r-1.0;  in BnDrv()
3833       t20 = t19-20.0/23.0*t2+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0;  in BnDrv()
3837       drv[0] = (4.0*t2-60.0/23.0*t1+120.0/253.0*r-40.0/1771.0)*(t5*s+30.0/13.0*  in BnDrv()
3845       t1 = r*r;  in BnDrv()
3847       t3 = t1*r;  in BnDrv()
3851       t8 = s+r-1.0;  in BnDrv()
3855       t18 = t2*r-30.0/23.0*t2+150.0/253.0*t3-200.0/1771.0*t1+15.0/1771.0*r-6.0/  in BnDrv()
3860       drv[0] = (5.0*t2-120.0/23.0*t3+450.0/253.0*t1-400.0/1771.0*r+15.0/1771.0)  in BnDrv()
3866       t1 = r*r;  in BnDrv()
3868       t3 = t1*r;  in BnDrv()
3869       t5 = s+2.0/11.0*r-2.0/11.0;  in BnDrv()
3871       t9 = s+r-1.0;  in BnDrv()
3874       t17 = t2*r-30.0/23.0*t2+150.0/253.0*t3-200.0/1771.0*t1+15.0/1771.0*r-6.0/  in BnDrv()
3880       drv[0] = (5.0*t2-120.0/23.0*t3+450.0/253.0*t1-400.0/1771.0*r+15.0/1771.0)  in BnDrv()
3886       t1 = r*r;  in BnDrv()
3888       t3 = t1*r;  in BnDrv()
3890       t6 = r-1.0;  in BnDrv()
3894       t12 = s+r-1.0;  in BnDrv()
3897       t18 = t2*r-30.0/23.0*t2+150.0/253.0*t3-200.0/1771.0*t1+15.0/1771.0*r-6.0/  in BnDrv()
3900       t24 = t22*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
3901       drv[0] = (5.0*t2-120.0/23.0*t3+450.0/253.0*t1-400.0/1771.0*r+15.0/1771.0)  in BnDrv()
3902 *t9*t15+t18*(6.0/11.0*s+6.0/55.0*r-6.0/55.0)*t15+t24;  in BnDrv()
3903       drv[1] = t18*(2.0*s+6.0/11.0*r-6.0/11.0)*t15+t24;  in BnDrv()
3904       drv[2] = t22*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
3907       t1 = r*r;  in BnDrv()
3909       t3 = t1*r;  in BnDrv()
3911       t7 = r-1.0;  in BnDrv()
3914       t14 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
3915       t17 = t2*r-30.0/23.0*t2+150.0/253.0*t3-200.0/1771.0*t1+15.0/1771.0*r-6.0/  in BnDrv()
3919       drv[0] = (5.0*t2-120.0/23.0*t3+450.0/253.0*t1-400.0/1771.0*r+15.0/1771.0)  in BnDrv()
3925       t1 = r*r;  in BnDrv()
3927       t3 = t1*r;  in BnDrv()
3931       t8 = r-1.0;  in BnDrv()
3935       t18 = t2*r-30.0/23.0*t2+150.0/253.0*t3-200.0/1771.0*t1+15.0/1771.0*r-6.0/  in BnDrv()
3939       drv[0] = (5.0*t2-120.0/23.0*t3+450.0/253.0*t1-400.0/1771.0*r+15.0/1771.0)  in BnDrv()
3946       t1 = r*r;  in BnDrv()
3948       t3 = t2*r;  in BnDrv()
3949       t4 = t1*r;  in BnDrv()
3951       t8 = s+r-1.0;  in BnDrv()
3954 /4807.0*r+1.0/14421.0;  in BnDrv()
3958 r-18.0/4807.0)*(t6*t+12.0/7.0*t6*t8+6.0/7.0*t*t10+4.0/35.0*t10*t8)+t19;  in BnDrv()
3963       t1 = r*r;  in BnDrv()
3965       t3 = t2*r;  in BnDrv()
3966       t4 = t1*r;  in BnDrv()
3967       t6 = s+2.0/9.0*r-2.0/9.0;  in BnDrv()
3969       t9 = s+r-1.0;  in BnDrv()
3973 /4807.0*r+1.0/14421.0;  in BnDrv()
3976       t19 = t17*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
3978 r-18.0/4807.0)*t6*t12+2.0/9.0*t16+t19;  in BnDrv()
3980       drv[2] = t17*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
3983       t1 = r*r;  in BnDrv()
3985       t3 = t2*r;  in BnDrv()
3986       t4 = t1*r;  in BnDrv()
3988       t7 = r-1.0;  in BnDrv()
3991       t12 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
3993 /4807.0*r+1.0/14421.0;  in BnDrv()
3996 r-18.0/4807.0)*t10*t12+t15*(2.0/3.0*s+r/6-1.0/6.0)*t12+2.0/3.0*t19;  in BnDrv()
3997       drv[1] = t15*(2.0*s+2.0/3.0*r-2.0/3.0)*t12+2.0/3.0*t19;  in BnDrv()
4001       t1 = r*r;  in BnDrv()
4003       t3 = t2*r;  in BnDrv()
4004       t4 = t1*r;  in BnDrv()
4006       t8 = r-1.0;  in BnDrv()
4009 /4807.0*r+1.0/14421.0;  in BnDrv()
4012 r-18.0/4807.0)*(t6*s+4.0/3.0*t6*t8+s*t10/2+t10*t8/21)+t16*(4.0/3.0*t6+t17+t10/7  in BnDrv()
4018       t1 = r*r;  in BnDrv()
4021       t4 = t2*r;  in BnDrv()
4022       t5 = t1*r;  in BnDrv()
4024       t8 = s+r-1.0;  in BnDrv()
4027 -168.0/4807.0*t1+28.0/14421.0*r-8.0/245157.0;  in BnDrv()
4028       t16 = t14*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
4030 253.0*t1-336.0/4807.0*r+28.0/14421.0)*(t7+6.0/5.0*t*t8+3.0/10.0*t10)+t16;  in BnDrv()
4032       drv[2] = t14*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
4035       t1 = r*r;  in BnDrv()
4038       t4 = t2*r;  in BnDrv()
4039       t5 = t1*r;  in BnDrv()
4040       t7 = s+2.0/7.0*r-2.0/7.0;  in BnDrv()
4041       t9 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
4043 -168.0/4807.0*t1+28.0/14421.0*r-8.0/245157.0;  in BnDrv()
4047 253.0*t1-336.0/4807.0*r+28.0/14421.0)*t7*t9+2.0/7.0*t13+2.0/3.0*t14;  in BnDrv()
4052       t1 = r*r;  in BnDrv()
4055       t4 = t2*r;  in BnDrv()
4056       t5 = t1*r;  in BnDrv()
4058       t8 = r-1.0;  in BnDrv()
4061 -168.0/4807.0*t1+28.0/14421.0*r-8.0/245157.0;  in BnDrv()
4063 253.0*t1-336.0/4807.0*r+28.0/14421.0)*(t7+6.0/7.0*s*t8+t10/7)+t14*(6.0/7.0*s+  in BnDrv()
4064 2.0/7.0*r-2.0/7.0);  in BnDrv()
4065       drv[1] = t14*(2.0*s+6.0/7.0*r-6.0/7.0);  in BnDrv()
4069       t1 = r*r;  in BnDrv()
4070       t2 = t1*r;  in BnDrv()
4074       t6 = t3*r;  in BnDrv()
4077 253.0*t2-3024.0/4807.0*t1+224.0/4807.0*r-96.0/81719.0)*(t+2.0/3.0*s+2.0/3.0*r  in BnDrv()
4079 -672.0/4807.0*t2+224.0/14421.0*t1-64.0/81719.0*r+1.0/81719.0;  in BnDrv()
4081 253.0*t3-672.0/4807.0*t2+224.0/14421.0*t1-64.0/81719.0*r+1.0/81719.0;  in BnDrv()
4083 -1008.0/4807.0*t2+112.0/4807.0*t1-96.0/81719.0*r+3.0/163438.0;  in BnDrv()
4086       t1 = r*r;  in BnDrv()
4087       t2 = t1*r;  in BnDrv()
4091       t6 = t3*r;  in BnDrv()
4094 253.0*t2-3024.0/4807.0*t1+224.0/4807.0*r-96.0/81719.0)*(s+2.0/5.0*r-2.0/5.0)+  in BnDrv()
4096 -2016.0/24035.0*t2+224.0/24035.0*t1-192.0/408595.0*r+3.0/408595.0;  in BnDrv()
4098 -1008.0/4807.0*t2+112.0/4807.0*t1-96.0/81719.0*r+3.0/163438.0;  in BnDrv()
4102       t1 = r*r;  in BnDrv()
4105       t4 = t1*r;  in BnDrv()
4106       drv[0] = 9.0*t3-720.0/23.0*t4*t2+11340.0/253.0*t2*t1-8640.0/253.0*t2*r+  in BnDrv()
4107 3780.0/253.0*t2-18144.0/4807.0*t4+2520.0/4807.0*t1-2880.0/81719.0*r+135.0/  in BnDrv()
4117       t5 = s+r-1.0;  in BnDrv()
4146       t5 = s+r-1.0;  in BnDrv()
4167       t26 = s+2.0/23.0*r-2.0/23.0;  in BnDrv()
4196       t7 = s+r-1.0;  in BnDrv()
4209       t27 = r-1.0;  in BnDrv()
4220       drv[0] = (6.0/23.0*s+6.0/253.0*r-6.0/253.0)*t24+t38;  in BnDrv()
4221       drv[1] = (2.0*s+6.0/23.0*r-6.0/23.0)*t24+t38;  in BnDrv()
4227       t2 = r-1.0;  in BnDrv()
4234       t11 = s+r-1.0;  in BnDrv()
4259       t3 = r-1.0;  in BnDrv()
4267       t13 = s+r-1.0;  in BnDrv()
4294       t4 = r-1.0;  in BnDrv()
4303       t15 = s+r-1.0;  in BnDrv()
4328       t4 = r-1.0;  in BnDrv()
4341       t18 = s+r-1.0;  in BnDrv()
4362       t5 = r-1.0;  in BnDrv()
4375       t20 = s+r-1.0;  in BnDrv()
4394       t6 = r-1.0;  in BnDrv()
4409       t22 = s+r-1.0;  in BnDrv()
4417       t38 = t36*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
4422       drv[2] = t36*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
4430       t6 = r-1.0;  in BnDrv()
4447       t24 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
4479       t5 = r-1.0;  in BnDrv()
4506       t5 = s+r-1.0;  in BnDrv()
4518       t26 = r-2.0/25.0;  in BnDrv()
4539       t1 = s+2.0/21.0*r-2.0/21.0;  in BnDrv()
4545       t7 = s+r-1.0;  in BnDrv()
4557       t26 = r-2.0/25.0;  in BnDrv()
4575       t2 = r-1.0;  in BnDrv()
4582       t11 = s+r-1.0;  in BnDrv()
4591       t27 = r-2.0/25.0;  in BnDrv()
4600       drv[0] = t5*t25+t27*(2.0/7.0*s+r/35-1.0/35.0)*t25+t38;  in BnDrv()
4601       drv[1] = t27*(2.0*s+2.0/7.0*r-2.0/7.0)*t25+t38;  in BnDrv()
4607       t3 = r-1.0;  in BnDrv()
4613       t13 = s+r-1.0;  in BnDrv()
4621       t27 = r-2.0/25.0;  in BnDrv()
4639       t4 = r-1.0;  in BnDrv()
4646       t15 = s+r-1.0;  in BnDrv()
4653       t27 = r-2.0/25.0;  in BnDrv()
4671       t4 = r-1.0;  in BnDrv()
4681       t18 = s+r-1.0;  in BnDrv()
4686       t27 = r-2.0/25.0;  in BnDrv()
4704       t5 = r-1.0;  in BnDrv()
4713       t20 = s+r-1.0;  in BnDrv()
4716       t27 = r-2.0/25.0;  in BnDrv()
4735       t6 = r-1.0;  in BnDrv()
4745       t22 = s+r-1.0;  in BnDrv()
4748       t27 = r-2.0/25.0;  in BnDrv()
4755       t38 = t36*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
4760       drv[2] = t36*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
4768       t6 = r-1.0;  in BnDrv()
4780       t24 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
4781       t26 = r-2.0/25.0;  in BnDrv()
4800       t5 = r-1.0;  in BnDrv()
4812       t26 = r-2.0/25.0;  in BnDrv()
4836       t7 = s+r-1.0;  in BnDrv()
4846       t26 = r*r;  in BnDrv()
4847       t27 = t26-6.0/25.0*r+1.0/100.0;  in BnDrv()
4856       drv[0] = (2.0*r-6.0/25.0)*(t4+72.0/17.0*t6*t7+126.0/17.0*t9*t10+588.0/  in BnDrv()
4864       t2 = s+2.0/19.0*r-2.0/19.0;  in BnDrv()
4869       t9 = s+r-1.0;  in BnDrv()
4878       t25 = r*r;  in BnDrv()
4879       t26 = t25-6.0/25.0*r+1.0/100.0;  in BnDrv()
4889       drv[0] = t2*(2.0*r-6.0/25.0)*t23+2.0/19.0*t27+t35;  in BnDrv()
4896       t3 = r-1.0;  in BnDrv()
4902       t12 = s+r-1.0;  in BnDrv()
4910       t26 = r*r;  in BnDrv()
4911       t27 = t26-6.0/25.0*r+1.0/100.0;  in BnDrv()
4919       drv[0] = (2.0*r-6.0/25.0)*t6*t24+t27*(6.0/19.0*s+2.0/57.0*r-2.0/57.0)*t24  in BnDrv()
4921       drv[1] = t27*(2.0*s+6.0/19.0*r-6.0/19.0)*t24+t37;  in BnDrv()
4927       t4 = r-1.0;  in BnDrv()
4932       t14 = s+r-1.0;  in BnDrv()
4939       t26 = r*r;  in BnDrv()
4940       t27 = t26-6.0/25.0*r+1.0/100.0;  in BnDrv()
4948       drv[0] = (2.0*r-6.0/25.0)*t9*t24+t27*(12.0/19.0*t2+4.0/19.0*t28+4.0/323.0  in BnDrv()
4958       t5 = r-1.0;  in BnDrv()
4966       t17 = s+r-1.0;  in BnDrv()
4971       t26 = r*r;  in BnDrv()
4972       t27 = t26-6.0/25.0*r+1.0/100.0;  in BnDrv()
4979       drv[0] = (2.0*r-6.0/25.0)*t12*t24+t27*(20.0/19.0*t4+40.0/57.0*t28+40.0/  in BnDrv()
4987       t5 = r-1.0;  in BnDrv()
4995       t19 = s+r-1.0;  in BnDrv()
4998       t26 = r*r;  in BnDrv()
4999       t27 = t26-6.0/25.0*r+1.0/100.0;  in BnDrv()
5006       drv[0] = (2.0*r-6.0/25.0)*t15*t24+t27*(30.0/19.0*t3+100.0/57.0*t28+200.0/  in BnDrv()
5016       t6 = r-1.0;  in BnDrv()
5025       t21 = s+r-1.0;  in BnDrv()
5028       t26 = r*r;  in BnDrv()
5029       t27 = t26-6.0/25.0*r+1.0/100.0;  in BnDrv()
5035       t37 = t35*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
5036       drv[0] = (2.0*r-6.0/25.0)*t18*t24+t27*(42.0/19.0*t5+70.0/19.0*t28+700.0/  in BnDrv()
5040       drv[2] = t35*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
5047       t7 = r-1.0;  in BnDrv()
5056       t23 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
5057       t25 = r*r;  in BnDrv()
5058       t26 = t25-6.0/25.0*r+1.0/100.0;  in BnDrv()
5065       drv[0] = (2.0*r-6.0/25.0)*t21*t23+t26*(56.0/19.0*t6+392.0/57.0*t27+1960.0  in BnDrv()
5078       t7 = r-1.0;  in BnDrv()
5088       t26 = r*r;  in BnDrv()
5089       t27 = t26-6.0/25.0*r+1.0/100.0;  in BnDrv()
5096       drv[0] = (2.0*r-6.0/25.0)*(t4+72.0/19.0*t6*t7+112.0/19.0*t9*t10+1568.0/  in BnDrv()
5105       t1 = r*r;  in BnDrv()
5110       t8 = s+r-1.0;  in BnDrv()
5117       t25 = t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0;  in BnDrv()
5125       drv[0] = (3.0*t1-24.0/25.0*r+3.0/50.0)*(t5*t4+56.0/15.0*t7*t8+28.0/5.0*  in BnDrv()
5133       t1 = r*r;  in BnDrv()
5134       t3 = s+2.0/17.0*r-2.0/17.0;  in BnDrv()
5138       t9 = s+r-1.0;  in BnDrv()
5146       t24 = t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0;  in BnDrv()
5155       drv[0] = t3*(3.0*t1-24.0/25.0*r+3.0/50.0)*t21+2.0/17.0*t25+t32;  in BnDrv()
5161       t1 = r*r;  in BnDrv()
5163       t4 = r-1.0;  in BnDrv()
5168       t12 = s+r-1.0;  in BnDrv()
5175       t25 = t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0;  in BnDrv()
5182       drv[0] = (3.0*t1-24.0/25.0*r+3.0/50.0)*t7*t22+t25*(6.0/17.0*s+3.0/68.0*r  in BnDrv()
5184       drv[1] = t25*(2.0*s+6.0/17.0*r-6.0/17.0)*t22+t34;  in BnDrv()
5189       t1 = r*r;  in BnDrv()
5191       t5 = r-1.0;  in BnDrv()
5197       t15 = s+r-1.0;  in BnDrv()
5202       t25 = t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0;  in BnDrv()
5208       drv[0] = (3.0*t1-24.0/25.0*r+3.0/50.0)*t10*t22+t25*(12.0/17.0*t3+9.0/34.0  in BnDrv()
5214       t1 = r*r;  in BnDrv()
5218       t6 = r-1.0;  in BnDrv()
5224       t17 = s+r-1.0;  in BnDrv()
5227       t25 = t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0;  in BnDrv()
5233       drv[0] = (3.0*t1-24.0/25.0*r+3.0/50.0)*t13*t22+t25*(20.0/17.0*t5+15.0/  in BnDrv()
5239       t1 = r*r;  in BnDrv()
5242       t6 = r-1.0;  in BnDrv()
5250       t19 = s+r-1.0;  in BnDrv()
5253       t25 = t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0;  in BnDrv()
5258       t34 = t32*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
5259       drv[0] = (3.0*t1-24.0/25.0*r+3.0/50.0)*t16*t22+t25*(30.0/17.0*t4+75.0/  in BnDrv()
5263       drv[2] = t32*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
5266       t1 = r*r;  in BnDrv()
5270       t7 = r-1.0;  in BnDrv()
5278       t21 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
5279       t24 = t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0;  in BnDrv()
5285       drv[0] = (3.0*t1-24.0/25.0*r+3.0/50.0)*t19*t21+t24*(42.0/17.0*t6+315.0/  in BnDrv()
5293       t1 = r*r;  in BnDrv()
5298       t8 = r-1.0;  in BnDrv()
5305       t25 = t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0;  in BnDrv()
5311       drv[0] = (3.0*t1-24.0/25.0*r+3.0/50.0)*(t5*t4+56.0/17.0*t7*t8+147.0/34.0*  in BnDrv()
5320       t1 = r*r;  in BnDrv()
5321       t2 = t1*r;  in BnDrv()
5325       t8 = s+r-1.0;  in BnDrv()
5332       t23 = t22-4.0/5.0*t2+t1/5-2.0/115.0*r+1.0/2530.0;  in BnDrv()
5339       drv[0] = (4.0*t2-12.0/5.0*t1+2.0/5.0*r-2.0/115.0)*(t5*t4+42.0/13.0*t7*t8+  in BnDrv()
5347       t1 = r*r;  in BnDrv()
5348       t2 = t1*r;  in BnDrv()
5349       t4 = s+2.0/15.0*r-2.0/15.0;  in BnDrv()
5352       t9 = s+r-1.0;  in BnDrv()
5360       t22 = t21-4.0/5.0*t2+t1/5-2.0/115.0*r+1.0/2530.0;  in BnDrv()
5368       drv[0] = t4*(4.0*t2-12.0/5.0*t1+2.0/5.0*r-2.0/115.0)*t19+2.0/15.0*t23+t29  in BnDrv()
5375       t1 = r*r;  in BnDrv()
5376       t2 = t1*r;  in BnDrv()
5378       t5 = r-1.0;  in BnDrv()
5384       t13 = s+r-1.0;  in BnDrv()
5390       t23 = t22-4.0/5.0*t2+t1/5-2.0/115.0*r+1.0/2530.0;  in BnDrv()
5395       drv[0] = (4.0*t2-12.0/5.0*t1+2.0/5.0*r-2.0/115.0)*t8*t20+t23*(2.0/5.0*s+  in BnDrv()
5396 2.0/35.0*r-2.0/35.0)*t20+t31;  in BnDrv()
5397       drv[1] = t23*(2.0*s+2.0/5.0*r-2.0/5.0)*t20+t31;  in BnDrv()
5401       t1 = r*r;  in BnDrv()
5402       t2 = t1*r;  in BnDrv()
5404       t6 = r-1.0;  in BnDrv()
5408       t15 = s+r-1.0;  in BnDrv()
5412       t23 = t22-4.0/5.0*t2+t1/5-2.0/115.0*r+1.0/2530.0;  in BnDrv()
5417       drv[0] = (4.0*t2-12.0/5.0*t1+2.0/5.0*r-2.0/115.0)*t11*t20+t23*(4.0/5.0*t4  in BnDrv()
5423       t1 = r*r;  in BnDrv()
5424       t2 = t1*r;  in BnDrv()
5428       t7 = r-1.0;  in BnDrv()
5434       t17 = s+r-1.0;  in BnDrv()
5438       t23 = t22-4.0/5.0*t2+t1/5-2.0/115.0*r+1.0/2530.0;  in BnDrv()
5442       t31 = t29*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
5443       drv[0] = (4.0*t2-12.0/5.0*t1+2.0/5.0*r-2.0/115.0)*t14*t20+t23*(4.0/3.0*t6  in BnDrv()
5446       drv[2] = t29*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
5449       t1 = r*r;  in BnDrv()
5450       t2 = t1*r;  in BnDrv()
5453       t7 = r-1.0;  in BnDrv()
5460       t19 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
5462       t22 = t21-4.0/5.0*t2+t1/5-2.0/115.0*r+1.0/2530.0;  in BnDrv()
5467       drv[0] = (4.0*t2-12.0/5.0*t1+2.0/5.0*r-2.0/115.0)*t17*t19+t22*(2.0*t5+  in BnDrv()
5474       t1 = r*r;  in BnDrv()
5475       t2 = t1*r;  in BnDrv()
5479       t8 = r-1.0;  in BnDrv()
5486       t23 = t22-4.0/5.0*t2+t1/5-2.0/115.0*r+1.0/2530.0;  in BnDrv()
5491       drv[0] = (4.0*t2-12.0/5.0*t1+2.0/5.0*r-2.0/115.0)*(t5*t4+14.0/5.0*t7*t8+  in BnDrv()
5499       t1 = r*r;  in BnDrv()
5501       t3 = t1*r;  in BnDrv()
5504       t8 = s+r-1.0;  in BnDrv()
5509       t21 = t2*r-6.0/5.0*t2+t3/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0;  in BnDrv()
5515       drv[0] = (5.0*t2-24.0/5.0*t3+3.0/2.0*t1-4.0/23.0*r+3.0/506.0)*(t6*t+30.0/  in BnDrv()
5522       t1 = r*r;  in BnDrv()
5524       t3 = t1*r;  in BnDrv()
5525       t5 = s+2.0/13.0*r-2.0/13.0;  in BnDrv()
5529       t10 = s+r-1.0;  in BnDrv()
5534       t20 = t2*r-6.0/5.0*t2+t3/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0;  in BnDrv()
5540       drv[0] = (5.0*t2-24.0/5.0*t3+3.0/2.0*t1-4.0/23.0*r+3.0/506.0)*t5*t17+2.0/  in BnDrv()
5546       t1 = r*r;  in BnDrv()
5548       t3 = t1*r;  in BnDrv()
5550       t6 = r-1.0;  in BnDrv()
5554       t13 = s+r-1.0;  in BnDrv()
5557       t21 = t2*r-6.0/5.0*t2+t3/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0;  in BnDrv()
5561       drv[0] = (5.0*t2-24.0/5.0*t3+3.0/2.0*t1-4.0/23.0*r+3.0/506.0)*t9*t18+t21*  in BnDrv()
5562 (6.0/13.0*s+r/13-1.0/13.0)*t18+t28;  in BnDrv()
5563       drv[1] = t21*(2.0*s+6.0/13.0*r-6.0/13.0)*t18+t28;  in BnDrv()
5567       t1 = r*r;  in BnDrv()
5569       t3 = t1*r;  in BnDrv()
5571       t7 = r-1.0;  in BnDrv()
5575       t15 = s+r-1.0;  in BnDrv()
5578       t21 = t2*r-6.0/5.0*t2+t3/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0;  in BnDrv()
5581       t28 = t26*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
5582       drv[0] = (5.0*t2-24.0/5.0*t3+3.0/2.0*t1-4.0/23.0*r+3.0/506.0)*t12*t18+t21  in BnDrv()
5585       drv[2] = t26*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
5588       t1 = r*r;  in BnDrv()
5590       t3 = t1*r;  in BnDrv()
5594       t8 = r-1.0;  in BnDrv()
5599       t17 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
5600       t20 = t2*r-6.0/5.0*t2+t3/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0;  in BnDrv()
5604       drv[0] = (5.0*t2-24.0/5.0*t3+3.0/2.0*t1-4.0/23.0*r+3.0/506.0)*t15*t17+t20  in BnDrv()
5611       t1 = r*r;  in BnDrv()
5613       t3 = t1*r;  in BnDrv()
5616       t8 = r-1.0;  in BnDrv()
5621       t21 = t2*r-6.0/5.0*t2+t3/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0;  in BnDrv()
5625       drv[0] = (5.0*t2-24.0/5.0*t3+3.0/2.0*t1-4.0/23.0*r+3.0/506.0)*(t6*s+30.0/  in BnDrv()
5634       t1 = r*r;  in BnDrv()
5636       t3 = t2*r;  in BnDrv()
5637       t4 = t1*r;  in BnDrv()
5641       t9 = s+r-1.0;  in BnDrv()
5646 1265.0*r+1.0/25300.0;  in BnDrv()
5650       drv[0] = (6.0*t3-42.0/5.0*t2+21.0/5.0*t4-21.0/23.0*t1+21.0/253.0*r-3.0/  in BnDrv()
5656       t1 = r*r;  in BnDrv()
5658       t3 = t2*r;  in BnDrv()
5659       t4 = t1*r;  in BnDrv()
5660       t6 = s+2.0/11.0*r-2.0/11.0;  in BnDrv()
5662       t10 = s+r-1.0;  in BnDrv()
5666 1265.0*r+1.0/25300.0;  in BnDrv()
5671       drv[0] = (6.0*t3-42.0/5.0*t2+21.0/5.0*t4-21.0/23.0*t1+21.0/253.0*r-3.0/  in BnDrv()
5677       t1 = r*r;  in BnDrv()
5679       t3 = t2*r;  in BnDrv()
5680       t4 = t1*r;  in BnDrv()
5682       t7 = r-1.0;  in BnDrv()
5686       t13 = s+r-1.0;  in BnDrv()
5690 1265.0*r+1.0/25300.0;  in BnDrv()
5692       t25 = t23*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
5693       drv[0] = (6.0*t3-42.0/5.0*t2+21.0/5.0*t4-21.0/23.0*t1+21.0/253.0*r-3.0/  in BnDrv()
5694 1265.0)*t10*t16+t19*(6.0/11.0*s+6.0/55.0*r-6.0/55.0)*t16+t25;  in BnDrv()
5695       drv[1] = t19*(2.0*s+6.0/11.0*r-6.0/11.0)*t16+t25;  in BnDrv()
5696       drv[2] = t23*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
5699       t1 = r*r;  in BnDrv()
5701       t3 = t2*r;  in BnDrv()
5702       t4 = t1*r;  in BnDrv()
5704       t8 = r-1.0;  in BnDrv()
5707       t15 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
5709 1265.0*r+1.0/25300.0;  in BnDrv()
5712       drv[0] = (6.0*t3-42.0/5.0*t2+21.0/5.0*t4-21.0/23.0*t1+21.0/253.0*r-3.0/  in BnDrv()
5718       t1 = r*r;  in BnDrv()
5720       t3 = t2*r;  in BnDrv()
5721       t4 = t1*r;  in BnDrv()
5725       t9 = r-1.0;  in BnDrv()
5730 1265.0*r+1.0/25300.0;  in BnDrv()
5733       drv[0] = (6.0*t3-42.0/5.0*t2+21.0/5.0*t4-21.0/23.0*t1+21.0/253.0*r-3.0/  in BnDrv()
5740       t1 = r*r;  in BnDrv()
5743       t4 = t2*r;  in BnDrv()
5744       t5 = t1*r;  in BnDrv()
5746       t9 = s+r-1.0;  in BnDrv()
5749 1265.0*t1+7.0/6325.0*r-2.0/120175.0;  in BnDrv()
5753 -56.0/1265.0*r+7.0/6325.0)*(t7*t+12.0/7.0*t7*t9+6.0/7.0*t*t11+4.0/35.0*t11*t9)+  in BnDrv()
5759       t1 = r*r;  in BnDrv()
5762       t4 = t2*r;  in BnDrv()
5763       t5 = t1*r;  in BnDrv()
5764       t7 = s+2.0/9.0*r-2.0/9.0;  in BnDrv()
5766       t10 = s+r-1.0;  in BnDrv()
5770 1265.0*t1+7.0/6325.0*r-2.0/120175.0;  in BnDrv()
5773       t20 = t18*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
5775 -56.0/1265.0*r+7.0/6325.0)*t7*t13+2.0/9.0*t17+t20;  in BnDrv()
5777       drv[2] = t18*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
5780       t1 = r*r;  in BnDrv()
5783       t4 = t2*r;  in BnDrv()
5784       t5 = t1*r;  in BnDrv()
5786       t8 = r-1.0;  in BnDrv()
5789       t13 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
5791 1265.0*t1+7.0/6325.0*r-2.0/120175.0;  in BnDrv()
5794 -56.0/1265.0*r+7.0/6325.0)*t11*t13+t16*(2.0/3.0*s+r/6-1.0/6.0)*t13+2.0/3.0*t20;  in BnDrv()
5795       drv[1] = t16*(2.0*s+2.0/3.0*r-2.0/3.0)*t13+2.0/3.0*t20;  in BnDrv()
5799       t1 = r*r;  in BnDrv()
5802       t4 = t2*r;  in BnDrv()
5803       t5 = t1*r;  in BnDrv()
5805       t9 = r-1.0;  in BnDrv()
5808 1265.0*t1+7.0/6325.0*r-2.0/120175.0;  in BnDrv()
5811 -56.0/1265.0*r+7.0/6325.0)*(t7*s+4.0/3.0*t7*t9+s*t11/2+t11*t9/21)+t17*(4.0/3.0*  in BnDrv()
5817       t1 = r*r;  in BnDrv()
5818       t2 = t1*r;  in BnDrv()
5822       t6 = t3*r;  in BnDrv()
5824       t9 = s+r-1.0;  in BnDrv()
5828 /1265.0*t2+84.0/6325.0*t1-72.0/120175.0*r+1.0/120175.0;  in BnDrv()
5829       t17 = t15*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
5831 1265.0*t2-504.0/1265.0*t1+168.0/6325.0*r-72.0/120175.0)*(t8+6.0/5.0*t*t9+3.0/  in BnDrv()
5834       drv[2] = t15*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
5837       t1 = r*r;  in BnDrv()
5838       t2 = t1*r;  in BnDrv()
5842       t6 = t3*r;  in BnDrv()
5843       t8 = s+2.0/7.0*r-2.0/7.0;  in BnDrv()
5844       t10 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
5847 /1265.0*t2+84.0/6325.0*t1-72.0/120175.0*r+1.0/120175.0;  in BnDrv()
5851 1265.0*t2-504.0/1265.0*t1+168.0/6325.0*r-72.0/120175.0)*t8*t10+2.0/7.0*t14+2.0/  in BnDrv()
5857       t1 = r*r;  in BnDrv()
5858       t2 = t1*r;  in BnDrv()
5862       t6 = t3*r;  in BnDrv()
5864       t9 = r-1.0;  in BnDrv()
5868 /1265.0*t2+84.0/6325.0*t1-72.0/120175.0*r+1.0/120175.0;  in BnDrv()
5870 1265.0*t2-504.0/1265.0*t1+168.0/6325.0*r-72.0/120175.0)*(t8+6.0/7.0*s*t9+t11/7)  in BnDrv()
5871 +t15*(6.0/7.0*s+2.0/7.0*r-2.0/7.0);  in BnDrv()
5872       drv[1] = t15*(2.0*s+6.0/7.0*r-6.0/7.0);  in BnDrv()
5876       t1 = r*r;  in BnDrv()
5879       t4 = t1*r;  in BnDrv()
5882       t7 = t2*r;  in BnDrv()
5883       t11 = t3*r;  in BnDrv()
5885 -3024.0/1265.0*t4+378.0/1265.0*t1-432.0/24035.0*r+9.0/24035.0)*(t+2.0/3.0*s+2.0  in BnDrv()
5886 /3.0*r-2.0/3.0)+2.0/3.0*t11-12.0/5.0*t3+18.0/5.0*t5-336.0/115.0*t6+1764.0/  in BnDrv()
5887 1265.0*t7-504.0/1265.0*t2+84.0/1265.0*t4-144.0/24035.0*t1+6.0/24035.0*r-4.0/  in BnDrv()
5891 *t7-504.0/1265.0*t2+84.0/1265.0*t4-144.0/24035.0*t1+6.0/24035.0*r-4.0/1225785.0  in BnDrv()
5894 -756.0/1265.0*t2+126.0/1265.0*t4-216.0/24035.0*t1+9.0/24035.0*r-2.0/408595.0;  in BnDrv()
5897       t1 = r*r;  in BnDrv()
5900       t4 = t1*r;  in BnDrv()
5903       t7 = t2*r;  in BnDrv()
5904       t11 = t3*r;  in BnDrv()
5906 -3024.0/1265.0*t4+378.0/1265.0*t1-432.0/24035.0*r+9.0/24035.0)*(s+2.0/5.0*r-2.0  in BnDrv()
5908 -1512.0/6325.0*t2+252.0/6325.0*t4-432.0/120175.0*t1+18.0/120175.0*r-4.0/  in BnDrv()
5912 -756.0/1265.0*t2+126.0/1265.0*t4-216.0/24035.0*t1+9.0/24035.0*r-2.0/408595.0;  in BnDrv()
5916       t1 = r*r;  in BnDrv()
5919       t5 = t1*r;  in BnDrv()
5920       drv[0] = 10.0*t3*r-198.0/5.0*t3+66.0*t2*t5-1386.0/23.0*t2*t1+756.0/23.0*  in BnDrv()
5921 t2*r-252.0/23.0*t2+252.0/115.0*t5-108.0/437.0*t1+6.0/437.0*r-2.0/7429.0;  in BnDrv()
5930       t6 = s+r-1.0;  in BnDrv()
5962       t6 = s+r-1.0;  in BnDrv()
5985       t29 = s+2.0/25.0*r-2.0/25.0;  in BnDrv()
6013       t6 = s+r-1.0;  in BnDrv()
6029       t30 = r-1.0;  in BnDrv()
6042       drv[0] = (6.0/25.0*s+r/50-1.0/50.0)*t27+t42;  in BnDrv()
6043       drv[1] = (2.0*s+6.0/25.0*r-6.0/25.0)*t27+t42;  in BnDrv()
6050       t2 = r-1.0;  in BnDrv()
6058       t11 = s+r-1.0;  in BnDrv()
6088       t3 = r-1.0;  in BnDrv()
6097       t14 = s+r-1.0;  in BnDrv()
6125       t4 = r-1.0;  in BnDrv()
6135       t16 = s+r-1.0;  in BnDrv()
6161       t4 = r-1.0;  in BnDrv()
6173       t18 = s+r-1.0;  in BnDrv()
6199       t5 = r-1.0;  in BnDrv()
6214       t21 = s+r-1.0;  in BnDrv()
6236       t6 = r-1.0;  in BnDrv()
6251       t23 = s+r-1.0;  in BnDrv()
6273       t6 = r-1.0;  in BnDrv()
6291       t25 = s+r-1.0;  in BnDrv()
6297       t42 = t40*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
6304       drv[2] = t40*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
6311       t5 = r-1.0;  in BnDrv()
6332       t27 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
6365       t6 = r-1.0;  in BnDrv()
6396       t6 = s+r-1.0;  in BnDrv()
6409       t29 = r-2.0/27.0;  in BnDrv()
6431       t1 = s+2.0/23.0*r-2.0/23.0;  in BnDrv()
6435       t6 = s+r-1.0;  in BnDrv()
6450       t29 = r-2.0/27.0;  in BnDrv()
6471       t2 = r-1.0;  in BnDrv()
6479       t11 = s+r-1.0;  in BnDrv()
6492       t30 = r-2.0/27.0;  in BnDrv()
6502       drv[0] = t5*t28+t30*(6.0/23.0*s+6.0/253.0*r-6.0/253.0)*t28+t42;  in BnDrv()
6503       drv[1] = t30*(2.0*s+6.0/23.0*r-6.0/23.0)*t28+t42;  in BnDrv()
6509       t3 = r-1.0;  in BnDrv()
6516       t14 = s+r-1.0;  in BnDrv()
6525       t30 = r-2.0/27.0;  in BnDrv()
6544       t4 = r-1.0;  in BnDrv()
6552       t16 = s+r-1.0;  in BnDrv()
6560       t30 = r-2.0/27.0;  in BnDrv()
6580       t4 = r-1.0;  in BnDrv()
6589       t18 = s+r-1.0;  in BnDrv()
6596       t30 = r-2.0/27.0;  in BnDrv()
6617       t5 = r-1.0;  in BnDrv()
6628       t21 = s+r-1.0;  in BnDrv()
6633       t30 = r-2.0/27.0;  in BnDrv()
6653       t6 = r-1.0;  in BnDrv()
6663       t23 = s+r-1.0;  in BnDrv()
6666       t30 = r-2.0/27.0;  in BnDrv()
6687       t6 = r-1.0;  in BnDrv()
6701       t25 = s+r-1.0;  in BnDrv()
6704       t30 = r-2.0/27.0;  in BnDrv()
6712       t42 = t40*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
6718       drv[2] = t40*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
6724       t5 = r-1.0;  in BnDrv()
6739       t27 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
6740       t29 = r-2.0/27.0;  in BnDrv()
6762       t6 = r-1.0;  in BnDrv()
6775       t29 = r-2.0/27.0;  in BnDrv()
6799       t6 = s+r-1.0;  in BnDrv()
6811       t29 = r*r;  in BnDrv()
6812       t30 = t29-2.0/9.0*r+1.0/117.0;  in BnDrv()
6823       drv[0] = (2.0*r-2.0/9.0)*(t4*t+90.0/19.0*t4*t6+180.0/19.0*t9*t10+3360.0/  in BnDrv()
6832       t2 = s+2.0/21.0*r-2.0/21.0;  in BnDrv()
6838       t9 = s+r-1.0;  in BnDrv()
6850       t28 = r*r;  in BnDrv()
6851       t29 = t28-2.0/9.0*r+1.0/117.0;  in BnDrv()
6862       drv[0] = t2*(2.0*r-2.0/9.0)*t26+2.0/21.0*t30+t39;  in BnDrv()
6869       t3 = r-1.0;  in BnDrv()
6876       t13 = s+r-1.0;  in BnDrv()
6885       t29 = r*r;  in BnDrv()
6886       t30 = t29-2.0/9.0*r+1.0/117.0;  in BnDrv()
6895       drv[0] = (2.0*r-2.0/9.0)*t6*t27+t30*(2.0/7.0*s+r/35-1.0/35.0)*t27+t41;  in BnDrv()
6896       drv[1] = t30*(2.0*s+2.0/7.0*r-2.0/7.0)*t27+t41;  in BnDrv()
6902       t4 = r-1.0;  in BnDrv()
6908       t15 = s+r-1.0;  in BnDrv()
6916       t29 = r*r;  in BnDrv()
6917       t30 = t29-2.0/9.0*r+1.0/117.0;  in BnDrv()
6926       drv[0] = (2.0*r-2.0/9.0)*t9*t27+t30*(4.0/7.0*t2+6.0/35.0*t31+6.0/665.0*t6  in BnDrv()
6936       t5 = r-1.0;  in BnDrv()
6943       t17 = s+r-1.0;  in BnDrv()
6950       t29 = r*r;  in BnDrv()
6951       t30 = t29-2.0/9.0*r+1.0/117.0;  in BnDrv()
6960       drv[0] = (2.0*r-2.0/9.0)*t12*t27+t30*(20.0/21.0*t4+4.0/7.0*t31+12.0/133.0  in BnDrv()
6969       t5 = r-1.0;  in BnDrv()
6979       t20 = s+r-1.0;  in BnDrv()
6984       t29 = r*r;  in BnDrv()
6985       t30 = t29-2.0/9.0*r+1.0/117.0;  in BnDrv()
6993       drv[0] = (2.0*r-2.0/9.0)*t15*t27+t30*(10.0/7.0*t3+10.0/7.0*t31+60.0/133.0  in BnDrv()
7003       t6 = r-1.0;  in BnDrv()
7012       t22 = s+r-1.0;  in BnDrv()
7015       t29 = r*r;  in BnDrv()
7016       t30 = t29-2.0/9.0*r+1.0/117.0;  in BnDrv()
7024       drv[0] = (2.0*r-2.0/9.0)*t18*t27+t30*(2.0*t5+3.0*t31+30.0/19.0*t32+20.0/  in BnDrv()
7035       t7 = r-1.0;  in BnDrv()
7045       t24 = s+r-1.0;  in BnDrv()
7048       t29 = r*r;  in BnDrv()
7049       t30 = t29-2.0/9.0*r+1.0/117.0;  in BnDrv()
7056       t41 = t39*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
7057       drv[0] = (2.0*r-2.0/9.0)*t21*t27+t30*(8.0/3.0*t6+28.0/5.0*t31+84.0/19.0*  in BnDrv()
7061       drv[2] = t39*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
7069       t7 = r-1.0;  in BnDrv()
7081       t26 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
7082       t28 = r*r;  in BnDrv()
7083       t29 = t28-2.0/9.0*r+1.0/117.0;  in BnDrv()
7091       drv[0] = (2.0*r-2.0/9.0)*t24*t26+t29*(24.0/7.0*t6+48.0/5.0*t30+1008.0/  in BnDrv()
7102       t6 = r-1.0;  in BnDrv()
7114       t29 = r*r;  in BnDrv()
7115       t30 = t29-2.0/9.0*r+1.0/117.0;  in BnDrv()
7123       drv[0] = (2.0*r-2.0/9.0)*(t4*s+30.0/7.0*t4*t6+54.0/7.0*t9*t10+144.0/19.0*  in BnDrv()
7133       t1 = r*r;  in BnDrv()
7139       t8 = s+r-1.0;  in BnDrv()
7149       t28 = t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0;  in BnDrv()
7158       drv[0] = (3.0*t1-8.0/9.0*r+2.0/39.0)*(t5+72.0/17.0*t7*t8+126.0/17.0*t10*  in BnDrv()
7166       t1 = r*r;  in BnDrv()
7167       t3 = s+2.0/19.0*r-2.0/19.0;  in BnDrv()
7172       t10 = s+r-1.0;  in BnDrv()
7181       t27 = t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0;  in BnDrv()
7191       drv[0] = t3*(3.0*t1-8.0/9.0*r+2.0/39.0)*t24+2.0/19.0*t28+t36;  in BnDrv()
7197       t1 = r*r;  in BnDrv()
7199       t4 = r-1.0;  in BnDrv()
7205       t13 = s+r-1.0;  in BnDrv()
7213       t28 = t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0;  in BnDrv()
7221       drv[0] = (3.0*t1-8.0/9.0*r+2.0/39.0)*t7*t25+t28*(6.0/19.0*s+2.0/57.0*r  in BnDrv()
7223       drv[1] = t28*(2.0*s+6.0/19.0*r-6.0/19.0)*t25+t38;  in BnDrv()
7228       t1 = r*r;  in BnDrv()
7230       t5 = r-1.0;  in BnDrv()
7235       t15 = s+r-1.0;  in BnDrv()
7242       t28 = t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0;  in BnDrv()
7250       drv[0] = (3.0*t1-8.0/9.0*r+2.0/39.0)*t10*t25+t28*(12.0/19.0*t3+4.0/19.0*  in BnDrv()
7257       t1 = r*r;  in BnDrv()
7261       t6 = r-1.0;  in BnDrv()
7269       t18 = s+r-1.0;  in BnDrv()
7274       t28 = t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0;  in BnDrv()
7281       drv[0] = (3.0*t1-8.0/9.0*r+2.0/39.0)*t13*t25+t28*(20.0/19.0*t5+40.0/57.0*  in BnDrv()
7287       t1 = r*r;  in BnDrv()
7290       t6 = r-1.0;  in BnDrv()
7298       t20 = s+r-1.0;  in BnDrv()
7301       t28 = t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0;  in BnDrv()
7308       drv[0] = (3.0*t1-8.0/9.0*r+2.0/39.0)*t16*t25+t28*(30.0/19.0*t4+100.0/57.0  in BnDrv()
7315       t1 = r*r;  in BnDrv()
7319       t7 = r-1.0;  in BnDrv()
7328       t22 = s+r-1.0;  in BnDrv()
7331       t28 = t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0;  in BnDrv()
7337       t38 = t36*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
7338       drv[0] = (3.0*t1-8.0/9.0*r+2.0/39.0)*t19*t25+t28*(42.0/19.0*t6+70.0/19.0*  in BnDrv()
7342       drv[2] = t36*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
7345       t1 = r*r;  in BnDrv()
7350       t8 = r-1.0;  in BnDrv()
7359       t24 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
7360       t27 = t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0;  in BnDrv()
7367       drv[0] = (3.0*t1-8.0/9.0*r+2.0/39.0)*t22*t24+t27*(56.0/19.0*t7+392.0/57.0  in BnDrv()
7375       t1 = r*r;  in BnDrv()
7381       t8 = r-1.0;  in BnDrv()
7391       t28 = t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0;  in BnDrv()
7398       drv[0] = (3.0*t1-8.0/9.0*r+2.0/39.0)*(t5+72.0/19.0*t7*t8+112.0/19.0*t10*  in BnDrv()
7408       t1 = r*r;  in BnDrv()
7409       t2 = t1*r;  in BnDrv()
7414       t9 = s+r-1.0;  in BnDrv()
7422       t26 = t25-20.0/27.0*t2+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0;  in BnDrv()
7430       drv[0] = (4.0*t2-20.0/9.0*t1+40.0/117.0*r-8.0/585.0)*(t6*t5+56.0/15.0*t8*  in BnDrv()
7438       t1 = r*r;  in BnDrv()
7439       t2 = t1*r;  in BnDrv()
7440       t4 = s+2.0/17.0*r-2.0/17.0;  in BnDrv()
7444       t10 = s+r-1.0;  in BnDrv()
7453       t25 = t24-20.0/27.0*t2+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0;  in BnDrv()
7462       drv[0] = (4.0*t2-20.0/9.0*t1+40.0/117.0*r-8.0/585.0)*t4*t22+2.0/17.0*t26+  in BnDrv()
7469       t1 = r*r;  in BnDrv()
7470       t2 = t1*r;  in BnDrv()
7472       t5 = r-1.0;  in BnDrv()
7477       t13 = s+r-1.0;  in BnDrv()
7485       t26 = t25-20.0/27.0*t2+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0;  in BnDrv()
7492       drv[0] = (4.0*t2-20.0/9.0*t1+40.0/117.0*r-8.0/585.0)*t8*t23+t26*(6.0/17.0  in BnDrv()
7493 *s+3.0/68.0*r-3.0/68.0)*t23+t35;  in BnDrv()
7494       drv[1] = t26*(2.0*s+6.0/17.0*r-6.0/17.0)*t23+t35;  in BnDrv()
7499       t1 = r*r;  in BnDrv()
7500       t2 = t1*r;  in BnDrv()
7502       t6 = r-1.0;  in BnDrv()
7508       t16 = s+r-1.0;  in BnDrv()
7514       t26 = t25-20.0/27.0*t2+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0;  in BnDrv()
7520       drv[0] = (4.0*t2-20.0/9.0*t1+40.0/117.0*r-8.0/585.0)*t11*t23+t26*(12.0/  in BnDrv()
7526       t1 = r*r;  in BnDrv()
7527       t2 = t1*r;  in BnDrv()
7531       t7 = r-1.0;  in BnDrv()
7537       t18 = s+r-1.0;  in BnDrv()
7541       t26 = t25-20.0/27.0*t2+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0;  in BnDrv()
7547       drv[0] = (4.0*t2-20.0/9.0*t1+40.0/117.0*r-8.0/585.0)*t14*t23+t26*(20.0/  in BnDrv()
7553       t1 = r*r;  in BnDrv()
7554       t2 = t1*r;  in BnDrv()
7557       t7 = r-1.0;  in BnDrv()
7565       t20 = s+r-1.0;  in BnDrv()
7569       t26 = t25-20.0/27.0*t2+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0;  in BnDrv()
7574       t35 = t33*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
7575       drv[0] = (4.0*t2-20.0/9.0*t1+40.0/117.0*r-8.0/585.0)*t17*t23+t26*(30.0/  in BnDrv()
7579       drv[2] = t33*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
7582       t1 = r*r;  in BnDrv()
7583       t2 = t1*r;  in BnDrv()
7587       t8 = r-1.0;  in BnDrv()
7595       t22 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
7597       t25 = t24-20.0/27.0*t2+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0;  in BnDrv()
7603       drv[0] = (4.0*t2-20.0/9.0*t1+40.0/117.0*r-8.0/585.0)*t20*t22+t25*(42.0/  in BnDrv()
7611       t1 = r*r;  in BnDrv()
7612       t2 = t1*r;  in BnDrv()
7617       t9 = r-1.0;  in BnDrv()
7625       t26 = t25-20.0/27.0*t2+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0;  in BnDrv()
7631       drv[0] = (4.0*t2-20.0/9.0*t1+40.0/117.0*r-8.0/585.0)*(t6*t5+56.0/17.0*t8*  in BnDrv()
7640       t1 = r*r;  in BnDrv()
7642       t3 = t1*r;  in BnDrv()
7646       t9 = s+r-1.0;  in BnDrv()
7652       t24 = t2*r-10.0/9.0*t2+50.0/117.0*t3-8.0/117.0*t1+r/234-1.0/13455.0;  in BnDrv()
7659       drv[0] = (5.0*t2-40.0/9.0*t3+50.0/39.0*t1-16.0/117.0*r+1.0/234.0)*(t6*t5+  in BnDrv()
7667       t1 = r*r;  in BnDrv()
7669       t3 = t1*r;  in BnDrv()
7670       t5 = s+2.0/15.0*r-2.0/15.0;  in BnDrv()
7673       t10 = s+r-1.0;  in BnDrv()
7680       t23 = t2*r-10.0/9.0*t2+50.0/117.0*t3-8.0/117.0*t1+r/234-1.0/13455.0;  in BnDrv()
7688       drv[0] = (5.0*t2-40.0/9.0*t3+50.0/39.0*t1-16.0/117.0*r+1.0/234.0)*t5*t20+  in BnDrv()
7695       t1 = r*r;  in BnDrv()
7697       t3 = t1*r;  in BnDrv()
7699       t6 = r-1.0;  in BnDrv()
7705       t14 = s+r-1.0;  in BnDrv()
7710       t24 = t2*r-10.0/9.0*t2+50.0/117.0*t3-8.0/117.0*t1+r/234-1.0/13455.0;  in BnDrv()
7715       drv[0] = (5.0*t2-40.0/9.0*t3+50.0/39.0*t1-16.0/117.0*r+1.0/234.0)*t9*t21+  in BnDrv()
7716 t24*(2.0/5.0*s+2.0/35.0*r-2.0/35.0)*t21+t32;  in BnDrv()
7717       drv[1] = t24*(2.0*s+2.0/5.0*r-2.0/5.0)*t21+t32;  in BnDrv()
7721       t1 = r*r;  in BnDrv()
7723       t3 = t1*r;  in BnDrv()
7725       t7 = r-1.0;  in BnDrv()
7729       t16 = s+r-1.0;  in BnDrv()
7732       t24 = t2*r-10.0/9.0*t2+50.0/117.0*t3-8.0/117.0*t1+r/234-1.0/13455.0;  in BnDrv()
7737       drv[0] = (5.0*t2-40.0/9.0*t3+50.0/39.0*t1-16.0/117.0*r+1.0/234.0)*t12*t21  in BnDrv()
7743       t1 = r*r;  in BnDrv()
7745       t3 = t1*r;  in BnDrv()
7749       t8 = r-1.0;  in BnDrv()
7755       t18 = s+r-1.0;  in BnDrv()
7758       t24 = t2*r-10.0/9.0*t2+50.0/117.0*t3-8.0/117.0*t1+r/234-1.0/13455.0;  in BnDrv()
7762       t32 = t30*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
7763       drv[0] = (5.0*t2-40.0/9.0*t3+50.0/39.0*t1-16.0/117.0*r+1.0/234.0)*t15*t21  in BnDrv()
7766       drv[2] = t30*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
7769       t1 = r*r;  in BnDrv()
7771       t3 = t1*r;  in BnDrv()
7774       t8 = r-1.0;  in BnDrv()
7781       t20 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
7782       t23 = t2*r-10.0/9.0*t2+50.0/117.0*t3-8.0/117.0*t1+r/234-1.0/13455.0;  in BnDrv()
7787       drv[0] = (5.0*t2-40.0/9.0*t3+50.0/39.0*t1-16.0/117.0*r+1.0/234.0)*t18*t20  in BnDrv()
7795       t1 = r*r;  in BnDrv()
7797       t3 = t1*r;  in BnDrv()
7801       t9 = r-1.0;  in BnDrv()
7807       t24 = t2*r-10.0/9.0*t2+50.0/117.0*t3-8.0/117.0*t1+r/234-1.0/13455.0;  in BnDrv()
7812       drv[0] = (5.0*t2-40.0/9.0*t3+50.0/39.0*t1-16.0/117.0*r+1.0/234.0)*(t6*t5+  in BnDrv()
7821       t1 = r*r;  in BnDrv()
7823       t3 = t2*r;  in BnDrv()
7824       t4 = t1*r;  in BnDrv()
7827       t9 = s+r-1.0;  in BnDrv()
7833 4485.0*r+7.0/296010.0;  in BnDrv()
7839       drv[0] = (6.0*t3-70.0/9.0*t2+140.0/39.0*t4-28.0/39.0*t1+7.0/117.0*r-7.0/  in BnDrv()
7847       t1 = r*r;  in BnDrv()
7849       t3 = t2*r;  in BnDrv()
7850       t4 = t1*r;  in BnDrv()
7851       t6 = s+2.0/13.0*r-2.0/13.0;  in BnDrv()
7855       t11 = s+r-1.0;  in BnDrv()
7861 4485.0*r+7.0/296010.0;  in BnDrv()
7867       drv[0] = (6.0*t3-70.0/9.0*t2+140.0/39.0*t4-28.0/39.0*t1+7.0/117.0*r-7.0/  in BnDrv()
7873       t1 = r*r;  in BnDrv()
7875       t3 = t2*r;  in BnDrv()
7876       t4 = t1*r;  in BnDrv()
7878       t7 = r-1.0;  in BnDrv()
7882       t14 = s+r-1.0;  in BnDrv()
7886 4485.0*r+7.0/296010.0;  in BnDrv()
7890       drv[0] = (6.0*t3-70.0/9.0*t2+140.0/39.0*t4-28.0/39.0*t1+7.0/117.0*r-7.0/  in BnDrv()
7891 4485.0)*t10*t19+t22*(6.0/13.0*s+r/13-1.0/13.0)*t19+t29;  in BnDrv()
7892       drv[1] = t22*(2.0*s+6.0/13.0*r-6.0/13.0)*t19+t29;  in BnDrv()
7896       t1 = r*r;  in BnDrv()
7898       t3 = t2*r;  in BnDrv()
7899       t4 = t1*r;  in BnDrv()
7901       t8 = r-1.0;  in BnDrv()
7905       t16 = s+r-1.0;  in BnDrv()
7909 4485.0*r+7.0/296010.0;  in BnDrv()
7912       t29 = t27*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
7913       drv[0] = (6.0*t3-70.0/9.0*t2+140.0/39.0*t4-28.0/39.0*t1+7.0/117.0*r-7.0/  in BnDrv()
7916       drv[2] = t27*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
7919       t1 = r*r;  in BnDrv()
7921       t3 = t2*r;  in BnDrv()
7922       t4 = t1*r;  in BnDrv()
7926       t9 = r-1.0;  in BnDrv()
7931       t18 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
7933 4485.0*r+7.0/296010.0;  in BnDrv()
7937       drv[0] = (6.0*t3-70.0/9.0*t2+140.0/39.0*t4-28.0/39.0*t1+7.0/117.0*r-7.0/  in BnDrv()
7945       t1 = r*r;  in BnDrv()
7947       t3 = t2*r;  in BnDrv()
7948       t4 = t1*r;  in BnDrv()
7951       t9 = r-1.0;  in BnDrv()
7957 4485.0*r+7.0/296010.0;  in BnDrv()
7961       drv[0] = (6.0*t3-70.0/9.0*t2+140.0/39.0*t4-28.0/39.0*t1+7.0/117.0*r-7.0/  in BnDrv()
7970       t1 = r*r;  in BnDrv()
7973       t4 = t2*r;  in BnDrv()
7974       t5 = t1*r;  in BnDrv()
7978       t10 = s+r-1.0;  in BnDrv()
7983 -196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0;  in BnDrv()
7988 t1-392.0/13455.0*r+98.0/148005.0)*(t8+20.0/9.0*t9*t10+5.0/3.0*t7*t12+10.0/21.0*  in BnDrv()
7994       t1 = r*r;  in BnDrv()
7997       t4 = t2*r;  in BnDrv()
7998       t5 = t1*r;  in BnDrv()
7999       t7 = s+2.0/11.0*r-2.0/11.0;  in BnDrv()
8001       t11 = s+r-1.0;  in BnDrv()
8005 -196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0;  in BnDrv()
8011 t1-392.0/13455.0*r+98.0/148005.0)*t7*t16+2.0/11.0*t20+t24;  in BnDrv()
8016       t1 = r*r;  in BnDrv()
8019       t4 = t2*r;  in BnDrv()
8020       t5 = t1*r;  in BnDrv()
8022       t8 = r-1.0;  in BnDrv()
8026       t14 = s+r-1.0;  in BnDrv()
8030 -196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0;  in BnDrv()
8032       t26 = t24*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
8034 t1-392.0/13455.0*r+98.0/148005.0)*t11*t17+t20*(6.0/11.0*s+6.0/55.0*r-6.0/55.0)*  in BnDrv()
8036       drv[1] = t20*(2.0*s+6.0/11.0*r-6.0/11.0)*t17+t26;  in BnDrv()
8037       drv[2] = t24*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
8040       t1 = r*r;  in BnDrv()
8043       t4 = t2*r;  in BnDrv()
8044       t5 = t1*r;  in BnDrv()
8046       t9 = r-1.0;  in BnDrv()
8049       t16 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
8051 -196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0;  in BnDrv()
8055 t1-392.0/13455.0*r+98.0/148005.0)*t14*t16+t19*(12.0/11.0*t7+36.0/55.0*t20+4.0/  in BnDrv()
8061       t1 = r*r;  in BnDrv()
8064       t4 = t2*r;  in BnDrv()
8065       t5 = t1*r;  in BnDrv()
8069       t10 = r-1.0;  in BnDrv()
8074 -196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0;  in BnDrv()
8078 t1-392.0/13455.0*r+98.0/148005.0)*(t8+20.0/11.0*t9*t10+12.0/11.0*t7*t12+8.0/  in BnDrv()
8084       t1 = r*r;  in BnDrv()
8085       t2 = t1*r;  in BnDrv()
8089       t6 = t3*r;  in BnDrv()
8091       t10 = s+r-1.0;  in BnDrv()
8095 4485.0*t2+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0;  in BnDrv()
8099 -392.0/1495.0*t1+784.0/49335.0*r-16.0/49335.0)*(t8*t+12.0/7.0*t8*t10+6.0/7.0*t*  in BnDrv()
8105       t1 = r*r;  in BnDrv()
8106       t2 = t1*r;  in BnDrv()
8110       t6 = t3*r;  in BnDrv()
8111       t8 = s+2.0/9.0*r-2.0/9.0;  in BnDrv()
8113       t11 = s+r-1.0;  in BnDrv()
8118 4485.0*t2+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0;  in BnDrv()
8121       t21 = t19*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
8123 -392.0/1495.0*t1+784.0/49335.0*r-16.0/49335.0)*t8*t14+2.0/9.0*t18+t21;  in BnDrv()
8125       drv[2] = t19*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
8128       t1 = r*r;  in BnDrv()
8129       t2 = t1*r;  in BnDrv()
8133       t6 = t3*r;  in BnDrv()
8135       t9 = r-1.0;  in BnDrv()
8138       t14 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
8141 4485.0*t2+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0;  in BnDrv()
8144 -392.0/1495.0*t1+784.0/49335.0*r-16.0/49335.0)*t12*t14+t17*(2.0/3.0*s+r/6-1.0/  in BnDrv()
8146       drv[1] = t17*(2.0*s+2.0/3.0*r-2.0/3.0)*t14+2.0/3.0*t21;  in BnDrv()
8150       t1 = r*r;  in BnDrv()
8151       t2 = t1*r;  in BnDrv()
8155       t6 = t3*r;  in BnDrv()
8157       t10 = r-1.0;  in BnDrv()
8161 4485.0*t2+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0;  in BnDrv()
8164 -392.0/1495.0*t1+784.0/49335.0*r-16.0/49335.0)*(t8*s+4.0/3.0*t8*t10+s*t12/2+t12  in BnDrv()
8170       t1 = r*r;  in BnDrv()
8173       t4 = t1*r;  in BnDrv()
8176       t7 = t2*r;  in BnDrv()
8178       t10 = s+r-1.0;  in BnDrv()
8180       t16 = t3*r-10.0/3.0*t3+60.0/13.0*t5-224.0/65.0*t6+98.0/65.0*t7-588.0/  in BnDrv()
8181 1495.0*t2+196.0/3289.0*t4-16.0/3289.0*t1+3.0/16445.0*r-2.0/937365.0;  in BnDrv()
8182       t18 = t16*(6.0/5.0*t+3.0/5.0*s+3.0/5.0*r-3.0/5.0);  in BnDrv()
8184 -2352.0/1495.0*t4+588.0/3289.0*t1-32.0/3289.0*r+3.0/16445.0)*(t9+6.0/5.0*t*t10+  in BnDrv()
8187       drv[2] = t16*(2.0*t+6.0/5.0*s+6.0/5.0*r-6.0/5.0);  in BnDrv()
8190       t1 = r*r;  in BnDrv()
8193       t4 = t1*r;  in BnDrv()
8196       t7 = t2*r;  in BnDrv()
8197       t9 = s+2.0/7.0*r-2.0/7.0;  in BnDrv()
8198       t11 = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0;  in BnDrv()
8199       t14 = t3*r-10.0/3.0*t3+60.0/13.0*t5-224.0/65.0*t6+98.0/65.0*t7-588.0/  in BnDrv()
8200 1495.0*t2+196.0/3289.0*t4-16.0/3289.0*t1+3.0/16445.0*r-2.0/937365.0;  in BnDrv()
8204 -2352.0/1495.0*t4+588.0/3289.0*t1-32.0/3289.0*r+3.0/16445.0)*t9*t11+2.0/7.0*t15  in BnDrv()
8210       t1 = r*r;  in BnDrv()
8213       t4 = t1*r;  in BnDrv()
8216       t7 = t2*r;  in BnDrv()
8218       t10 = r-1.0;  in BnDrv()
8220       t16 = t3*r-10.0/3.0*t3+60.0/13.0*t5-224.0/65.0*t6+98.0/65.0*t7-588.0/  in BnDrv()
8221 1495.0*t2+196.0/3289.0*t4-16.0/3289.0*t1+3.0/16445.0*r-2.0/937365.0;  in BnDrv()
8223 -2352.0/1495.0*t4+588.0/3289.0*t1-32.0/3289.0*r+3.0/16445.0)*(t9+6.0/7.0*s*t10+  in BnDrv()
8224 t12/7)+t16*(6.0/7.0*s+2.0/7.0*r-2.0/7.0);  in BnDrv()
8225       drv[1] = t16*(2.0*s+6.0/7.0*r-6.0/7.0);  in BnDrv()
8229       t1 = r*r;  in BnDrv()
8232       t4 = t3*r;  in BnDrv()
8233       t5 = t1*r;  in BnDrv()
8236       t8 = t2*r;  in BnDrv()
8239 -2156.0/299.0*t2+392.0/299.0*t5-40.0/299.0*t1+2.0/299.0*r-2.0/17043.0)*(t+2.0/  in BnDrv()
8240 3.0*s+2.0/3.0*r-2.0/3.0)+2.0/3.0*t12-220.0/81.0*t4+550.0/117.0*t3-176.0/39.0*t6  in BnDrv()
8242 /51129.0*r+2.0/2300805.0;  in BnDrv()
8244 t7-4312.0/4485.0*t8+196.0/897.0*t2-80.0/2691.0*t5+2.0/897.0*t1-4.0/51129.0*r+  in BnDrv()
8247 1495.0*t8+98.0/299.0*t2-40.0/897.0*t5+t1/299-2.0/17043.0*r+1.0/766935.0;  in BnDrv()
8253       t1 = r*r;  in BnDrv()
8256       t4 = t3*r;  in BnDrv()
8257       t5 = t1*r;  in BnDrv()
8260       t8 = t2*r;  in BnDrv()
8263 -2156.0/299.0*t2+392.0/299.0*t5-40.0/299.0*t1+2.0/299.0*r-2.0/17043.0)*(s+2.0/  in BnDrv()
8264 5.0*r-2.0/5.0)+2.0/5.0*t12-44.0/27.0*t4+110.0/39.0*t3-176.0/65.0*t6+308.0/195.0  in BnDrv()
8265 *t7-4312.0/7475.0*t8+196.0/1495.0*t2-16.0/897.0*t5+2.0/1495.0*t1-4.0/85215.0*r+  in BnDrv()
8268 1495.0*t8+98.0/299.0*t2-40.0/897.0*t5+t1/299-2.0/17043.0*r+1.0/766935.0;  in BnDrv()
8274       t1 = r*r;  in BnDrv()
8277       t6 = t1*r;  in BnDrv()
8278       t10 = 11.0*t3*t1-440.0/9.0*t3*r+1210.0/13.0*t3-3872.0/39.0*t2*t6+847.0/  in BnDrv()
8279 13.0*t2*t1-40656.0/1495.0*t2*r+2156.0/299.0*t2-352.0/299.0*t6+33.0/299.0*t1  in BnDrv()
8280 -88.0/17043.0*r+22.0/255645.0;  in BnDrv()