21    double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18;  in Bn()  local
480       t14 = t12*t12;  in Bn()
482 t11+6.0/5.0*t*t12+3.0/10.0*t14);  in Bn()
663       t14 = t12*t12;  in Bn()
664       t18 = t14*t14;  in Bn()
666 t12+5.0/3.0*t9*t14+10.0/21.0*t*t14*t12+5.0/126.0*t18);  in Bn()
675       t14 = s+r-1.0;  in Bn()
676       t16 = t14*t14;  in Bn()
678 12.0/7.0*t14*t12+6.0/7.0*t*t16+4.0/35.0*t16*t14);  in Bn()
727       t14 = t10*t10;  in Bn()
729 11.0*t4*t*t10+40.0/33.0*t4*t10*t7+5.0/22.0*t*t14+t14*t7/77);  in Bn()
759       t14 = t*t;  in Bn()
763 /273)*(t14+6.0/5.0*t*t15+3.0/10.0*t17);  in Bn()
799       t14 = t10*t10;  in Bn()
801 t8+5.0/3.0*t5*t10+10.0/21.0*t*t10*t8+5.0/126.0*t14);  in Bn()
979       t14 = t9*t9;  in Bn()
980       t22 = t14*t14;  in Bn()
982 17.0*t2*t14+252.0/221.0*t4*t14*t6+42.0/221.0*t1*t14*t9+36.0/2431.0*t*t14*t12+  in Bn()
992       t14 = t10*t10;  in Bn()
994 +56.0/13.0*t4*t12+70.0/39.0*t3*t14+56.0/143.0*t2*t14*t7+28.0/715.0*t*t14*t10+  in Bn()
995 8.0/6435.0*t14*t12);  in Bn()
1084       t14 = t9*t9;  in Bn()
1085       t22 = t14*t14;  in Bn()
1087 735.0/323.0*t2*t14+196.0/323.0*t4*t14*t6+28.0/323.0*t1*t14*t9+24.0/4199.0*s*t14  in Bn()
1097       t14 = t10*t10;  in Bn()
1099 t4*t12+70.0/39.0*t3*t14+56.0/143.0*t2*t14*t7+28.0/715.0*t*t14*t10+8.0/6435.0*  in Bn()
1100 t14*t12);  in Bn()
1119       t14 = t11*t11;  in Bn()
1120       t18 = t14*t14;  in Bn()
1122 +30.0/11.0*t8*t*t14+40.0/33.0*t8*t14*t11+5.0/22.0*t*t18+t18*t11/77);  in Bn()
1130       t14 = s+r-1.0;  in Bn()
1131       t16 = t14*t14;  in Bn()
1134 20.0/9.0*t11*t*t14+5.0/3.0*t11*t16+10.0/21.0*t*t16*t14+5.0/126.0*t20);  in Bn()
1142       t14 = t*t;  in Bn()
1146 /476)*(t14*t+12.0/7.0*t16*t14+6.0/7.0*t*t18+4.0/35.0*t18*t16);  in Bn()
1177       t14 = t10*t10;  in Bn()
1179 17.0*t4*t12+35.0/34.0*t3*t14+42.0/221.0*t2*t14*t7+7.0/442.0*s*t14*t10+t14*t12  in Bn()
1188       t14 = t9*t9;  in Bn()
1190 350.0/143.0*t3*t*t9*t7+105.0/143.0*t3*t14+14.0/143.0*t*t14*t7+7.0/1716.0*t14*t9  in Bn()
1211       t14 = t12*t12;  in Bn()
1212       t18 = t14*t14;  in Bn()
1214 t12+5.0/3.0*t14*t9+10.0/21.0*t*t14*t12+5.0/126.0*t18);  in Bn()
1222       t14 = s+r-1.0;  in Bn()
1223       t16 = t14*t14;  in Bn()
1225 *t7*t5)*(t12*t+12.0/7.0*t14*t12+6.0/7.0*t*t16+4.0/35.0*t16*t14);  in Bn()
1256       t14 = t9*t9;  in Bn()
1258 *t3*s*t9*t7+5.0/13.0*t3*t14+6.0/143.0*s*t14*t7+t14*t9/715);  in Bn()
1266       t14 = t10*t10;  in Bn()
1268 11.0*t4*t*t10+40.0/33.0*t4*t10*t7+5.0/22.0*t*t14+t14*t7/77);  in Bn()
1287       t14 = t12*t12;  in Bn()
1289 t10*t+12.0/7.0*t10*t12+6.0/7.0*t*t14+4.0/35.0*t14*t12);  in Bn()
1297       t14 = s+r-1.0;  in Bn()
1298       t16 = t14*t14;  in Bn()
1300 13.0*s*t8+2.0/143.0*t8*t6)*(t13+6.0/5.0*t*t14+3.0/10.0*t16);  in Bn()
1318       t14 = t10*t10;  in Bn()
1320 13.0*t4*s*t10+100.0/143.0*t4*t10*t7+15.0/143.0*s*t14+2.0/429.0*t14*t7);  in Bn()
1329       t14 = t10*t10;  in Bn()
1331 t5*t*t8+5.0/3.0*t5*t10+10.0/21.0*t*t10*t8+5.0/126.0*t14);  in Bn()
1350       t14 = t12*t12;  in Bn()
1352 s*t6+3.0/55.0*t8)*(t11+6.0/5.0*t*t12+3.0/10.0*t14);  in Bn()
1370       t14 = t10*t10;  in Bn()
1372 *t5*s*t8+12.0/11.0*t5*t10+8.0/33.0*s*t10*t8+t14/66);  in Bn()
1492       t14 = t11*t11;  in Bn()
1493       t16 = t14*t11;  in Bn()
1494       t18 = t14*t14;  in Bn()
1496 t14+56.0/13.0*t8*t16+70.0/39.0*t7*t18+56.0/143.0*t6*t18*t11+28.0/715.0*t*t18*  in Bn()
1497 t14+8.0/6435.0*t18*t16);  in Bn()
1578       t14 = t9*t9;  in Bn()
1579       t22 = t14*t14;  in Bn()
1581 19.0*t2*t14+112.0/323.0*t4*t14*t6+14.0/323.0*t1*t14*t9+4.0/1615.0*s*t14*t12+t22  in Bn()
1618       t14 = t12*t9;  in Bn()
1621 /5.0*t6*t*t12+56.0/13.0*t6*t14+70.0/39.0*t5*t16+56.0/143.0*t4*t16*t9+28.0/715.0  in Bn()
1622 *t*t16*t12+8.0/6435.0*t16*t14);  in Bn()
1631       t14 = t12*t12;  in Bn()
1632       t19 = t14*t14;  in Bn()
1634 105.0/26.0*t9*t14+350.0/143.0*t8*t*t14*t12+105.0/143.0*t8*t19+14.0/143.0*t*t19*  in Bn()
1635 t12+7.0/1716.0*t19*t14);  in Bn()
1643       t14 = s+r-1.0;  in Bn()
1644       t17 = t14*t14;  in Bn()
1647 t12*t+30.0/11.0*t14*t12+30.0/11.0*t11*t*t17+40.0/33.0*t11*t17*t14+5.0/22.0*t*  in Bn()
1648 t21+t21*t14/77);  in Bn()
1656       t14 = t*t;  in Bn()
1657       t15 = t14*t14;  in Bn()
1662 t5+5.0/3876.0*t11)*(t15+20.0/9.0*t14*t*t17+5.0/3.0*t19*t14+10.0/21.0*t*t19*t17+  in Bn()
1698       t14 = t10*t10;  in Bn()
1700 969.0*t4*t12+1225.0/1938.0*t3*t14+98.0/969.0*t2*t14*t7+7.0/969.0*s*t14*t10+2.0/  in Bn()
1701 12597.0*t14*t12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);  in Bn()
1802       t14 = t9*t9;  in Bn()
1804 t9+35.0/34.0*t3*s*t9*t7+15.0/68.0*t3*t14+9.0/442.0*s*t14*t7+t14*t9/1768)*(t+2.0  in Bn()
1886       t14 = t10*t10;  in Bn()
1888 7.0*t4*s*t10+40.0/91.0*t4*t10*t7+5.0/91.0*s*t14+2.0/1001.0*t14*t7)*(t+2.0/3.0*s  in Bn()
1944       t14 = t*t;  in Bn()
1948 12.0/13.0*t5*t7+3.0/13.0*s*t9+2.0/143.0*t9*t7)*(t14+6.0/5.0*t*t15+3.0/10.0*t17)  in Bn()
1958       t14 = t10*t10;  in Bn()
1960 /13.0*t5*s*t8+10.0/13.0*t5*t10+20.0/143.0*s*t10*t8+t14/143)*(t+2.0/3.0*s+2.0/  in Bn()
2162       t14 = t11*t11;  in Bn()
2163       t17 = t14*t11;  in Bn()
2164       t19 = t14*t14;  in Bn()
2167 t7*t6*t14+588.0/85.0*t7*t*t17+63.0/17.0*t7*t19+252.0/221.0*t9*t19*t11+42.0/  in Bn()
2168 221.0*t6*t19*t14+36.0/2431.0*t*t19*t17+9.0/24310.0*t27);  in Bn()
2177       t14 = s+r-1.0;  in Bn()
2178       t17 = t14*t14;  in Bn()
2179       t19 = t17*t14;  in Bn()
2182 56.0/15.0*t11*t9*t14+28.0/5.0*t11*t*t17+56.0/13.0*t19*t11+70.0/39.0*t10*t21+  in Bn()
2183 56.0/143.0*t9*t21*t14+28.0/715.0*t*t21*t17+8.0/6435.0*t21*t19);  in Bn()
2253       t14 = t9*t9;  in Bn()
2254       t22 = t14*t14;  in Bn()
2259 252.0/253.0*t2*t14+1008.0/4807.0*t4*t14*t6+112.0/4807.0*t1*t14*t9+96.0/81719.0*  in Bn()
2260 s*t14*t12+3.0/163438.0*t22)*(t24+6.0/5.0*t*t25+3.0/10.0*t27);  in Bn()
2353       t14 = t*t;  in Bn()
2354       t15 = t14*t14;  in Bn()
2359 t11/1197)*(t15*t+30.0/11.0*t17*t15+30.0/11.0*t14*t*t20+40.0/33.0*t14*t20*t17+  in Bn()
2397       t14 = t10*t10;  in Bn()
2402 *t12+70.0/171.0*t3*t14+56.0/969.0*t2*t14*t7+7.0/1938.0*s*t14*t10+t14*t12/14535)  in Bn()
2441       t14 = t11*t8;  in Bn()
2445 t11+588.0/85.0*t4*t*t14+63.0/17.0*t4*t16+252.0/221.0*t6*t16*t8+42.0/221.0*t3*  in Bn()
2446 t16*t11+36.0/2431.0*t*t16*t14+9.0/24310.0*t24);  in Bn()
2525       t14 = t9*t9;  in Bn()
2530 700.0/969.0*t3*s*t9*t7+175.0/1292.0*t3*t14+7.0/646.0*s*t14*t7+t14*t9/3876)*(t21  in Bn()
2554       t14 = t11*t8;  in Bn()
2558 t11+1568.0/323.0*t4*s*t14+735.0/323.0*t4*t16+196.0/323.0*t6*t16*t8+28.0/323.0*  in Bn()
2559 t3*t16*t11+24.0/4199.0*s*t16*t14+t24/8398);  in Bn()
2568       t14 = t12*t9;  in Bn()
2571 28.0/5.0*t6*t*t12+56.0/13.0*t6*t14+70.0/39.0*t5*t16+56.0/143.0*t4*t16*t9+28.0/  in Bn()
2572 715.0*t*t16*t12+8.0/6435.0*t16*t14);  in Bn()
2605       t14 = t13*t13;  in Bn()
2610 68.0*s*t8+t8*t6/170)*(t14+20.0/9.0*t13*t*t16+5.0/3.0*t13*t18+10.0/21.0*t*t18*  in Bn()
2633       t14 = t10*t10;  in Bn()
2638 /68.0*t4*s*t10+5.0/17.0*t4*t10*t7+15.0/476.0*s*t14+3.0/3094.0*t14*t7)*(t19+6.0/  in Bn()
2659       t14 = t12*t9;  in Bn()
2662 147.0/34.0*t6*s*t12+49.0/17.0*t6*t14+35.0/34.0*t5*t16+42.0/221.0*t4*t16*t9+7.0/  in Bn()
2663 442.0*s*t16*t12+t16*t14/2431);  in Bn()
2697       t14 = s+r-1.0;  in Bn()
2698       t16 = t14*t14;  in Bn()
2701 /35)*(t12+20.0/9.0*t11*t*t14+5.0/3.0*t16*t11+10.0/21.0*t*t16*t14+5.0/126.0*t20)  in Bn()
2710       t14 = t*t;  in Bn()
2714 6.0/35.0*s*t9+4.0/455.0*t9*t7)*(t14*t+12.0/7.0*t16*t14+6.0/7.0*t*t18+4.0/35.0*  in Bn()
2724       t14 = t10*t10;  in Bn()
2729 4.0/7.0*t5*t10+8.0/91.0*s*t8*t10+t14/273)*(t17+6.0/5.0*t*t18+3.0/10.0*t20);  in Bn()
2786       t14 = s+r-1.0;  in Bn()
2787       t16 = t14*t14;  in Bn()
2789 /13.0*s*t7+t9/26)*(t12*t+12.0/7.0*t12*t14+6.0/7.0*t*t16+4.0/35.0*t16*t14);  in Bn()
2857       t14 = s+r-1.0;  in Bn()
2858       t16 = t14*t14;  in Bn()
2860 1265.0*r+1.0/25300.0)*(t7+6.0/11.0*s*t8+3.0/55.0*t10)*(t13+6.0/5.0*t*t14+3.0/  in Bn()
3029       t14 = t10*t10;  in Bn()
3030       t17 = t14*t10;  in Bn()
3031       t20 = t14*t14;  in Bn()
3034 t14+3360.0/323.0*t7*t6*t17+2205.0/323.0*t7*t*t20+882.0/323.0*t7*t20*t10+210.0/  in Bn()
3035 323.0*t12*t20*t14+360.0/4199.0*t6*t20*t17+45.0/8398.0*t*t28+5.0/46189.0*t28*t10  in Bn()
3046       t14 = s+r-1.0;  in Bn()
3047       t17 = t14*t14;  in Bn()
3048       t20 = t17*t14;  in Bn()
3052 t12*t14+126.0/17.0*t10*t9*t17+588.0/85.0*t10*t*t20+63.0/17.0*t10*t22+252.0/  in Bn()
3053 221.0*t12*t22*t14+42.0/221.0*t9*t22*t17+36.0/2431.0*t*t22*t20+9.0/24310.0*t30);  in Bn()
3063       t14 = t12*t12;  in Bn()
3068       f = (t2+4.0/5.0*t1*s*t4+t1*t6/5+2.0/115.0*s*t6*t4+t10/2530)*(t13*t14+  in Bn()
3069 56.0/15.0*t14*t12*t17+28.0/5.0*t14*t*t20+56.0/13.0*t22*t14+70.0/39.0*t13*t24+  in Bn()
3127       t14 = t9*t9;  in Bn()
3128       t22 = t14*t14;  in Bn()
3133 882.0/1265.0*t2*t14+168.0/1265.0*t4*t14*t6+84.0/6325.0*t1*t14*t9+72.0/120175.0*  in Bn()
3134 s*t14*t12+t22/120175)*(t24*t+12.0/7.0*t26*t24+6.0/7.0*t*t28+4.0/35.0*t28*t26);  in Bn()
3253       t14 = t*t;  in Bn()
3254       t15 = t14*t14;  in Bn()
3259 t7*t5+t11/1771)*(t14*t15+42.0/13.0*t15*t*t18+105.0/26.0*t15*t20+350.0/143.0*t14  in Bn()
3260 *t*t20*t18+105.0/143.0*t14*t25+14.0/143.0*t*t25*t18+7.0/1716.0*t25*t20);  in Bn()
3299       t14 = t10*t10;  in Bn()
3304 253.0*t4*t12+70.0/253.0*t3*t14+168.0/4807.0*t2*t14*t7+28.0/14421.0*s*t14*t10+  in Bn()
3305 8.0/245157.0*t14*t12)*(t23*t+12.0/7.0*t25*t23+6.0/7.0*t*t27+4.0/35.0*t27*t25);  in Bn()
3362       t14 = t11*t7;  in Bn()
3366 +3360.0/323.0*t4*t3*t14+2205.0/323.0*t4*t*t17+882.0/323.0*t4*t17*t7+210.0/323.0  in Bn()
3367 *t9*t17*t11+360.0/4199.0*t3*t17*t14+45.0/8398.0*t*t25+5.0/46189.0*t25*t7);  in Bn()
3392       t14 = s+r-1.0;  in Bn()
3393       t17 = t14*t14;  in Bn()
3394       t19 = t17*t14;  in Bn()
3397 *t11*t9*t14+28.0/5.0*t11*t*t17+56.0/13.0*t11*t19+70.0/39.0*t10*t21+56.0/143.0*  in Bn()
3398 t9*t21*t14+28.0/715.0*t*t21*t17+8.0/6435.0*t21*t19);  in Bn()
3452       t14 = t9*t9;  in Bn()
3457 *t3*s*t9*t7+5.0/57.0*t3*t14+2.0/323.0*s*t14*t7+t14*t9/7752)*(t21*t+12.0/7.0*t23  in Bn()
3484       t14 = t11*t8;  in Bn()
3488 336.0/95.0*t4*s*t14+28.0/19.0*t4*t16+112.0/323.0*t6*t16*t8+14.0/323.0*t3*t16*  in Bn()
3489 t11+4.0/1615.0*s*t16*t14+t24/22610)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);  in Bn()
3499       t14 = t11*t7;  in Bn()
3503 144.0/19.0*t4*t3*t14+84.0/19.0*t4*s*t17+504.0/323.0*t4*t17*t7+105.0/323.0*t9*  in Bn()
3504 t17*t11+12.0/323.0*t3*t17*t14+9.0/4522.0*s*t25+t25*t7/29393);  in Bn()
3527       t14 = t11*t11;  in Bn()
3528       t16 = t14*t11;  in Bn()
3529       t18 = t14*t14;  in Bn()
3531 *t8+56.0/15.0*t8*t6*t11+28.0/5.0*t8*t*t14+56.0/13.0*t16*t8+70.0/39.0*t7*t18+  in Bn()
3532 56.0/143.0*t6*t18*t11+28.0/715.0*t*t18*t14+8.0/6435.0*t18*t16);  in Bn()
3541       t14 = s+r-1.0;  in Bn()
3542       t16 = t14*t14;  in Bn()
3545 t10*t11+42.0/13.0*t11*t*t14+105.0/26.0*t16*t11+350.0/143.0*t10*t*t16*t14+105.0/  in Bn()
3546 143.0*t10*t21+14.0/143.0*t*t21*t14+7.0/1716.0*t21*t16);  in Bn()
3554       t14 = t13*t13;  in Bn()
3559 19.0*s*t8+4.0/969.0*t8*t6)*(t14*t+30.0/11.0*t16*t14+30.0/11.0*t13*t*t19+40.0/  in Bn()
3584       t14 = t10*t10;  in Bn()
3589 57.0*t4*s*t10+200.0/969.0*t4*t10*t7+25.0/1292.0*s*t14+t14*t7/1938)*(t19*t+12.0/  in Bn()
3613       t14 = t12*t9;  in Bn()
3616 196.0/57.0*t6*s*t12+1960.0/969.0*t6*t14+1225.0/1938.0*t5*t16+98.0/969.0*t4*t16*  in Bn()
3617 t9+7.0/969.0*s*t16*t12+2.0/12597.0*t16*t14)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0);  in Bn()
3668       t14 = s+r-1.0;  in Bn()
3669       t17 = t14*t14;  in Bn()
3672 17.0*s*t6+3.0/136.0*t8)*(t12*t+30.0/11.0*t12*t14+30.0/11.0*t11*t*t17+40.0/33.0*  in Bn()
3673 t11*t17*t14+5.0/22.0*t*t21+t21*t14/77);  in Bn()
3681       t14 = t*t;  in Bn()
3682       t15 = t14*t14;  in Bn()
3687 /17.0*t5*t7+9.0/68.0*s*t9+t9*t7/170)*(t15+20.0/9.0*t14*t*t17+5.0/3.0*t14*t19+  in Bn()
3697       t14 = t10*t10;  in Bn()
3702 17.0*t8*t5*s+15.0/34.0*t5*t10+s*t10*t8/17+t14/476)*(t17*t+12.0/7.0*t19*t17+6.0/  in Bn()
3764       t14 = t11*t11;  in Bn()
3765       t18 = t14*t14;  in Bn()
3767 *t1+r/234-1.0/13455.0)*(t9*t+30.0/11.0*t11*t9+30.0/11.0*t8*t*t14+40.0/33.0*t8*  in Bn()
3768 t14*t11+5.0/22.0*t*t18+t18*t11/77);  in Bn()
3855       t14 = t12*t12;  in Bn()
3856       t18 = t14*t14;  in Bn()
3859 3.0*t9*t14+10.0/21.0*t*t14*t12+5.0/126.0*t18);  in Bn()
3930       t14 = t12*t12;  in Bn()
3933 t+12.0/7.0*t12*t10+6.0/7.0*t*t14+4.0/35.0*t14*t12);  in Bn()
3942       t14 = t*t;  in Bn()
3947 t11)*(t14+6.0/5.0*t*t15+3.0/10.0*t17);  in Bn()
3992       t14 = t12*t12;  in Bn()
3995 /9.0)*(t11+6.0/5.0*t*t12+3.0/10.0*t14);  in Bn()