17 double En(int ip, double r, double s) {  in En()  argument
27       f = s-r;  in En()
31       t3 = r*r;  in En()
32       f = t1-3.0*r*s+t3;  in En()
36       t4 = r*r;  in En()
37       f = t1*s-6.0*r*t1+6.0*t4*s-t4*r;  in En()
42       t5 = r*r;  in En()
44       f = t2-10.0*r*t1*s+20.0*t5*t1-10.0*t5*r*s+t9;  in En()
49       t5 = r*r;  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()
56       t6 = r*r;  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()
59 *r*s+t11*t6;  in En()
65       t7 = r*r;  in En()
66       t10 = t7*r;  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()
69 t12*r*t1+28.0*t12*t7*s-t12*t10;  in En()
76       t7 = r*r;  in En()
77       t10 = t7*r;  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()
81 *t13*r*t4+336.0*t13*t7*t1-36.0*t13*t10*s+t21;  in En()
87       t6 = r*r;  in En()
89       t10 = t6*r;  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()
99       t7 = r*r;  in En()
100       t9 = t7*r;  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()
113       t8 = r*r;  in En()
114       t11 = t8*r;  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()
119 -1210.0*t24*r*t1+66.0*t24*t8*s-t24*t11;  in En()
126       t8 = r*r;  in En()
127       t11 = t8*r;  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()
139       t8 = r*r;  in En()
141       t12 = t8*r;  in En()
143       t18 = t15*r;  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()
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()