19    double f = 0.0;  in Fn()  local
26       f = 1.0;  in Fn()
29       f = s-1.0/3.0;  in Fn()
32       f = r-1.0/3.0;  in Fn()
36       f = t1-3.0/4.0*s+3.0/28.0;  in Fn()
39       f = r*s-r/4-s/4+1.0/14.0;  in Fn()
43       f = t1-3.0/4.0*r+3.0/28.0;  in Fn()
47       f = t1*s-6.0/5.0*t1+2.0/5.0*s-1.0/30.0;  in Fn()
51       f = r*t1-3.0/5.0*r*s+r/15-t1/5+2.0/15.0*s-1.0/60.0;  in Fn()
55       f = t1*s-t1/5-3.0/5.0*r*s+2.0/15.0*r+s/15-1.0/60.0;  in Fn()
59       f = r*t1-6.0/5.0*t1+2.0/5.0*r-1.0/30.0;  in Fn()
64       f = t2-5.0/3.0*t1*s+10.0/11.0*t1-2.0/11.0*s+1.0/99.0;  in Fn()
69       f = r*t2-r*t1+3.0/11.0*r*s-r/55-t2/6+2.0/11.0*t1-3.0/55.0*s+2.0/495.0;  in Fn()
74       f = t1*t2-t1*s/2+t1/22-r*t2/2+3.0/11.0*r*s-3.0/110.0*r+t2/22-3.0/110.0*s  in Fn()
80       f = t2*s-t2/6-t1*s+2.0/11.0*t1+3.0/11.0*r*s-3.0/55.0*r-s/55+2.0/495.0;  in Fn()
85       f = t2-5.0/3.0*r*t1+10.0/11.0*t1-2.0/11.0*r+1.0/99.0;  in Fn()
90       f = t2*s-15.0/7.0*t2+150.0/91.0*t1*s-50.0/91.0*t1+75.0/1001.0*s-3.0/  in Fn()
97       f = r*t2-10.0/7.0*r*t4+60.0/91.0*r*t1-10.0/91.0*r*s+5.0/1001.0*r-t2/7+  in Fn()
104       f = t1*t3-6.0/7.0*t1*t2+18.0/91.0*t1*s-t1/91-3.0/7.0*r*t3+36.0/91.0*r*  in Fn()
111       f = t2*t3-3.0/7.0*t2*s+3.0/91.0*t2-6.0/7.0*t1*t3+36.0/91.0*t1*s-3.0/  in Fn()
118       f = t2*s-t2/7-10.0/7.0*t4*s+20.0/91.0*t4+60.0/91.0*t1*s-10.0/91.0*t1  in Fn()
124       f = r*t2-15.0/7.0*t2+150.0/91.0*r*t1-50.0/91.0*t1+75.0/1001.0*r-3.0/  in Fn()
130       f = t1*t2-21.0/8.0*t2*s+21.0/8.0*t2-5.0/4.0*t1*s+15.0/52.0*t1-3.0/104.0*  in Fn()
138       f = r*t3-15.0/8.0*r*t2+5.0/4.0*r*t6-5.0/14.0*r*t1+15.0/364.0*r*s-r/728-  in Fn()
146       f = t1*t3-5.0/4.0*t1*t5+t1*t2/2-t1*s/14+t1/364-3.0/8.0*r*t3+r*t5/2-3.0/  in Fn()
154       f = t2*t4-3.0/4.0*t2*t3+3.0/20.0*t2*s-t2/140-3.0/4.0*t1*t4+3.0/5.0*t1*  in Fn()
163       f = t2*t3-3.0/8.0*t2*s+t2/40-5.0/4.0*t6*t3+t6*s/2-t6/28+t1*t3/2-3.0/  in Fn()
171       f = t3*s-t3/8-15.0/8.0*t2*s+t2/4+5.0/4.0*t6*s-5.0/28.0*t6-5.0/14.0*t1*s  in Fn()
177       f = t1*t2-21.0/8.0*r*t2+21.0/8.0*t2-5.0/4.0*r*t1+15.0/52.0*t1-3.0/104.0*  in Fn()
184       f = t2*t3-28.0/9.0*t1*t3+196.0/51.0*t3*s-245.0/102.0*t3+245.0/306.0*t2  in Fn()
193       f = r*t3-7.0/3.0*r*t5+35.0/17.0*r*t2-175.0/204.0*r*t8+35.0/204.0*r*t1-r  in Fn()
203       f = t1*t4-5.0/3.0*t1*t3+50.0/51.0*t1*t7-25.0/102.0*t1*t2+5.0/204.0*t1*s  in Fn()
213       f = t2*t4-10.0/9.0*t2*t6+20.0/51.0*t2*t3-5.0/102.0*t2*s+t2/612-2.0/3.0*  in Fn()
223       f = t2*t4-2.0/3.0*t2*t3+2.0/17.0*t2*s-t2/204-10.0/9.0*t8*t4+40.0/51.0*  in Fn()
233       f = t3*t4-t3*s/3+t3/51-5.0/3.0*t2*t4+10.0/17.0*t2*s-5.0/136.0*t2+50.0/  in Fn()
243       f = t3*s-t3/9-7.0/3.0*t5*s+14.0/51.0*t5+35.0/17.0*t2*s-35.0/136.0*t2  in Fn()
251       f = t2*t3-28.0/9.0*t1*t3+196.0/51.0*r*t3-245.0/102.0*t3+245.0/306.0*t2  in Fn()
259       f = t3-18.0/5.0*t2*t4+504.0/95.0*t1*t2-392.0/95.0*t2*s+588.0/323.0*t2  in Fn()
269       f = r*t4-14.0/5.0*r*t6+294.0/95.0*r*t8-98.0/57.0*r*t3+490.0/969.0*r*t2  in Fn()
283       f = t1*t4-21.0/10.0*t1*t6+63.0/38.0*t1*t3-35.0/57.0*t1*t9+35.0/323.0*t1  in Fn()
300       f = t17+t23;  in Fn()
314       f = t12+t24;  in Fn()
329       f = t16+t23;  in Fn()
341       f = t3*t4-3.0/10.0*t3*s+3.0/190.0*t3-21.0/10.0*t7*t4+63.0/95.0*t7*s-7.0  in Fn()
351       f = t4*s-t4/10-14.0/5.0*t6*s+28.0/95.0*t6+294.0/95.0*t8*s-98.0/285.0*t8  in Fn()
360       f = t3-18.0/5.0*t2*t4+504.0/95.0*t2*t1-392.0/95.0*r*t2+588.0/323.0*t2  in Fn()
368       f = t3*s-45.0/11.0*t3+540.0/77.0*t2*t5-72.0/11.0*t2*t1+756.0/209.0*t2*s  in Fn()
379       f = r*t3-36.0/11.0*r*t6+48.0/11.0*r*t8-168.0/55.0*r*t10+252.0/209.0*r*  in Fn()
398       f = t19+t23;  in Fn()
414       f = t21+t27;  in Fn()
431       f = t22+t29;  in Fn()
448       f = t21+t29;  in Fn()
464       f = t15+t27;  in Fn()
480       f = t16+t23;  in Fn()
490       f = t3*s-t3/11-36.0/11.0*t6*s+24.0/77.0*t6+48.0/11.0*t8*s-24.0/55.0*t8  in Fn()
508       f = t3*t1-55.0/12.0*t3*s+825.0/92.0*t3-225.0/23.0*t2*t6+150.0/23.0*t2*  in Fn()
521       f = r*t4-15.0/4.0*r*t3+135.0/23.0*r*t8-1260.0/253.0*r*t10+630.0/253.0*r  in Fn()
541       f = t19+t26;  in Fn()
560       f = t14+t31;  in Fn()
579       f = t15+t34;  in Fn()
600       f = t26+t35;  in Fn()
619       f = t25+t34;  in Fn()
637       f = t16+t31;  in Fn()
654       f = t17+t26;  in Fn()
665       f = t4*s-t4/12-15.0/4.0*t3*s+15.0/46.0*t3+135.0/23.0*t8*s-135.0/253.0*  in Fn()
675       f = t3*t1-55.0/12.0*t3*r+825.0/92.0*t3-225.0/23.0*t2*t6+150.0/23.0*t2*  in Fn()
684       f = t4*t2-66.0/13.0*t4*t1+726.0/65.0*t4*s-363.0/26.0*t4+3267.0/299.0*t3  in Fn()
704       f = t20+t21;  in Fn()
724       f = t22+t29;  in Fn()
746       f = t25+t35;  in Fn()
768       f = t24+t39;  in Fn()
786       f = 3.0/124982.0*t2-3.0/32890.0*t5+3.0/16445.0*t6-9.0/65780.0*t4-9.0/  in Fn()
808       f = -315.0/6578.0*t3*t1-15.0/3289.0*t1*t2-63.0/6578.0*t7*s-63.0/32890.0  in Fn()
834       f = t24+t39;  in Fn()
856       f = t22+t35;  in Fn()
876       f = t24+t29;  in Fn()
894       f = t16+t21;  in Fn()
901       f = t4*t2-66.0/13.0*t4*t1+726.0/65.0*r*t4-363.0/26.0*t4+3267.0/299.0*t3  in Fn()
910       f = t3*t2-39.0/7.0*t3*t5+286.0/21.0*t3*t1-1210.0/63.0*t3*s+121.0/7.0*t3  in Fn()
931       f = t22+t23;  in Fn()
952       f = t25+t32;  in Fn()
975       f = t28+t39;  in Fn()
1000       f = t17+t31+t37+t47;  in Fn()
1025       f = t22+t33+t38+t50;  in Fn()
1050       f = t21+t34+t42+t51;  in Fn()
1075       f = t20+t33+t47+t50;  in Fn()
1100       f = t13+t22+t34+t47;  in Fn()
1122       f = t30+t39;  in Fn()
1142       f = t15+t32;  in Fn()
1161       f = t16+t23;  in Fn()
1168       f = t3*t2-39.0/7.0*t3*t5+286.0/21.0*t3*t1-1210.0/63.0*t3*r+121.0/7.0*t3  in Fn()
1172     return f;  in Fn()