19    double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t18;  in En()  local
41       t2 = t1*t1;  in En()
44       f = t2-10.0*r*t1*s+20.0*t5*t1-10.0*t5*r*s+t9;  in En()
48       t2 = t1*t1;  in En()
51       f = t2*s-15.0*t2*r+50.0*t5*t1*s-50.0*t5*r*t1+15.0*t10*s-t10*r;  in En()
55       t2 = t1*t1;  in En()
58       f = t2*t1-21.0*r*t2*s+105.0*t6*t2-175.0*t6*r*t1*s+105.0*t11*t1-21.0*t11  in En()
63       t2 = t1*s;  in En()
68       f = t3*t2-28.0*r*t3*t1+196.0*t7*t3*s-490.0*t10*t3+490.0*t12*t2-196.0*  in En()
73       t2 = t1*t1;  in En()
74       t3 = t2*t2;  in En()
80       f = t3-36.0*r*t2*t4+336.0*t7*t2*t1-1176.0*t10*t2*s+1764.0*t13*t2-1176.0  in En()
85       t2 = t1*t1;  in En()
86       t3 = t2*t2;  in En()
92       f = t3*s-45.0*t3*r+540.0*t6*t2*t7-2520.0*t10*t2*t1+5292.0*t13*t2*s  in En()
93 -5292.0*t13*r*t2+2520.0*t13*t6*t7-540.0*t13*t10*t1+45.0*t22*s-t22*r;  in En()
97       t2 = t1*t1;  in En()
98       t3 = t2*t2;  in En()
104       f = t3*t1-55.0*r*t3*s+825.0*t7*t3-4950.0*t9*t2*t10+13860.0*t13*t2*t1  in En()
105 -19404.0*t13*r*t2*s+13860.0*t13*t7*t2-4950.0*t13*t9*t10+825.0*t23*t1-55.0*t23*r  in En()
110       t2 = t1*s;  in En()
117       f = t2*t4-66.0*r*t4*t1+1210.0*t8*t4*s-9075.0*t11*t4+32670.0*t13*t3*t2  in En()
118 -60984.0*t13*r*t3*t1+60984.0*t13*t8*t3*s-32670.0*t13*t11*t3+9075.0*t24*t2  in En()
123       t2 = t1*t1;  in En()
124       t3 = t2*t2;  in En()
130       f = t3*t2-78.0*r*t3*t5+1716.0*t8*t3*t1-15730.0*t11*t3*s+70785.0*t14*t3  in En()
131 -169884.0*t14*r*t2*t5+226512.0*t14*t8*t2*t1-169884.0*t14*t11*t2*s+70785.0*t25*  in En()
132 t2-15730.0*t25*r*t5+1716.0*t25*t8*t1-78.0*t25*t11*s+t25*t14;  in En()
136       t2 = t1*t1;  in En()
137       t3 = t2*s;  in En()
138       t4 = t2*t2;  in En()
145       f = t4*t3-91.0*r*t2*t4+2366.0*t8*t4*t9-26026.0*t12*t4*t1+143143.0*t15*  in En()
146 t4*s-429429.0*t18*t4+736164.0*t15*t8*t2*t9-736164.0*t15*t12*t2*t1+429429.0*t26*  in En()
147 t3-143143.0*t26*r*t2+26026.0*t26*t8*t9-2366.0*t26*t12*t1+91.0*t26*t15*s-t26*t18  in En()