17 int EnDrv(int ip, double r, double s, double drv[2]) {  in EnDrv()
23       drv[0] = 0.0;  in EnDrv()
24       drv[1] = 0.0;  in EnDrv()
27       drv[0] = -1.0;  in EnDrv()
28       drv[1] = 1.0;  in EnDrv()
31       drv[0] = -3.0*s+2.0*r;  in EnDrv()
32       drv[1] = 2.0*s-3.0*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()
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()
59       drv[0] = -15.0*t2+100.0*t4-150.0*t6+60.0*t8-5.0*t9;  in EnDrv()
60       drv[1] = 5.0*t2-60.0*t4+150.0*t6-100.0*t8+15.0*t9;  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()
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()
108       drv[0] = -36.0*t4+672.0*t6-3528.0*t9+7056.0*t11-5880.0*t13+2016.0*t15  in EnDrv()
110       drv[1] = 8.0*t4-252.0*t6+2016.0*t9-5880.0*t11+7056.0*t13-3528.0*t15+672.0  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()
153       drv[0] = -55.0*t4+1650.0*t5-14850.0*t9+55440.0*t12-97020.0*t15+83160.0*  in EnDrv()
155       drv[1] = 10.0*t4-495.0*t5+6600.0*t9-34650.0*t12+83160.0*t15-97020.0*t17+  in EnDrv()
182       drv[0] = t28;  in EnDrv()
183       drv[1] = t29;  in EnDrv()
210       drv[0] = t31;  in EnDrv()
211       drv[1] = t32;  in EnDrv()
241       drv[0] = t34;  in EnDrv()
242       drv[1] = t35;  in EnDrv()