18    double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18;  in EnDrv()  local
37       t3 = r*r;  in EnDrv()
38       drv[0] = -6.0*t1+12.0*t2-3.0*t3;  in EnDrv()
39       drv[1] = 3.0*t1-12.0*t2+6.0*t3;  in EnDrv()
44       t3 = r*t1;  in EnDrv()
48       drv[0] = -10.0*t2+40.0*t3-30.0*t5+4.0*t6;  in EnDrv()
49       drv[1] = 4.0*t2-30.0*t3+40.0*t5-10.0*t6;  in EnDrv()
65       t3 = t2*s;  in EnDrv()
73       drv[0] = -21.0*t3+210.0*t4-525.0*t7+420.0*t9-105.0*t11+6.0*t12;  in EnDrv()
74       drv[1] = 6.0*t3-105.0*t4+420.0*t7-525.0*t9+210.0*t11-21.0*t12;  in EnDrv()
79       t3 = t2*t1;  in EnDrv()
88       drv[0] = -28.0*t3+392.0*t5-1470.0*t7+1960.0*t10-980.0*t12+168.0*t14-7.0*  in EnDrv()
90       drv[1] = 7.0*t3-168.0*t5+980.0*t7-1960.0*t10+1470.0*t12-392.0*t14+28.0*  in EnDrv()
96       t3 = t1*t1;  in EnDrv()
97       t4 = t3*t2;  in EnDrv()
98       t6 = r*t3*t1;  in EnDrv()
100       t9 = t7*t3*s;  in EnDrv()
102       t11 = t10*t3;  in EnDrv()
116       t3 = t2*t2;  in EnDrv()
129       drv[0] = -45.0*t3+1080.0*t6-7560.0*t9+21168.0*t12-26460.0*t14+15120.0*t16  in EnDrv()
131       drv[1] = 9.0*t3-360.0*t6+3780.0*t9-15120.0*t12+26460.0*t14-21168.0*t16+  in EnDrv()
137       t3 = t2*t2;  in EnDrv()
138       t4 = t3*s;  in EnDrv()
139       t5 = r*t3;  in EnDrv()
161       t3 = t2*t2;  in EnDrv()
162       t4 = t3*t1;  in EnDrv()
163       t6 = r*t3*s;  in EnDrv()
165       t8 = t7*t3;  in EnDrv()
188       t3 = t1*t1;  in EnDrv()
189       t4 = t3*t3;  in EnDrv()
197       t15 = t13*t3*t2;  in EnDrv()
198       t18 = t13*r*t3*t1;  in EnDrv()
199       t21 = t13*t8*t3*s;  in EnDrv()
200       t23 = t13*t11*t3;  in EnDrv()
216       t3 = t2*t2;  in EnDrv()
217       t4 = t3*t2;  in EnDrv()
219       t7 = r*t3*t5;  in EnDrv()
221       t10 = t8*t3*t1;  in EnDrv()
223       t13 = t11*t3*s;  in EnDrv()
225       t15 = t14*t3;  in EnDrv()