1 /*------------------------------------------------------------------------- 2 Scientific Computation Research Center, RPI, Troy NY 3 (C) Copyright 1995, RPI-SCOREC 4 5 Project : shapeFuntions 6 Author(s): Saikat Dey 7 Creation : Oct., 95 8 Modifi. : 9 Function : 10 return the "body" mode for a tetrahedron. i,j,k are the highest 11 monomial orders of r,s,t in the resulting polynomial order 12 NOTE: code generated by Maple. 13 -------------------------------------------------------------------------*/ 14 15 #ifdef __cplusplus 16 extern "C" { 17 #endif 18 19 double Bn(int i, int j, int k, double r, double s, double t) { 20 double f = 0.0; 21 double t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18; 22 double t19,t20,t21,t22,t23,t24,t25,t26,t27,t28,t29,t30,t31,t32,t33,t34; 23 double t35,t36,t37,t38,t39,t40,t41,t42,t43,t44,t45; 24 25 /* p=4 0 0 0 */ 26 if( i==0 && j==0 && k==0 ) { 27 f = 1.0; 28 /* p=5 0 0 1 */ 29 } else if( i==0 && j==0 && k==1 ) { 30 f = t+2.0/3.0*s+2.0/3.0*r-2.0/3.0; 31 /* p=5 0 1 0 */ 32 } else if( i==0 && j==1 && k==0 ) { 33 f = s+2.0/5.0*r-2.0/5.0; 34 /* p=5 1 0 0 */ 35 } else if( i==1 && j==0 && k==0 ) { 36 f = r-2.0/7.0; 37 /* p=6 0 0 2 */ 38 } else if( i==0 && j==0 && k==2 ) { 39 t1 = t*t; 40 t2 = s+r-1.0; 41 t4 = t2*t2; 42 f = t1+6.0/5.0*t*t2+3.0/10.0*t4; 43 /* p=6 0 1 1 */ 44 } else if( i==0 && j==1 && k==1 ) { 45 f = (s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 46 /* p=6 0 2 0 */ 47 } else if( i==0 && j==2 && k==0 ) { 48 t1 = s*s; 49 t2 = r-1.0; 50 t4 = t2*t2; 51 f = t1+6.0/7.0*t2*s+t4/7; 52 /* p=6 1 0 1 */ 53 } else if( i==1 && j==0 && k==1 ) { 54 f = (r-2.0/9.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 55 /* p=6 1 1 0 */ 56 } else if( i==1 && j==1 && k==0 ) { 57 f = (r-2.0/9.0)*(s+2.0/5.0*r-2.0/5.0); 58 /* p=6 2 0 0 */ 59 } else if( i==2 && j==0 && k==0 ) { 60 t1 = r*r; 61 f = t1-2.0/3.0*r+1.0/12.0; 62 /* p=7 0 0 3 */ 63 } else if( i==0 && j==0 && k==3 ) { 64 t1 = t*t; 65 t3 = s+r-1.0; 66 t5 = t3*t3; 67 f = t1*t+12.0/7.0*t3*t1+6.0/7.0*t*t5+4.0/35.0*t3*t5; 68 /* p=7 0 1 2 */ 69 } else if( i==0 && j==1 && k==2 ) { 70 t2 = t*t; 71 t3 = s+r-1.0; 72 t5 = t3*t3; 73 f = (s+2.0/9.0*r-2.0/9.0)*(t2+6.0/5.0*t*t3+3.0/10.0*t5); 74 /* p=7 0 2 1 */ 75 } else if( i==0 && j==2 && k==1 ) { 76 t1 = s*s; 77 t2 = r-1.0; 78 t4 = t2*t2; 79 f = (t1+2.0/3.0*s*t2+t4/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 80 /* p=7 0 3 0 */ 81 } else if( i==0 && j==3 && k==0 ) { 82 t1 = s*s; 83 t3 = r-1.0; 84 t5 = t3*t3; 85 f = t1*s+4.0/3.0*t1*t3+s*t5/2+t5*t3/21; 86 /* p=7 1 0 2 */ 87 } else if( i==1 && j==0 && k==2 ) { 88 t2 = t*t; 89 t3 = s+r-1.0; 90 t5 = t3*t3; 91 f = (r-2.0/11.0)*(t2+6.0/5.0*t*t3+3.0/10.0*t5); 92 /* p=7 1 1 1 */ 93 } else if( i==1 && j==1 && k==1 ) { 94 f = (r-2.0/11.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 95 /* p=7 1 2 0 */ 96 } else if( i==1 && j==2 && k==0 ) { 97 t2 = s*s; 98 t3 = r-1.0; 99 t5 = t3*t3; 100 f = (r-2.0/11.0)*(t2+6.0/7.0*s*t3+t5/7); 101 /* p=7 2 0 1 */ 102 } else if( i==2 && j==0 && k==1 ) { 103 t1 = r*r; 104 f = (t1-6.0/11.0*r+3.0/55.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 105 /* p=7 2 1 0 */ 106 } else if( i==2 && j==1 && k==0 ) { 107 t1 = r*r; 108 f = (t1-6.0/11.0*r+3.0/55.0)*(s+2.0/5.0*r-2.0/5.0); 109 /* p=7 3 0 0 */ 110 } else if( i==3 && j==0 && k==0 ) { 111 t1 = r*r; 112 f = t1*r-12.0/11.0*t1+18.0/55.0*r-4.0/165.0; 113 /* p=8 0 0 4 */ 114 } else if( i==0 && j==0 && k==4 ) { 115 t1 = t*t; 116 t2 = t1*t1; 117 t4 = s+r-1.0; 118 t6 = t4*t4; 119 t10 = t6*t6; 120 f = t2+20.0/9.0*t1*t*t4+5.0/3.0*t1*t6+10.0/21.0*t*t6*t4+5.0/126.0*t10; 121 /* p=8 0 1 3 */ 122 } else if( i==0 && j==1 && k==3 ) { 123 t2 = t*t; 124 t4 = s+r-1.0; 125 t6 = t4*t4; 126 f = (s+2.0/11.0*r-2.0/11.0)*(t2*t+12.0/7.0*t2*t4+6.0/7.0*t*t6+4.0/35.0* 127 t6*t4); 128 /* p=8 0 2 2 */ 129 } else if( i==0 && j==2 && k==2 ) { 130 t1 = s*s; 131 t2 = r-1.0; 132 t4 = t2*t2; 133 t6 = t*t; 134 t7 = s+r-1.0; 135 t9 = t7*t7; 136 f = (t1+6.0/11.0*s*t2+3.0/55.0*t4)*(t6+6.0/5.0*t*t7+3.0/10.0*t9); 137 /* p=8 0 3 1 */ 138 } else if( i==0 && j==3 && k==1 ) { 139 t1 = s*s; 140 t3 = r-1.0; 141 t5 = t3*t3; 142 f = (t1*s+12.0/11.0*t1*t3+18.0/55.0*s*t5+4.0/165.0*t5*t3)*(t+2.0/3.0*s+ 143 2.0/3.0*r-2.0/3.0); 144 /* p=8 0 4 0 */ 145 } else if( i==0 && j==4 && k==0 ) { 146 t1 = s*s; 147 t2 = t1*t1; 148 t4 = r-1.0; 149 t6 = t4*t4; 150 t10 = t6*t6; 151 f = t2+20.0/11.0*t1*s*t4+12.0/11.0*t1*t6+8.0/33.0*s*t6*t4+t10/66; 152 /* p=8 1 0 3 */ 153 } else if( i==1 && j==0 && k==3 ) { 154 t2 = t*t; 155 t4 = s+r-1.0; 156 t6 = t4*t4; 157 f = (r-2.0/13.0)*(t2*t+12.0/7.0*t2*t4+6.0/7.0*t*t6+4.0/35.0*t6*t4); 158 /* p=8 1 1 2 */ 159 } else if( i==1 && j==1 && k==2 ) { 160 t4 = t*t; 161 t5 = s+r-1.0; 162 t7 = t5*t5; 163 f = (r-2.0/13.0)*(s+2.0/9.0*r-2.0/9.0)*(t4+6.0/5.0*t*t5+3.0/10.0*t7); 164 /* p=8 1 2 1 */ 165 } else if( i==1 && j==2 && k==1 ) { 166 t2 = s*s; 167 t3 = r-1.0; 168 t5 = t3*t3; 169 f = (r-2.0/13.0)*(t2+2.0/3.0*s*t3+t5/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 170 /* p=8 1 3 0 */ 171 } else if( i==1 && j==3 && k==0 ) { 172 t2 = s*s; 173 t4 = r-1.0; 174 t6 = t4*t4; 175 f = (r-2.0/13.0)*(t2*s+4.0/3.0*t2*t4+s*t6/2+t6*t4/21); 176 /* p=8 2 0 2 */ 177 } else if( i==2 && j==0 && k==2 ) { 178 t1 = r*r; 179 t3 = t*t; 180 t4 = s+r-1.0; 181 t6 = t4*t4; 182 f = (t1-6.0/13.0*r+1.0/26.0)*(t3+6.0/5.0*t*t4+3.0/10.0*t6); 183 /* p=8 2 1 1 */ 184 } else if( i==2 && j==1 && k==1 ) { 185 t1 = r*r; 186 f = (t1-6.0/13.0*r+1.0/26.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0* 187 r-2.0/3.0); 188 /* p=8 2 2 0 */ 189 } else if( i==2 && j==2 && k==0 ) { 190 t1 = r*r; 191 t3 = s*s; 192 t4 = r-1.0; 193 t6 = t4*t4; 194 f = (t1-6.0/13.0*r+1.0/26.0)*(t3+6.0/7.0*s*t4+t6/7); 195 /* p=8 3 0 1 */ 196 } else if( i==3 && j==0 && k==1 ) { 197 t1 = r*r; 198 f = (t1*r-12.0/13.0*t1+3.0/13.0*r-2.0/143.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/ 199 3.0); 200 /* p=8 3 1 0 */ 201 } else if( i==3 && j==1 && k==0 ) { 202 t1 = r*r; 203 f = (t1*r-12.0/13.0*t1+3.0/13.0*r-2.0/143.0)*(s+2.0/5.0*r-2.0/5.0); 204 /* p=8 4 0 0 */ 205 } else if( i==4 && j==0 && k==0 ) { 206 t1 = r*r; 207 t2 = t1*t1; 208 f = t2-20.0/13.0*t1*r+10.0/13.0*t1-20.0/143.0*r+1.0/143.0; 209 /* p=9 0 0 5 */ 210 } else if( i==0 && j==0 && k==5 ) { 211 t1 = t*t; 212 t2 = t1*t1; 213 t4 = s+r-1.0; 214 t7 = t4*t4; 215 t11 = t7*t7; 216 f = t2*t+30.0/11.0*t2*t4+30.0/11.0*t1*t*t7+40.0/33.0*t1*t7*t4+5.0/22.0* 217 t*t11+t11*t4/77; 218 /* p=9 0 1 4 */ 219 } else if( i==0 && j==1 && k==4 ) { 220 t2 = t*t; 221 t3 = t2*t2; 222 t5 = s+r-1.0; 223 t7 = t5*t5; 224 t11 = t7*t7; 225 f = (s+2.0/13.0*r-2.0/13.0)*(t3+20.0/9.0*t2*t*t5+5.0/3.0*t2*t7+10.0/ 226 21.0*t*t7*t5+5.0/126.0*t11); 227 /* p=9 0 2 3 */ 228 } else if( i==0 && j==2 && k==3 ) { 229 t1 = s*s; 230 t2 = r-1.0; 231 t4 = t2*t2; 232 t6 = t*t; 233 t8 = s+r-1.0; 234 t10 = t8*t8; 235 f = (t1+6.0/13.0*t2*s+t4/26)*(t6*t+12.0/7.0*t6*t8+6.0/7.0*t*t10+4.0/ 236 35.0*t10*t8); 237 /* p=9 0 3 2 */ 238 } else if( i==0 && j==3 && k==2 ) { 239 t1 = s*s; 240 t3 = r-1.0; 241 t5 = t3*t3; 242 t9 = t*t; 243 t10 = s+r-1.0; 244 t12 = t10*t10; 245 f = (t1*s+12.0/13.0*t1*t3+3.0/13.0*s*t5+2.0/143.0*t5*t3)*(t9+6.0/5.0*t* 246 t10+3.0/10.0*t12); 247 /* p=9 0 4 1 */ 248 } else if( i==0 && j==4 && k==1 ) { 249 t1 = s*s; 250 t2 = t1*t1; 251 t4 = r-1.0; 252 t6 = t4*t4; 253 t10 = t6*t6; 254 f = (t2+20.0/13.0*t1*s*t4+10.0/13.0*t1*t6+20.0/143.0*s*t6*t4+t10/143)*( 255 t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 256 /* p=9 0 5 0 */ 257 } else if( i==0 && j==5 && k==0 ) { 258 t1 = s*s; 259 t2 = t1*t1; 260 t4 = r-1.0; 261 t7 = t4*t4; 262 t11 = t7*t7; 263 f = t2*s+30.0/13.0*t2*t4+25.0/13.0*t1*s*t7+100.0/143.0*t1*t7*t4+15.0/ 264 143.0*s*t11+2.0/429.0*t11*t4; 265 /* p=9 1 0 4 */ 266 } else if( i==1 && j==0 && k==4 ) { 267 t2 = t*t; 268 t3 = t2*t2; 269 t5 = s+r-1.0; 270 t7 = t5*t5; 271 t11 = t7*t7; 272 f = (r-2.0/15.0)*(t3+20.0/9.0*t2*t*t5+5.0/3.0*t2*t7+10.0/21.0*t*t7*t5+ 273 5.0/126.0*t11); 274 /* p=9 1 1 3 */ 275 } else if( i==1 && j==1 && k==3 ) { 276 t4 = t*t; 277 t6 = s+r-1.0; 278 t8 = t6*t6; 279 f = (r-2.0/15.0)*(s+2.0/11.0*r-2.0/11.0)*(t4*t+12.0/7.0*t4*t6+6.0/7.0*t 280 *t8+4.0/35.0*t8*t6); 281 /* p=9 1 2 2 */ 282 } else if( i==1 && j==2 && k==2 ) { 283 t2 = s*s; 284 t3 = r-1.0; 285 t5 = t3*t3; 286 t8 = t*t; 287 t9 = s+r-1.0; 288 t11 = t9*t9; 289 f = (r-2.0/15.0)*(t2+6.0/11.0*s*t3+3.0/55.0*t5)*(t8+6.0/5.0*t*t9+3.0/ 290 10.0*t11); 291 /* p=9 1 3 1 */ 292 } else if( i==1 && j==3 && k==1 ) { 293 t2 = s*s; 294 t4 = r-1.0; 295 t6 = t4*t4; 296 f = (r-2.0/15.0)*(t2*s+12.0/11.0*t2*t4+18.0/55.0*s*t6+4.0/165.0*t4*t6)* 297 (t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 298 /* p=9 1 4 0 */ 299 } else if( i==1 && j==4 && k==0 ) { 300 t2 = s*s; 301 t3 = t2*t2; 302 t5 = r-1.0; 303 t7 = t5*t5; 304 t11 = t7*t7; 305 f = (r-2.0/15.0)*(t3+20.0/11.0*t2*s*t5+12.0/11.0*t2*t7+8.0/33.0*s*t7*t5 306 +t11/66); 307 /* p=9 2 0 3 */ 308 } else if( i==2 && j==0 && k==3 ) { 309 t1 = r*r; 310 t3 = t*t; 311 t5 = s+r-1.0; 312 t7 = t5*t5; 313 f = (t1-2.0/5.0*r+1.0/35.0)*(t3*t+12.0/7.0*t3*t5+6.0/7.0*t*t7+4.0/35.0* 314 t7*t5); 315 /* p=9 2 1 2 */ 316 } else if( i==2 && j==1 && k==2 ) { 317 t1 = r*r; 318 t5 = t*t; 319 t6 = s+r-1.0; 320 t8 = t6*t6; 321 f = (t1-2.0/5.0*r+1.0/35.0)*(s+2.0/9.0*r-2.0/9.0)*(t5+6.0/5.0*t*t6+3.0/ 322 10.0*t8); 323 /* p=9 2 2 1 */ 324 } else if( i==2 && j==2 && k==1 ) { 325 t1 = r*r; 326 t3 = s*s; 327 t4 = r-1.0; 328 t6 = t4*t4; 329 f = (t1-2.0/5.0*r+1.0/35.0)*(t3+2.0/3.0*s*t4+t6/12)*(t+2.0/3.0*s+2.0/ 330 3.0*r-2.0/3.0); 331 /* p=9 2 3 0 */ 332 } else if( i==2 && j==3 && k==0 ) { 333 t1 = r*r; 334 t3 = s*s; 335 t5 = r-1.0; 336 t7 = t5*t5; 337 f = (t1-2.0/5.0*r+1.0/35.0)*(s*t3+4.0/3.0*t3*t5+s*t7/2+t7*t5/21); 338 /* p=9 3 0 2 */ 339 } else if( i==3 && j==0 && k==2 ) { 340 t1 = r*r; 341 t4 = t*t; 342 t5 = s+r-1.0; 343 t7 = t5*t5; 344 f = (t1*r-4.0/5.0*t1+6.0/35.0*r-4.0/455.0)*(t4+6.0/5.0*t*t5+3.0/10.0*t7) 345 ; 346 /* p=9 3 1 1 */ 347 } else if( i==3 && j==1 && k==1 ) { 348 t1 = r*r; 349 f = (t1*r-4.0/5.0*t1+6.0/35.0*r-4.0/455.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/ 350 3.0*s+2.0/3.0*r-2.0/3.0); 351 /* p=9 3 2 0 */ 352 } else if( i==3 && j==2 && k==0 ) { 353 t1 = r*r; 354 t4 = s*s; 355 t5 = r-1.0; 356 t7 = t5*t5; 357 f = (t1*r-4.0/5.0*t1+6.0/35.0*r-4.0/455.0)*(t4+6.0/7.0*s*t5+t7/7); 358 /* p=9 4 0 1 */ 359 } else if( i==4 && j==0 && k==1 ) { 360 t1 = r*r; 361 t2 = t1*t1; 362 f = (t2-4.0/3.0*t1*r+4.0/7.0*t1-8.0/91.0*r+1.0/273.0)*(t+2.0/3.0*s+2.0/ 363 3.0*r-2.0/3.0); 364 /* p=9 4 1 0 */ 365 } else if( i==4 && j==1 && k==0 ) { 366 t1 = r*r; 367 t2 = t1*t1; 368 f = (t2-4.0/3.0*t1*r+4.0/7.0*t1-8.0/91.0*r+1.0/273.0)*(s+2.0/5.0*r-2.0/ 369 5.0); 370 /* p=9 5 0 0 */ 371 } else if( i==5 && j==0 && k==0 ) { 372 t1 = r*r; 373 t2 = t1*t1; 374 f = t2*r-2.0*t2+10.0/7.0*t1*r-40.0/91.0*t1+5.0/91.0*r-2.0/1001.0; 375 /* p=10 0 0 6 */ 376 } else if( i==0 && j==0 && k==6 ) { 377 t1 = t*t; 378 t2 = t1*t1; 379 t5 = s+r-1.0; 380 t7 = t5*t5; 381 t12 = t7*t7; 382 f = t2*t1+42.0/13.0*t2*t*t5+105.0/26.0*t2*t7+350.0/143.0*t1*t*t7*t5+ 383 105.0/143.0*t1*t12+14.0/143.0*t*t12*t5+7.0/1716.0*t12*t7; 384 /* p=10 0 1 5 */ 385 } else if( i==0 && j==1 && k==5 ) { 386 t2 = t*t; 387 t3 = t2*t2; 388 t5 = s+r-1.0; 389 t8 = t5*t5; 390 t12 = t8*t8; 391 f = (s+2.0/15.0*r-2.0/15.0)*(t3*t+30.0/11.0*t3*t5+30.0/11.0*t2*t*t8+ 392 40.0/33.0*t2*t8*t5+5.0/22.0*t*t12+t12*t5/77); 393 /* p=10 0 2 4 */ 394 } else if( i==0 && j==2 && k==4 ) { 395 t1 = s*s; 396 t2 = r-1.0; 397 t4 = t2*t2; 398 t6 = t*t; 399 t7 = t6*t6; 400 t9 = s+r-1.0; 401 t11 = t9*t9; 402 t15 = t11*t11; 403 f = (t1+2.0/5.0*s*t2+t4/35)*(t7+20.0/9.0*t6*t*t9+5.0/3.0*t6*t11+10.0/ 404 21.0*t*t11*t9+5.0/126.0*t15); 405 /* p=10 0 3 3 */ 406 } else if( i==0 && j==3 && k==3 ) { 407 t1 = s*s; 408 t3 = r-1.0; 409 t5 = t3*t3; 410 t9 = t*t; 411 t11 = s+r-1.0; 412 t13 = t11*t11; 413 f = (t1*s+4.0/5.0*t1*t3+6.0/35.0*s*t5+4.0/455.0*t3*t5)*(t9*t+12.0/7.0* 414 t11*t9+6.0/7.0*t*t13+4.0/35.0*t13*t11); 415 /* p=10 0 4 2 */ 416 } else if( i==0 && j==4 && k==2 ) { 417 t1 = s*s; 418 t2 = t1*t1; 419 t4 = r-1.0; 420 t6 = t4*t4; 421 t10 = t6*t6; 422 t12 = t*t; 423 t13 = s+r-1.0; 424 t15 = t13*t13; 425 f = (t2+4.0/3.0*t1*s*t4+4.0/7.0*t1*t6+8.0/91.0*s*t6*t4+t10/273)*(t12+ 426 6.0/5.0*t*t13+3.0/10.0*t15); 427 /* p=10 0 5 1 */ 428 } else if( i==0 && j==5 && k==1 ) { 429 t1 = s*s; 430 t2 = t1*t1; 431 t4 = r-1.0; 432 t7 = t4*t4; 433 t11 = t7*t7; 434 f = (s*t2+2.0*t2*t4+10.0/7.0*t1*s*t7+40.0/91.0*t1*t7*t4+5.0/91.0*s*t11+ 435 2.0/1001.0*t11*t4)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 436 /* p=10 0 6 0 */ 437 } else if( i==0 && j==6 && k==0 ) { 438 t1 = s*s; 439 t2 = t1*t1; 440 t5 = r-1.0; 441 t7 = t5*t5; 442 t12 = t7*t7; 443 f = t2*t1+14.0/5.0*t2*s*t5+3.0*t2*t7+20.0/13.0*t1*s*t7*t5+5.0/13.0*t1* 444 t12+6.0/143.0*s*t12*t5+t12*t7/715; 445 /* p=10 1 0 5 */ 446 } else if( i==1 && j==0 && k==5 ) { 447 t2 = t*t; 448 t3 = t2*t2; 449 t5 = s+r-1.0; 450 t8 = t5*t5; 451 t12 = t8*t8; 452 f = (r-2.0/17.0)*(t3*t+30.0/11.0*t3*t5+30.0/11.0*t2*t*t8+40.0/33.0*t2* 453 t8*t5+5.0/22.0*t*t12+t12*t5/77); 454 /* p=10 1 1 4 */ 455 } else if( i==1 && j==1 && k==4 ) { 456 t4 = t*t; 457 t5 = t4*t4; 458 t7 = s+r-1.0; 459 t9 = t7*t7; 460 t13 = t9*t9; 461 f = (r-2.0/17.0)*(s+2.0/13.0*r-2.0/13.0)*(t5+20.0/9.0*t4*t*t7+5.0/3.0* 462 t4*t9+10.0/21.0*t*t9*t7+5.0/126.0*t13); 463 /* p=10 1 2 3 */ 464 } else if( i==1 && j==2 && k==3 ) { 465 t2 = s*s; 466 t3 = r-1.0; 467 t5 = t3*t3; 468 t8 = t*t; 469 t10 = s+r-1.0; 470 t12 = t10*t10; 471 f = (r-2.0/17.0)*(t2+6.0/13.0*s*t3+t5/26)*(t8*t+12.0/7.0*t8*t10+6.0/7.0 472 *t*t12+4.0/35.0*t12*t10); 473 /* p=10 1 3 2 */ 474 } else if( i==1 && j==3 && k==2 ) { 475 t2 = s*s; 476 t4 = r-1.0; 477 t6 = t4*t4; 478 t11 = t*t; 479 t12 = s+r-1.0; 480 t14 = t12*t12; 481 f = (r-2.0/17.0)*(t2*s+12.0/13.0*t2*t4+3.0/13.0*s*t6+2.0/143.0*t6*t4)*( 482 t11+6.0/5.0*t*t12+3.0/10.0*t14); 483 /* p=10 1 4 1 */ 484 } else if( i==1 && j==4 && k==1 ) { 485 t2 = s*s; 486 t3 = t2*t2; 487 t5 = r-1.0; 488 t7 = t5*t5; 489 t11 = t7*t7; 490 f = (r-2.0/17.0)*(t3+20.0/13.0*t2*s*t5+10.0/13.0*t2*t7+20.0/143.0*s*t7* 491 t5+t11/143)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 492 /* p=10 1 5 0 */ 493 } else if( i==1 && j==5 && k==0 ) { 494 t2 = s*s; 495 t3 = t2*t2; 496 t5 = r-1.0; 497 t8 = t5*t5; 498 t12 = t8*t8; 499 f = (r-2.0/17.0)*(s*t3+30.0/13.0*t3*t5+25.0/13.0*t2*s*t8+100.0/143.0*t2 500 *t8*t5+15.0/143.0*s*t12+2.0/429.0*t12*t5); 501 /* p=10 2 0 4 */ 502 } else if( i==2 && j==0 && k==4 ) { 503 t1 = r*r; 504 t3 = t*t; 505 t4 = t3*t3; 506 t6 = s+r-1.0; 507 t8 = t6*t6; 508 t12 = t8*t8; 509 f = (t1-6.0/17.0*r+3.0/136.0)*(t4+20.0/9.0*t3*t*t6+5.0/3.0*t3*t8+10.0/ 510 21.0*t*t8*t6+5.0/126.0*t12); 511 /* p=10 2 1 3 */ 512 } else if( i==2 && j==1 && k==3 ) { 513 t1 = r*r; 514 t5 = t*t; 515 t7 = s+r-1.0; 516 t9 = t7*t7; 517 f = (t1-6.0/17.0*r+3.0/136.0)*(s+2.0/11.0*r-2.0/11.0)*(t5*t+12.0/7.0*t7 518 *t5+6.0/7.0*t*t9+4.0/35.0*t9*t7); 519 /* p=10 2 2 2 */ 520 } else if( i==2 && j==2 && k==2 ) { 521 t1 = r*r; 522 t3 = s*s; 523 t4 = r-1.0; 524 t6 = t4*t4; 525 t9 = t*t; 526 t10 = s+r-1.0; 527 t12 = t10*t10; 528 f = (t1-6.0/17.0*r+3.0/136.0)*(t3+6.0/11.0*s*t4+3.0/55.0*t6)*(t9+6.0/ 529 5.0*t*t10+3.0/10.0*t12); 530 /* p=10 2 3 1 */ 531 } else if( i==2 && j==3 && k==1 ) { 532 t1 = r*r; 533 t3 = s*s; 534 t5 = r-1.0; 535 t7 = t5*t5; 536 f = (t1-6.0/17.0*r+3.0/136.0)*(s*t3+12.0/11.0*t3*t5+18.0/55.0*s*t7+4.0/ 537 165.0*t7*t5)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 538 /* p=10 2 4 0 */ 539 } else if( i==2 && j==4 && k==0 ) { 540 t1 = r*r; 541 t3 = s*s; 542 t4 = t3*t3; 543 t6 = r-1.0; 544 t8 = t6*t6; 545 t12 = t8*t8; 546 f = (t1-6.0/17.0*r+3.0/136.0)*(t4+20.0/11.0*t3*s*t6+12.0/11.0*t3*t8+8.0 547 /33.0*s*t8*t6+t12/66); 548 /* p=10 3 0 3 */ 549 } else if( i==3 && j==0 && k==3 ) { 550 t1 = r*r; 551 t4 = t*t; 552 t6 = s+r-1.0; 553 t8 = t6*t6; 554 f = (t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0)*(t4*t+12.0/7.0*t4*t6+6.0/ 555 7.0*t*t8+4.0/35.0*t8*t6); 556 /* p=10 3 1 2 */ 557 } else if( i==3 && j==1 && k==2 ) { 558 t1 = r*r; 559 t6 = t*t; 560 t7 = s+r-1.0; 561 t9 = t7*t7; 562 f = (t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0)*(s+2.0/9.0*r-2.0/9.0)*(t6+ 563 6.0/5.0*t*t7+3.0/10.0*t9); 564 /* p=10 3 2 1 */ 565 } else if( i==3 && j==2 && k==1 ) { 566 t1 = r*r; 567 t4 = s*s; 568 t5 = r-1.0; 569 t7 = t5*t5; 570 f = (t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0)*(t4+2.0/3.0*s*t5+t7/12)*(t 571 +2.0/3.0*s+2.0/3.0*r-2.0/3.0); 572 /* p=10 3 3 0 */ 573 } else if( i==3 && j==3 && k==0 ) { 574 t1 = r*r; 575 t4 = s*s; 576 t6 = r-1.0; 577 t8 = t6*t6; 578 f = (t1*r-12.0/17.0*t1+9.0/68.0*r-1.0/170.0)*(t4*s+4.0/3.0*t4*t6+s*t8/2 579 +t8*t6/21); 580 /* p=10 4 0 2 */ 581 } else if( i==4 && j==0 && k==2 ) { 582 t1 = r*r; 583 t2 = t1*t1; 584 t5 = t*t; 585 t6 = s+r-1.0; 586 t8 = t6*t6; 587 f = (t2-20.0/17.0*t1*r+15.0/34.0*t1-r/17+1.0/476.0)*(t5+6.0/5.0*t*t6+ 588 3.0/10.0*t8); 589 /* p=10 4 1 1 */ 590 } else if( i==4 && j==1 && k==1 ) { 591 t1 = r*r; 592 t2 = t1*t1; 593 f = (t2-20.0/17.0*t1*r+15.0/34.0*t1-r/17+1.0/476.0)*(s+2.0/7.0*r-2.0/7.0 594 )*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 595 /* p=10 4 2 0 */ 596 } else if( i==4 && j==2 && k==0 ) { 597 t1 = r*r; 598 t2 = t1*t1; 599 t5 = s*s; 600 t6 = r-1.0; 601 t8 = t6*t6; 602 f = (t2-20.0/17.0*t1*r+15.0/34.0*t1-r/17+1.0/476.0)*(t5+6.0/7.0*s*t6+t8 603 /7); 604 /* p=10 5 0 1 */ 605 } else if( i==5 && j==0 && k==1 ) { 606 t1 = r*r; 607 t2 = t1*t1; 608 f = (t2*r-30.0/17.0*t2+75.0/68.0*t1*r-5.0/17.0*t1+15.0/476.0*r-3.0/ 609 3094.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 610 /* p=10 5 1 0 */ 611 } else if( i==5 && j==1 && k==0 ) { 612 t1 = r*r; 613 t2 = t1*t1; 614 f = (t2*r-30.0/17.0*t2+75.0/68.0*t1*r-5.0/17.0*t1+15.0/476.0*r-3.0/ 615 3094.0)*(s+2.0/5.0*r-2.0/5.0); 616 /* p=10 6 0 0 */ 617 } else if( i==6 && j==0 && k==0 ) { 618 t1 = r*r; 619 t2 = t1*t1; 620 f = t2*t1-42.0/17.0*t2*r+315.0/136.0*t2-35.0/34.0*t1*r+15.0/68.0*t1-9.0/ 621 442.0*r+1.0/1768.0; 622 /* p=11 0 0 7 */ 623 } else if( i==0 && j==0 && k==7 ) { 624 t1 = t*t; 625 t2 = t1*t; 626 t3 = t1*t1; 627 t6 = s+r-1.0; 628 t9 = t6*t6; 629 t11 = t9*t6; 630 t13 = t9*t9; 631 f = t3*t2+56.0/15.0*t3*t1*t6+28.0/5.0*t3*t*t9+56.0/13.0*t3*t11+70.0/ 632 39.0*t2*t13+56.0/143.0*t1*t13*t6+28.0/715.0*t*t13*t9+8.0/6435.0*t13*t11; 633 /* p=11 0 1 6 */ 634 } else if( i==0 && j==1 && k==6 ) { 635 t2 = t*t; 636 t3 = t2*t2; 637 t6 = s+r-1.0; 638 t8 = t6*t6; 639 t13 = t8*t8; 640 f = (s+2.0/17.0*r-2.0/17.0)*(t3*t2+42.0/13.0*t3*t*t6+105.0/26.0*t3*t8+ 641 350.0/143.0*t2*t*t8*t6+105.0/143.0*t2*t13+14.0/143.0*t*t13*t6+7.0/1716.0*t13*t8 642 ); 643 /* p=11 0 2 5 */ 644 } else if( i==0 && j==2 && k==5 ) { 645 t1 = s*s; 646 t2 = r-1.0; 647 t4 = t2*t2; 648 t6 = t*t; 649 t7 = t6*t6; 650 t9 = s+r-1.0; 651 t12 = t9*t9; 652 t16 = t12*t12; 653 f = (t1+6.0/17.0*s*t2+3.0/136.0*t4)*(t7*t+30.0/11.0*t7*t9+30.0/11.0*t6* 654 t*t12+40.0/33.0*t6*t12*t9+5.0/22.0*t*t16+t16*t9/77); 655 /* p=11 0 3 4 */ 656 } else if( i==0 && j==3 && k==4 ) { 657 t1 = s*s; 658 t3 = r-1.0; 659 t5 = t3*t3; 660 t9 = t*t; 661 t10 = t9*t9; 662 t12 = s+r-1.0; 663 t14 = t12*t12; 664 t18 = t14*t14; 665 f = (t1*s+12.0/17.0*t3*t1+9.0/68.0*s*t5+t5*t3/170)*(t10+20.0/9.0*t9*t* 666 t12+5.0/3.0*t9*t14+10.0/21.0*t*t14*t12+5.0/126.0*t18); 667 /* p=11 0 4 3 */ 668 } else if( i==0 && j==4 && k==3 ) { 669 t1 = s*s; 670 t2 = t1*t1; 671 t4 = r-1.0; 672 t6 = t4*t4; 673 t10 = t6*t6; 674 t12 = t*t; 675 t14 = s+r-1.0; 676 t16 = t14*t14; 677 f = (t2+20.0/17.0*t1*s*t4+15.0/34.0*t1*t6+s*t6*t4/17+t10/476)*(t12*t+ 678 12.0/7.0*t14*t12+6.0/7.0*t*t16+4.0/35.0*t16*t14); 679 /* p=11 0 5 2 */ 680 } else if( i==0 && j==5 && k==2 ) { 681 t1 = s*s; 682 t2 = t1*t1; 683 t4 = r-1.0; 684 t7 = t4*t4; 685 t11 = t7*t7; 686 t15 = t*t; 687 t16 = s+r-1.0; 688 t18 = t16*t16; 689 f = (s*t2+30.0/17.0*t2*t4+75.0/68.0*t1*s*t7+5.0/17.0*t1*t7*t4+15.0/ 690 476.0*s*t11+3.0/3094.0*t11*t4)*(t15+6.0/5.0*t*t16+3.0/10.0*t18); 691 /* p=11 0 6 1 */ 692 } else if( i==0 && j==6 && k==1 ) { 693 t1 = s*s; 694 t2 = t1*t1; 695 t5 = r-1.0; 696 t7 = t5*t5; 697 t12 = t7*t7; 698 f = (t2*t1+42.0/17.0*t2*s*t5+315.0/136.0*t2*t7+35.0/34.0*t1*s*t7*t5+ 699 15.0/68.0*t1*t12+9.0/442.0*s*t12*t5+t12*t7/1768)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0 700 ); 701 /* p=11 0 7 0 */ 702 } else if( i==0 && j==7 && k==0 ) { 703 t1 = s*s; 704 t2 = t1*s; 705 t3 = t1*t1; 706 t6 = r-1.0; 707 t9 = t6*t6; 708 t11 = t9*t6; 709 t13 = t9*t9; 710 f = t3*t2+56.0/17.0*t3*t1*t6+147.0/34.0*t3*s*t9+49.0/17.0*t3*t11+35.0/ 711 34.0*t2*t13+42.0/221.0*t1*t13*t6+7.0/442.0*s*t13*t9+t13*t11/2431; 712 /* p=11 1 0 6 */ 713 } else if( i==1 && j==0 && k==6 ) { 714 t2 = t*t; 715 t3 = t2*t2; 716 t6 = s+r-1.0; 717 t8 = t6*t6; 718 t13 = t8*t8; 719 f = (r-2.0/19.0)*(t3*t2+42.0/13.0*t3*t*t6+105.0/26.0*t3*t8+350.0/143.0* 720 t2*t*t8*t6+105.0/143.0*t2*t13+14.0/143.0*t*t13*t6+7.0/1716.0*t13*t8); 721 /* p=11 1 1 5 */ 722 } else if( i==1 && j==1 && k==5 ) { 723 t4 = t*t; 724 t5 = t4*t4; 725 t7 = s+r-1.0; 726 t10 = t7*t7; 727 t14 = t10*t10; 728 f = (r-2.0/19.0)*(s+2.0/15.0*r-2.0/15.0)*(t5*t+30.0/11.0*t5*t7+30.0/ 729 11.0*t4*t*t10+40.0/33.0*t4*t10*t7+5.0/22.0*t*t14+t14*t7/77); 730 /* p=11 1 2 4 */ 731 } else if( i==1 && j==2 && k==4 ) { 732 t2 = s*s; 733 t3 = r-1.0; 734 t5 = t3*t3; 735 t8 = t*t; 736 t9 = t8*t8; 737 t11 = s+r-1.0; 738 t13 = t11*t11; 739 t17 = t13*t13; 740 f = (r-2.0/19.0)*(t2+2.0/5.0*s*t3+t5/35)*(t9+20.0/9.0*t8*t*t11+5.0/3.0* 741 t13*t8+10.0/21.0*t*t13*t11+5.0/126.0*t17); 742 /* p=11 1 3 3 */ 743 } else if( i==1 && j==3 && k==3 ) { 744 t2 = s*s; 745 t4 = r-1.0; 746 t6 = t4*t4; 747 t11 = t*t; 748 t13 = s+r-1.0; 749 t15 = t13*t13; 750 f = (r-2.0/19.0)*(t2*s+4.0/5.0*t2*t4+6.0/35.0*s*t6+4.0/455.0*t6*t4)*( 751 t11*t+12.0/7.0*t13*t11+6.0/7.0*t*t15+4.0/35.0*t15*t13); 752 /* p=11 1 4 2 */ 753 } else if( i==1 && j==4 && k==2 ) { 754 t2 = s*s; 755 t3 = t2*t2; 756 t5 = r-1.0; 757 t7 = t5*t5; 758 t11 = t7*t7; 759 t14 = t*t; 760 t15 = s+r-1.0; 761 t17 = t15*t15; 762 f = (r-2.0/19.0)*(t3+4.0/3.0*t2*s*t5+4.0/7.0*t2*t7+8.0/91.0*s*t5*t7+t11 763 /273)*(t14+6.0/5.0*t*t15+3.0/10.0*t17); 764 /* p=11 1 5 1 */ 765 } else if( i==1 && j==5 && k==1 ) { 766 t2 = s*s; 767 t3 = t2*t2; 768 t5 = r-1.0; 769 t8 = t5*t5; 770 t12 = t8*t8; 771 f = (r-2.0/19.0)*(s*t3+2.0*t3*t5+10.0/7.0*t2*s*t8+40.0/91.0*t2*t8*t5+ 772 5.0/91.0*s*t12+2.0/1001.0*t12*t5)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 773 /* p=11 1 6 0 */ 774 } else if( i==1 && j==6 && k==0 ) { 775 t2 = s*s; 776 t3 = t2*t2; 777 t6 = r-1.0; 778 t8 = t6*t6; 779 t13 = t8*t8; 780 f = (r-2.0/19.0)*(t3*t2+14.0/5.0*t3*s*t6+3.0*t3*t8+20.0/13.0*t2*s*t8*t6 781 +5.0/13.0*t2*t13+6.0/143.0*s*t13*t6+t13*t8/715); 782 /* p=11 2 0 5 */ 783 } else if( i==2 && j==0 && k==5 ) { 784 t1 = r*r; 785 t3 = t*t; 786 t4 = t3*t3; 787 t6 = s+r-1.0; 788 t9 = t6*t6; 789 t13 = t9*t9; 790 f = (t1-6.0/19.0*r+1.0/57.0)*(t4*t+30.0/11.0*t6*t4+30.0/11.0*t3*t*t9+ 791 40.0/33.0*t3*t9*t6+5.0/22.0*t*t13+t13*t6/77); 792 /* p=11 2 1 4 */ 793 } else if( i==2 && j==1 && k==4 ) { 794 t1 = r*r; 795 t5 = t*t; 796 t6 = t5*t5; 797 t8 = s+r-1.0; 798 t10 = t8*t8; 799 t14 = t10*t10; 800 f = (s+2.0/13.0*r-2.0/13.0)*(t1-6.0/19.0*r+1.0/57.0)*(t6+20.0/9.0*t5*t* 801 t8+5.0/3.0*t5*t10+10.0/21.0*t*t10*t8+5.0/126.0*t14); 802 /* p=11 2 2 3 */ 803 } else if( i==2 && j==2 && k==3 ) { 804 t1 = r*r; 805 t3 = s*s; 806 t4 = r-1.0; 807 t6 = t4*t4; 808 t9 = t*t; 809 t11 = s+r-1.0; 810 t13 = t11*t11; 811 f = (t1-6.0/19.0*r+1.0/57.0)*(t3+6.0/13.0*s*t4+t6/26)*(t9*t+12.0/7.0*t9 812 *t11+6.0/7.0*t*t13+4.0/35.0*t13*t11); 813 /* p=11 2 3 2 */ 814 } else if( i==2 && j==3 && k==2 ) { 815 t1 = r*r; 816 t3 = s*s; 817 t5 = r-1.0; 818 t7 = t5*t5; 819 t12 = t*t; 820 t13 = s+r-1.0; 821 t15 = t13*t13; 822 f = (t1-6.0/19.0*r+1.0/57.0)*(t3*s+12.0/13.0*t3*t5+3.0/13.0*s*t7+2.0/ 823 143.0*t7*t5)*(t12+6.0/5.0*t*t13+3.0/10.0*t15); 824 /* p=11 2 4 1 */ 825 } else if( i==2 && j==4 && k==1 ) { 826 t1 = r*r; 827 t3 = s*s; 828 t4 = t3*t3; 829 t6 = r-1.0; 830 t8 = t6*t6; 831 t12 = t8*t8; 832 f = (t1-6.0/19.0*r+1.0/57.0)*(t4+20.0/13.0*t3*s*t6+10.0/13.0*t3*t8+20.0 833 /143.0*s*t8*t6+t12/143)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 834 /* p=11 2 5 0 */ 835 } else if( i==2 && j==5 && k==0 ) { 836 t1 = r*r; 837 t3 = s*s; 838 t4 = t3*t3; 839 t6 = r-1.0; 840 t9 = t6*t6; 841 t13 = t9*t9; 842 f = (t1-6.0/19.0*r+1.0/57.0)*(s*t4+30.0/13.0*t4*t6+25.0/13.0*t3*s*t9+ 843 100.0/143.0*t3*t9*t6+15.0/143.0*s*t13+2.0/429.0*t13*t6); 844 /* p=11 3 0 4 */ 845 } else if( i==3 && j==0 && k==4 ) { 846 t1 = r*r; 847 t4 = t*t; 848 t5 = t4*t4; 849 t7 = s+r-1.0; 850 t9 = t7*t7; 851 t13 = t9*t9; 852 f = (t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0)*(t5+20.0/9.0*t4*t*t7+5.0/ 853 3.0*t4*t9+10.0/21.0*t*t9*t7+5.0/126.0*t13); 854 /* p=11 3 1 3 */ 855 } else if( i==3 && j==1 && k==3 ) { 856 t1 = r*r; 857 t6 = t*t; 858 t8 = s+r-1.0; 859 t10 = t8*t8; 860 f = (t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0)*(s+2.0/11.0*r-2.0/11.0)*( 861 t6*t+12.0/7.0*t8*t6+6.0/7.0*t*t10+4.0/35.0*t10*t8); 862 /* p=11 3 2 2 */ 863 } else if( i==3 && j==2 && k==2 ) { 864 t1 = r*r; 865 t4 = s*s; 866 t5 = r-1.0; 867 t7 = t5*t5; 868 t10 = t*t; 869 t11 = s+r-1.0; 870 t13 = t11*t11; 871 f = (t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0)*(t4+6.0/11.0*s*t5+3.0/55.0 872 *t7)*(t10+6.0/5.0*t*t11+3.0/10.0*t13); 873 /* p=11 3 3 1 */ 874 } else if( i==3 && j==3 && k==1 ) { 875 t1 = r*r; 876 t4 = s*s; 877 t6 = r-1.0; 878 t8 = t6*t6; 879 f = (t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0)*(s*t4+12.0/11.0*t4*t6+18.0 880 /55.0*s*t8+4.0/165.0*t8*t6)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 881 /* p=11 3 4 0 */ 882 } else if( i==3 && j==4 && k==0 ) { 883 t1 = r*r; 884 t4 = s*s; 885 t5 = t4*t4; 886 t7 = r-1.0; 887 t9 = t7*t7; 888 t13 = t9*t9; 889 f = (t1*r-12.0/19.0*t1+2.0/19.0*r-4.0/969.0)*(t5+20.0/11.0*t4*s*t7+12.0 890 /11.0*t4*t9+8.0/33.0*s*t9*t7+t13/66); 891 /* p=11 4 0 3 */ 892 } else if( i==4 && j==0 && k==3 ) { 893 t1 = r*r; 894 t2 = t1*t1; 895 t5 = t*t; 896 t7 = s+r-1.0; 897 t9 = t7*t7; 898 f = (t2-20.0/19.0*t1*r+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0)*(t5*t+12.0 899 /7.0*t7*t5+6.0/7.0*t9*t+4.0/35.0*t9*t7); 900 /* p=11 4 1 2 */ 901 } else if( i==4 && j==1 && k==2 ) { 902 t1 = r*r; 903 t2 = t1*t1; 904 t7 = t*t; 905 t8 = s+r-1.0; 906 t10 = t8*t8; 907 f = (t2-20.0/19.0*t1*r+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0)*(s+2.0/9.0 908 *r-2.0/9.0)*(t7+6.0/5.0*t*t8+3.0/10.0*t10); 909 /* p=11 4 2 1 */ 910 } else if( i==4 && j==2 && k==1 ) { 911 t1 = r*r; 912 t2 = t1*t1; 913 t5 = s*s; 914 t6 = r-1.0; 915 t8 = t6*t6; 916 f = (t2-20.0/19.0*t1*r+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0)*(t5+2.0/ 917 3.0*s*t6+t8/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 918 /* p=11 4 3 0 */ 919 } else if( i==4 && j==3 && k==0 ) { 920 t1 = r*r; 921 t2 = t1*t1; 922 t5 = s*s; 923 t7 = r-1.0; 924 t9 = t7*t7; 925 f = (t2-20.0/19.0*t1*r+20.0/57.0*t1-40.0/969.0*r+5.0/3876.0)*(t5*s+4.0/ 926 3.0*t7*t5+s*t9/2+t9*t7/21); 927 /* p=11 5 0 2 */ 928 } else if( i==5 && j==0 && k==2 ) { 929 t1 = r*r; 930 t2 = t1*t1; 931 t6 = t*t; 932 t7 = s+r-1.0; 933 t9 = t7*t7; 934 f = (t2*r-30.0/19.0*t2+50.0/57.0*t1*r-200.0/969.0*t1+25.0/1292.0*r-1.0/ 935 1938.0)*(t6+6.0/5.0*t*t7+3.0/10.0*t9); 936 /* p=11 5 1 1 */ 937 } else if( i==5 && j==1 && k==1 ) { 938 t1 = r*r; 939 t2 = t1*t1; 940 f = (t2*r-30.0/19.0*t2+50.0/57.0*t1*r-200.0/969.0*t1+25.0/1292.0*r-1.0/ 941 1938.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 942 /* p=11 5 2 0 */ 943 } else if( i==5 && j==2 && k==0 ) { 944 t1 = r*r; 945 t2 = t1*t1; 946 t6 = s*s; 947 t7 = r-1.0; 948 t9 = t7*t7; 949 f = (t2*r-30.0/19.0*t2+50.0/57.0*t1*r-200.0/969.0*t1+25.0/1292.0*r-1.0/ 950 1938.0)*(t6+6.0/7.0*s*t7+t9/7); 951 /* p=11 6 0 1 */ 952 } else if( i==6 && j==0 && k==1 ) { 953 t1 = r*r; 954 t2 = t1*t1; 955 f = (t2*t1-42.0/19.0*t2*r+35.0/19.0*t2-700.0/969.0*t1*r+175.0/1292.0*t1 956 -7.0/646.0*r+1.0/3876.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 957 /* p=11 6 1 0 */ 958 } else if( i==6 && j==1 && k==0 ) { 959 t1 = r*r; 960 t2 = t1*t1; 961 f = (t2*t1-42.0/19.0*t2*r+35.0/19.0*t2-700.0/969.0*t1*r+175.0/1292.0*t1 962 -7.0/646.0*r+1.0/3876.0)*(s+2.0/5.0*r-2.0/5.0); 963 /* p=11 7 0 0 */ 964 } else if( i==7 && j==0 && k==0 ) { 965 t1 = r*r; 966 t2 = t1*r; 967 t3 = t1*t1; 968 f = t3*t2-56.0/19.0*t3*t1+196.0/57.0*t3*r-1960.0/969.0*t3+1225.0/1938.0* 969 t2-98.0/969.0*t1+7.0/969.0*r-2.0/12597.0; 970 /* p=12 0 0 8 */ 971 } else if( i==0 && j==0 && k==8 ) { 972 t1 = t*t; 973 t2 = t1*t1; 974 t3 = t2*t2; 975 t4 = t1*t; 976 t6 = s+r-1.0; 977 t9 = t6*t6; 978 t12 = t9*t6; 979 t14 = t9*t9; 980 t22 = t14*t14; 981 f = t3+72.0/17.0*t2*t4*t6+126.0/17.0*t2*t1*t9+588.0/85.0*t2*t*t12+63.0/ 982 17.0*t2*t14+252.0/221.0*t4*t14*t6+42.0/221.0*t1*t14*t9+36.0/2431.0*t*t14*t12+ 983 9.0/24310.0*t22; 984 /* p=12 0 1 7 */ 985 } else if( i==0 && j==1 && k==7 ) { 986 t2 = t*t; 987 t3 = t2*t; 988 t4 = t2*t2; 989 t7 = s+r-1.0; 990 t10 = t7*t7; 991 t12 = t10*t7; 992 t14 = t10*t10; 993 f = (s+2.0/19.0*r-2.0/19.0)*(t4*t3+56.0/15.0*t4*t2*t7+28.0/5.0*t4*t*t10 994 +56.0/13.0*t4*t12+70.0/39.0*t3*t14+56.0/143.0*t2*t14*t7+28.0/715.0*t*t14*t10+ 995 8.0/6435.0*t14*t12); 996 /* p=12 0 2 6 */ 997 } else if( i==0 && j==2 && k==6 ) { 998 t1 = s*s; 999 t2 = r-1.0; 1000 t4 = t2*t2; 1001 t6 = t*t; 1002 t7 = t6*t6; 1003 t10 = s+r-1.0; 1004 t12 = t10*t10; 1005 t17 = t12*t12; 1006 f = (t1+6.0/19.0*s*t2+t4/57)*(t6*t7+42.0/13.0*t7*t*t10+105.0/26.0*t7* 1007 t12+350.0/143.0*t6*t*t12*t10+105.0/143.0*t6*t17+14.0/143.0*t*t17*t10+7.0/1716.0 1008 *t17*t12); 1009 /* p=12 0 3 5 */ 1010 } else if( i==0 && j==3 && k==5 ) { 1011 t1 = s*s; 1012 t3 = r-1.0; 1013 t5 = t3*t3; 1014 t9 = t*t; 1015 t10 = t9*t9; 1016 t12 = s+r-1.0; 1017 t15 = t12*t12; 1018 t19 = t15*t15; 1019 f = (t1*s+12.0/19.0*t1*t3+2.0/19.0*s*t5+4.0/969.0*t5*t3)*(t10*t+30.0/ 1020 11.0*t12*t10+30.0/11.0*t9*t*t15+40.0/33.0*t9*t15*t12+5.0/22.0*t*t19+t19*t12/77) 1021 ; 1022 /* p=12 0 4 4 */ 1023 } else if( i==0 && j==4 && k==4 ) { 1024 t1 = s*s; 1025 t2 = t1*t1; 1026 t4 = r-1.0; 1027 t6 = t4*t4; 1028 t10 = t6*t6; 1029 t12 = t*t; 1030 t13 = t12*t12; 1031 t15 = s+r-1.0; 1032 t17 = t15*t15; 1033 t21 = t17*t17; 1034 f = (t2+20.0/19.0*t4*t1*s+20.0/57.0*t1*t6+40.0/969.0*s*t6*t4+5.0/3876.0 1035 *t10)*(t13+20.0/9.0*t12*t*t15+5.0/3.0*t17*t12+10.0/21.0*t*t17*t15+5.0/126.0*t21 1036 ); 1037 /* p=12 0 5 3 */ 1038 } else if( i==0 && j==5 && k==3 ) { 1039 t1 = s*s; 1040 t2 = t1*t1; 1041 t4 = r-1.0; 1042 t7 = t4*t4; 1043 t11 = t7*t7; 1044 t15 = t*t; 1045 t17 = s+r-1.0; 1046 t19 = t17*t17; 1047 f = (s*t2+30.0/19.0*t4*t2+50.0/57.0*t1*s*t7+200.0/969.0*t1*t7*t4+25.0/ 1048 1292.0*s*t11+t11*t4/1938)*(t15*t+12.0/7.0*t17*t15+6.0/7.0*t*t19+4.0/35.0*t19* 1049 t17); 1050 /* p=12 0 6 2 */ 1051 } else if( i==0 && j==6 && k==2 ) { 1052 t1 = s*s; 1053 t2 = t1*t1; 1054 t5 = r-1.0; 1055 t7 = t5*t5; 1056 t12 = t7*t7; 1057 t18 = t*t; 1058 t19 = s+r-1.0; 1059 t21 = t19*t19; 1060 f = (t2*t1+42.0/19.0*t2*s*t5+35.0/19.0*t2*t7+700.0/969.0*t1*s*t7*t5+ 1061 175.0/1292.0*t1*t12+7.0/646.0*s*t12*t5+t7*t12/3876)*(t18+6.0/5.0*t*t19+3.0/10.0 1062 *t21); 1063 /* p=12 0 7 1 */ 1064 } else if( i==0 && j==7 && k==1 ) { 1065 t1 = s*s; 1066 t2 = t1*s; 1067 t3 = t1*t1; 1068 t6 = r-1.0; 1069 t9 = t6*t6; 1070 t11 = t9*t6; 1071 t13 = t9*t9; 1072 f = (t3*t2+56.0/19.0*t1*t3*t6+196.0/57.0*t3*s*t9+1960.0/969.0*t3*t11+ 1073 1225.0/1938.0*t2*t13+98.0/969.0*t1*t13*t6+7.0/969.0*s*t13*t9+2.0/12597.0*t13* 1074 t11)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 1075 /* p=12 0 8 0 */ 1076 } else if( i==0 && j==8 && k==0 ) { 1077 t1 = s*s; 1078 t2 = t1*t1; 1079 t3 = t2*t2; 1080 t4 = t1*s; 1081 t6 = r-1.0; 1082 t9 = t6*t6; 1083 t12 = t9*t6; 1084 t14 = t9*t9; 1085 t22 = t14*t14; 1086 f = t3+72.0/19.0*t2*t4*t6+112.0/19.0*t2*t1*t9+1568.0/323.0*t2*s*t12+ 1087 735.0/323.0*t2*t14+196.0/323.0*t4*t14*t6+28.0/323.0*t1*t14*t9+24.0/4199.0*s*t14 1088 *t12+t22/8398; 1089 /* p=12 1 0 7 */ 1090 } else if( i==1 && j==0 && k==7 ) { 1091 t2 = t*t; 1092 t3 = t2*t; 1093 t4 = t2*t2; 1094 t7 = s+r-1.0; 1095 t10 = t7*t7; 1096 t12 = t10*t7; 1097 t14 = t10*t10; 1098 f = (r-2.0/21.0)*(t4*t3+56.0/15.0*t2*t4*t7+28.0/5.0*t4*t*t10+56.0/13.0* 1099 t4*t12+70.0/39.0*t3*t14+56.0/143.0*t2*t14*t7+28.0/715.0*t*t14*t10+8.0/6435.0* 1100 t14*t12); 1101 /* p=12 1 1 6 */ 1102 } else if( i==1 && j==1 && k==6 ) { 1103 t4 = t*t; 1104 t5 = t4*t4; 1105 t8 = s+r-1.0; 1106 t10 = t8*t8; 1107 t15 = t10*t10; 1108 f = (s+2.0/17.0*r-2.0/17.0)*(r-2.0/21.0)*(t5*t4+42.0/13.0*t5*t*t8+105.0 1109 /26.0*t5*t10+350.0/143.0*t4*t*t10*t8+105.0/143.0*t4*t15+14.0/143.0*t*t15*t8+7.0 1110 /1716.0*t15*t10); 1111 /* p=12 1 2 5 */ 1112 } else if( i==1 && j==2 && k==5 ) { 1113 t2 = s*s; 1114 t3 = r-1.0; 1115 t5 = t3*t3; 1116 t8 = t*t; 1117 t9 = t8*t8; 1118 t11 = s+r-1.0; 1119 t14 = t11*t11; 1120 t18 = t14*t14; 1121 f = (r-2.0/21.0)*(t2+6.0/17.0*s*t3+3.0/136.0*t5)*(t9*t+30.0/11.0*t9*t11 1122 +30.0/11.0*t8*t*t14+40.0/33.0*t8*t14*t11+5.0/22.0*t*t18+t18*t11/77); 1123 /* p=12 1 3 4 */ 1124 } else if( i==1 && j==3 && k==4 ) { 1125 t2 = s*s; 1126 t4 = r-1.0; 1127 t6 = t4*t4; 1128 t11 = t*t; 1129 t12 = t11*t11; 1130 t14 = s+r-1.0; 1131 t16 = t14*t14; 1132 t20 = t16*t16; 1133 f = (r-2.0/21.0)*(t2*s+12.0/17.0*t2*t4+9.0/68.0*s*t6+t6*t4/170)*(t12+ 1134 20.0/9.0*t11*t*t14+5.0/3.0*t11*t16+10.0/21.0*t*t16*t14+5.0/126.0*t20); 1135 /* p=12 1 4 3 */ 1136 } else if( i==1 && j==4 && k==3 ) { 1137 t2 = s*s; 1138 t3 = t2*t2; 1139 t5 = r-1.0; 1140 t7 = t5*t5; 1141 t11 = t7*t7; 1142 t14 = t*t; 1143 t16 = s+r-1.0; 1144 t18 = t16*t16; 1145 f = (r-2.0/21.0)*(t3+20.0/17.0*t5*t2*s+15.0/34.0*t2*t7+s*t7*t5/17+t11 1146 /476)*(t14*t+12.0/7.0*t16*t14+6.0/7.0*t*t18+4.0/35.0*t18*t16); 1147 /* p=12 1 5 2 */ 1148 } else if( i==1 && j==5 && k==2 ) { 1149 t2 = s*s; 1150 t3 = t2*t2; 1151 t5 = r-1.0; 1152 t8 = t5*t5; 1153 t12 = t8*t8; 1154 t17 = t*t; 1155 t18 = s+r-1.0; 1156 t20 = t18*t18; 1157 f = (r-2.0/21.0)*(s*t3+30.0/17.0*t3*t5+75.0/68.0*t2*s*t8+5.0/17.0*t2*t8 1158 *t5+15.0/476.0*s*t12+3.0/3094.0*t12*t5)*(t17+6.0/5.0*t*t18+3.0/10.0*t20); 1159 /* p=12 1 6 1 */ 1160 } else if( i==1 && j==6 && k==1 ) { 1161 t2 = s*s; 1162 t3 = t2*t2; 1163 t6 = r-1.0; 1164 t8 = t6*t6; 1165 t13 = t8*t8; 1166 f = (r-2.0/21.0)*(t3*t2+42.0/17.0*t3*s*t6+315.0/136.0*t3*t8+35.0/34.0* 1167 t2*s*t8*t6+15.0/68.0*t2*t13+9.0/442.0*s*t13*t6+t13*t8/1768)*(t+2.0/3.0*s+2.0/ 1168 3.0*r-2.0/3.0); 1169 /* p=12 1 7 0 */ 1170 } else if( i==1 && j==7 && k==0 ) { 1171 t2 = s*s; 1172 t3 = t2*s; 1173 t4 = t2*t2; 1174 t7 = r-1.0; 1175 t10 = t7*t7; 1176 t12 = t10*t7; 1177 t14 = t10*t10; 1178 f = (r-2.0/21.0)*(t4*t3+56.0/17.0*t4*t2*t7+147.0/34.0*t4*s*t10+49.0/ 1179 17.0*t4*t12+35.0/34.0*t3*t14+42.0/221.0*t2*t14*t7+7.0/442.0*s*t14*t10+t14*t12 1180 /2431); 1181 /* p=12 2 0 6 */ 1182 } else if( i==2 && j==0 && k==6 ) { 1183 t1 = r*r; 1184 t3 = t*t; 1185 t4 = t3*t3; 1186 t7 = s+r-1.0; 1187 t9 = t7*t7; 1188 t14 = t9*t9; 1189 f = (t1-2.0/7.0*r+1.0/70.0)*(t4*t3+42.0/13.0*t4*t*t7+105.0/26.0*t4*t9+ 1190 350.0/143.0*t3*t*t9*t7+105.0/143.0*t3*t14+14.0/143.0*t*t14*t7+7.0/1716.0*t14*t9 1191 ); 1192 /* p=12 2 1 5 */ 1193 } else if( i==2 && j==1 && k==5 ) { 1194 t1 = r*r; 1195 t5 = t*t; 1196 t6 = t5*t5; 1197 t8 = s+r-1.0; 1198 t11 = t8*t8; 1199 t15 = t11*t11; 1200 f = (s+2.0/15.0*r-2.0/15.0)*(t1-2.0/7.0*r+1.0/70.0)*(t6*t+30.0/11.0*t6* 1201 t8+30.0/11.0*t5*t*t11+40.0/33.0*t5*t11*t8+5.0/22.0*t*t15+t15*t8/77); 1202 /* p=12 2 2 4 */ 1203 } else if( i==2 && j==2 && k==4 ) { 1204 t1 = r*r; 1205 t3 = s*s; 1206 t4 = r-1.0; 1207 t6 = t4*t4; 1208 t9 = t*t; 1209 t10 = t9*t9; 1210 t12 = s+r-1.0; 1211 t14 = t12*t12; 1212 t18 = t14*t14; 1213 f = (t1-2.0/7.0*r+1.0/70.0)*(t3+2.0/5.0*s*t4+t6/35)*(t10+20.0/9.0*t9*t* 1214 t12+5.0/3.0*t14*t9+10.0/21.0*t*t14*t12+5.0/126.0*t18); 1215 /* p=12 2 3 3 */ 1216 } else if( i==2 && j==3 && k==3 ) { 1217 t1 = r*r; 1218 t3 = s*s; 1219 t5 = r-1.0; 1220 t7 = t5*t5; 1221 t12 = t*t; 1222 t14 = s+r-1.0; 1223 t16 = t14*t14; 1224 f = (t1-2.0/7.0*r+1.0/70.0)*(t3*s+4.0/5.0*t3*t5+6.0/35.0*s*t7+4.0/455.0 1225 *t7*t5)*(t12*t+12.0/7.0*t14*t12+6.0/7.0*t*t16+4.0/35.0*t16*t14); 1226 /* p=12 2 4 2 */ 1227 } else if( i==2 && j==4 && k==2 ) { 1228 t1 = r*r; 1229 t3 = s*s; 1230 t4 = t3*t3; 1231 t6 = r-1.0; 1232 t8 = t6*t6; 1233 t12 = t8*t8; 1234 t15 = t*t; 1235 t16 = s+r-1.0; 1236 t18 = t16*t16; 1237 f = (t1-2.0/7.0*r+1.0/70.0)*(t4+4.0/3.0*t3*s*t6+4.0/7.0*t3*t8+8.0/91.0* 1238 s*t6*t8+t12/273)*(t15+6.0/5.0*t*t16+3.0/10.0*t18); 1239 /* p=12 2 5 1 */ 1240 } else if( i==2 && j==5 && k==1 ) { 1241 t1 = r*r; 1242 t3 = s*s; 1243 t4 = t3*t3; 1244 t6 = r-1.0; 1245 t9 = t6*t6; 1246 t13 = t9*t9; 1247 f = (t1-2.0/7.0*r+1.0/70.0)*(s*t4+2.0*t4*t6+10.0/7.0*t3*s*t9+40.0/91.0* 1248 t3*t9*t6+5.0/91.0*s*t13+2.0/1001.0*t13*t6)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 1249 /* p=12 2 6 0 */ 1250 } else if( i==2 && j==6 && k==0 ) { 1251 t1 = r*r; 1252 t3 = s*s; 1253 t4 = t3*t3; 1254 t7 = r-1.0; 1255 t9 = t7*t7; 1256 t14 = t9*t9; 1257 f = (t1-2.0/7.0*r+1.0/70.0)*(t4*t3+14.0/5.0*t4*s*t7+3.0*t4*t9+20.0/13.0 1258 *t3*s*t9*t7+5.0/13.0*t3*t14+6.0/143.0*s*t14*t7+t14*t9/715); 1259 /* p=12 3 0 5 */ 1260 } else if( i==3 && j==0 && k==5 ) { 1261 t1 = r*r; 1262 t4 = t*t; 1263 t5 = t4*t4; 1264 t7 = s+r-1.0; 1265 t10 = t7*t7; 1266 t14 = t10*t10; 1267 f = (t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(t5*t+30.0/11.0*t5*t7+30.0/ 1268 11.0*t4*t*t10+40.0/33.0*t4*t10*t7+5.0/22.0*t*t14+t14*t7/77); 1269 /* p=12 3 1 4 */ 1270 } else if( i==3 && j==1 && k==4 ) { 1271 t1 = r*r; 1272 t6 = t*t; 1273 t7 = t6*t6; 1274 t9 = s+r-1.0; 1275 t11 = t9*t9; 1276 t15 = t11*t11; 1277 f = (s+2.0/13.0*r-2.0/13.0)*(t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(t7+ 1278 20.0/9.0*t6*t*t9+5.0/3.0*t6*t11+10.0/21.0*t*t11*t9+5.0/126.0*t15); 1279 /* p=12 3 2 3 */ 1280 } else if( i==3 && j==2 && k==3 ) { 1281 t1 = r*r; 1282 t4 = s*s; 1283 t5 = r-1.0; 1284 t7 = t5*t5; 1285 t10 = t*t; 1286 t12 = s+r-1.0; 1287 t14 = t12*t12; 1288 f = (t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(t4+6.0/13.0*s*t5+t7/26)*( 1289 t10*t+12.0/7.0*t10*t12+6.0/7.0*t*t14+4.0/35.0*t14*t12); 1290 /* p=12 3 3 2 */ 1291 } else if( i==3 && j==3 && k==2 ) { 1292 t1 = r*r; 1293 t4 = s*s; 1294 t6 = r-1.0; 1295 t8 = t6*t6; 1296 t13 = t*t; 1297 t14 = s+r-1.0; 1298 t16 = t14*t14; 1299 f = (t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(t4*s+12.0/13.0*t4*t6+3.0/ 1300 13.0*s*t8+2.0/143.0*t8*t6)*(t13+6.0/5.0*t*t14+3.0/10.0*t16); 1301 /* p=12 3 4 1 */ 1302 } else if( i==3 && j==4 && k==1 ) { 1303 t1 = r*r; 1304 t4 = s*s; 1305 t5 = t4*t4; 1306 t7 = r-1.0; 1307 t9 = t7*t7; 1308 t13 = t9*t9; 1309 f = (t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(t5+20.0/13.0*t4*s*t7+10.0/ 1310 13.0*t4*t9+20.0/143.0*s*t9*t7+t13/143)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 1311 /* p=12 3 5 0 */ 1312 } else if( i==3 && j==5 && k==0 ) { 1313 t1 = r*r; 1314 t4 = s*s; 1315 t5 = t4*t4; 1316 t7 = r-1.0; 1317 t10 = t7*t7; 1318 t14 = t10*t10; 1319 f = (t1*r-4.0/7.0*t1+3.0/35.0*r-2.0/665.0)*(s*t5+30.0/13.0*t5*t7+25.0/ 1320 13.0*t4*s*t10+100.0/143.0*t4*t10*t7+15.0/143.0*s*t14+2.0/429.0*t14*t7); 1321 /* p=12 4 0 4 */ 1322 } else if( i==4 && j==0 && k==4 ) { 1323 t1 = r*r; 1324 t2 = t1*t1; 1325 t5 = t*t; 1326 t6 = t5*t5; 1327 t8 = s+r-1.0; 1328 t10 = t8*t8; 1329 t14 = t10*t10; 1330 f = (t2-20.0/21.0*t1*r+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0)*(t6+20.0/9.0* 1331 t5*t*t8+5.0/3.0*t5*t10+10.0/21.0*t*t10*t8+5.0/126.0*t14); 1332 /* p=12 4 1 3 */ 1333 } else if( i==4 && j==1 && k==3 ) { 1334 t1 = r*r; 1335 t2 = t1*t1; 1336 t7 = t*t; 1337 t9 = s+r-1.0; 1338 t11 = t9*t9; 1339 f = (t2-20.0/21.0*t1*r+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0)*(s+2.0/11.0*r 1340 -2.0/11.0)*(t7*t+12.0/7.0*t9*t7+6.0/7.0*t*t11+4.0/35.0*t11*t9); 1341 /* p=12 4 2 2 */ 1342 } else if( i==4 && j==2 && k==2 ) { 1343 t1 = r*r; 1344 t2 = t1*t1; 1345 t5 = s*s; 1346 t6 = r-1.0; 1347 t8 = t6*t6; 1348 t11 = t*t; 1349 t12 = s+r-1.0; 1350 t14 = t12*t12; 1351 f = (t2-20.0/21.0*t1*r+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0)*(t5+6.0/11.0* 1352 s*t6+3.0/55.0*t8)*(t11+6.0/5.0*t*t12+3.0/10.0*t14); 1353 /* p=12 4 3 1 */ 1354 } else if( i==4 && j==3 && k==1 ) { 1355 t1 = r*r; 1356 t2 = t1*t1; 1357 t5 = s*s; 1358 t7 = r-1.0; 1359 t9 = t7*t7; 1360 f = (t2-20.0/21.0*t1*r+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0)*(t5*s+12.0/ 1361 11.0*t5*t7+18.0/55.0*s*t9+4.0/165.0*t9*t7)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 1362 /* p=12 4 4 0 */ 1363 } else if( i==4 && j==4 && k==0 ) { 1364 t1 = r*r; 1365 t2 = t1*t1; 1366 t5 = s*s; 1367 t6 = t5*t5; 1368 t8 = r-1.0; 1369 t10 = t8*t8; 1370 t14 = t10*t10; 1371 f = (t2-20.0/21.0*t1*r+2.0/7.0*t1-4.0/133.0*r+1.0/1197.0)*(t6+20.0/11.0 1372 *t5*s*t8+12.0/11.0*t5*t10+8.0/33.0*s*t10*t8+t14/66); 1373 /* p=12 5 0 3 */ 1374 } else if( i==5 && j==0 && k==3 ) { 1375 t1 = r*r; 1376 t2 = t1*t1; 1377 t6 = t*t; 1378 t8 = s+r-1.0; 1379 t10 = t8*t8; 1380 f = (t2*r-10.0/7.0*t2+5.0/7.0*t1*r-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0 1381 )*(t6*t+12.0/7.0*t6*t8+6.0/7.0*t*t10+4.0/35.0*t10*t8); 1382 /* p=12 5 1 2 */ 1383 } else if( i==5 && j==1 && k==2 ) { 1384 t1 = r*r; 1385 t2 = t1*t1; 1386 t8 = t*t; 1387 t9 = s+r-1.0; 1388 t11 = t9*t9; 1389 f = (t2*r-10.0/7.0*t2+5.0/7.0*t1*r-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0 1390 )*(s+2.0/9.0*r-2.0/9.0)*(t8+6.0/5.0*t*t9+3.0/10.0*t11); 1391 /* p=12 5 2 1 */ 1392 } else if( i==5 && j==2 && k==1 ) { 1393 t1 = r*r; 1394 t2 = t1*t1; 1395 t6 = s*s; 1396 t7 = r-1.0; 1397 t9 = t7*t7; 1398 f = (t2*r-10.0/7.0*t2+5.0/7.0*t1*r-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0 1399 )*(t6+2.0/3.0*s*t7+t9/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 1400 /* p=12 5 3 0 */ 1401 } else if( i==5 && j==3 && k==0 ) { 1402 t1 = r*r; 1403 t2 = t1*t1; 1404 t6 = s*s; 1405 t8 = r-1.0; 1406 t10 = t8*t8; 1407 f = (t2*r-10.0/7.0*t2+5.0/7.0*t1*r-20.0/133.0*t1+5.0/399.0*r-2.0/6783.0 1408 )*(s*t6+4.0/3.0*t6*t8+s*t10/2+t10*t8/21); 1409 /* p=12 6 0 2 */ 1410 } else if( i==6 && j==0 && k==2 ) { 1411 t1 = r*r; 1412 t2 = t1*t1; 1413 t7 = t*t; 1414 t8 = s+r-1.0; 1415 t10 = t8*t8; 1416 f = (t2*t1-2.0*t2*r+3.0/2.0*t2-10.0/19.0*t1*r+5.0/57.0*t1-2.0/323.0*r+ 1417 1.0/7752.0)*(t7+6.0/5.0*t*t8+3.0/10.0*t10); 1418 /* p=12 6 1 1 */ 1419 } else if( i==6 && j==1 && k==1 ) { 1420 t1 = r*r; 1421 t2 = t1*t1; 1422 f = (t2*t1-2.0*t2*r+3.0/2.0*t2-10.0/19.0*t1*r+5.0/57.0*t1-2.0/323.0*r+ 1423 1.0/7752.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 1424 /* p=12 6 2 0 */ 1425 } else if( i==6 && j==2 && k==0 ) { 1426 t1 = r*r; 1427 t2 = t1*t1; 1428 t7 = s*s; 1429 t8 = r-1.0; 1430 t10 = t8*t8; 1431 f = (t2*t1-2.0*t2*r+3.0/2.0*t2-10.0/19.0*t1*r+5.0/57.0*t1-2.0/323.0*r+ 1432 1.0/7752.0)*(t7+6.0/7.0*s*t8+t10/7); 1433 /* p=12 7 0 1 */ 1434 } else if( i==7 && j==0 && k==1 ) { 1435 t1 = r*r; 1436 t2 = t1*r; 1437 t3 = t1*t1; 1438 f = (t3*t2-8.0/3.0*t3*t1+14.0/5.0*t3*r-28.0/19.0*t3+70.0/171.0*t2-56.0/ 1439 969.0*t1+7.0/1938.0*r-1.0/14535.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 1440 /* p=12 7 1 0 */ 1441 } else if( i==7 && j==1 && k==0 ) { 1442 t1 = r*r; 1443 t2 = t1*r; 1444 t3 = t1*t1; 1445 f = (t3*t2-8.0/3.0*t3*t1+14.0/5.0*t3*r-28.0/19.0*t3+70.0/171.0*t2-56.0/ 1446 969.0*t1+7.0/1938.0*r-1.0/14535.0)*(s+2.0/5.0*r-2.0/5.0); 1447 /* p=12 8 0 0 */ 1448 } else if( i==8 && j==0 && k==0 ) { 1449 t1 = r*r; 1450 t2 = t1*t1; 1451 t3 = t2*t2; 1452 t4 = t1*r; 1453 f = t3-24.0/7.0*t2*t4+24.0/5.0*t2*t1-336.0/95.0*t2*r+28.0/19.0*t2-112.0/ 1454 323.0*t4+14.0/323.0*t1-4.0/1615.0*r+1.0/22610.0; 1455 /* p=13 0 0 9 */ 1456 } else if( i==0 && j==0 && k==9 ) { 1457 t1 = t*t; 1458 t2 = t1*t1; 1459 t3 = t2*t2; 1460 t5 = s+r-1.0; 1461 t7 = t1*t; 1462 t9 = t5*t5; 1463 t12 = t9*t5; 1464 t15 = t9*t9; 1465 t23 = t15*t15; 1466 f = t3*t+90.0/19.0*t3*t5+180.0/19.0*t2*t7*t9+3360.0/323.0*t2*t1*t12+ 1467 2205.0/323.0*t2*t*t15+882.0/323.0*t2*t15*t5+210.0/323.0*t7*t15*t9+360.0/4199.0* 1468 t1*t15*t12+45.0/8398.0*t*t23+5.0/46189.0*t23*t5; 1469 /* p=13 0 1 8 */ 1470 } else if( i==0 && j==1 && k==8 ) { 1471 t2 = t*t; 1472 t3 = t2*t2; 1473 t4 = t3*t3; 1474 t5 = t2*t; 1475 t7 = s+r-1.0; 1476 t10 = t7*t7; 1477 t13 = t10*t7; 1478 t15 = t10*t10; 1479 t23 = t15*t15; 1480 f = (s+2.0/21.0*r-2.0/21.0)*(t4+72.0/17.0*t3*t5*t7+126.0/17.0*t3*t2*t10 1481 +588.0/85.0*t3*t*t13+63.0/17.0*t3*t15+252.0/221.0*t5*t15*t7+42.0/221.0*t2*t15* 1482 t10+36.0/2431.0*t*t15*t13+9.0/24310.0*t23); 1483 /* p=13 0 2 7 */ 1484 } else if( i==0 && j==2 && k==7 ) { 1485 t1 = s*s; 1486 t2 = r-1.0; 1487 t4 = t2*t2; 1488 t6 = t*t; 1489 t7 = t6*t; 1490 t8 = t6*t6; 1491 t11 = s+r-1.0; 1492 t14 = t11*t11; 1493 t16 = t14*t11; 1494 t18 = t14*t14; 1495 f = (t1+2.0/7.0*s*t2+t4/70)*(t7*t8+56.0/15.0*t8*t6*t11+28.0/5.0*t8*t* 1496 t14+56.0/13.0*t8*t16+70.0/39.0*t7*t18+56.0/143.0*t6*t18*t11+28.0/715.0*t*t18* 1497 t14+8.0/6435.0*t18*t16); 1498 /* p=13 0 3 6 */ 1499 } else if( i==0 && j==3 && k==6 ) { 1500 t1 = s*s; 1501 t3 = r-1.0; 1502 t5 = t3*t3; 1503 t9 = t*t; 1504 t10 = t9*t9; 1505 t13 = s+r-1.0; 1506 t15 = t13*t13; 1507 t20 = t15*t15; 1508 f = (t1*s+4.0/7.0*t3*t1+3.0/35.0*s*t5+2.0/665.0*t3*t5)*(t9*t10+42.0/ 1509 13.0*t10*t*t13+105.0/26.0*t15*t10+350.0/143.0*t9*t*t15*t13+105.0/143.0*t9*t20+ 1510 14.0/143.0*t*t20*t13+7.0/1716.0*t20*t15); 1511 /* p=13 0 4 5 */ 1512 } else if( i==0 && j==4 && k==5 ) { 1513 t1 = s*s; 1514 t2 = t1*t1; 1515 t4 = r-1.0; 1516 t6 = t4*t4; 1517 t10 = t6*t6; 1518 t12 = t*t; 1519 t13 = t12*t12; 1520 t15 = s+r-1.0; 1521 t18 = t15*t15; 1522 t22 = t18*t18; 1523 f = (t2+20.0/21.0*t1*s*t4+2.0/7.0*t1*t6+4.0/133.0*s*t6*t4+t10/1197)*( 1524 t13*t+30.0/11.0*t13*t15+30.0/11.0*t12*t*t18+40.0/33.0*t12*t18*t15+5.0/22.0*t* 1525 t22+t22*t15/77); 1526 /* p=13 0 5 4 */ 1527 } else if( i==0 && j==5 && k==4 ) { 1528 t1 = s*s; 1529 t2 = t1*t1; 1530 t4 = r-1.0; 1531 t7 = t4*t4; 1532 t11 = t7*t7; 1533 t15 = t*t; 1534 t16 = t15*t15; 1535 t18 = s+r-1.0; 1536 t20 = t18*t18; 1537 t24 = t20*t20; 1538 f = (t2*s+10.0/7.0*t2*t4+5.0/7.0*t1*s*t7+20.0/133.0*t1*t7*t4+5.0/399.0* 1539 s*t11+2.0/6783.0*t11*t4)*(t16+20.0/9.0*t15*t*t18+5.0/3.0*t20*t15+10.0/21.0*t* 1540 t20*t18+5.0/126.0*t24); 1541 /* p=13 0 6 3 */ 1542 } else if( i==0 && j==6 && k==3 ) { 1543 t1 = s*s; 1544 t2 = t1*t1; 1545 t5 = r-1.0; 1546 t7 = t5*t5; 1547 t12 = t7*t7; 1548 t18 = t*t; 1549 t20 = s+r-1.0; 1550 t22 = t20*t20; 1551 f = (t2*t1+2.0*t2*s*t5+3.0/2.0*t2*t7+10.0/19.0*t1*s*t7*t5+5.0/57.0*t1* 1552 t12+2.0/323.0*s*t12*t5+t12*t7/7752)*(t18*t+12.0/7.0*t20*t18+6.0/7.0*t*t22+4.0/ 1553 35.0*t22*t20); 1554 /* p=13 0 7 2 */ 1555 } else if( i==0 && j==7 && k==2 ) { 1556 t1 = s*s; 1557 t2 = t1*s; 1558 t3 = t1*t1; 1559 t6 = r-1.0; 1560 t9 = t6*t6; 1561 t11 = t9*t6; 1562 t13 = t9*t9; 1563 t21 = t*t; 1564 t22 = s+r-1.0; 1565 t24 = t22*t22; 1566 f = (t3*t2+8.0/3.0*t3*t1*t6+14.0/5.0*t3*s*t9+28.0/19.0*t3*t11+70.0/ 1567 171.0*t2*t13+56.0/969.0*t1*t13*t6+7.0/1938.0*s*t13*t9+t13*t11/14535)*(t21+6.0/ 1568 5.0*t*t22+3.0/10.0*t24); 1569 /* p=13 0 8 1 */ 1570 } else if( i==0 && j==8 && k==1 ) { 1571 t1 = s*s; 1572 t2 = t1*t1; 1573 t3 = t2*t2; 1574 t4 = t1*s; 1575 t6 = r-1.0; 1576 t9 = t6*t6; 1577 t12 = t9*t6; 1578 t14 = t9*t9; 1579 t22 = t14*t14; 1580 f = (t3+24.0/7.0*t2*t4*t6+24.0/5.0*t2*t1*t9+336.0/95.0*t2*s*t12+28.0/ 1581 19.0*t2*t14+112.0/323.0*t4*t14*t6+14.0/323.0*t1*t14*t9+4.0/1615.0*s*t14*t12+t22 1582 /22610)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 1583 /* p=13 0 9 0 */ 1584 } else if( i==0 && j==9 && k==0 ) { 1585 t1 = s*s; 1586 t2 = t1*t1; 1587 t3 = t2*t2; 1588 t5 = r-1.0; 1589 t7 = t1*s; 1590 t9 = t5*t5; 1591 t12 = t9*t5; 1592 t15 = t9*t9; 1593 t23 = t15*t15; 1594 f = t3*s+30.0/7.0*t3*t5+54.0/7.0*t2*t7*t9+144.0/19.0*t2*t1*t12+84.0/ 1595 19.0*t2*s*t15+504.0/323.0*t2*t15*t5+105.0/323.0*t7*t15*t9+12.0/323.0*t1*t15*t12 1596 +9.0/4522.0*s*t23+t23*t5/29393; 1597 /* p=13 1 0 8 */ 1598 } else if( i==1 && j==0 && k==8 ) { 1599 t2 = t*t; 1600 t3 = t2*t2; 1601 t4 = t3*t3; 1602 t5 = t2*t; 1603 t7 = s+r-1.0; 1604 t10 = t7*t7; 1605 t13 = t10*t7; 1606 t15 = t10*t10; 1607 t23 = t15*t15; 1608 f = (r-2.0/23.0)*(t4+72.0/17.0*t3*t5*t7+126.0/17.0*t3*t2*t10+588.0/85.0 1609 *t3*t*t13+63.0/17.0*t3*t15+252.0/221.0*t5*t15*t7+42.0/221.0*t2*t15*t10+36.0/ 1610 2431.0*t*t13*t15+9.0/24310.0*t23); 1611 /* p=13 1 1 7 */ 1612 } else if( i==1 && j==1 && k==7 ) { 1613 t4 = t*t; 1614 t5 = t4*t; 1615 t6 = t4*t4; 1616 t9 = s+r-1.0; 1617 t12 = t9*t9; 1618 t14 = t12*t9; 1619 t16 = t12*t12; 1620 f = (r-2.0/23.0)*(s+2.0/19.0*r-2.0/19.0)*(t6*t5+56.0/15.0*t6*t4*t9+28.0 1621 /5.0*t6*t*t12+56.0/13.0*t6*t14+70.0/39.0*t5*t16+56.0/143.0*t4*t16*t9+28.0/715.0 1622 *t*t16*t12+8.0/6435.0*t16*t14); 1623 /* p=13 1 2 6 */ 1624 } else if( i==1 && j==2 && k==6 ) { 1625 t2 = s*s; 1626 t3 = r-1.0; 1627 t5 = t3*t3; 1628 t8 = t*t; 1629 t9 = t8*t8; 1630 t12 = s+r-1.0; 1631 t14 = t12*t12; 1632 t19 = t14*t14; 1633 f = (r-2.0/23.0)*(t2+6.0/19.0*s*t3+t5/57)*(t8*t9+42.0/13.0*t9*t*t12+ 1634 105.0/26.0*t9*t14+350.0/143.0*t8*t*t14*t12+105.0/143.0*t8*t19+14.0/143.0*t*t19* 1635 t12+7.0/1716.0*t19*t14); 1636 /* p=13 1 3 5 */ 1637 } else if( i==1 && j==3 && k==5 ) { 1638 t2 = s*s; 1639 t4 = r-1.0; 1640 t6 = t4*t4; 1641 t11 = t*t; 1642 t12 = t11*t11; 1643 t14 = s+r-1.0; 1644 t17 = t14*t14; 1645 t21 = t17*t17; 1646 f = (r-2.0/23.0)*(t2*s+12.0/19.0*t2*t4+2.0/19.0*s*t6+4.0/969.0*t6*t4)*( 1647 t12*t+30.0/11.0*t14*t12+30.0/11.0*t11*t*t17+40.0/33.0*t11*t17*t14+5.0/22.0*t* 1648 t21+t21*t14/77); 1649 /* p=13 1 4 4 */ 1650 } else if( i==1 && j==4 && k==4 ) { 1651 t2 = s*s; 1652 t3 = t2*t2; 1653 t5 = r-1.0; 1654 t7 = t5*t5; 1655 t11 = t7*t7; 1656 t14 = t*t; 1657 t15 = t14*t14; 1658 t17 = s+r-1.0; 1659 t19 = t17*t17; 1660 t23 = t19*t19; 1661 f = (r-2.0/23.0)*(t3+20.0/19.0*t2*s*t5+20.0/57.0*t2*t7+40.0/969.0*s*t7* 1662 t5+5.0/3876.0*t11)*(t15+20.0/9.0*t14*t*t17+5.0/3.0*t19*t14+10.0/21.0*t*t19*t17+ 1663 5.0/126.0*t23); 1664 /* p=13 1 5 3 */ 1665 } else if( i==1 && j==5 && k==3 ) { 1666 t2 = s*s; 1667 t3 = t2*t2; 1668 t5 = r-1.0; 1669 t8 = t5*t5; 1670 t12 = t8*t8; 1671 t17 = t*t; 1672 t19 = s+r-1.0; 1673 t21 = t19*t19; 1674 f = (r-2.0/23.0)*(s*t3+30.0/19.0*t3*t5+50.0/57.0*t2*s*t8+200.0/969.0*t2 1675 *t8*t5+25.0/1292.0*s*t12+t12*t5/1938)*(t17*t+12.0/7.0*t19*t17+6.0/7.0*t*t21+4.0 1676 /35.0*t21*t19); 1677 /* p=13 1 6 2 */ 1678 } else if( i==1 && j==6 && k==2 ) { 1679 t2 = s*s; 1680 t3 = t2*t2; 1681 t6 = r-1.0; 1682 t8 = t6*t6; 1683 t13 = t8*t8; 1684 t20 = t*t; 1685 t21 = s+r-1.0; 1686 t23 = t21*t21; 1687 f = (r-2.0/23.0)*(t3*t2+42.0/19.0*t6*t3*s+35.0/19.0*t3*t8+700.0/969.0* 1688 t2*s*t8*t6+175.0/1292.0*t2*t13+7.0/646.0*s*t13*t6+t13*t8/3876)*(t20+6.0/5.0*t* 1689 t21+3.0/10.0*t23); 1690 /* p=13 1 7 1 */ 1691 } else if( i==1 && j==7 && k==1 ) { 1692 t2 = s*s; 1693 t3 = t2*s; 1694 t4 = t2*t2; 1695 t7 = r-1.0; 1696 t10 = t7*t7; 1697 t12 = t10*t7; 1698 t14 = t10*t10; 1699 f = (r-2.0/23.0)*(t4*t3+56.0/19.0*t4*t2*t7+196.0/57.0*t4*s*t10+1960.0/ 1700 969.0*t4*t12+1225.0/1938.0*t3*t14+98.0/969.0*t2*t14*t7+7.0/969.0*s*t14*t10+2.0/ 1701 12597.0*t14*t12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 1702 /* p=13 1 8 0 */ 1703 } else if( i==1 && j==8 && k==0 ) { 1704 t2 = s*s; 1705 t3 = t2*t2; 1706 t4 = t3*t3; 1707 t5 = t2*s; 1708 t7 = r-1.0; 1709 t10 = t7*t7; 1710 t13 = t10*t7; 1711 t15 = t10*t10; 1712 t23 = t15*t15; 1713 f = (r-2.0/23.0)*(t4+72.0/19.0*t3*t5*t7+112.0/19.0*t3*t2*t10+1568.0/ 1714 323.0*t3*s*t13+735.0/323.0*t3*t15+196.0/323.0*t5*t15*t7+28.0/323.0*t2*t15*t10+ 1715 24.0/4199.0*s*t15*t13+t23/8398); 1716 /* p=13 2 0 7 */ 1717 } else if( i==2 && j==0 && k==7 ) { 1718 t1 = r*r; 1719 t3 = t*t; 1720 t4 = t3*t; 1721 t5 = t3*t3; 1722 t8 = s+r-1.0; 1723 t11 = t8*t8; 1724 t13 = t11*t8; 1725 t15 = t11*t11; 1726 f = (t1-6.0/23.0*r+3.0/253.0)*(t5*t4+56.0/15.0*t3*t5*t8+28.0/5.0*t5*t* 1727 t11+56.0/13.0*t5*t13+70.0/39.0*t4*t15+56.0/143.0*t3*t15*t8+28.0/715.0*t*t15*t11 1728 +8.0/6435.0*t15*t13); 1729 /* p=13 2 1 6 */ 1730 } else if( i==2 && j==1 && k==6 ) { 1731 t1 = r*r; 1732 t5 = t*t; 1733 t6 = t5*t5; 1734 t9 = s+r-1.0; 1735 t11 = t9*t9; 1736 t16 = t11*t11; 1737 f = (s+2.0/17.0*r-2.0/17.0)*(t1-6.0/23.0*r+3.0/253.0)*(t6*t5+42.0/13.0* 1738 t6*t*t9+105.0/26.0*t6*t11+350.0/143.0*t5*t*t11*t9+105.0/143.0*t5*t16+14.0/143.0 1739 *t*t16*t9+7.0/1716.0*t16*t11); 1740 /* p=13 2 2 5 */ 1741 } else if( i==2 && j==2 && k==5 ) { 1742 t1 = r*r; 1743 t3 = s*s; 1744 t4 = r-1.0; 1745 t6 = t4*t4; 1746 t9 = t*t; 1747 t10 = t9*t9; 1748 t12 = s+r-1.0; 1749 t15 = t12*t12; 1750 t19 = t15*t15; 1751 f = (t1-6.0/23.0*r+3.0/253.0)*(t3+6.0/17.0*s*t4+3.0/136.0*t6)*(t10*t+ 1752 30.0/11.0*t10*t12+30.0/11.0*t9*t*t15+40.0/33.0*t9*t15*t12+5.0/22.0*t*t19+t19* 1753 t12/77); 1754 /* p=13 2 3 4 */ 1755 } else if( i==2 && j==3 && k==4 ) { 1756 t1 = r*r; 1757 t3 = s*s; 1758 t5 = r-1.0; 1759 t7 = t5*t5; 1760 t12 = t*t; 1761 t13 = t12*t12; 1762 t15 = s+r-1.0; 1763 t17 = t15*t15; 1764 t21 = t17*t17; 1765 f = (t1-6.0/23.0*r+3.0/253.0)*(t3*s+12.0/17.0*t3*t5+9.0/68.0*s*t7+t7*t5 1766 /170)*(t13+20.0/9.0*t12*t*t15+5.0/3.0*t12*t17+10.0/21.0*t*t17*t15+5.0/126.0*t21 1767 ); 1768 /* p=13 2 4 3 */ 1769 } else if( i==2 && j==4 && k==3 ) { 1770 t1 = r*r; 1771 t3 = s*s; 1772 t4 = t3*t3; 1773 t6 = r-1.0; 1774 t8 = t6*t6; 1775 t12 = t8*t8; 1776 t15 = t*t; 1777 t17 = s+r-1.0; 1778 t19 = t17*t17; 1779 f = (t1-6.0/23.0*r+3.0/253.0)*(t4+20.0/17.0*t6*t3*s+15.0/34.0*t3*t8+s* 1780 t8*t6/17+t12/476)*(t15*t+12.0/7.0*t17*t15+6.0/7.0*t*t19+4.0/35.0*t19*t17); 1781 /* p=13 2 5 2 */ 1782 } else if( i==2 && j==5 && k==2 ) { 1783 t1 = r*r; 1784 t3 = s*s; 1785 t4 = t3*t3; 1786 t6 = r-1.0; 1787 t9 = t6*t6; 1788 t13 = t9*t9; 1789 t18 = t*t; 1790 t19 = s+r-1.0; 1791 t21 = t19*t19; 1792 f = (t1-6.0/23.0*r+3.0/253.0)*(t4*s+30.0/17.0*t4*t6+75.0/68.0*t3*s*t9+ 1793 5.0/17.0*t3*t9*t6+15.0/476.0*s*t13+3.0/3094.0*t13*t6)*(t18+6.0/5.0*t*t19+3.0/ 1794 10.0*t21); 1795 /* p=13 2 6 1 */ 1796 } else if( i==2 && j==6 && k==1 ) { 1797 t1 = r*r; 1798 t3 = s*s; 1799 t4 = t3*t3; 1800 t7 = r-1.0; 1801 t9 = t7*t7; 1802 t14 = t9*t9; 1803 f = (t1-6.0/23.0*r+3.0/253.0)*(t4*t3+42.0/17.0*t4*s*t7+315.0/136.0*t4* 1804 t9+35.0/34.0*t3*s*t9*t7+15.0/68.0*t3*t14+9.0/442.0*s*t14*t7+t14*t9/1768)*(t+2.0 1805 /3.0*s+2.0/3.0*r-2.0/3.0); 1806 /* p=13 2 7 0 */ 1807 } else if( i==2 && j==7 && k==0 ) { 1808 t1 = r*r; 1809 t3 = s*s; 1810 t4 = t3*s; 1811 t5 = t3*t3; 1812 t8 = r-1.0; 1813 t11 = t8*t8; 1814 t13 = t11*t8; 1815 t15 = t11*t11; 1816 f = (t1-6.0/23.0*r+3.0/253.0)*(t5*t4+56.0/17.0*t5*t3*t8+147.0/34.0*t5*s 1817 *t11+49.0/17.0*t5*t13+35.0/34.0*t4*t15+42.0/221.0*t3*t15*t8+7.0/442.0*s*t15*t11 1818 +t15*t13/2431); 1819 /* p=13 3 0 6 */ 1820 } else if( i==3 && j==0 && k==6 ) { 1821 t1 = r*r; 1822 t4 = t*t; 1823 t5 = t4*t4; 1824 t8 = s+r-1.0; 1825 t10 = t8*t8; 1826 t15 = t10*t10; 1827 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t5*t4+42.0/13.0*t5*t* 1828 t8+105.0/26.0*t5*t10+350.0/143.0*t4*t*t10*t8+105.0/143.0*t4*t15+14.0/143.0*t* 1829 t15*t8+7.0/1716.0*t15*t10); 1830 /* p=13 3 1 5 */ 1831 } else if( i==3 && j==1 && k==5 ) { 1832 t1 = r*r; 1833 t6 = t*t; 1834 t7 = t6*t6; 1835 t9 = s+r-1.0; 1836 t12 = t9*t9; 1837 t16 = t12*t12; 1838 f = (s+2.0/15.0*r-2.0/15.0)*(t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0) 1839 *(t7*t+30.0/11.0*t9*t7+30.0/11.0*t6*t*t12+40.0/33.0*t6*t12*t9+5.0/22.0*t*t16+ 1840 t16*t9/77); 1841 /* p=13 3 2 4 */ 1842 } else if( i==3 && j==2 && k==4 ) { 1843 t1 = r*r; 1844 t4 = s*s; 1845 t5 = r-1.0; 1846 t7 = t5*t5; 1847 t10 = t*t; 1848 t11 = t10*t10; 1849 t13 = s+r-1.0; 1850 t15 = t13*t13; 1851 t19 = t15*t15; 1852 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t4+2.0/5.0*s*t5+t7/35) 1853 *(t11+20.0/9.0*t10*t*t13+5.0/3.0*t15*t10+10.0/21.0*t*t15*t13+5.0/126.0*t19); 1854 /* p=13 3 3 3 */ 1855 } else if( i==3 && j==3 && k==3 ) { 1856 t1 = r*r; 1857 t4 = s*s; 1858 t6 = r-1.0; 1859 t8 = t6*t6; 1860 t13 = t*t; 1861 t15 = s+r-1.0; 1862 t17 = t15*t15; 1863 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t4*s+4.0/5.0*t4*t6+6.0 1864 /35.0*s*t8+4.0/455.0*t8*t6)*(t13*t+12.0/7.0*t15*t13+6.0/7.0*t*t17+4.0/35.0*t17* 1865 t15); 1866 /* p=13 3 4 2 */ 1867 } else if( i==3 && j==4 && k==2 ) { 1868 t1 = r*r; 1869 t4 = s*s; 1870 t5 = t4*t4; 1871 t7 = r-1.0; 1872 t9 = t7*t7; 1873 t13 = t9*t9; 1874 t16 = t*t; 1875 t17 = s+r-1.0; 1876 t19 = t17*t17; 1877 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t5+4.0/3.0*t4*s*t7+4.0 1878 /7.0*t4*t9+8.0/91.0*s*t9*t7+t13/273)*(t16+6.0/5.0*t*t17+3.0/10.0*t19); 1879 /* p=13 3 5 1 */ 1880 } else if( i==3 && j==5 && k==1 ) { 1881 t1 = r*r; 1882 t4 = s*s; 1883 t5 = t4*t4; 1884 t7 = r-1.0; 1885 t10 = t7*t7; 1886 t14 = t10*t10; 1887 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t5*s+2.0*t5*t7+10.0/ 1888 7.0*t4*s*t10+40.0/91.0*t4*t10*t7+5.0/91.0*s*t14+2.0/1001.0*t14*t7)*(t+2.0/3.0*s 1889 +2.0/3.0*r-2.0/3.0); 1890 /* p=13 3 6 0 */ 1891 } else if( i==3 && j==6 && k==0 ) { 1892 t1 = r*r; 1893 t4 = s*s; 1894 t5 = t4*t4; 1895 t8 = r-1.0; 1896 t10 = t8*t8; 1897 t15 = t10*t10; 1898 f = (t1*r-12.0/23.0*t1+18.0/253.0*r-4.0/1771.0)*(t5*t4+14.0/5.0*t5*s*t8 1899 +3.0*t5*t10+20.0/13.0*t4*s*t10*t8+5.0/13.0*t4*t15+6.0/143.0*s*t15*t8+t15*t10 1900 /715); 1901 /* p=13 4 0 5 */ 1902 } else if( i==4 && j==0 && k==5 ) { 1903 t1 = r*r; 1904 t2 = t1*t1; 1905 t5 = t*t; 1906 t6 = t5*t5; 1907 t8 = s+r-1.0; 1908 t11 = t8*t8; 1909 t15 = t11*t11; 1910 f = (t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0)*(t6*t+ 1911 30.0/11.0*t6*t8+30.0/11.0*t5*t*t11+40.0/33.0*t5*t11*t8+5.0/22.0*t*t15+t15*t8/77 1912 ); 1913 /* p=13 4 1 4 */ 1914 } else if( i==4 && j==1 && k==4 ) { 1915 t1 = r*r; 1916 t2 = t1*t1; 1917 t7 = t*t; 1918 t8 = t7*t7; 1919 t10 = s+r-1.0; 1920 t12 = t10*t10; 1921 t16 = t12*t12; 1922 f = (s+2.0/13.0*r-2.0/13.0)*(t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/ 1923 1771.0*r+1.0/1771.0)*(t8+20.0/9.0*t7*t*t10+5.0/3.0*t7*t12+10.0/21.0*t*t12*t10+ 1924 5.0/126.0*t16); 1925 /* p=13 4 2 3 */ 1926 } else if( i==4 && j==2 && k==3 ) { 1927 t1 = r*r; 1928 t2 = t1*t1; 1929 t5 = s*s; 1930 t6 = r-1.0; 1931 t8 = t6*t6; 1932 t11 = t*t; 1933 t13 = s+r-1.0; 1934 t15 = t13*t13; 1935 f = (t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0)*(t5+6.0/ 1936 13.0*s*t6+t8/26)*(t11*t+12.0/7.0*t11*t13+6.0/7.0*t*t15+4.0/35.0*t15*t13); 1937 /* p=13 4 3 2 */ 1938 } else if( i==4 && j==3 && k==2 ) { 1939 t1 = r*r; 1940 t2 = t1*t1; 1941 t5 = s*s; 1942 t7 = r-1.0; 1943 t9 = t7*t7; 1944 t14 = t*t; 1945 t15 = s+r-1.0; 1946 t17 = t15*t15; 1947 f = (t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0)*(t5*s+ 1948 12.0/13.0*t5*t7+3.0/13.0*s*t9+2.0/143.0*t9*t7)*(t14+6.0/5.0*t*t15+3.0/10.0*t17) 1949 ; 1950 /* p=13 4 4 1 */ 1951 } else if( i==4 && j==4 && k==1 ) { 1952 t1 = r*r; 1953 t2 = t1*t1; 1954 t5 = s*s; 1955 t6 = t5*t5; 1956 t8 = r-1.0; 1957 t10 = t8*t8; 1958 t14 = t10*t10; 1959 f = (t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0)*(t6+20.0 1960 /13.0*t5*s*t8+10.0/13.0*t5*t10+20.0/143.0*s*t10*t8+t14/143)*(t+2.0/3.0*s+2.0/ 1961 3.0*r-2.0/3.0); 1962 /* p=13 4 5 0 */ 1963 } else if( i==4 && j==5 && k==0 ) { 1964 t1 = r*r; 1965 t2 = t1*t1; 1966 t5 = s*s; 1967 t6 = t5*t5; 1968 t8 = r-1.0; 1969 t11 = t8*t8; 1970 t15 = t11*t11; 1971 f = (t2-20.0/23.0*t1*r+60.0/253.0*t1-40.0/1771.0*r+1.0/1771.0)*(t6*s+ 1972 30.0/13.0*t6*t8+25.0/13.0*t5*s*t11+100.0/143.0*t5*t11*t8+15.0/143.0*s*t15+2.0/ 1973 429.0*t15*t8); 1974 /* p=13 5 0 4 */ 1975 } else if( i==5 && j==0 && k==4 ) { 1976 t1 = r*r; 1977 t2 = t1*t1; 1978 t6 = t*t; 1979 t7 = t6*t6; 1980 t9 = s+r-1.0; 1981 t11 = t9*t9; 1982 t15 = t11*t11; 1983 f = (t2*r-30.0/23.0*t2+150.0/253.0*t1*r-200.0/1771.0*t1+15.0/1771.0*r 1984 -6.0/33649.0)*(t7+20.0/9.0*t6*t*t9+5.0/3.0*t6*t11+10.0/21.0*t*t11*t9+5.0/126.0* 1985 t15); 1986 /* p=13 5 1 3 */ 1987 } else if( i==5 && j==1 && k==3 ) { 1988 t1 = r*r; 1989 t2 = t1*t1; 1990 t8 = t*t; 1991 t10 = s+r-1.0; 1992 t12 = t10*t10; 1993 f = (t2*r-30.0/23.0*t2+150.0/253.0*t1*r-200.0/1771.0*t1+15.0/1771.0*r 1994 -6.0/33649.0)*(s+2.0/11.0*r-2.0/11.0)*(t8*t+12.0/7.0*t8*t10+6.0/7.0*t*t12+4.0/ 1995 35.0*t12*t10); 1996 /* p=13 5 2 2 */ 1997 } else if( i==5 && j==2 && k==2 ) { 1998 t1 = r*r; 1999 t2 = t1*t1; 2000 t6 = s*s; 2001 t7 = r-1.0; 2002 t9 = t7*t7; 2003 t12 = t*t; 2004 t13 = s+r-1.0; 2005 t15 = t13*t13; 2006 f = (t2*r-30.0/23.0*t2+150.0/253.0*t1*r-200.0/1771.0*t1+15.0/1771.0*r 2007 -6.0/33649.0)*(t6+6.0/11.0*s*t7+3.0/55.0*t9)*(t12+6.0/5.0*t*t13+3.0/10.0*t15); 2008 /* p=13 5 3 1 */ 2009 } else if( i==5 && j==3 && k==1 ) { 2010 t1 = r*r; 2011 t2 = t1*t1; 2012 t6 = s*s; 2013 t8 = r-1.0; 2014 t10 = t8*t8; 2015 f = (t2*r-30.0/23.0*t2+150.0/253.0*t1*r-200.0/1771.0*t1+15.0/1771.0*r 2016 -6.0/33649.0)*(t6*s+12.0/11.0*t6*t8+18.0/55.0*s*t10+4.0/165.0*t8*t10)*(t+2.0/ 2017 3.0*s+2.0/3.0*r-2.0/3.0); 2018 /* p=13 5 4 0 */ 2019 } else if( i==5 && j==4 && k==0 ) { 2020 t1 = r*r; 2021 t2 = t1*t1; 2022 t6 = s*s; 2023 t7 = t6*t6; 2024 t9 = r-1.0; 2025 t11 = t9*t9; 2026 t15 = t11*t11; 2027 f = (t2*r-30.0/23.0*t2+150.0/253.0*t1*r-200.0/1771.0*t1+15.0/1771.0*r 2028 -6.0/33649.0)*(t7+20.0/11.0*t6*s*t9+12.0/11.0*t6*t11+8.0/33.0*s*t11*t9+t15/66); 2029 /* p=13 6 0 3 */ 2030 } else if( i==6 && j==0 && k==3 ) { 2031 t1 = r*r; 2032 t2 = t1*t1; 2033 t7 = t*t; 2034 t9 = s+r-1.0; 2035 t11 = t9*t9; 2036 f = (t2*t1-42.0/23.0*t2*r+315.0/253.0*t2-100.0/253.0*t1*r+15.0/253.0*t1 2037 -18.0/4807.0*r+1.0/14421.0)*(t7*t+12.0/7.0*t7*t9+6.0/7.0*t*t11+4.0/35.0*t11*t9) 2038 ; 2039 /* p=13 6 1 2 */ 2040 } else if( i==6 && j==1 && k==2 ) { 2041 t1 = r*r; 2042 t2 = t1*t1; 2043 t9 = t*t; 2044 t10 = s+r-1.0; 2045 t12 = t10*t10; 2046 f = (t2*t1-42.0/23.0*t2*r+315.0/253.0*t2-100.0/253.0*t1*r+15.0/253.0*t1 2047 -18.0/4807.0*r+1.0/14421.0)*(s+2.0/9.0*r-2.0/9.0)*(t9+6.0/5.0*t*t10+3.0/10.0* 2048 t12); 2049 /* p=13 6 2 1 */ 2050 } else if( i==6 && j==2 && k==1 ) { 2051 t1 = r*r; 2052 t2 = t1*t1; 2053 t7 = s*s; 2054 t8 = r-1.0; 2055 t10 = t8*t8; 2056 f = (t2*t1-42.0/23.0*t2*r+315.0/253.0*t2-100.0/253.0*t1*r+15.0/253.0*t1 2057 -18.0/4807.0*r+1.0/14421.0)*(t7+2.0/3.0*s*t8+t10/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0 2058 /3.0); 2059 /* p=13 6 3 0 */ 2060 } else if( i==6 && j==3 && k==0 ) { 2061 t1 = r*r; 2062 t2 = t1*t1; 2063 t7 = s*s; 2064 t9 = r-1.0; 2065 t11 = t9*t9; 2066 f = (t2*t1-42.0/23.0*t2*r+315.0/253.0*t2-100.0/253.0*t1*r+15.0/253.0*t1 2067 -18.0/4807.0*r+1.0/14421.0)*(t7*s+4.0/3.0*t7*t9+s*t11/2+t11*t9/21); 2068 /* p=13 7 0 2 */ 2069 } else if( i==7 && j==0 && k==2 ) { 2070 t1 = r*r; 2071 t2 = t1*r; 2072 t3 = t1*t1; 2073 t8 = t*t; 2074 t9 = s+r-1.0; 2075 t11 = t9*t9; 2076 f = (t3*t2-56.0/23.0*t3*t1+588.0/253.0*t3*r-280.0/253.0*t3+70.0/253.0* 2077 t2-168.0/4807.0*t1+28.0/14421.0*r-8.0/245157.0)*(t8+6.0/5.0*t*t9+3.0/10.0*t11); 2078 /* p=13 7 1 1 */ 2079 } else if( i==7 && j==1 && k==1 ) { 2080 t1 = r*r; 2081 t2 = t1*r; 2082 t3 = t1*t1; 2083 f = (t3*t2-56.0/23.0*t3*t1+588.0/253.0*t3*r-280.0/253.0*t3+70.0/253.0* 2084 t2-168.0/4807.0*t1+28.0/14421.0*r-8.0/245157.0)*(s+2.0/7.0*r-2.0/7.0)*(t+2.0/ 2085 3.0*s+2.0/3.0*r-2.0/3.0); 2086 /* p=13 7 2 0 */ 2087 } else if( i==7 && j==2 && k==0 ) { 2088 t1 = r*r; 2089 t2 = t1*r; 2090 t3 = t1*t1; 2091 t8 = s*s; 2092 t9 = r-1.0; 2093 t11 = t9*t9; 2094 f = (t3*t2-56.0/23.0*t3*t1+588.0/253.0*t3*r-280.0/253.0*t3+70.0/253.0* 2095 t2-168.0/4807.0*t1+28.0/14421.0*r-8.0/245157.0)*(t8+6.0/7.0*s*t9+t11/7); 2096 /* p=13 8 0 1 */ 2097 } else if( i==8 && j==0 && k==1 ) { 2098 t1 = r*r; 2099 t2 = t1*t1; 2100 t3 = t2*t2; 2101 t4 = t1*r; 2102 f = (t3-72.0/23.0*t2*t4+1008.0/253.0*t2*t1-672.0/253.0*t2*r+252.0/253.0 2103 *t2-1008.0/4807.0*t4+112.0/4807.0*t1-96.0/81719.0*r+3.0/163438.0)*(t+2.0/3.0*s+ 2104 2.0/3.0*r-2.0/3.0); 2105 /* p=13 8 1 0 */ 2106 } else if( i==8 && j==1 && k==0 ) { 2107 t1 = r*r; 2108 t2 = t1*t1; 2109 t3 = t2*t2; 2110 t4 = t1*r; 2111 f = (t3-72.0/23.0*t2*t4+1008.0/253.0*t2*t1-672.0/253.0*t2*r+252.0/253.0 2112 *t2-1008.0/4807.0*t4+112.0/4807.0*t1-96.0/81719.0*r+3.0/163438.0)*(s+2.0/5.0*r 2113 -2.0/5.0); 2114 /* p=13 9 0 0 */ 2115 } else if( i==9 && j==0 && k==0 ) { 2116 t1 = r*r; 2117 t2 = t1*t1; 2118 t3 = t2*t2; 2119 t5 = t1*r; 2120 f = t3*r-90.0/23.0*t3+1620.0/253.0*t2*t5-1440.0/253.0*t2*t1+756.0/253.0* 2121 t2*r-4536.0/4807.0*t2+840.0/4807.0*t5-1440.0/81719.0*t1+135.0/163438.0*r-1.0/ 2122 81719.0; 2123 /* p=14 0 0 10 */ 2124 } else if( i==0 && j==0 && k==10 ) { 2125 t1 = t*t; 2126 t2 = t1*t1; 2127 t3 = t2*t2; 2128 t6 = s+r-1.0; 2129 t8 = t6*t6; 2130 t10 = t1*t; 2131 t12 = t8*t6; 2132 t15 = t8*t8; 2133 t24 = t15*t15; 2134 f = t3*t1+110.0/21.0*t3*t*t6+165.0/14.0*t3*t8+1980.0/133.0*t2*t10*t12+ 2135 220.0/19.0*t2*t1*t15+1848.0/323.0*t2*t*t15*t6+1155.0/646.0*t2*t15*t8+110.0/ 2136 323.0*t10*t15*t12+165.0/4522.0*t1*t24+55.0/29393.0*t*t24*t6+11.0/352716.0*t24* 2137 t8; 2138 /* p=14 0 1 9 */ 2139 } else if( i==0 && j==1 && k==9 ) { 2140 t2 = t*t; 2141 t3 = t2*t2; 2142 t4 = t3*t3; 2143 t6 = s+r-1.0; 2144 t8 = t2*t; 2145 t10 = t6*t6; 2146 t13 = t10*t6; 2147 t16 = t10*t10; 2148 t24 = t16*t16; 2149 f = (s+2.0/23.0*r-2.0/23.0)*(t4*t+90.0/19.0*t4*t6+180.0/19.0*t3*t8*t10+ 2150 3360.0/323.0*t3*t2*t13+2205.0/323.0*t3*t*t16+882.0/323.0*t3*t16*t6+210.0/323.0* 2151 t8*t16*t10+360.0/4199.0*t2*t16*t13+45.0/8398.0*t*t24+5.0/46189.0*t24*t6); 2152 /* p=14 0 2 8 */ 2153 } else if( i==0 && j==2 && k==8 ) { 2154 t1 = s*s; 2155 t2 = r-1.0; 2156 t4 = t2*t2; 2157 t6 = t*t; 2158 t7 = t6*t6; 2159 t8 = t7*t7; 2160 t9 = t6*t; 2161 t11 = s+r-1.0; 2162 t14 = t11*t11; 2163 t17 = t14*t11; 2164 t19 = t14*t14; 2165 t27 = t19*t19; 2166 f = (t1+6.0/23.0*s*t2+3.0/253.0*t4)*(t8+72.0/17.0*t7*t9*t11+126.0/17.0* 2167 t7*t6*t14+588.0/85.0*t7*t*t17+63.0/17.0*t7*t19+252.0/221.0*t9*t19*t11+42.0/ 2168 221.0*t6*t19*t14+36.0/2431.0*t*t19*t17+9.0/24310.0*t27); 2169 /* p=14 0 3 7 */ 2170 } else if( i==0 && j==3 && k==7 ) { 2171 t1 = s*s; 2172 t3 = r-1.0; 2173 t5 = t3*t3; 2174 t9 = t*t; 2175 t10 = t9*t; 2176 t11 = t9*t9; 2177 t14 = s+r-1.0; 2178 t17 = t14*t14; 2179 t19 = t17*t14; 2180 t21 = t17*t17; 2181 f = (t1*s+12.0/23.0*t3*t1+18.0/253.0*s*t5+4.0/1771.0*t5*t3)*(t10*t11+ 2182 56.0/15.0*t11*t9*t14+28.0/5.0*t11*t*t17+56.0/13.0*t19*t11+70.0/39.0*t10*t21+ 2183 56.0/143.0*t9*t21*t14+28.0/715.0*t*t21*t17+8.0/6435.0*t21*t19); 2184 /* p=14 0 4 6 */ 2185 } else if( i==0 && j==4 && k==6 ) { 2186 t1 = s*s; 2187 t2 = t1*t1; 2188 t4 = r-1.0; 2189 t6 = t4*t4; 2190 t10 = t6*t6; 2191 t12 = t*t; 2192 t13 = t12*t12; 2193 t16 = s+r-1.0; 2194 t18 = t16*t16; 2195 t23 = t18*t18; 2196 f = (t2+20.0/23.0*t1*s*t4+60.0/253.0*t1*t6+40.0/1771.0*s*t4*t6+t10/1771 2197 )*(t12*t13+42.0/13.0*t13*t*t16+105.0/26.0*t13*t18+350.0/143.0*t12*t*t18*t16+ 2198 105.0/143.0*t12*t23+14.0/143.0*t*t23*t16+7.0/1716.0*t23*t18); 2199 /* p=14 0 5 5 */ 2200 } else if( i==0 && j==5 && k==5 ) { 2201 t1 = s*s; 2202 t2 = t1*t1; 2203 t4 = r-1.0; 2204 t7 = t4*t4; 2205 t11 = t7*t7; 2206 t15 = t*t; 2207 t16 = t15*t15; 2208 t18 = s+r-1.0; 2209 t21 = t18*t18; 2210 t25 = t21*t21; 2211 f = (s*t2+30.0/23.0*t2*t4+150.0/253.0*t7*t1*s+200.0/1771.0*t1*t7*t4+ 2212 15.0/1771.0*s*t11+6.0/33649.0*t11*t4)*(t16*t+30.0/11.0*t18*t16+30.0/11.0*t15*t* 2213 t21+40.0/33.0*t15*t21*t18+5.0/22.0*t*t25+t25*t18/77); 2214 /* p=14 0 6 4 */ 2215 } else if( i==0 && j==6 && k==4 ) { 2216 t1 = s*s; 2217 t2 = t1*t1; 2218 t5 = r-1.0; 2219 t7 = t5*t5; 2220 t12 = t7*t7; 2221 t18 = t*t; 2222 t19 = t18*t18; 2223 t21 = s+r-1.0; 2224 t23 = t21*t21; 2225 t27 = t23*t23; 2226 f = (t2*t1+42.0/23.0*t2*s*t5+315.0/253.0*t2*t7+100.0/253.0*t1*s*t7*t5+ 2227 15.0/253.0*t1*t12+18.0/4807.0*s*t12*t5+t12*t7/14421)*(t19+20.0/9.0*t18*t*t21+ 2228 5.0/3.0*t23*t18+10.0/21.0*t*t23*t21+5.0/126.0*t27); 2229 /* p=14 0 7 3 */ 2230 } else if( i==0 && j==7 && k==3 ) { 2231 t1 = s*s; 2232 t2 = t1*s; 2233 t3 = t1*t1; 2234 t6 = r-1.0; 2235 t9 = t6*t6; 2236 t11 = t9*t6; 2237 t13 = t9*t9; 2238 t21 = t*t; 2239 t23 = s+r-1.0; 2240 t25 = t23*t23; 2241 f = (t3*t2+56.0/23.0*t3*t1*t6+588.0/253.0*t3*s*t9+280.0/253.0*t3*t11+ 2242 70.0/253.0*t2*t13+168.0/4807.0*t1*t13*t6+28.0/14421.0*s*t13*t9+8.0/245157.0*t13 2243 *t11)*(t21*t+12.0/7.0*t21*t23+6.0/7.0*t*t25+4.0/35.0*t25*t23); 2244 /* p=14 0 8 2 */ 2245 } else if( i==0 && j==8 && k==2 ) { 2246 t1 = s*s; 2247 t2 = t1*t1; 2248 t3 = t2*t2; 2249 t4 = t1*s; 2250 t6 = r-1.0; 2251 t9 = t6*t6; 2252 t12 = t9*t6; 2253 t14 = t9*t9; 2254 t22 = t14*t14; 2255 t24 = t*t; 2256 t25 = s+r-1.0; 2257 t27 = t25*t25; 2258 f = (t3+72.0/23.0*t2*t4*t6+1008.0/253.0*t2*t1*t9+672.0/253.0*t2*s*t12+ 2259 252.0/253.0*t2*t14+1008.0/4807.0*t4*t14*t6+112.0/4807.0*t1*t14*t9+96.0/81719.0* 2260 s*t14*t12+3.0/163438.0*t22)*(t24+6.0/5.0*t*t25+3.0/10.0*t27); 2261 /* p=14 0 9 1 */ 2262 } else if( i==0 && j==9 && k==1 ) { 2263 t1 = s*s; 2264 t2 = t1*t1; 2265 t3 = t2*t2; 2266 t5 = r-1.0; 2267 t7 = t1*s; 2268 t9 = t5*t5; 2269 t12 = t9*t5; 2270 t15 = t9*t9; 2271 t23 = t15*t15; 2272 f = (t3*s+90.0/23.0*t3*t5+1620.0/253.0*t2*t7*t9+1440.0/253.0*t2*t1*t12+ 2273 756.0/253.0*t2*s*t15+4536.0/4807.0*t2*t15*t5+840.0/4807.0*t7*t15*t9+1440.0/ 2274 81719.0*t1*t15*t12+135.0/163438.0*s*t23+t23*t5/81719)*(t+2.0/3.0*s+2.0/3.0*r 2275 -2.0/3.0); 2276 /* p=14 0 10 0 */ 2277 } else if( i==0 && j==10 && k==0 ) { 2278 t1 = s*s; 2279 t2 = t1*t1; 2280 t3 = t2*t2; 2281 t6 = r-1.0; 2282 t8 = t6*t6; 2283 t10 = t1*s; 2284 t12 = t8*t6; 2285 t15 = t8*t8; 2286 t24 = t15*t15; 2287 f = t3*t1+110.0/23.0*t3*s*t6+225.0/23.0*t3*t8+1800.0/161.0*t2*t10*t12+ 2288 180.0/23.0*t2*t1*t15+1512.0/437.0*t2*s*t15*t6+420.0/437.0*t2*t15*t8+1200.0/ 2289 7429.0*t10*t15*t12+225.0/14858.0*t1*t24+5.0/7429.0*s*t24*t6+t24*t8/104006; 2290 /* p=14 1 0 9 */ 2291 } else if( i==1 && j==0 && k==9 ) { 2292 t2 = t*t; 2293 t3 = t2*t2; 2294 t4 = t3*t3; 2295 t6 = s+r-1.0; 2296 t8 = t2*t; 2297 t10 = t6*t6; 2298 t13 = t10*t6; 2299 t16 = t10*t10; 2300 t24 = t16*t16; 2301 f = (r-2.0/25.0)*(t4*t+90.0/19.0*t4*t6+180.0/19.0*t3*t8*t10+3360.0/ 2302 323.0*t3*t2*t13+2205.0/323.0*t3*t*t16+882.0/323.0*t3*t16*t6+210.0/323.0*t8*t16* 2303 t10+360.0/4199.0*t2*t16*t13+45.0/8398.0*t*t24+5.0/46189.0*t24*t6); 2304 /* p=14 1 1 8 */ 2305 } else if( i==1 && j==1 && k==8 ) { 2306 t4 = t*t; 2307 t5 = t4*t4; 2308 t6 = t5*t5; 2309 t7 = t4*t; 2310 t9 = s+r-1.0; 2311 t12 = t9*t9; 2312 t15 = t12*t9; 2313 t17 = t12*t12; 2314 t25 = t17*t17; 2315 f = (s+2.0/21.0*r-2.0/21.0)*(r-2.0/25.0)*(t6+72.0/17.0*t5*t7*t9+126.0/ 2316 17.0*t5*t4*t12+588.0/85.0*t5*t*t15+63.0/17.0*t5*t17+252.0/221.0*t7*t17*t9+42.0/ 2317 221.0*t4*t17*t12+36.0/2431.0*t*t17*t15+9.0/24310.0*t25); 2318 /* p=14 1 2 7 */ 2319 } else if( i==1 && j==2 && k==7 ) { 2320 t2 = s*s; 2321 t3 = r-1.0; 2322 t5 = t3*t3; 2323 t8 = t*t; 2324 t9 = t8*t; 2325 t10 = t8*t8; 2326 t13 = s+r-1.0; 2327 t16 = t13*t13; 2328 t18 = t16*t13; 2329 t20 = t16*t16; 2330 f = (r-2.0/25.0)*(t2+2.0/7.0*s*t3+t5/70)*(t10*t9+56.0/15.0*t10*t8*t13+ 2331 28.0/5.0*t10*t*t16+56.0/13.0*t10*t18+70.0/39.0*t9*t20+56.0/143.0*t8*t20*t13+ 2332 28.0/715.0*t*t20*t16+8.0/6435.0*t20*t18); 2333 /* p=14 1 3 6 */ 2334 } else if( i==1 && j==3 && k==6 ) { 2335 t2 = s*s; 2336 t4 = r-1.0; 2337 t6 = t4*t4; 2338 t11 = t*t; 2339 t12 = t11*t11; 2340 t15 = s+r-1.0; 2341 t17 = t15*t15; 2342 t22 = t17*t17; 2343 f = (r-2.0/25.0)*(t2*s+4.0/7.0*t2*t4+3.0/35.0*s*t6+2.0/665.0*t6*t4)*( 2344 t11*t12+42.0/13.0*t12*t*t15+105.0/26.0*t17*t12+350.0/143.0*t11*t*t17*t15+105.0/ 2345 143.0*t11*t22+14.0/143.0*t*t22*t15+7.0/1716.0*t22*t17); 2346 /* p=14 1 4 5 */ 2347 } else if( i==1 && j==4 && k==5 ) { 2348 t2 = s*s; 2349 t3 = t2*t2; 2350 t5 = r-1.0; 2351 t7 = t5*t5; 2352 t11 = t7*t7; 2353 t14 = t*t; 2354 t15 = t14*t14; 2355 t17 = s+r-1.0; 2356 t20 = t17*t17; 2357 t24 = t20*t20; 2358 f = (r-2.0/25.0)*(t3+20.0/21.0*t5*t2*s+2.0/7.0*t2*t7+4.0/133.0*s*t5*t7+ 2359 t11/1197)*(t15*t+30.0/11.0*t17*t15+30.0/11.0*t14*t*t20+40.0/33.0*t14*t20*t17+ 2360 5.0/22.0*t*t24+t24*t17/77); 2361 /* p=14 1 5 4 */ 2362 } else if( i==1 && j==5 && k==4 ) { 2363 t2 = s*s; 2364 t3 = t2*t2; 2365 t5 = r-1.0; 2366 t8 = t5*t5; 2367 t12 = t8*t8; 2368 t17 = t*t; 2369 t18 = t17*t17; 2370 t20 = s+r-1.0; 2371 t22 = t20*t20; 2372 t26 = t22*t22; 2373 f = (r-2.0/25.0)*(s*t3+10.0/7.0*t3*t5+5.0/7.0*t2*s*t8+20.0/133.0*t2*t8* 2374 t5+5.0/399.0*s*t12+2.0/6783.0*t12*t5)*(t18+20.0/9.0*t17*t*t20+5.0/3.0*t22*t17+ 2375 10.0/21.0*t*t22*t20+5.0/126.0*t26); 2376 /* p=14 1 6 3 */ 2377 } else if( i==1 && j==6 && k==3 ) { 2378 t2 = s*s; 2379 t3 = t2*t2; 2380 t6 = r-1.0; 2381 t8 = t6*t6; 2382 t13 = t8*t8; 2383 t20 = t*t; 2384 t22 = s+r-1.0; 2385 t24 = t22*t22; 2386 f = (r-2.0/25.0)*(t3*t2+2.0*t3*s*t6+3.0/2.0*t3*t8+10.0/19.0*t2*s*t8*t6+ 2387 5.0/57.0*t2*t13+2.0/323.0*s*t13*t6+t13*t8/7752)*(t20*t+12.0/7.0*t22*t20+6.0/7.0 2388 *t*t24+4.0/35.0*t24*t22); 2389 /* p=14 1 7 2 */ 2390 } else if( i==1 && j==7 && k==2 ) { 2391 t2 = s*s; 2392 t3 = t2*s; 2393 t4 = t2*t2; 2394 t7 = r-1.0; 2395 t10 = t7*t7; 2396 t12 = t10*t7; 2397 t14 = t10*t10; 2398 t23 = t*t; 2399 t24 = s+r-1.0; 2400 t26 = t24*t24; 2401 f = (r-2.0/25.0)*(t4*t3+8.0/3.0*t2*t4*t7+14.0/5.0*t10*t4*s+28.0/19.0*t4 2402 *t12+70.0/171.0*t3*t14+56.0/969.0*t2*t14*t7+7.0/1938.0*s*t14*t10+t14*t12/14535) 2403 *(t23+6.0/5.0*t*t24+3.0/10.0*t26); 2404 /* p=14 1 8 1 */ 2405 } else if( i==1 && j==8 && k==1 ) { 2406 t2 = s*s; 2407 t3 = t2*t2; 2408 t4 = t3*t3; 2409 t5 = t2*s; 2410 t7 = r-1.0; 2411 t10 = t7*t7; 2412 t13 = t10*t7; 2413 t15 = t10*t10; 2414 t23 = t15*t15; 2415 f = (r-2.0/25.0)*(t4+24.0/7.0*t3*t5*t7+24.0/5.0*t3*t2*t10+336.0/95.0*t3 2416 *s*t13+28.0/19.0*t3*t15+112.0/323.0*t5*t15*t7+14.0/323.0*t2*t15*t10+4.0/1615.0* 2417 s*t15*t13+t23/22610)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 2418 /* p=14 1 9 0 */ 2419 } else if( i==1 && j==9 && k==0 ) { 2420 t2 = s*s; 2421 t3 = t2*t2; 2422 t4 = t3*t3; 2423 t6 = r-1.0; 2424 t8 = t2*s; 2425 t10 = t6*t6; 2426 t13 = t10*t6; 2427 t16 = t10*t10; 2428 t24 = t16*t16; 2429 f = (r-2.0/25.0)*(t4*s+30.0/7.0*t4*t6+54.0/7.0*t3*t8*t10+144.0/19.0*t3* 2430 t2*t13+84.0/19.0*t3*s*t16+504.0/323.0*t3*t16*t6+105.0/323.0*t8*t16*t10+12.0/ 2431 323.0*t2*t16*t13+9.0/4522.0*s*t24+t24*t6/29393); 2432 /* p=14 2 0 8 */ 2433 } else if( i==2 && j==0 && k==8 ) { 2434 t1 = r*r; 2435 t3 = t*t; 2436 t4 = t3*t3; 2437 t5 = t4*t4; 2438 t6 = t3*t; 2439 t8 = s+r-1.0; 2440 t11 = t8*t8; 2441 t14 = t11*t8; 2442 t16 = t11*t11; 2443 t24 = t16*t16; 2444 f = (t1-6.0/25.0*r+1.0/100.0)*(t5+72.0/17.0*t4*t6*t8+126.0/17.0*t4*t3* 2445 t11+588.0/85.0*t4*t*t14+63.0/17.0*t4*t16+252.0/221.0*t6*t16*t8+42.0/221.0*t3* 2446 t16*t11+36.0/2431.0*t*t16*t14+9.0/24310.0*t24); 2447 /* p=14 2 1 7 */ 2448 } else if( i==2 && j==1 && k==7 ) { 2449 t1 = r*r; 2450 t5 = t*t; 2451 t6 = t5*t; 2452 t7 = t5*t5; 2453 t10 = s+r-1.0; 2454 t13 = t10*t10; 2455 t15 = t13*t10; 2456 t17 = t13*t13; 2457 f = (s+2.0/19.0*r-2.0/19.0)*(t1-6.0/25.0*r+1.0/100.0)*(t6*t7+56.0/15.0* 2458 t7*t5*t10+28.0/5.0*t7*t*t13+56.0/13.0*t15*t7+70.0/39.0*t6*t17+56.0/143.0*t5*t17 2459 *t10+28.0/715.0*t*t17*t13+8.0/6435.0*t17*t15); 2460 /* p=14 2 2 6 */ 2461 } else if( i==2 && j==2 && k==6 ) { 2462 t1 = r*r; 2463 t3 = s*s; 2464 t4 = r-1.0; 2465 t6 = t4*t4; 2466 t9 = t*t; 2467 t10 = t9*t9; 2468 t13 = s+r-1.0; 2469 t15 = t13*t13; 2470 t20 = t15*t15; 2471 f = (t1-6.0/25.0*r+1.0/100.0)*(t3+6.0/19.0*t4*s+t6/57)*(t9*t10+42.0/ 2472 13.0*t10*t*t13+105.0/26.0*t15*t10+350.0/143.0*t9*t*t15*t13+105.0/143.0*t9*t20+ 2473 14.0/143.0*t*t20*t13+7.0/1716.0*t20*t15); 2474 /* p=14 2 3 5 */ 2475 } else if( i==2 && j==3 && k==5 ) { 2476 t1 = r*r; 2477 t3 = s*s; 2478 t5 = r-1.0; 2479 t7 = t5*t5; 2480 t12 = t*t; 2481 t13 = t12*t12; 2482 t15 = s+r-1.0; 2483 t18 = t15*t15; 2484 t22 = t18*t18; 2485 f = (t1-6.0/25.0*r+1.0/100.0)*(t3*s+12.0/19.0*t3*t5+2.0/19.0*s*t7+4.0/ 2486 969.0*t7*t5)*(t13*t+30.0/11.0*t15*t13+30.0/11.0*t12*t*t18+40.0/33.0*t12*t18*t15 2487 +5.0/22.0*t*t22+t22*t15/77); 2488 /* p=14 2 4 4 */ 2489 } else if( i==2 && j==4 && k==4 ) { 2490 t1 = r*r; 2491 t3 = s*s; 2492 t4 = t3*t3; 2493 t6 = r-1.0; 2494 t8 = t6*t6; 2495 t12 = t8*t8; 2496 t15 = t*t; 2497 t16 = t15*t15; 2498 t18 = s+r-1.0; 2499 t20 = t18*t18; 2500 t24 = t20*t20; 2501 f = (t1-6.0/25.0*r+1.0/100.0)*(t4+20.0/19.0*t3*s*t6+20.0/57.0*t3*t8+ 2502 40.0/969.0*s*t8*t6+5.0/3876.0*t12)*(t16+20.0/9.0*t15*t*t18+5.0/3.0*t15*t20+10.0 2503 /21.0*t*t20*t18+5.0/126.0*t24); 2504 /* p=14 2 5 3 */ 2505 } else if( i==2 && j==5 && k==3 ) { 2506 t1 = r*r; 2507 t3 = s*s; 2508 t4 = t3*t3; 2509 t6 = r-1.0; 2510 t9 = t6*t6; 2511 t13 = t9*t9; 2512 t18 = t*t; 2513 t20 = s+r-1.0; 2514 t22 = t20*t20; 2515 f = (t1-6.0/25.0*r+1.0/100.0)*(t4*s+30.0/19.0*t4*t6+50.0/57.0*t3*s*t9+ 2516 200.0/969.0*t3*t9*t6+25.0/1292.0*s*t13+t13*t6/1938)*(t18*t+12.0/7.0*t20*t18+6.0 2517 /7.0*t*t22+4.0/35.0*t22*t20); 2518 /* p=14 2 6 2 */ 2519 } else if( i==2 && j==6 && k==2 ) { 2520 t1 = r*r; 2521 t3 = s*s; 2522 t4 = t3*t3; 2523 t7 = r-1.0; 2524 t9 = t7*t7; 2525 t14 = t9*t9; 2526 t21 = t*t; 2527 t22 = s+r-1.0; 2528 t24 = t22*t22; 2529 f = (t1-6.0/25.0*r+1.0/100.0)*(t4*t3+42.0/19.0*t4*s*t7+35.0/19.0*t4*t9+ 2530 700.0/969.0*t3*s*t9*t7+175.0/1292.0*t3*t14+7.0/646.0*s*t14*t7+t14*t9/3876)*(t21 2531 +6.0/5.0*t*t22+3.0/10.0*t24); 2532 /* p=14 2 7 1 */ 2533 } else if( i==2 && j==7 && k==1 ) { 2534 t1 = r*r; 2535 t3 = s*s; 2536 t4 = t3*s; 2537 t5 = t3*t3; 2538 t8 = r-1.0; 2539 t11 = t8*t8; 2540 t13 = t11*t8; 2541 t15 = t11*t11; 2542 f = (t1-6.0/25.0*r+1.0/100.0)*(t5*t4+56.0/19.0*t5*t3*t8+196.0/57.0*t5*s 2543 *t11+1960.0/969.0*t5*t13+1225.0/1938.0*t4*t15+98.0/969.0*t3*t15*t8+7.0/969.0*s* 2544 t15*t11+2.0/12597.0*t15*t13)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 2545 /* p=14 2 8 0 */ 2546 } else if( i==2 && j==8 && k==0 ) { 2547 t1 = r*r; 2548 t3 = s*s; 2549 t4 = t3*t3; 2550 t5 = t4*t4; 2551 t6 = t3*s; 2552 t8 = r-1.0; 2553 t11 = t8*t8; 2554 t14 = t11*t8; 2555 t16 = t11*t11; 2556 t24 = t16*t16; 2557 f = (t1-6.0/25.0*r+1.0/100.0)*(t5+72.0/19.0*t4*t6*t8+112.0/19.0*t4*t3* 2558 t11+1568.0/323.0*t4*s*t14+735.0/323.0*t4*t16+196.0/323.0*t6*t16*t8+28.0/323.0* 2559 t3*t16*t11+24.0/4199.0*s*t16*t14+t24/8398); 2560 /* p=14 3 0 7 */ 2561 } else if( i==3 && j==0 && k==7 ) { 2562 t1 = r*r; 2563 t4 = t*t; 2564 t5 = t4*t; 2565 t6 = t4*t4; 2566 t9 = s+r-1.0; 2567 t12 = t9*t9; 2568 t14 = t12*t9; 2569 t16 = t12*t12; 2570 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t5*t6+56.0/15.0*t4*t6*t9+ 2571 28.0/5.0*t6*t*t12+56.0/13.0*t6*t14+70.0/39.0*t5*t16+56.0/143.0*t4*t16*t9+28.0/ 2572 715.0*t*t16*t12+8.0/6435.0*t16*t14); 2573 /* p=14 3 1 6 */ 2574 } else if( i==3 && j==1 && k==6 ) { 2575 t1 = r*r; 2576 t6 = t*t; 2577 t7 = t6*t6; 2578 t10 = s+r-1.0; 2579 t12 = t10*t10; 2580 t17 = t12*t12; 2581 f = (s+2.0/17.0*r-2.0/17.0)*(t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*( 2582 t6*t7+42.0/13.0*t7*t*t10+105.0/26.0*t7*t12+350.0/143.0*t6*t*t12*t10+105.0/143.0 2583 *t6*t17+14.0/143.0*t*t17*t10+7.0/1716.0*t17*t12); 2584 /* p=14 3 2 5 */ 2585 } else if( i==3 && j==2 && k==5 ) { 2586 t1 = r*r; 2587 t4 = s*s; 2588 t5 = r-1.0; 2589 t7 = t5*t5; 2590 t10 = t*t; 2591 t11 = t10*t10; 2592 t13 = s+r-1.0; 2593 t16 = t13*t13; 2594 t20 = t16*t16; 2595 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t4+6.0/17.0*s*t5+3.0/ 2596 136.0*t7)*(t11*t+30.0/11.0*t11*t13+30.0/11.0*t10*t*t16+40.0/33.0*t10*t16*t13+ 2597 5.0/22.0*t*t20+t20*t13/77); 2598 /* p=14 3 3 4 */ 2599 } else if( i==3 && j==3 && k==4 ) { 2600 t1 = r*r; 2601 t4 = s*s; 2602 t6 = r-1.0; 2603 t8 = t6*t6; 2604 t13 = t*t; 2605 t14 = t13*t13; 2606 t16 = s+r-1.0; 2607 t18 = t16*t16; 2608 t22 = t18*t18; 2609 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t4*s+12.0/17.0*t4*t6+9.0/ 2610 68.0*s*t8+t8*t6/170)*(t14+20.0/9.0*t13*t*t16+5.0/3.0*t13*t18+10.0/21.0*t*t18* 2611 t16+5.0/126.0*t22); 2612 /* p=14 3 4 3 */ 2613 } else if( i==3 && j==4 && k==3 ) { 2614 t1 = r*r; 2615 t4 = s*s; 2616 t5 = t4*t4; 2617 t7 = r-1.0; 2618 t9 = t7*t7; 2619 t13 = t9*t9; 2620 t16 = t*t; 2621 t18 = s+r-1.0; 2622 t20 = t18*t18; 2623 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t5+20.0/17.0*t4*s*t7+15.0 2624 /34.0*t4*t9+s*t9*t7/17+t13/476)*(t16*t+12.0/7.0*t18*t16+6.0/7.0*t*t20+4.0/35.0* 2625 t20*t18); 2626 /* p=14 3 5 2 */ 2627 } else if( i==3 && j==5 && k==2 ) { 2628 t1 = r*r; 2629 t4 = s*s; 2630 t5 = t4*t4; 2631 t7 = r-1.0; 2632 t10 = t7*t7; 2633 t14 = t10*t10; 2634 t19 = t*t; 2635 t20 = s+r-1.0; 2636 t22 = t20*t20; 2637 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(s*t5+30.0/17.0*t5*t7+75.0 2638 /68.0*t4*s*t10+5.0/17.0*t4*t10*t7+15.0/476.0*s*t14+3.0/3094.0*t14*t7)*(t19+6.0/ 2639 5.0*t*t20+3.0/10.0*t22); 2640 /* p=14 3 6 1 */ 2641 } else if( i==3 && j==6 && k==1 ) { 2642 t1 = r*r; 2643 t4 = s*s; 2644 t5 = t4*t4; 2645 t8 = r-1.0; 2646 t10 = t8*t8; 2647 t15 = t10*t10; 2648 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t5*t4+42.0/17.0*t5*s*t8+ 2649 315.0/136.0*t5*t10+35.0/34.0*t4*s*t10*t8+15.0/68.0*t4*t15+9.0/442.0*s*t15*t8+ 2650 t15*t10/1768)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 2651 /* p=14 3 7 0 */ 2652 } else if( i==3 && j==7 && k==0 ) { 2653 t1 = r*r; 2654 t4 = s*s; 2655 t5 = t4*s; 2656 t6 = t4*t4; 2657 t9 = r-1.0; 2658 t12 = t9*t9; 2659 t14 = t12*t9; 2660 t16 = t12*t12; 2661 f = (t1*r-12.0/25.0*t1+3.0/50.0*r-1.0/575.0)*(t6*t5+56.0/17.0*t4*t6*t9+ 2662 147.0/34.0*t6*s*t12+49.0/17.0*t6*t14+35.0/34.0*t5*t16+42.0/221.0*t4*t16*t9+7.0/ 2663 442.0*s*t16*t12+t16*t14/2431); 2664 /* p=14 4 0 6 */ 2665 } else if( i==4 && j==0 && k==6 ) { 2666 t1 = r*r; 2667 t2 = t1*t1; 2668 t5 = t*t; 2669 t6 = t5*t5; 2670 t9 = s+r-1.0; 2671 t11 = t9*t9; 2672 t16 = t11*t11; 2673 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(t6*t5+42.0/13.0*t6*t 2674 *t9+105.0/26.0*t6*t11+350.0/143.0*t5*t*t11*t9+105.0/143.0*t5*t16+14.0/143.0*t* 2675 t16*t9+7.0/1716.0*t16*t11); 2676 /* p=14 4 1 5 */ 2677 } else if( i==4 && j==1 && k==5 ) { 2678 t1 = r*r; 2679 t2 = t1*t1; 2680 t7 = t*t; 2681 t8 = t7*t7; 2682 t10 = s+r-1.0; 2683 t13 = t10*t10; 2684 t17 = t13*t13; 2685 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(s+2.0/15.0*r-2.0/ 2686 15.0)*(t8*t+30.0/11.0*t8*t10+30.0/11.0*t7*t*t13+40.0/33.0*t7*t13*t10+5.0/22.0*t 2687 *t17+t17*t10/77); 2688 /* p=14 4 2 4 */ 2689 } else if( i==4 && j==2 && k==4 ) { 2690 t1 = r*r; 2691 t2 = t1*t1; 2692 t5 = s*s; 2693 t6 = r-1.0; 2694 t8 = t6*t6; 2695 t11 = t*t; 2696 t12 = t11*t11; 2697 t14 = s+r-1.0; 2698 t16 = t14*t14; 2699 t20 = t16*t16; 2700 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(t5+2.0/5.0*s*t6+t8 2701 /35)*(t12+20.0/9.0*t11*t*t14+5.0/3.0*t16*t11+10.0/21.0*t*t16*t14+5.0/126.0*t20) 2702 ; 2703 /* p=14 4 3 3 */ 2704 } else if( i==4 && j==3 && k==3 ) { 2705 t1 = r*r; 2706 t2 = t1*t1; 2707 t5 = s*s; 2708 t7 = r-1.0; 2709 t9 = t7*t7; 2710 t14 = t*t; 2711 t16 = s+r-1.0; 2712 t18 = t16*t16; 2713 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(t5*s+4.0/5.0*t5*t7+ 2714 6.0/35.0*s*t9+4.0/455.0*t9*t7)*(t14*t+12.0/7.0*t16*t14+6.0/7.0*t*t18+4.0/35.0* 2715 t18*t16); 2716 /* p=14 4 4 2 */ 2717 } else if( i==4 && j==4 && k==2 ) { 2718 t1 = r*r; 2719 t2 = t1*t1; 2720 t5 = s*s; 2721 t6 = t5*t5; 2722 t8 = r-1.0; 2723 t10 = t8*t8; 2724 t14 = t10*t10; 2725 t17 = t*t; 2726 t18 = s+r-1.0; 2727 t20 = t18*t18; 2728 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(t6+4.0/3.0*t5*s*t8+ 2729 4.0/7.0*t5*t10+8.0/91.0*s*t8*t10+t14/273)*(t17+6.0/5.0*t*t18+3.0/10.0*t20); 2730 /* p=14 4 5 1 */ 2731 } else if( i==4 && j==5 && k==1 ) { 2732 t1 = r*r; 2733 t2 = t1*t1; 2734 t5 = s*s; 2735 t6 = t5*t5; 2736 t8 = r-1.0; 2737 t11 = t8*t8; 2738 t15 = t11*t11; 2739 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(s*t6+2.0*t6*t8+10.0/ 2740 7.0*t5*s*t11+40.0/91.0*t5*t11*t8+5.0/91.0*s*t15+2.0/1001.0*t15*t8)*(t+2.0/3.0*s 2741 +2.0/3.0*r-2.0/3.0); 2742 /* p=14 4 6 0 */ 2743 } else if( i==4 && j==6 && k==0 ) { 2744 t1 = r*r; 2745 t2 = t1*t1; 2746 t5 = s*s; 2747 t6 = t5*t5; 2748 t9 = r-1.0; 2749 t11 = t9*t9; 2750 t16 = t11*t11; 2751 f = (t2-4.0/5.0*t1*r+t1/5-2.0/115.0*r+1.0/2530.0)*(t6*t5+14.0/5.0*t6*s* 2752 t9+3.0*t6*t11+20.0/13.0*t5*s*t11*t9+5.0/13.0*t5*t16+6.0/143.0*s*t16*t9+t16*t11 2753 /715); 2754 /* p=14 5 0 5 */ 2755 } else if( i==5 && j==0 && k==5 ) { 2756 t1 = r*r; 2757 t2 = t1*t1; 2758 t6 = t*t; 2759 t7 = t6*t6; 2760 t9 = s+r-1.0; 2761 t12 = t9*t9; 2762 t16 = t12*t12; 2763 f = (t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0)*(t7*t+ 2764 30.0/11.0*t7*t9+30.0/11.0*t6*t*t12+40.0/33.0*t6*t12*t9+5.0/22.0*t*t16+t16*t9/77 2765 ); 2766 /* p=14 5 1 4 */ 2767 } else if( i==5 && j==1 && k==4 ) { 2768 t1 = r*r; 2769 t2 = t1*t1; 2770 t8 = t*t; 2771 t9 = t8*t8; 2772 t11 = s+r-1.0; 2773 t13 = t11*t11; 2774 t17 = t13*t13; 2775 f = (s+2.0/13.0*r-2.0/13.0)*(t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/ 2776 506.0*r-1.0/8855.0)*(t9+20.0/9.0*t8*t*t11+5.0/3.0*t8*t13+10.0/21.0*t*t13*t11+ 2777 5.0/126.0*t17); 2778 /* p=14 5 2 3 */ 2779 } else if( i==5 && j==2 && k==3 ) { 2780 t1 = r*r; 2781 t2 = t1*t1; 2782 t6 = s*s; 2783 t7 = r-1.0; 2784 t9 = t7*t7; 2785 t12 = t*t; 2786 t14 = s+r-1.0; 2787 t16 = t14*t14; 2788 f = (t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0)*(t6+6.0 2789 /13.0*s*t7+t9/26)*(t12*t+12.0/7.0*t12*t14+6.0/7.0*t*t16+4.0/35.0*t16*t14); 2790 /* p=14 5 3 2 */ 2791 } else if( i==5 && j==3 && k==2 ) { 2792 t1 = r*r; 2793 t2 = t1*t1; 2794 t6 = s*s; 2795 t8 = r-1.0; 2796 t10 = t8*t8; 2797 t15 = t*t; 2798 t16 = s+r-1.0; 2799 t18 = t16*t16; 2800 f = (t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0)*(t6*s+ 2801 12.0/13.0*t6*t8+3.0/13.0*s*t10+2.0/143.0*t10*t8)*(t15+6.0/5.0*t*t16+3.0/10.0* 2802 t18); 2803 /* p=14 5 4 1 */ 2804 } else if( i==5 && j==4 && k==1 ) { 2805 t1 = r*r; 2806 t2 = t1*t1; 2807 t6 = s*s; 2808 t7 = t6*t6; 2809 t9 = r-1.0; 2810 t11 = t9*t9; 2811 t15 = t11*t11; 2812 f = (t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0)*(t7+ 2813 20.0/13.0*t6*s*t9+10.0/13.0*t6*t11+20.0/143.0*s*t11*t9+t15/143)*(t+2.0/3.0*s+ 2814 2.0/3.0*r-2.0/3.0); 2815 /* p=14 5 5 0 */ 2816 } else if( i==5 && j==5 && k==0 ) { 2817 t1 = r*r; 2818 t2 = t1*t1; 2819 t6 = s*s; 2820 t7 = t6*t6; 2821 t9 = r-1.0; 2822 t12 = t9*t9; 2823 t16 = t12*t12; 2824 f = (t2*r-6.0/5.0*t2+t1*r/2-2.0/23.0*t1+3.0/506.0*r-1.0/8855.0)*(s*t7+ 2825 30.0/13.0*t7*t9+25.0/13.0*t6*s*t12+100.0/143.0*t6*t12*t9+15.0/143.0*s*t16+2.0/ 2826 429.0*t16*t9); 2827 /* p=14 6 0 4 */ 2828 } else if( i==6 && j==0 && k==4 ) { 2829 t1 = r*r; 2830 t2 = t1*t1; 2831 t7 = t*t; 2832 t8 = t7*t7; 2833 t10 = s+r-1.0; 2834 t12 = t10*t10; 2835 t16 = t12*t12; 2836 f = (t2*t1-42.0/25.0*t2*r+21.0/20.0*t2-7.0/23.0*t1*r+21.0/506.0*t1-3.0/ 2837 1265.0*r+1.0/25300.0)*(t8+20.0/9.0*t7*t*t10+5.0/3.0*t7*t12+10.0/21.0*t*t12*t10+ 2838 5.0/126.0*t16); 2839 /* p=14 6 1 3 */ 2840 } else if( i==6 && j==1 && k==3 ) { 2841 t1 = r*r; 2842 t2 = t1*t1; 2843 t9 = t*t; 2844 t11 = s+r-1.0; 2845 t13 = t11*t11; 2846 f = (t2*t1-42.0/25.0*t2*r+21.0/20.0*t2-7.0/23.0*t1*r+21.0/506.0*t1-3.0/ 2847 1265.0*r+1.0/25300.0)*(s+2.0/11.0*r-2.0/11.0)*(t9*t+12.0/7.0*t11*t9+6.0/7.0*t* 2848 t13+4.0/35.0*t13*t11); 2849 /* p=14 6 2 2 */ 2850 } else if( i==6 && j==2 && k==2 ) { 2851 t1 = r*r; 2852 t2 = t1*t1; 2853 t7 = s*s; 2854 t8 = r-1.0; 2855 t10 = t8*t8; 2856 t13 = t*t; 2857 t14 = s+r-1.0; 2858 t16 = t14*t14; 2859 f = (t2*t1-42.0/25.0*t2*r+21.0/20.0*t2-7.0/23.0*t1*r+21.0/506.0*t1-3.0/ 2860 1265.0*r+1.0/25300.0)*(t7+6.0/11.0*s*t8+3.0/55.0*t10)*(t13+6.0/5.0*t*t14+3.0/ 2861 10.0*t16); 2862 /* p=14 6 3 1 */ 2863 } else if( i==6 && j==3 && k==1 ) { 2864 t1 = r*r; 2865 t2 = t1*t1; 2866 t7 = s*s; 2867 t9 = r-1.0; 2868 t11 = t9*t9; 2869 f = (t2*t1-42.0/25.0*t2*r+21.0/20.0*t2-7.0/23.0*t1*r+21.0/506.0*t1-3.0/ 2870 1265.0*r+1.0/25300.0)*(t7*s+12.0/11.0*t7*t9+18.0/55.0*s*t11+4.0/165.0*t11*t9)*( 2871 t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 2872 /* p=14 6 4 0 */ 2873 } else if( i==6 && j==4 && k==0 ) { 2874 t1 = r*r; 2875 t2 = t1*t1; 2876 t7 = s*s; 2877 t8 = t7*t7; 2878 t10 = r-1.0; 2879 t12 = t10*t10; 2880 t16 = t12*t12; 2881 f = (t2*t1-42.0/25.0*t2*r+21.0/20.0*t2-7.0/23.0*t1*r+21.0/506.0*t1-3.0/ 2882 1265.0*r+1.0/25300.0)*(t8+20.0/11.0*t7*s*t10+12.0/11.0*t7*t12+8.0/33.0*s*t12* 2883 t10+t16/66); 2884 /* p=14 7 0 3 */ 2885 } else if( i==7 && j==0 && k==3 ) { 2886 t1 = r*r; 2887 t2 = t1*r; 2888 t3 = t1*t1; 2889 t8 = t*t; 2890 t10 = s+r-1.0; 2891 t12 = t10*t10; 2892 f = (t3*t2-56.0/25.0*t3*t1+49.0/25.0*t3*r-98.0/115.0*t3+49.0/253.0*t2 2893 -28.0/1265.0*t1+7.0/6325.0*r-2.0/120175.0)*(t8*t+12.0/7.0*t8*t10+6.0/7.0*t*t12+ 2894 4.0/35.0*t12*t10); 2895 /* p=14 7 1 2 */ 2896 } else if( i==7 && j==1 && k==2 ) { 2897 t1 = r*r; 2898 t2 = t1*r; 2899 t3 = t1*t1; 2900 t10 = t*t; 2901 t11 = s+r-1.0; 2902 t13 = t11*t11; 2903 f = (t3*t2-56.0/25.0*t3*t1+49.0/25.0*t3*r-98.0/115.0*t3+49.0/253.0*t2 2904 -28.0/1265.0*t1+7.0/6325.0*r-2.0/120175.0)*(s+2.0/9.0*r-2.0/9.0)*(t10+6.0/5.0*t 2905 *t11+3.0/10.0*t13); 2906 /* p=14 7 2 1 */ 2907 } else if( i==7 && j==2 && k==1 ) { 2908 t1 = r*r; 2909 t2 = t1*r; 2910 t3 = t1*t1; 2911 t8 = s*s; 2912 t9 = r-1.0; 2913 t11 = t9*t9; 2914 f = (t3*t2-56.0/25.0*t3*t1+49.0/25.0*t3*r-98.0/115.0*t3+49.0/253.0*t2 2915 -28.0/1265.0*t1+7.0/6325.0*r-2.0/120175.0)*(t8+2.0/3.0*s*t9+t11/12)*(t+2.0/3.0* 2916 s+2.0/3.0*r-2.0/3.0); 2917 /* p=14 7 3 0 */ 2918 } else if( i==7 && j==3 && k==0 ) { 2919 t1 = r*r; 2920 t2 = t1*r; 2921 t3 = t1*t1; 2922 t8 = s*s; 2923 t10 = r-1.0; 2924 t12 = t10*t10; 2925 f = (t3*t2-56.0/25.0*t3*t1+49.0/25.0*t3*r-98.0/115.0*t3+49.0/253.0*t2 2926 -28.0/1265.0*t1+7.0/6325.0*r-2.0/120175.0)*(s*t8+4.0/3.0*t8*t10+s*t12/2+t12*t10 2927 /21); 2928 /* p=14 8 0 2 */ 2929 } else if( i==8 && j==0 && k==2 ) { 2930 t1 = r*r; 2931 t2 = t1*t1; 2932 t3 = t2*t2; 2933 t4 = t1*r; 2934 t9 = t*t; 2935 t10 = s+r-1.0; 2936 t12 = t10*t10; 2937 f = (t3-72.0/25.0*t2*t4+84.0/25.0*t2*t1-1176.0/575.0*t2*r+882.0/1265.0* 2938 t2-168.0/1265.0*t4+84.0/6325.0*t1-72.0/120175.0*r+1.0/120175.0)*(t9+6.0/5.0*t* 2939 t10+3.0/10.0*t12); 2940 /* p=14 8 1 1 */ 2941 } else if( i==8 && j==1 && k==1 ) { 2942 t1 = r*r; 2943 t2 = t1*t1; 2944 t3 = t2*t2; 2945 t4 = t1*r; 2946 f = (t3-72.0/25.0*t2*t4+84.0/25.0*t2*t1-1176.0/575.0*t2*r+882.0/1265.0* 2947 t2-168.0/1265.0*t4+84.0/6325.0*t1-72.0/120175.0*r+1.0/120175.0)*(s+2.0/7.0*r 2948 -2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 2949 /* p=14 8 2 0 */ 2950 } else if( i==8 && j==2 && k==0 ) { 2951 t1 = r*r; 2952 t2 = t1*t1; 2953 t3 = t2*t2; 2954 t4 = t1*r; 2955 t9 = s*s; 2956 t10 = r-1.0; 2957 t12 = t10*t10; 2958 f = (t3-72.0/25.0*t2*t4+84.0/25.0*t2*t1-1176.0/575.0*t2*r+882.0/1265.0* 2959 t2-168.0/1265.0*t4+84.0/6325.0*t1-72.0/120175.0*r+1.0/120175.0)*(t9+6.0/7.0*s* 2960 t10+t12/7); 2961 /* p=14 9 0 1 */ 2962 } else if( i==9 && j==0 && k==1 ) { 2963 t1 = r*r; 2964 t2 = t1*t1; 2965 t3 = t2*t2; 2966 t5 = t1*r; 2967 f = (t3*r-18.0/5.0*t3+27.0/5.0*t2*t5-504.0/115.0*t2*t1+2646.0/1265.0*t2 2968 *r-756.0/1265.0*t2+126.0/1265.0*t5-216.0/24035.0*t1+9.0/24035.0*r-2.0/408595.0) 2969 *(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 2970 /* p=14 9 1 0 */ 2971 } else if( i==9 && j==1 && k==0 ) { 2972 t1 = r*r; 2973 t2 = t1*t1; 2974 t3 = t2*t2; 2975 t5 = t1*r; 2976 f = (t3*r-18.0/5.0*t3+27.0/5.0*t2*t5-504.0/115.0*t2*t1+2646.0/1265.0*t2 2977 *r-756.0/1265.0*t2+126.0/1265.0*t5-216.0/24035.0*t1+9.0/24035.0*r-2.0/408595.0) 2978 *(s+2.0/5.0*r-2.0/5.0); 2979 /* p=14 10 0 0 */ 2980 } else if( i==10 && j==0 && k==0 ) { 2981 t1 = r*r; 2982 t2 = t1*t1; 2983 t3 = t2*t2; 2984 t6 = t1*r; 2985 f = t3*t1-22.0/5.0*t3*r+33.0/4.0*t3-198.0/23.0*t2*t6+126.0/23.0*t2*t1 2986 -252.0/115.0*t2*r+63.0/115.0*t2-36.0/437.0*t6+3.0/437.0*t1-2.0/7429.0*r+1.0/ 2987 297160.0; 2988 /* p=15 0 0 11 */ 2989 } else if( i==0 && j==0 && k==11 ) { 2990 t1 = t*t; 2991 t2 = t1*t; 2992 t3 = t1*t1; 2993 t4 = t3*t3; 2994 t7 = s+r-1.0; 2995 t10 = t7*t7; 2996 t12 = t10*t7; 2997 t15 = t10*t10; 2998 t25 = t15*t15; 2999 f = t2*t4+132.0/23.0*t4*t1*t7+330.0/23.0*t4*t*t10+3300.0/161.0*t4*t12+ 3000 2970.0/161.0*t3*t2*t15+4752.0/437.0*t3*t1*t15*t7+1848.0/437.0*t3*t*t15*t10+ 3001 7920.0/7429.0*t3*t15*t12+2475.0/14858.0*t2*t25+110.0/7429.0*t1*t25*t7+33.0/ 3002 52003.0*t*t25*t10+6.0/676039.0*t25*t12; 3003 /* p=15 0 1 10 */ 3004 } else if( i==0 && j==1 && k==10 ) { 3005 t2 = t*t; 3006 t3 = t2*t2; 3007 t4 = t3*t3; 3008 t7 = s+r-1.0; 3009 t9 = t7*t7; 3010 t11 = t2*t; 3011 t13 = t9*t7; 3012 t16 = t9*t9; 3013 t25 = t16*t16; 3014 t30 = t2*t4+110.0/21.0*t4*t*t7+165.0/14.0*t4*t9+1980.0/133.0*t3*t11*t13+ 3015 220.0/19.0*t3*t2*t16+1848.0/323.0*t3*t*t16*t7+1155.0/646.0*t3*t16*t9+110.0/ 3016 323.0*t11*t16*t13+165.0/4522.0*t2*t25+55.0/29393.0*t*t25*t7+11.0/352716.0*t25* 3017 t9; 3018 f = (s+2.0/25.0*r-2.0/25.0)*t30; 3019 /* p=15 0 2 9 */ 3020 } else if( i==0 && j==2 && k==9 ) { 3021 t1 = s*s; 3022 t2 = r-1.0; 3023 t4 = t2*t2; 3024 t6 = t*t; 3025 t7 = t6*t6; 3026 t8 = t7*t7; 3027 t10 = s+r-1.0; 3028 t12 = t6*t; 3029 t14 = t10*t10; 3030 t17 = t14*t10; 3031 t20 = t14*t14; 3032 t28 = t20*t20; 3033 f = (t1+6.0/25.0*s*t2+t4/100)*(t8*t+90.0/19.0*t8*t10+180.0/19.0*t7*t12* 3034 t14+3360.0/323.0*t7*t6*t17+2205.0/323.0*t7*t*t20+882.0/323.0*t7*t20*t10+210.0/ 3035 323.0*t12*t20*t14+360.0/4199.0*t6*t20*t17+45.0/8398.0*t*t28+5.0/46189.0*t28*t10 3036 ); 3037 /* p=15 0 3 8 */ 3038 } else if( i==0 && j==3 && k==8 ) { 3039 t1 = s*s; 3040 t3 = r-1.0; 3041 t5 = t3*t3; 3042 t9 = t*t; 3043 t10 = t9*t9; 3044 t11 = t10*t10; 3045 t12 = t9*t; 3046 t14 = s+r-1.0; 3047 t17 = t14*t14; 3048 t20 = t17*t14; 3049 t22 = t17*t17; 3050 t30 = t22*t22; 3051 f = (t1*s+12.0/25.0*t3*t1+3.0/50.0*s*t5+t5*t3/575)*(t11+72.0/17.0*t10* 3052 t12*t14+126.0/17.0*t10*t9*t17+588.0/85.0*t10*t*t20+63.0/17.0*t10*t22+252.0/ 3053 221.0*t12*t22*t14+42.0/221.0*t9*t22*t17+36.0/2431.0*t*t22*t20+9.0/24310.0*t30); 3054 /* p=15 0 4 7 */ 3055 } else if( i==0 && j==4 && k==7 ) { 3056 t1 = s*s; 3057 t2 = t1*t1; 3058 t4 = r-1.0; 3059 t6 = t4*t4; 3060 t10 = t6*t6; 3061 t12 = t*t; 3062 t13 = t12*t; 3063 t14 = t12*t12; 3064 t17 = s+r-1.0; 3065 t20 = t17*t17; 3066 t22 = t20*t17; 3067 t24 = t20*t20; 3068 f = (t2+4.0/5.0*t1*s*t4+t1*t6/5+2.0/115.0*s*t6*t4+t10/2530)*(t13*t14+ 3069 56.0/15.0*t14*t12*t17+28.0/5.0*t14*t*t20+56.0/13.0*t22*t14+70.0/39.0*t13*t24+ 3070 56.0/143.0*t12*t24*t17+28.0/715.0*t*t24*t20+8.0/6435.0*t24*t22); 3071 /* p=15 0 5 6 */ 3072 } else if( i==0 && j==5 && k==6 ) { 3073 t1 = s*s; 3074 t2 = t1*t1; 3075 t4 = r-1.0; 3076 t7 = t4*t4; 3077 t11 = t7*t7; 3078 t15 = t*t; 3079 t16 = t15*t15; 3080 t19 = s+r-1.0; 3081 t21 = t19*t19; 3082 t26 = t21*t21; 3083 f = (s*t2+6.0/5.0*t2*t4+t7*t1*s/2+2.0/23.0*t1*t7*t4+3.0/506.0*s*t11+t11 3084 *t4/8855)*(t16*t15+42.0/13.0*t16*t*t19+105.0/26.0*t16*t21+350.0/143.0*t15*t*t21 3085 *t19+105.0/143.0*t15*t26+14.0/143.0*t*t26*t19+7.0/1716.0*t26*t21); 3086 /* p=15 0 6 5 */ 3087 } else if( i==0 && j==6 && k==5 ) { 3088 t1 = s*s; 3089 t2 = t1*t1; 3090 t5 = r-1.0; 3091 t7 = t5*t5; 3092 t12 = t7*t7; 3093 t18 = t*t; 3094 t19 = t18*t18; 3095 t21 = s+r-1.0; 3096 t24 = t21*t21; 3097 t28 = t24*t24; 3098 f = (t2*t1+42.0/25.0*t2*s*t5+21.0/20.0*t2*t7+7.0/23.0*t7*t5*t1*s+21.0/ 3099 506.0*t1*t12+3.0/1265.0*s*t12*t5+t7*t12/25300)*(t19*t+30.0/11.0*t21*t19+30.0/ 3100 11.0*t18*t*t24+40.0/33.0*t18*t24*t21+5.0/22.0*t*t28+t28*t21/77); 3101 /* p=15 0 7 4 */ 3102 } else if( i==0 && j==7 && k==4 ) { 3103 t1 = s*s; 3104 t2 = t1*s; 3105 t3 = t1*t1; 3106 t6 = r-1.0; 3107 t9 = t6*t6; 3108 t11 = t9*t6; 3109 t13 = t9*t9; 3110 t21 = t*t; 3111 t22 = t21*t21; 3112 t24 = s+r-1.0; 3113 t26 = t24*t24; 3114 t30 = t26*t26; 3115 f = (t3*t2+56.0/25.0*t3*t1*t6+49.0/25.0*t3*s*t9+98.0/115.0*t3*t11+49.0/ 3116 253.0*t2*t13+28.0/1265.0*t1*t13*t6+7.0/6325.0*s*t13*t9+2.0/120175.0*t13*t11)*( 3117 t22+20.0/9.0*t21*t*t24+5.0/3.0*t21*t26+10.0/21.0*t*t26*t24+5.0/126.0*t30); 3118 /* p=15 0 8 3 */ 3119 } else if( i==0 && j==8 && k==3 ) { 3120 t1 = s*s; 3121 t2 = t1*t1; 3122 t3 = t2*t2; 3123 t4 = t1*s; 3124 t6 = r-1.0; 3125 t9 = t6*t6; 3126 t12 = t9*t6; 3127 t14 = t9*t9; 3128 t22 = t14*t14; 3129 t24 = t*t; 3130 t26 = s+r-1.0; 3131 t28 = t26*t26; 3132 f = (t3+72.0/25.0*t2*t4*t6+84.0/25.0*t2*t1*t9+1176.0/575.0*t2*s*t12+ 3133 882.0/1265.0*t2*t14+168.0/1265.0*t4*t14*t6+84.0/6325.0*t1*t14*t9+72.0/120175.0* 3134 s*t14*t12+t22/120175)*(t24*t+12.0/7.0*t26*t24+6.0/7.0*t*t28+4.0/35.0*t28*t26); 3135 /* p=15 0 9 2 */ 3136 } else if( i==0 && j==9 && k==2 ) { 3137 t1 = s*s; 3138 t2 = t1*t1; 3139 t3 = t2*t2; 3140 t5 = r-1.0; 3141 t7 = t1*s; 3142 t9 = t5*t5; 3143 t12 = t9*t5; 3144 t15 = t9*t9; 3145 t23 = t15*t15; 3146 t27 = t*t; 3147 t28 = s+r-1.0; 3148 t30 = t28*t28; 3149 f = (t3*s+18.0/5.0*t3*t5+27.0/5.0*t2*t7*t9+504.0/115.0*t2*t1*t12+2646.0 3150 /1265.0*t2*s*t15+756.0/1265.0*t2*t15*t5+126.0/1265.0*t7*t15*t9+216.0/24035.0*t1 3151 *t15*t12+9.0/24035.0*s*t23+2.0/408595.0*t23*t5)*(t27+6.0/5.0*t*t28+3.0/10.0*t30 3152 ); 3153 /* p=15 0 10 1 */ 3154 } else if( i==0 && j==10 && k==1 ) { 3155 t1 = s*s; 3156 t2 = t1*t1; 3157 t3 = t2*t2; 3158 t6 = r-1.0; 3159 t8 = t6*t6; 3160 t10 = t1*s; 3161 t12 = t8*t6; 3162 t15 = t8*t8; 3163 t24 = t15*t15; 3164 t29 = t3*t1+22.0/5.0*t3*s*t6+33.0/4.0*t3*t8+198.0/23.0*t2*t10*t12+126.0/ 3165 23.0*t2*t1*t15+252.0/115.0*t2*s*t15*t6+63.0/115.0*t2*t15*t8+36.0/437.0*t10*t15* 3166 t12+3.0/437.0*t1*t24+2.0/7429.0*s*t24*t6+t24*t8/297160; 3167 f = t29*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 3168 /* p=15 0 11 0 */ 3169 } else if( i==0 && j==11 && k==0 ) { 3170 t1 = s*s; 3171 t2 = t1*s; 3172 t3 = t1*t1; 3173 t4 = t3*t3; 3174 t7 = r-1.0; 3175 t10 = t7*t7; 3176 t12 = t10*t7; 3177 t15 = t10*t10; 3178 t25 = t15*t15; 3179 f = t2*t4+132.0/25.0*t4*t1*t7+121.0/10.0*t4*s*t10+363.0/23.0*t4*t12+ 3180 297.0/23.0*t3*t2*t15+792.0/115.0*t3*t1*t15*t7+1386.0/575.0*t3*s*t15*t10+1188.0/ 3181 2185.0*t3*t15*t12+33.0/437.0*t2*t25+44.0/7429.0*t1*t25*t7+33.0/148580.0*s*t25* 3182 t10+t25*t12/371450; 3183 /* p=15 1 0 10 */ 3184 } else if( i==1 && j==0 && k==10 ) { 3185 t2 = t*t; 3186 t3 = t2*t2; 3187 t4 = t3*t3; 3188 t7 = s+r-1.0; 3189 t9 = t7*t7; 3190 t11 = t2*t; 3191 t13 = t9*t7; 3192 t16 = t9*t9; 3193 t25 = t16*t16; 3194 t30 = t2*t4+110.0/21.0*t4*t*t7+165.0/14.0*t4*t9+1980.0/133.0*t3*t11*t13+ 3195 220.0/19.0*t3*t2*t16+1848.0/323.0*t3*t*t16*t7+1155.0/646.0*t3*t16*t9+110.0/ 3196 323.0*t11*t16*t13+165.0/4522.0*t2*t25+55.0/29393.0*t*t25*t7+11.0/352716.0*t25* 3197 t9; 3198 f = (r-2.0/27.0)*t30; 3199 /* p=15 1 1 9 */ 3200 } else if( i==1 && j==1 && k==9 ) { 3201 t4 = t*t; 3202 t5 = t4*t4; 3203 t6 = t5*t5; 3204 t8 = s+r-1.0; 3205 t10 = t4*t; 3206 t12 = t8*t8; 3207 t15 = t12*t8; 3208 t18 = t12*t12; 3209 t26 = t18*t18; 3210 f = (s+2.0/23.0*r-2.0/23.0)*(r-2.0/27.0)*(t6*t+90.0/19.0*t6*t8+180.0/ 3211 19.0*t5*t10*t12+3360.0/323.0*t5*t4*t15+2205.0/323.0*t5*t*t18+882.0/323.0*t5*t18 3212 *t8+210.0/323.0*t10*t18*t12+360.0/4199.0*t4*t18*t15+45.0/8398.0*t*t26+5.0/ 3213 46189.0*t26*t8); 3214 /* p=15 1 2 8 */ 3215 } else if( i==1 && j==2 && k==8 ) { 3216 t2 = s*s; 3217 t3 = r-1.0; 3218 t5 = t3*t3; 3219 t8 = t*t; 3220 t9 = t8*t8; 3221 t10 = t9*t9; 3222 t11 = t8*t; 3223 t13 = s+r-1.0; 3224 t16 = t13*t13; 3225 t19 = t16*t13; 3226 t21 = t16*t16; 3227 t29 = t21*t21; 3228 f = (r-2.0/27.0)*(t2+6.0/23.0*s*t3+3.0/253.0*t5)*(t10+72.0/17.0*t9*t11* 3229 t13+126.0/17.0*t9*t8*t16+588.0/85.0*t9*t*t19+63.0/17.0*t9*t21+252.0/221.0*t11* 3230 t21*t13+42.0/221.0*t8*t21*t16+36.0/2431.0*t*t21*t19+9.0/24310.0*t29); 3231 /* p=15 1 3 7 */ 3232 } else if( i==1 && j==3 && k==7 ) { 3233 t2 = s*s; 3234 t4 = r-1.0; 3235 t6 = t4*t4; 3236 t11 = t*t; 3237 t12 = t11*t; 3238 t13 = t11*t11; 3239 t16 = s+r-1.0; 3240 t19 = t16*t16; 3241 t21 = t19*t16; 3242 t23 = t19*t19; 3243 f = (r-2.0/27.0)*(t2*s+12.0/23.0*t2*t4+18.0/253.0*s*t6+4.0/1771.0*t6*t4 3244 )*(t12*t13+56.0/15.0*t11*t16*t13+28.0/5.0*t13*t*t19+56.0/13.0*t21*t13+70.0/39.0 3245 *t12*t23+56.0/143.0*t11*t23*t16+28.0/715.0*t*t23*t19+8.0/6435.0*t23*t21); 3246 /* p=15 1 4 6 */ 3247 } else if( i==1 && j==4 && k==6 ) { 3248 t2 = s*s; 3249 t3 = t2*t2; 3250 t5 = r-1.0; 3251 t7 = t5*t5; 3252 t11 = t7*t7; 3253 t14 = t*t; 3254 t15 = t14*t14; 3255 t18 = s+r-1.0; 3256 t20 = t18*t18; 3257 t25 = t20*t20; 3258 f = (r-2.0/27.0)*(t3+20.0/23.0*t5*t2*s+60.0/253.0*t2*t7+40.0/1771.0*s* 3259 t7*t5+t11/1771)*(t14*t15+42.0/13.0*t15*t*t18+105.0/26.0*t15*t20+350.0/143.0*t14 3260 *t*t20*t18+105.0/143.0*t14*t25+14.0/143.0*t*t25*t18+7.0/1716.0*t25*t20); 3261 /* p=15 1 5 5 */ 3262 } else if( i==1 && j==5 && k==5 ) { 3263 t2 = s*s; 3264 t3 = t2*t2; 3265 t5 = r-1.0; 3266 t8 = t5*t5; 3267 t12 = t8*t8; 3268 t17 = t*t; 3269 t18 = t17*t17; 3270 t20 = s+r-1.0; 3271 t23 = t20*t20; 3272 t27 = t23*t23; 3273 f = (r-2.0/27.0)*(s*t3+30.0/23.0*t3*t5+150.0/253.0*t2*s*t8+200.0/1771.0 3274 *t2*t8*t5+15.0/1771.0*s*t12+6.0/33649.0*t12*t5)*(t18*t+30.0/11.0*t20*t18+30.0/ 3275 11.0*t17*t*t23+40.0/33.0*t17*t23*t20+5.0/22.0*t*t27+t27*t20/77); 3276 /* p=15 1 6 4 */ 3277 } else if( i==1 && j==6 && k==4 ) { 3278 t2 = s*s; 3279 t3 = t2*t2; 3280 t6 = r-1.0; 3281 t8 = t6*t6; 3282 t13 = t8*t8; 3283 t20 = t*t; 3284 t21 = t20*t20; 3285 t23 = s+r-1.0; 3286 t25 = t23*t23; 3287 t29 = t25*t25; 3288 f = (r-2.0/27.0)*(t3*t2+42.0/23.0*t3*s*t6+315.0/253.0*t3*t8+100.0/253.0 3289 *t2*s*t8*t6+15.0/253.0*t2*t13+18.0/4807.0*s*t13*t6+t13*t8/14421)*(t21+20.0/9.0* 3290 t20*t*t23+5.0/3.0*t20*t25+10.0/21.0*t*t25*t23+5.0/126.0*t29); 3291 /* p=15 1 7 3 */ 3292 } else if( i==1 && j==7 && k==3 ) { 3293 t2 = s*s; 3294 t3 = t2*s; 3295 t4 = t2*t2; 3296 t7 = r-1.0; 3297 t10 = t7*t7; 3298 t12 = t10*t7; 3299 t14 = t10*t10; 3300 t23 = t*t; 3301 t25 = s+r-1.0; 3302 t27 = t25*t25; 3303 f = (r-2.0/27.0)*(t4*t3+56.0/23.0*t4*t2*t7+588.0/253.0*t4*s*t10+280.0/ 3304 253.0*t4*t12+70.0/253.0*t3*t14+168.0/4807.0*t2*t14*t7+28.0/14421.0*s*t14*t10+ 3305 8.0/245157.0*t14*t12)*(t23*t+12.0/7.0*t25*t23+6.0/7.0*t*t27+4.0/35.0*t27*t25); 3306 /* p=15 1 8 2 */ 3307 } else if( i==1 && j==8 && k==2 ) { 3308 t2 = s*s; 3309 t3 = t2*t2; 3310 t4 = t3*t3; 3311 t5 = t2*s; 3312 t7 = r-1.0; 3313 t10 = t7*t7; 3314 t13 = t10*t7; 3315 t15 = t10*t10; 3316 t23 = t15*t15; 3317 t26 = t*t; 3318 t27 = s+r-1.0; 3319 t29 = t27*t27; 3320 f = (r-2.0/27.0)*(t4+72.0/23.0*t3*t5*t7+1008.0/253.0*t3*t2*t10+672.0/ 3321 253.0*t3*s*t13+252.0/253.0*t3*t15+1008.0/4807.0*t5*t15*t7+112.0/4807.0*t2*t15* 3322 t10+96.0/81719.0*s*t15*t13+3.0/163438.0*t23)*(t26+6.0/5.0*t*t27+3.0/10.0*t29); 3323 /* p=15 1 9 1 */ 3324 } else if( i==1 && j==9 && k==1 ) { 3325 t2 = s*s; 3326 t3 = t2*t2; 3327 t4 = t3*t3; 3328 t6 = r-1.0; 3329 t8 = t2*s; 3330 t10 = t6*t6; 3331 t13 = t10*t6; 3332 t16 = t10*t10; 3333 t24 = t16*t16; 3334 f = (r-2.0/27.0)*(t4*s+90.0/23.0*t4*t6+1620.0/253.0*t3*t8*t10+1440.0/ 3335 253.0*t3*t2*t13+756.0/253.0*t3*s*t16+4536.0/4807.0*t3*t16*t6+840.0/4807.0*t8* 3336 t16*t10+1440.0/81719.0*t2*t16*t13+135.0/163438.0*s*t24+t24*t6/81719)*(t+2.0/3.0 3337 *s+2.0/3.0*r-2.0/3.0); 3338 /* p=15 1 10 0 */ 3339 } else if( i==1 && j==10 && k==0 ) { 3340 t2 = s*s; 3341 t3 = t2*t2; 3342 t4 = t3*t3; 3343 t7 = r-1.0; 3344 t9 = t7*t7; 3345 t11 = t2*s; 3346 t13 = t9*t7; 3347 t16 = t9*t9; 3348 t25 = t16*t16; 3349 t30 = t4*t2+110.0/23.0*t4*s*t7+225.0/23.0*t4*t9+1800.0/161.0*t3*t11*t13+ 3350 180.0/23.0*t3*t2*t16+1512.0/437.0*t3*s*t16*t7+420.0/437.0*t3*t16*t9+1200.0/ 3351 7429.0*t11*t16*t13+225.0/14858.0*t2*t25+5.0/7429.0*s*t25*t7+t25*t9/104006; 3352 f = (r-2.0/27.0)*t30; 3353 /* p=15 2 0 9 */ 3354 } else if( i==2 && j==0 && k==9 ) { 3355 t1 = r*r; 3356 t3 = t*t; 3357 t4 = t3*t3; 3358 t5 = t4*t4; 3359 t7 = s+r-1.0; 3360 t9 = t3*t; 3361 t11 = t7*t7; 3362 t14 = t11*t7; 3363 t17 = t11*t11; 3364 t25 = t17*t17; 3365 f = (t1-2.0/9.0*r+1.0/117.0)*(t5*t+90.0/19.0*t5*t7+180.0/19.0*t4*t9*t11 3366 +3360.0/323.0*t4*t3*t14+2205.0/323.0*t4*t*t17+882.0/323.0*t4*t17*t7+210.0/323.0 3367 *t9*t17*t11+360.0/4199.0*t3*t17*t14+45.0/8398.0*t*t25+5.0/46189.0*t25*t7); 3368 /* p=15 2 1 8 */ 3369 } else if( i==2 && j==1 && k==8 ) { 3370 t1 = r*r; 3371 t5 = t*t; 3372 t6 = t5*t5; 3373 t7 = t6*t6; 3374 t8 = t5*t; 3375 t10 = s+r-1.0; 3376 t13 = t10*t10; 3377 t16 = t13*t10; 3378 t18 = t13*t13; 3379 t26 = t18*t18; 3380 f = (t1-2.0/9.0*r+1.0/117.0)*(s+2.0/21.0*r-2.0/21.0)*(t7+72.0/17.0*t6* 3381 t8*t10+126.0/17.0*t6*t5*t13+588.0/85.0*t6*t*t16+63.0/17.0*t6*t18+252.0/221.0*t8 3382 *t18*t10+42.0/221.0*t5*t18*t13+36.0/2431.0*t*t18*t16+9.0/24310.0*t26); 3383 /* p=15 2 2 7 */ 3384 } else if( i==2 && j==2 && k==7 ) { 3385 t1 = r*r; 3386 t3 = s*s; 3387 t4 = r-1.0; 3388 t6 = t4*t4; 3389 t9 = t*t; 3390 t10 = t9*t; 3391 t11 = t9*t9; 3392 t14 = s+r-1.0; 3393 t17 = t14*t14; 3394 t19 = t17*t14; 3395 t21 = t17*t17; 3396 f = (t1-2.0/9.0*r+1.0/117.0)*(t3+2.0/7.0*s*t4+t6/70)*(t11*t10+56.0/15.0 3397 *t11*t9*t14+28.0/5.0*t11*t*t17+56.0/13.0*t11*t19+70.0/39.0*t10*t21+56.0/143.0* 3398 t9*t21*t14+28.0/715.0*t*t21*t17+8.0/6435.0*t21*t19); 3399 /* p=15 2 3 6 */ 3400 } else if( i==2 && j==3 && k==6 ) { 3401 t1 = r*r; 3402 t3 = s*s; 3403 t5 = r-1.0; 3404 t7 = t5*t5; 3405 t12 = t*t; 3406 t13 = t12*t12; 3407 t16 = s+r-1.0; 3408 t18 = t16*t16; 3409 t23 = t18*t18; 3410 f = (t1-2.0/9.0*r+1.0/117.0)*(t3*s+4.0/7.0*t3*t5+3.0/35.0*s*t7+2.0/ 3411 665.0*t5*t7)*(t12*t13+42.0/13.0*t13*t*t16+105.0/26.0*t18*t13+350.0/143.0*t12*t* 3412 t18*t16+105.0/143.0*t12*t23+14.0/143.0*t*t23*t16+7.0/1716.0*t23*t18); 3413 /* p=15 2 4 5 */ 3414 } else if( i==2 && j==4 && k==5 ) { 3415 t1 = r*r; 3416 t3 = s*s; 3417 t4 = t3*t3; 3418 t6 = r-1.0; 3419 t8 = t6*t6; 3420 t12 = t8*t8; 3421 t15 = t*t; 3422 t16 = t15*t15; 3423 t18 = s+r-1.0; 3424 t21 = t18*t18; 3425 t25 = t21*t21; 3426 f = (t1-2.0/9.0*r+1.0/117.0)*(t4+20.0/21.0*t6*t3*s+2.0/7.0*t3*t8+4.0/ 3427 133.0*s*t6*t8+t12/1197)*(t16*t+30.0/11.0*t18*t16+30.0/11.0*t15*t*t21+40.0/33.0* 3428 t15*t21*t18+5.0/22.0*t*t25+t25*t18/77); 3429 /* p=15 2 5 4 */ 3430 } else if( i==2 && j==5 && k==4 ) { 3431 t1 = r*r; 3432 t3 = s*s; 3433 t4 = t3*t3; 3434 t6 = r-1.0; 3435 t9 = t6*t6; 3436 t13 = t9*t9; 3437 t18 = t*t; 3438 t19 = t18*t18; 3439 t21 = s+r-1.0; 3440 t23 = t21*t21; 3441 t27 = t23*t23; 3442 f = (t1-2.0/9.0*r+1.0/117.0)*(s*t4+10.0/7.0*t4*t6+5.0/7.0*t3*s*t9+20.0/ 3443 133.0*t3*t9*t6+5.0/399.0*s*t13+2.0/6783.0*t13*t6)*(t19+20.0/9.0*t18*t*t21+5.0/ 3444 3.0*t23*t18+10.0/21.0*t*t23*t21+5.0/126.0*t27); 3445 /* p=15 2 6 3 */ 3446 } else if( i==2 && j==6 && k==3 ) { 3447 t1 = r*r; 3448 t3 = s*s; 3449 t4 = t3*t3; 3450 t7 = r-1.0; 3451 t9 = t7*t7; 3452 t14 = t9*t9; 3453 t21 = t*t; 3454 t23 = s+r-1.0; 3455 t25 = t23*t23; 3456 f = (t1-2.0/9.0*r+1.0/117.0)*(t4*t3+2.0*t4*s*t7+3.0/2.0*t4*t9+10.0/19.0 3457 *t3*s*t9*t7+5.0/57.0*t3*t14+2.0/323.0*s*t14*t7+t14*t9/7752)*(t21*t+12.0/7.0*t23 3458 *t21+6.0/7.0*t*t25+4.0/35.0*t25*t23); 3459 /* p=15 2 7 2 */ 3460 } else if( i==2 && j==7 && k==2 ) { 3461 t1 = r*r; 3462 t3 = s*s; 3463 t4 = t3*s; 3464 t5 = t3*t3; 3465 t8 = r-1.0; 3466 t11 = t8*t8; 3467 t13 = t11*t8; 3468 t15 = t11*t11; 3469 t24 = t*t; 3470 t25 = s+r-1.0; 3471 t27 = t25*t25; 3472 f = (t1-2.0/9.0*r+1.0/117.0)*(t5*t4+8.0/3.0*t5*t3*t8+14.0/5.0*t5*s*t11+ 3473 28.0/19.0*t5*t13+70.0/171.0*t4*t15+56.0/969.0*t3*t15*t8+7.0/1938.0*s*t15*t11+ 3474 t15*t13/14535)*(t24+6.0/5.0*t*t25+3.0/10.0*t27); 3475 /* p=15 2 8 1 */ 3476 } else if( i==2 && j==8 && k==1 ) { 3477 t1 = r*r; 3478 t3 = s*s; 3479 t4 = t3*t3; 3480 t5 = t4*t4; 3481 t6 = t3*s; 3482 t8 = r-1.0; 3483 t11 = t8*t8; 3484 t14 = t11*t8; 3485 t16 = t11*t11; 3486 t24 = t16*t16; 3487 f = (t1-2.0/9.0*r+1.0/117.0)*(t5+24.0/7.0*t4*t6*t8+24.0/5.0*t4*t3*t11+ 3488 336.0/95.0*t4*s*t14+28.0/19.0*t4*t16+112.0/323.0*t6*t16*t8+14.0/323.0*t3*t16* 3489 t11+4.0/1615.0*s*t16*t14+t24/22610)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 3490 /* p=15 2 9 0 */ 3491 } else if( i==2 && j==9 && k==0 ) { 3492 t1 = r*r; 3493 t3 = s*s; 3494 t4 = t3*t3; 3495 t5 = t4*t4; 3496 t7 = r-1.0; 3497 t9 = t3*s; 3498 t11 = t7*t7; 3499 t14 = t11*t7; 3500 t17 = t11*t11; 3501 t25 = t17*t17; 3502 f = (t1-2.0/9.0*r+1.0/117.0)*(t5*s+30.0/7.0*t5*t7+54.0/7.0*t4*t9*t11+ 3503 144.0/19.0*t4*t3*t14+84.0/19.0*t4*s*t17+504.0/323.0*t4*t17*t7+105.0/323.0*t9* 3504 t17*t11+12.0/323.0*t3*t17*t14+9.0/4522.0*s*t25+t25*t7/29393); 3505 /* p=15 3 0 8 */ 3506 } else if( i==3 && j==0 && k==8 ) { 3507 t1 = r*r; 3508 t4 = t*t; 3509 t5 = t4*t4; 3510 t6 = t5*t5; 3511 t7 = t4*t; 3512 t9 = s+r-1.0; 3513 t12 = t9*t9; 3514 t15 = t12*t9; 3515 t17 = t12*t12; 3516 t25 = t17*t17; 3517 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t6+72.0/17.0*t5*t7*t9+ 3518 126.0/17.0*t5*t4*t12+588.0/85.0*t5*t*t15+63.0/17.0*t5*t17+252.0/221.0*t7*t17*t9 3519 +42.0/221.0*t4*t17*t12+36.0/2431.0*t*t17*t15+9.0/24310.0*t25); 3520 /* p=15 3 1 7 */ 3521 } else if( i==3 && j==1 && k==7 ) { 3522 t1 = r*r; 3523 t6 = t*t; 3524 t7 = t6*t; 3525 t8 = t6*t6; 3526 t11 = s+r-1.0; 3527 t14 = t11*t11; 3528 t16 = t14*t11; 3529 t18 = t14*t14; 3530 f = (s+2.0/19.0*r-2.0/19.0)*(t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t7 3531 *t8+56.0/15.0*t8*t6*t11+28.0/5.0*t8*t*t14+56.0/13.0*t16*t8+70.0/39.0*t7*t18+ 3532 56.0/143.0*t6*t18*t11+28.0/715.0*t*t18*t14+8.0/6435.0*t18*t16); 3533 /* p=15 3 2 6 */ 3534 } else if( i==3 && j==2 && k==6 ) { 3535 t1 = r*r; 3536 t4 = s*s; 3537 t5 = r-1.0; 3538 t7 = t5*t5; 3539 t10 = t*t; 3540 t11 = t10*t10; 3541 t14 = s+r-1.0; 3542 t16 = t14*t14; 3543 t21 = t16*t16; 3544 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t4+6.0/19.0*s*t5+t7/57)*( 3545 t10*t11+42.0/13.0*t11*t*t14+105.0/26.0*t16*t11+350.0/143.0*t10*t*t16*t14+105.0/ 3546 143.0*t10*t21+14.0/143.0*t*t21*t14+7.0/1716.0*t21*t16); 3547 /* p=15 3 3 5 */ 3548 } else if( i==3 && j==3 && k==5 ) { 3549 t1 = r*r; 3550 t4 = s*s; 3551 t6 = r-1.0; 3552 t8 = t6*t6; 3553 t13 = t*t; 3554 t14 = t13*t13; 3555 t16 = s+r-1.0; 3556 t19 = t16*t16; 3557 t23 = t19*t19; 3558 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t4*s+12.0/19.0*t4*t6+2.0/ 3559 19.0*s*t8+4.0/969.0*t8*t6)*(t14*t+30.0/11.0*t16*t14+30.0/11.0*t13*t*t19+40.0/ 3560 33.0*t13*t19*t16+5.0/22.0*t*t23+t23*t16/77); 3561 /* p=15 3 4 4 */ 3562 } else if( i==3 && j==4 && k==4 ) { 3563 t1 = r*r; 3564 t4 = s*s; 3565 t5 = t4*t4; 3566 t7 = r-1.0; 3567 t9 = t7*t7; 3568 t13 = t9*t9; 3569 t16 = t*t; 3570 t17 = t16*t16; 3571 t19 = s+r-1.0; 3572 t21 = t19*t19; 3573 t25 = t21*t21; 3574 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t5+20.0/19.0*t4*s*t7+20.0/ 3575 57.0*t4*t9+40.0/969.0*s*t9*t7+5.0/3876.0*t13)*(t17+20.0/9.0*t16*t*t19+5.0/3.0* 3576 t21*t16+10.0/21.0*t*t21*t19+5.0/126.0*t25); 3577 /* p=15 3 5 3 */ 3578 } else if( i==3 && j==5 && k==3 ) { 3579 t1 = r*r; 3580 t4 = s*s; 3581 t5 = t4*t4; 3582 t7 = r-1.0; 3583 t10 = t7*t7; 3584 t14 = t10*t10; 3585 t19 = t*t; 3586 t21 = s+r-1.0; 3587 t23 = t21*t21; 3588 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(s*t5+30.0/19.0*t5*t7+50.0/ 3589 57.0*t4*s*t10+200.0/969.0*t4*t10*t7+25.0/1292.0*s*t14+t14*t7/1938)*(t19*t+12.0/ 3590 7.0*t21*t19+6.0/7.0*t*t23+4.0/35.0*t23*t21); 3591 /* p=15 3 6 2 */ 3592 } else if( i==3 && j==6 && k==2 ) { 3593 t1 = r*r; 3594 t4 = s*s; 3595 t5 = t4*t4; 3596 t8 = r-1.0; 3597 t10 = t8*t8; 3598 t15 = t10*t10; 3599 t22 = t*t; 3600 t23 = s+r-1.0; 3601 t25 = t23*t23; 3602 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t5*t4+42.0/19.0*t5*s*t8+ 3603 35.0/19.0*t5*t10+700.0/969.0*t4*s*t10*t8+175.0/1292.0*t4*t15+7.0/646.0*s*t15*t8 3604 +t15*t10/3876)*(t22+6.0/5.0*t*t23+3.0/10.0*t25); 3605 /* p=15 3 7 1 */ 3606 } else if( i==3 && j==7 && k==1 ) { 3607 t1 = r*r; 3608 t4 = s*s; 3609 t5 = t4*s; 3610 t6 = t4*t4; 3611 t9 = r-1.0; 3612 t12 = t9*t9; 3613 t14 = t12*t9; 3614 t16 = t12*t12; 3615 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t6*t5+56.0/19.0*t4*t6*t9+ 3616 196.0/57.0*t6*s*t12+1960.0/969.0*t6*t14+1225.0/1938.0*t5*t16+98.0/969.0*t4*t16* 3617 t9+7.0/969.0*s*t16*t12+2.0/12597.0*t16*t14)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 3618 /* p=15 3 8 0 */ 3619 } else if( i==3 && j==8 && k==0 ) { 3620 t1 = r*r; 3621 t4 = s*s; 3622 t5 = t4*t4; 3623 t6 = t5*t5; 3624 t7 = t4*s; 3625 t9 = r-1.0; 3626 t12 = t9*t9; 3627 t15 = t12*t9; 3628 t17 = t12*t12; 3629 t25 = t17*t17; 3630 f = (t1*r-4.0/9.0*t1+2.0/39.0*r-4.0/2925.0)*(t6+72.0/19.0*t5*t7*t9+ 3631 112.0/19.0*t5*t4*t12+1568.0/323.0*t5*s*t15+735.0/323.0*t5*t17+196.0/323.0*t7* 3632 t17*t9+28.0/323.0*t4*t17*t12+24.0/4199.0*s*t17*t15+t25/8398); 3633 /* p=15 4 0 7 */ 3634 } else if( i==4 && j==0 && k==7 ) { 3635 t1 = r*r; 3636 t2 = t1*t1; 3637 t5 = t*t; 3638 t6 = t5*t; 3639 t7 = t5*t5; 3640 t10 = s+r-1.0; 3641 t13 = t10*t10; 3642 t15 = t13*t10; 3643 t17 = t13*t13; 3644 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t7*t6+ 3645 56.0/15.0*t5*t7*t10+28.0/5.0*t7*t*t13+56.0/13.0*t7*t15+70.0/39.0*t6*t17+56.0/ 3646 143.0*t5*t17*t10+28.0/715.0*t*t17*t13+8.0/6435.0*t17*t15); 3647 /* p=15 4 1 6 */ 3648 } else if( i==4 && j==1 && k==6 ) { 3649 t1 = r*r; 3650 t2 = t1*t1; 3651 t7 = t*t; 3652 t8 = t7*t7; 3653 t11 = s+r-1.0; 3654 t13 = t11*t11; 3655 t18 = t13*t13; 3656 f = (s+2.0/17.0*r-2.0/17.0)*(t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0* 3657 r+1.0/3510.0)*(t8*t7+42.0/13.0*t8*t*t11+105.0/26.0*t8*t13+350.0/143.0*t7*t*t13* 3658 t11+105.0/143.0*t7*t18+14.0/143.0*t*t18*t11+7.0/1716.0*t18*t13); 3659 /* p=15 4 2 5 */ 3660 } else if( i==4 && j==2 && k==5 ) { 3661 t1 = r*r; 3662 t2 = t1*t1; 3663 t5 = s*s; 3664 t6 = r-1.0; 3665 t8 = t6*t6; 3666 t11 = t*t; 3667 t12 = t11*t11; 3668 t14 = s+r-1.0; 3669 t17 = t14*t14; 3670 t21 = t17*t17; 3671 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t5+6.0/ 3672 17.0*s*t6+3.0/136.0*t8)*(t12*t+30.0/11.0*t12*t14+30.0/11.0*t11*t*t17+40.0/33.0* 3673 t11*t17*t14+5.0/22.0*t*t21+t21*t14/77); 3674 /* p=15 4 3 4 */ 3675 } else if( i==4 && j==3 && k==4 ) { 3676 t1 = r*r; 3677 t2 = t1*t1; 3678 t5 = s*s; 3679 t7 = r-1.0; 3680 t9 = t7*t7; 3681 t14 = t*t; 3682 t15 = t14*t14; 3683 t17 = s+r-1.0; 3684 t19 = t17*t17; 3685 t23 = t19*t19; 3686 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t5*s+12.0 3687 /17.0*t5*t7+9.0/68.0*s*t9+t9*t7/170)*(t15+20.0/9.0*t14*t*t17+5.0/3.0*t14*t19+ 3688 10.0/21.0*t*t19*t17+5.0/126.0*t23); 3689 /* p=15 4 4 3 */ 3690 } else if( i==4 && j==4 && k==3 ) { 3691 t1 = r*r; 3692 t2 = t1*t1; 3693 t5 = s*s; 3694 t6 = t5*t5; 3695 t8 = r-1.0; 3696 t10 = t8*t8; 3697 t14 = t10*t10; 3698 t17 = t*t; 3699 t19 = s+r-1.0; 3700 t21 = t19*t19; 3701 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t6+20.0/ 3702 17.0*t8*t5*s+15.0/34.0*t5*t10+s*t10*t8/17+t14/476)*(t17*t+12.0/7.0*t19*t17+6.0/ 3703 7.0*t*t21+4.0/35.0*t21*t19); 3704 /* p=15 4 5 2 */ 3705 } else if( i==4 && j==5 && k==2 ) { 3706 t1 = r*r; 3707 t2 = t1*t1; 3708 t5 = s*s; 3709 t6 = t5*t5; 3710 t8 = r-1.0; 3711 t11 = t8*t8; 3712 t15 = t11*t11; 3713 t20 = t*t; 3714 t21 = s+r-1.0; 3715 t23 = t21*t21; 3716 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(s*t6+30.0 3717 /17.0*t6*t8+75.0/68.0*t5*s*t11+5.0/17.0*t5*t11*t8+15.0/476.0*s*t15+3.0/3094.0* 3718 t15*t8)*(t20+6.0/5.0*t*t21+3.0/10.0*t23); 3719 /* p=15 4 6 1 */ 3720 } else if( i==4 && j==6 && k==1 ) { 3721 t1 = r*r; 3722 t2 = t1*t1; 3723 t5 = s*s; 3724 t6 = t5*t5; 3725 t9 = r-1.0; 3726 t11 = t9*t9; 3727 t16 = t11*t11; 3728 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t6*t5+ 3729 42.0/17.0*t6*s*t9+315.0/136.0*t6*t11+35.0/34.0*t5*s*t11*t9+15.0/68.0*t5*t16+9.0 3730 /442.0*s*t16*t9+t16*t11/1768)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 3731 /* p=15 4 7 0 */ 3732 } else if( i==4 && j==7 && k==0 ) { 3733 t1 = r*r; 3734 t2 = t1*t1; 3735 t5 = s*s; 3736 t6 = t5*s; 3737 t7 = t5*t5; 3738 t10 = r-1.0; 3739 t13 = t10*t10; 3740 t15 = t13*t10; 3741 t17 = t13*t13; 3742 f = (t2-20.0/27.0*t1*r+20.0/117.0*t1-8.0/585.0*r+1.0/3510.0)*(t7*t6+ 3743 56.0/17.0*t7*t5*t10+147.0/34.0*t7*s*t13+49.0/17.0*t7*t15+35.0/34.0*t6*t17+42.0/ 3744 221.0*t5*t17*t10+7.0/442.0*s*t17*t13+t17*t15/2431); 3745 /* p=15 5 0 6 */ 3746 } else if( i==5 && j==0 && k==6 ) { 3747 t1 = r*r; 3748 t2 = t1*t1; 3749 t6 = t*t; 3750 t7 = t6*t6; 3751 t10 = s+r-1.0; 3752 t12 = t10*t10; 3753 t17 = t12*t12; 3754 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*( 3755 t7*t6+42.0/13.0*t7*t*t10+105.0/26.0*t7*t12+350.0/143.0*t6*t*t12*t10+105.0/143.0 3756 *t6*t17+14.0/143.0*t*t17*t10+7.0/1716.0*t17*t12); 3757 /* p=15 5 1 5 */ 3758 } else if( i==5 && j==1 && k==5 ) { 3759 t1 = r*r; 3760 t2 = t1*t1; 3761 t8 = t*t; 3762 t9 = t8*t8; 3763 t11 = s+r-1.0; 3764 t14 = t11*t11; 3765 t18 = t14*t14; 3766 f = (s+2.0/15.0*r-2.0/15.0)*(t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0 3767 *t1+r/234-1.0/13455.0)*(t9*t+30.0/11.0*t11*t9+30.0/11.0*t8*t*t14+40.0/33.0*t8* 3768 t14*t11+5.0/22.0*t*t18+t18*t11/77); 3769 /* p=15 5 2 4 */ 3770 } else if( i==5 && j==2 && k==4 ) { 3771 t1 = r*r; 3772 t2 = t1*t1; 3773 t6 = s*s; 3774 t7 = r-1.0; 3775 t9 = t7*t7; 3776 t12 = t*t; 3777 t13 = t12*t12; 3778 t15 = s+r-1.0; 3779 t17 = t15*t15; 3780 t21 = t17*t17; 3781 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*( 3782 t6+2.0/5.0*s*t7+t9/35)*(t13+20.0/9.0*t12*t*t15+5.0/3.0*t17*t12+10.0/21.0*t*t17* 3783 t15+5.0/126.0*t21); 3784 /* p=15 5 3 3 */ 3785 } else if( i==5 && j==3 && k==3 ) { 3786 t1 = r*r; 3787 t2 = t1*t1; 3788 t6 = s*s; 3789 t8 = r-1.0; 3790 t10 = t8*t8; 3791 t15 = t*t; 3792 t17 = s+r-1.0; 3793 t19 = t17*t17; 3794 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*( 3795 s*t6+4.0/5.0*t6*t8+6.0/35.0*s*t10+4.0/455.0*t10*t8)*(t15*t+12.0/7.0*t17*t15+6.0 3796 /7.0*t*t19+4.0/35.0*t19*t17); 3797 /* p=15 5 4 2 */ 3798 } else if( i==5 && j==4 && k==2 ) { 3799 t1 = r*r; 3800 t2 = t1*t1; 3801 t6 = s*s; 3802 t7 = t6*t6; 3803 t9 = r-1.0; 3804 t11 = t9*t9; 3805 t15 = t11*t11; 3806 t18 = t*t; 3807 t19 = s+r-1.0; 3808 t21 = t19*t19; 3809 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*( 3810 t7+4.0/3.0*t6*s*t9+4.0/7.0*t6*t11+8.0/91.0*s*t11*t9+t15/273)*(t18+6.0/5.0*t*t19 3811 +3.0/10.0*t21); 3812 /* p=15 5 5 1 */ 3813 } else if( i==5 && j==5 && k==1 ) { 3814 t1 = r*r; 3815 t2 = t1*t1; 3816 t6 = s*s; 3817 t7 = t6*t6; 3818 t9 = r-1.0; 3819 t12 = t9*t9; 3820 t16 = t12*t12; 3821 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*( 3822 t7*s+2.0*t7*t9+10.0/7.0*t6*s*t12+40.0/91.0*t6*t12*t9+5.0/91.0*s*t16+2.0/1001.0* 3823 t16*t9)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 3824 /* p=15 5 6 0 */ 3825 } else if( i==5 && j==6 && k==0 ) { 3826 t1 = r*r; 3827 t2 = t1*t1; 3828 t6 = s*s; 3829 t7 = t6*t6; 3830 t10 = r-1.0; 3831 t12 = t10*t10; 3832 t17 = t12*t12; 3833 f = (t2*r-10.0/9.0*t2+50.0/117.0*t1*r-8.0/117.0*t1+r/234-1.0/13455.0)*( 3834 t7*t6+14.0/5.0*t7*s*t10+3.0*t7*t12+20.0/13.0*t6*s*t12*t10+5.0/13.0*t6*t17+6.0/ 3835 143.0*s*t17*t10+t17*t12/715); 3836 /* p=15 6 0 5 */ 3837 } else if( i==6 && j==0 && k==5 ) { 3838 t1 = r*r; 3839 t2 = t1*t1; 3840 t7 = t*t; 3841 t8 = t7*t7; 3842 t10 = s+r-1.0; 3843 t13 = t10*t10; 3844 t17 = t13*t13; 3845 f = (t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/117.0*t1*r+7.0/234.0*t1-7.0/ 3846 4485.0*r+7.0/296010.0)*(t8*t+30.0/11.0*t8*t10+30.0/11.0*t7*t*t13+40.0/33.0*t7* 3847 t13*t10+5.0/22.0*t*t17+t17*t10/77); 3848 /* p=15 6 1 4 */ 3849 } else if( i==6 && j==1 && k==4 ) { 3850 t1 = r*r; 3851 t2 = t1*t1; 3852 t9 = t*t; 3853 t10 = t9*t9; 3854 t12 = s+r-1.0; 3855 t14 = t12*t12; 3856 t18 = t14*t14; 3857 f = (s+2.0/13.0*r-2.0/13.0)*(t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/ 3858 117.0*t1*r+7.0/234.0*t1-7.0/4485.0*r+7.0/296010.0)*(t10+20.0/9.0*t9*t*t12+5.0/ 3859 3.0*t9*t14+10.0/21.0*t*t14*t12+5.0/126.0*t18); 3860 /* p=15 6 2 3 */ 3861 } else if( i==6 && j==2 && k==3 ) { 3862 t1 = r*r; 3863 t2 = t1*t1; 3864 t7 = s*s; 3865 t8 = r-1.0; 3866 t10 = t8*t8; 3867 t13 = t*t; 3868 t15 = s+r-1.0; 3869 t17 = t15*t15; 3870 f = (t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/117.0*t1*r+7.0/234.0*t1-7.0/ 3871 4485.0*r+7.0/296010.0)*(t7+6.0/13.0*s*t8+t10/26)*(t13*t+12.0/7.0*t13*t15+6.0/ 3872 7.0*t*t17+4.0/35.0*t17*t15); 3873 /* p=15 6 3 2 */ 3874 } else if( i==6 && j==3 && k==2 ) { 3875 t1 = r*r; 3876 t2 = t1*t1; 3877 t7 = s*s; 3878 t9 = r-1.0; 3879 t11 = t9*t9; 3880 t16 = t*t; 3881 t17 = s+r-1.0; 3882 t19 = t17*t17; 3883 f = (t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/117.0*t1*r+7.0/234.0*t1-7.0/ 3884 4485.0*r+7.0/296010.0)*(t7*s+12.0/13.0*t7*t9+3.0/13.0*s*t11+2.0/143.0*t11*t9)*( 3885 t16+6.0/5.0*t*t17+3.0/10.0*t19); 3886 /* p=15 6 4 1 */ 3887 } else if( i==6 && j==4 && k==1 ) { 3888 t1 = r*r; 3889 t2 = t1*t1; 3890 t7 = s*s; 3891 t8 = t7*t7; 3892 t10 = r-1.0; 3893 t12 = t10*t10; 3894 t16 = t12*t12; 3895 f = (t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/117.0*t1*r+7.0/234.0*t1-7.0/ 3896 4485.0*r+7.0/296010.0)*(t8+20.0/13.0*t7*s*t10+10.0/13.0*t7*t12+20.0/143.0*s*t12 3897 *t10+t16/143)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 3898 /* p=15 6 5 0 */ 3899 } else if( i==6 && j==5 && k==0 ) { 3900 t1 = r*r; 3901 t2 = t1*t1; 3902 t7 = s*s; 3903 t8 = t7*t7; 3904 t10 = r-1.0; 3905 t13 = t10*t10; 3906 t17 = t13*t13; 3907 f = (t2*t1-14.0/9.0*t2*r+35.0/39.0*t2-28.0/117.0*t1*r+7.0/234.0*t1-7.0/ 3908 4485.0*r+7.0/296010.0)*(t8*s+30.0/13.0*t8*t10+25.0/13.0*t7*s*t13+100.0/143.0*t7 3909 *t13*t10+15.0/143.0*s*t17+2.0/429.0*t17*t10); 3910 /* p=15 7 0 4 */ 3911 } else if( i==7 && j==0 && k==4 ) { 3912 t1 = r*r; 3913 t2 = t1*r; 3914 t3 = t1*t1; 3915 t8 = t*t; 3916 t9 = t8*t8; 3917 t11 = s+r-1.0; 3918 t13 = t11*t11; 3919 t17 = t13*t13; 3920 f = (t3*t2-56.0/27.0*t3*t1+196.0/117.0*t3*r-392.0/585.0*t3+49.0/351.0* 3921 t2-196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0)*(t9+20.0/9.0*t8*t*t11+5.0/3.0 3922 *t8*t13+10.0/21.0*t*t13*t11+5.0/126.0*t17); 3923 /* p=15 7 1 3 */ 3924 } else if( i==7 && j==1 && k==3 ) { 3925 t1 = r*r; 3926 t2 = t1*r; 3927 t3 = t1*t1; 3928 t10 = t*t; 3929 t12 = s+r-1.0; 3930 t14 = t12*t12; 3931 f = (t3*t2-56.0/27.0*t3*t1+196.0/117.0*t3*r-392.0/585.0*t3+49.0/351.0* 3932 t2-196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0)*(s+2.0/11.0*r-2.0/11.0)*(t10* 3933 t+12.0/7.0*t12*t10+6.0/7.0*t*t14+4.0/35.0*t14*t12); 3934 /* p=15 7 2 2 */ 3935 } else if( i==7 && j==2 && k==2 ) { 3936 t1 = r*r; 3937 t2 = t1*r; 3938 t3 = t1*t1; 3939 t8 = s*s; 3940 t9 = r-1.0; 3941 t11 = t9*t9; 3942 t14 = t*t; 3943 t15 = s+r-1.0; 3944 t17 = t15*t15; 3945 f = (t3*t2-56.0/27.0*t3*t1+196.0/117.0*t3*r-392.0/585.0*t3+49.0/351.0* 3946 t2-196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0)*(t8+6.0/11.0*s*t9+3.0/55.0* 3947 t11)*(t14+6.0/5.0*t*t15+3.0/10.0*t17); 3948 /* p=15 7 3 1 */ 3949 } else if( i==7 && j==3 && k==1 ) { 3950 t1 = r*r; 3951 t2 = t1*r; 3952 t3 = t1*t1; 3953 t8 = s*s; 3954 t10 = r-1.0; 3955 t12 = t10*t10; 3956 f = (t3*t2-56.0/27.0*t3*t1+196.0/117.0*t3*r-392.0/585.0*t3+49.0/351.0* 3957 t2-196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0)*(t8*s+12.0/11.0*t8*t10+18.0/ 3958 55.0*s*t12+4.0/165.0*t12*t10)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 3959 /* p=15 7 4 0 */ 3960 } else if( i==7 && j==4 && k==0 ) { 3961 t1 = r*r; 3962 t2 = t1*r; 3963 t3 = t1*t1; 3964 t8 = s*s; 3965 t9 = t8*t8; 3966 t11 = r-1.0; 3967 t13 = t11*t11; 3968 t17 = t13*t13; 3969 f = (t3*t2-56.0/27.0*t3*t1+196.0/117.0*t3*r-392.0/585.0*t3+49.0/351.0* 3970 t2-196.0/13455.0*t1+98.0/148005.0*r-4.0/444015.0)*(t9+20.0/11.0*t8*s*t11+12.0/ 3971 11.0*t8*t13+8.0/33.0*s*t13*t11+t17/66); 3972 /* p=15 8 0 3 */ 3973 } else if( i==8 && j==0 && k==3 ) { 3974 t1 = r*r; 3975 t2 = t1*t1; 3976 t3 = t2*t2; 3977 t4 = t1*r; 3978 t9 = t*t; 3979 t11 = s+r-1.0; 3980 t13 = t11*t11; 3981 f = (t3-8.0/3.0*t2*t4+112.0/39.0*t2*t1-1568.0/975.0*t2*r+98.0/195.0*t2 3982 -392.0/4485.0*t4+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0)*(t9*t+12.0/7.0* 3983 t9*t11+6.0/7.0*t*t13+4.0/35.0*t13*t11); 3984 /* p=15 8 1 2 */ 3985 } else if( i==8 && j==1 && k==2 ) { 3986 t1 = r*r; 3987 t2 = t1*t1; 3988 t3 = t2*t2; 3989 t4 = t1*r; 3990 t11 = t*t; 3991 t12 = s+r-1.0; 3992 t14 = t12*t12; 3993 f = (t3-8.0/3.0*t2*t4+112.0/39.0*t2*t1-1568.0/975.0*t2*r+98.0/195.0*t2 3994 -392.0/4485.0*t4+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0)*(s+2.0/9.0*r-2.0 3995 /9.0)*(t11+6.0/5.0*t*t12+3.0/10.0*t14); 3996 /* p=15 8 2 1 */ 3997 } else if( i==8 && j==2 && k==1 ) { 3998 t1 = r*r; 3999 t2 = t1*t1; 4000 t3 = t2*t2; 4001 t4 = t1*r; 4002 t9 = s*s; 4003 t10 = r-1.0; 4004 t12 = t10*t10; 4005 f = (t3-8.0/3.0*t2*t4+112.0/39.0*t2*t1-1568.0/975.0*t2*r+98.0/195.0*t2 4006 -392.0/4485.0*t4+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0)*(t9+2.0/3.0*s* 4007 t10+t12/12)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 4008 /* p=15 8 3 0 */ 4009 } else if( i==8 && j==3 && k==0 ) { 4010 t1 = r*r; 4011 t2 = t1*t1; 4012 t3 = t2*t2; 4013 t4 = t1*r; 4014 t9 = s*s; 4015 t11 = r-1.0; 4016 t13 = t11*t11; 4017 f = (t3-8.0/3.0*t2*t4+112.0/39.0*t2*t1-1568.0/975.0*t2*r+98.0/195.0*t2 4018 -392.0/4485.0*t4+392.0/49335.0*t1-16.0/49335.0*r+1.0/246675.0)*(t9*s+4.0/3.0*t9 4019 *t11+s*t13/2+t13*t11/21); 4020 /* p=15 9 0 2 */ 4021 } else if( i==9 && j==0 && k==2 ) { 4022 t1 = r*r; 4023 t2 = t1*t1; 4024 t3 = t2*t2; 4025 t5 = t1*r; 4026 t10 = t*t; 4027 t11 = s+r-1.0; 4028 t13 = t11*t11; 4029 f = (t3*r-10.0/3.0*t3+60.0/13.0*t2*t5-224.0/65.0*t2*t1+98.0/65.0*t2*r 4030 -588.0/1495.0*t2+196.0/3289.0*t5-16.0/3289.0*t1+3.0/16445.0*r-2.0/937365.0)*( 4031 t10+6.0/5.0*t*t11+3.0/10.0*t13); 4032 /* p=15 9 1 1 */ 4033 } else if( i==9 && j==1 && k==1 ) { 4034 t1 = r*r; 4035 t2 = t1*t1; 4036 t3 = t2*t2; 4037 t5 = t1*r; 4038 f = (t3*r-10.0/3.0*t3+60.0/13.0*t2*t5-224.0/65.0*t2*t1+98.0/65.0*t2*r 4039 -588.0/1495.0*t2+196.0/3289.0*t5-16.0/3289.0*t1+3.0/16445.0*r-2.0/937365.0)*(s+ 4040 2.0/7.0*r-2.0/7.0)*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 4041 /* p=15 9 2 0 */ 4042 } else if( i==9 && j==2 && k==0 ) { 4043 t1 = r*r; 4044 t2 = t1*t1; 4045 t3 = t2*t2; 4046 t5 = t1*r; 4047 t10 = s*s; 4048 t11 = r-1.0; 4049 t13 = t11*t11; 4050 f = (t3*r-10.0/3.0*t3+60.0/13.0*t2*t5-224.0/65.0*t2*t1+98.0/65.0*t2*r 4051 -588.0/1495.0*t2+196.0/3289.0*t5-16.0/3289.0*t1+3.0/16445.0*r-2.0/937365.0)*( 4052 t10+6.0/7.0*s*t11+t13/7); 4053 /* p=15 10 0 1 */ 4054 } else if( i==10 && j==0 && k==1 ) { 4055 t1 = r*r; 4056 t2 = t1*t1; 4057 t3 = t2*t2; 4058 t6 = t1*r; 4059 t10 = t3*t1-110.0/27.0*t3*r+275.0/39.0*t3-88.0/13.0*t2*t6+154.0/39.0*t2* 4060 t1-2156.0/1495.0*t2*r+98.0/299.0*t2-40.0/897.0*t6+t1/299-2.0/17043.0*r+1.0/ 4061 766935.0; 4062 f = t10*(t+2.0/3.0*s+2.0/3.0*r-2.0/3.0); 4063 /* p=15 10 1 0 */ 4064 } else if( i==10 && j==1 && k==0 ) { 4065 t1 = r*r; 4066 t2 = t1*t1; 4067 t3 = t2*t2; 4068 t6 = t1*r; 4069 t10 = t3*t1-110.0/27.0*t3*r+275.0/39.0*t3-88.0/13.0*t2*t6+154.0/39.0*t2* 4070 t1-2156.0/1495.0*t2*r+98.0/299.0*t2-40.0/897.0*t6+t1/299-2.0/17043.0*r+1.0/ 4071 766935.0; 4072 f = t10*(s+2.0/5.0*r-2.0/5.0); 4073 /* p=15 11 0 0 */ 4074 } else if( i==11 && j==0 && k==0 ) { 4075 t1 = r*r; 4076 t2 = t1*r; 4077 t3 = t1*t1; 4078 t4 = t3*t3; 4079 f = t4*t2-44.0/9.0*t4*t1+1210.0/117.0*t4*r-484.0/39.0*t4+121.0/13.0*t3* 4080 t2-6776.0/1495.0*t3*t1+2156.0/1495.0*t3*r-88.0/299.0*t3+11.0/299.0*t2-44.0/ 4081 17043.0*t1+22.0/255645.0*r-4.0/4345965.0; 4082 } 4083 return f ; 4084 } 4085 4086 #ifdef __cplusplus 4087 } 4088 #endif 4089