1 Subroutine hfmkr(R,IJK,PQ,NAB,Lp,Lp3) 2c $Id$ 3 4 Implicit real*8 (a-h,o-z) 5 Implicit integer (i-n) 6 7 Dimension R(NAB,0:Lp,Lp3),IJK(0:Lp,0:Lp,0:Lp),PQ(NAB,3) 8 9 Dimension Nxyz(3) 10c 11c Recursively define the auxiliary integrals needed for ERIs and NAIs. 12c 13c Recursion Formula: 14c 15c R = PQx R + (N-1) R 16c N,L,M,j N-1,L,M,j+1 N-2,L,M,j+1 17c 18c R = PQy R + (L-1) R 19c N,L,M,j N,L-1,M,j+1 N,L-2,M,j+1 20c 21c R = PQz R + (M-1) R 22c N,L,M,j N,L,M-1,j+1 N,L,M-2,j+1 23c 24c****************************************************************************** 25 26#if defined(USE_SPECIAL_MAKER) 27 if (Lp.le.1) then 28 call hfmkr_1(R,PQ,IJK,NAB,Lp,Lp3) 29 goto 00099 30 elseif(Lp.eq.2) then 31 call hfmkr_2(R,PQ,IJK,NAB,Lp,Lp3) 32 goto 00099 33 elseif(Lp.eq.3) then 34 call hfmkr_3(R,PQ,IJK,NAB,Lp,Lp3) 35 goto 00099 36 elseif(Lp.eq.4) then 37 call hfmkr_4(R,PQ,IJK,NAB,Lp,Lp3) 38 goto 00099 39 elseif(Lp.eq.5) then 40 call hfmkr_5(R,PQ,IJK,NAB,Lp,Lp3) 41 goto 00099 42 elseif(Lp.eq.6) then 43 call hfmkr_6(R,PQ,IJK,NAB,Lp,Lp3) 44 goto 00099 45 elseif(Lp.eq.7) then 46 call hfmkr_7(R,PQ,IJK,NAB,Lp,Lp3) 47 goto 00099 48 elseif(Lp.eq.8) then 49 call hfmkr_8(R,PQ,IJK,NAB,Lp,Lp3) 50 goto 00099 51 elseif (Lp.eq.9) then 52 call hfmkr_9(R,PQ,IJK,NAB,Lp,Lp3) 53 goto 00099 54 endif 55#endif 56 57#if defined(INT_DEBUG) 58 write(6,*)'lp',lp,' lp3',lp3,' nab',nab 59 write(6,*)' PQ coords start' 60 call intintp(PQ,nab,3,'hfmkr pq') 61 write(6,*)' PQ coords end' 62 write(6,*)' r b4 recursion' 63 do ii = 1,NAB 64 do jj = 0,Lp 65 do kk = 1,lp3 66 if (abs(R(ii,jj,kk)).gt.1.0d-12) then 67 write(6,*)'R (',ii,jj,kk,') =',R(ii,jj,kk) 68 endif 69 enddo 70 enddo 71 enddo 72#endif 73 74 n0 = 1 75 76 IJK(0,0,0) = n0 77 78 do 180 l = 1,Lp 79 80 l2 = ((l+1)*(l+2))/2 81 82 do 170 i = 1,l2 83 84 call getNxyz(l,i,Nxyz) 85 86 Ip = Nxyz(1) 87 Jp = Nxyz(2) 88 Kp = Nxyz(3) 89 90 n0 = n0 + 1 91 92 IJK(Ip,Jp,Kp) = n0 93 94 if( Ip.gt.0 )then 95 96 if( Ip.eq.1 )then 97 98 n1 = IJK(Ip-1,Jp,Kp) 99 100 do 115 j = 0,Lp-l 101 do 110 m = 1,NAB 102 R(m,j,n0) = PQ(m,1)*R(m,j+1,n1) 103 110 continue 104 115 continue 105 106 else 107 108 n1 = IJK(Ip-1,Jp,Kp) 109 n2 = IJK(Ip-2,Jp,Kp) 110 111 do 125 j = 0,Lp-l 112 do 120 m = 1,NAB 113 R(m,j,n0) = PQ(m,1)*R(m,j+1,n1) + (Ip-1)*R(m,j+1,n2) 114 120 continue 115 125 continue 116 117 end if 118 119 else if( Jp.gt.0 )then 120 121 if( Jp.eq.1 )then 122 123 n1 = IJK(Ip,Jp-1,Kp) 124 125 do 135 j = 0,Lp-l 126 do 130 m = 1,NAB 127 R(m,j,n0) = PQ(m,2)*R(m,j+1,n1) 128 130 continue 129 135 continue 130 131 else 132 133 n1 = IJK(Ip,Jp-1,Kp) 134 n2 = IJK(Ip,Jp-2,Kp) 135 136 do 145 j = 0,Lp-l 137 do 140 m = 1,NAB 138 R(m,j,n0) = PQ(m,2)*R(m,j+1,n1) + (Jp-1)*R(m,j+1,n2) 139 140 continue 140 145 continue 141 142 end if 143 144 else if( Kp.gt.0 )then 145 146 if( Kp.eq.1 )then 147 148 n1 = IJK(Ip,Jp,Kp-1) 149 150 do 155 j = 0,Lp-l 151 do 150 m = 1,NAB 152 R(m,j,n0) = PQ(m,3)*R(m,j+1,n1) 153 150 continue 154 155 continue 155 156 else 157 158 n1 = IJK(Ip,Jp,Kp-1) 159 n2 = IJK(Ip,Jp,Kp-2) 160 161 do 165 j = 0,Lp-l 162 do 160 m = 1,NAB 163 R(m,j,n0) = PQ(m,3)*R(m,j+1,n1) + (Kp-1)*R(m,j+1,n2) 164 160 continue 165 165 continue 166 167 end if 168 169 end if 170 171 170 continue 172 173 180 continue 174 17500099 continue 176 177#if defined(INT_DEBUG) 178 write(6,*)' r after recursion' 179 do ii = 1,NAB 180 do jj = 0,Lp 181 do kk = 1,lp3 182 if (abs(R(ii,jj,kk)).gt.1.0d-12) then 183 write(6,*)'R (',ii,jj,kk,') =',R(ii,jj,kk) 184 endif 185 enddo 186 enddo 187 enddo 188#endif 189 end 190 191 192#if defined(USE_SPECIAL_MAKER) 193 subroutine hfmkr_1(R,PQ,IJK,NAB,Lr,Lr3) 194c $Id$ 195 implicit none 196c 197 integer NAB, Lr, Lr3 198 integer IJK(0:Lr,0:Lr,0:Lr) 199 double precision PQ(NAB,3) 200 double precision R(NAB,0:Lr,Lr3) 201c 202 integer m 203 double precision x,y,z 204c 205 IJK(0,0,0)=1 206 if(Lr.eq.0) goto 00001 207 IJK(1,0,0)=2 208 IJK(0,1,0)=3 209 IJK(0,0,1)=4 210c 21100001 continue 212 if(Lr.eq.0) goto 00002 213 do 00100 m=1,NAB 214 X=PQ(m,1) 215 Y=PQ(m,2) 216 Z=PQ(m,3) 217 R(M,0,2)=X*R(m,1,1) 218 R(M,0,3)=Y*R(m,1,1) 219 R(M,0,4)=Z*R(m,1,1) 22000100 continue 221c 22200002 continue 223c 224 end 225 subroutine hfmkr_2(R,PQ,IJK,NAB,Lr,Lr3) 226 implicit none 227c $Id$ 228c 229 integer NAB, Lr, Lr3 230 integer IJK(0:Lr,0:Lr,0:Lr) 231 double precision PQ(NAB,3) 232 double precision R(NAB,0:Lr,Lr3) 233c 234 integer m 235 double precision x,y,z 236 double precision R21, R11 237c 238 IJK(0,0,0)=1 239 IJK(1,0,0)=2 240 IJK(0,1,0)=3 241 IJK(0,0,1)=4 242 IJK(2,0,0)=5 243 IJK(1,1,0)=6 244 IJK(1,0,1)=7 245 IJK(0,2,0)=8 246 IJK(0,1,1)=9 247 IJK(0,0,2)=10 248c 249 do 00100 m=1,NAB 250 X=PQ(m,1) 251 Y=PQ(m,2) 252 Z=PQ(m,3) 253 R11=R(m,1,1) 254 R21=R(m,2,1) 255 R(M,0,2) =X*R11 256 R(M,0,3) =Y*R11 257 R(M,0,4) =Z*R11 258 R(M,0,5) =X*X*R21+R11 259 R(M,0,6) =X*Y*R21 260 R(M,0,7) =X*Z*R21 261 R(M,0,8) =Y*Y*R21+R11 262 R(M,0,9) =Y*Z*R21 263 R(M,0,10)=Z*Z*R21+R11 26400100 continue 265c 266 end 267 subroutine hfmkr_3(R,PQ,IJK,NAB,Lr,Lr3) 268 implicit none 269c $Id$ 270c 271 integer NAB, Lr, Lr3 272 integer IJK(0:Lr,0:Lr,0:Lr) 273 double precision PQ(NAB,3) 274 double precision R(NAB,0:Lr,Lr3) 275c 276 integer m 277 double precision x,y,z 278 double precision x_2, y_2, z_2 279 double precision x_3, y_3, z_3 280 double precision R11,R21,R31 281c 282 IJK(0,0,0)=1 283 IJK(1,0,0)=2 284 IJK(0,1,0)=3 285 IJK(0,0,1)=4 286 IJK(2,0,0)=5 287 IJK(1,1,0)=6 288 IJK(1,0,1)=7 289 IJK(0,2,0)=8 290 IJK(0,1,1)=9 291 IJK(0,0,2)=10 292 IJK(3,0,0)=11 293 IJK(2,1,0)=12 294 IJK(2,0,1)=13 295 IJK(1,2,0)=14 296 IJK(1,1,1)=15 297 IJK(1,0,2)=16 298 IJK(0,3,0)=17 299 IJK(0,2,1)=18 300 IJK(0,1,2)=19 301 IJK(0,0,3)=20 302c 303 do 00100 m=1,NAB 304 X=PQ(m,1) 305 Y=PQ(m,2) 306 Z=PQ(m,3) 307 X_2=X*X 308 Y_2=Y*Y 309 Z_2=Z*Z 310 X_3=X_2*X 311 Y_3=Y_2*Y 312 Z_3=Z_2*Z 313 R11=R(m,1,1) 314 R21=R(m,2,1) 315 R31=R(m,3,1) 316 R(M,0,2)=X*R11 317 R(M,0,3)=Y*R11 318 R(M,0,4)=Z*R11 319 R(M,0,5)=X_2*R21+R11 320 R(M,0,6)=X*Y*R21 321 R(M,0,7)=X*Z*R21 322 R(M,0,8)=Y_2*R21+R11 323 R(M,0,9)=Y*Z*R21 324 R(M,0,10)=Z_2*R21+R11 325 R(M,0,11)=X_3*R31+3.0d00*X*R21 326 R(M,0,12)=X_2*Y*R31+Y*R21 327 R(M,0,13)=X_2*Z*R31+Z*R21 328 R(M,0,14)=X*Y_2*R31+X*R21 329 R(M,0,15)=X*Y*Z*R31 330 R(M,0,16)=X*Z_2*R31+X*R21 331 R(M,0,17)=Y_3*R31+3.0d00*Y*R21 332 R(M,0,18)=Y_2*Z*R31+Z*R21 333 R(M,0,19)=Y*Z_2*R31+Y*R21 334 R(M,0,20)=Z_3*R31+3.0d00*Z*R21 33500100 continue 336c 337 end 338 subroutine hfmkr_4(R,PQ,IJK,NAB,Lr,Lr3) 339 implicit none 340c $Id$ 341c 342 integer NAB, Lr, Lr3 343 integer IJK(0:Lr,0:Lr,0:Lr) 344 double precision PQ(NAB,3) 345 double precision R(NAB,0:Lr,Lr3) 346c 347 integer m 348 double precision x,y,z 349 double precision x_2, y_2, z_2 350 double precision x_3, y_3, z_3 351 double precision x_4, y_4, z_4 352 double precision R11,R21,R31,R41 353c 354 IJK(0,0,0)=1 355 IJK(1,0,0)=2 356 IJK(0,1,0)=3 357 IJK(0,0,1)=4 358 IJK(2,0,0)=5 359 IJK(1,1,0)=6 360 IJK(1,0,1)=7 361 IJK(0,2,0)=8 362 IJK(0,1,1)=9 363 IJK(0,0,2)=10 364 IJK(3,0,0)=11 365 IJK(2,1,0)=12 366 IJK(2,0,1)=13 367 IJK(1,2,0)=14 368 IJK(1,1,1)=15 369 IJK(1,0,2)=16 370 IJK(0,3,0)=17 371 IJK(0,2,1)=18 372 IJK(0,1,2)=19 373 IJK(0,0,3)=20 374 IJK(4,0,0)=21 375 IJK(3,1,0)=22 376 IJK(3,0,1)=23 377 IJK(2,2,0)=24 378 IJK(2,1,1)=25 379 IJK(2,0,2)=26 380 IJK(1,3,0)=27 381 IJK(1,2,1)=28 382 IJK(1,1,2)=29 383 IJK(1,0,3)=30 384 IJK(0,4,0)=31 385 IJK(0,3,1)=32 386 IJK(0,2,2)=33 387 IJK(0,1,3)=34 388 IJK(0,0,4)=35 389c 390 do 00100 m=1,NAB 391 X=PQ(m,1) 392 Y=PQ(m,2) 393 Z=PQ(m,3) 394 X_2=X*X 395 Y_2=Y*Y 396 Z_2=Z*Z 397 X_3=X_2*X 398 Y_3=Y_2*Y 399 Z_3=Z_2*Z 400 X_4=X_3*X 401 Y_4=Y_3*Y 402 Z_4=Z_3*Z 403c 404 R11=R(m,1,1) 405 R21=R(m,2,1) 406 R31=R(m,3,1) 407 R41=R(m,4,1) 408c 409 R(M,0,2)=X*R11 410 R(M,0,3)=Y*R11 411 R(M,0,4)=Z*R11 412 R(M,0,5)=X_2*R21+R11 413 R(M,0,6)=X*Y*R21 414 R(M,0,7)=X*Z*R21 415 R(M,0,8)=Y_2*R21+R11 416 R(M,0,9)=Y*Z*R21 417 R(M,0,10)=Z_2*R21+R11 418 R(M,0,11)=X_3*R31+3.0d00*X*R21 419 R(M,0,12)=X_2*Y*R31+Y*R21 420 R(M,0,13)=X_2*Z*R31+Z*R21 421 R(M,0,14)=X*Y_2*R31+X*R21 422 R(M,0,15)=X*Y*Z*R31 423 R(M,0,16)=X*Z_2*R31+X*R21 424 R(M,0,17)=Y_3*R31+3.0d00*Y*R21 425 R(M,0,18)=Y_2*Z*R31+Z*R21 426 R(M,0,19)=Y*Z_2*R31+Y*R21 427 R(M,0,20)=Z_3*R31+3.0d00*Z*R21 428 R(M,0,21)=X_4*R41+6.0d00*X_2*R31+3.0d00*R21 429 R(M,0,22)=X_3*Y*R41+3.0d00*X*Y*R31 430 R(M,0,23)=X_3*Z*R41+3.0d00*X*Z*R31 431 R(M,0,24)=X_2*Y_2*R41+(X_2+Y_2)*R31+R21 432 R(M,0,25)=X_2*Y*Z*R41+Y*Z*R31 433 R(M,0,26)=X_2*Z_2*R41+(X_2+Z_2)*R31+R21 434 R(M,0,27)=X*Y_3*R41+3.0d00*X*Y*R31 435 R(M,0,28)=X*Y_2*Z*R41+X*Z*R31 436 R(M,0,29)=X*Y*Z_2*R41+X*Y*R31 437 R(M,0,30)=X*Z_3*R41+3.0d00*X*Z*R31 438 R(M,0,31)=Y_4*R41+6.0d00*Y_2*R31+3.0d00*R21 439 R(M,0,32)=Y_3*Z*R41+3.0d00*Y*Z*R31 440 R(M,0,33)=Y_2*Z_2*R41+(Y_2+Z_2)*R31+R21 441 R(M,0,34)=Y*Z_3*R41+3.0d00*Y*Z*R31 442 R(M,0,35)=Z_4*R41+6.0d00*Z_2*R31+3.0d00*R21 44300100 continue 444c 445 end 446 subroutine hfmkr_5(R,PQ,IJK,NAB,Lr,Lr3) 447 implicit none 448c $Id$ 449c 450 integer NAB, Lr, Lr3 451 integer IJK(0:Lr,0:Lr,0:Lr) 452 double precision PQ(NAB,3) 453 double precision R(NAB,0:Lr,Lr3) 454c 455 integer m 456 double precision x,y,z 457 double precision x_2, y_2, z_2 458 double precision x_3, y_3, z_3 459 double precision x_4, y_4, z_4 460 double precision x_5, y_5, z_5 461 double precision R11,R21,R31,R41,R51 462c 463 IJK(0,0,0)=1 464 IJK(1,0,0)=2 465 IJK(0,1,0)=3 466 IJK(0,0,1)=4 467 IJK(2,0,0)=5 468 IJK(1,1,0)=6 469 IJK(1,0,1)=7 470 IJK(0,2,0)=8 471 IJK(0,1,1)=9 472 IJK(0,0,2)=10 473 IJK(3,0,0)=11 474 IJK(2,1,0)=12 475 IJK(2,0,1)=13 476 IJK(1,2,0)=14 477 IJK(1,1,1)=15 478 IJK(1,0,2)=16 479 IJK(0,3,0)=17 480 IJK(0,2,1)=18 481 IJK(0,1,2)=19 482 IJK(0,0,3)=20 483 IJK(4,0,0)=21 484 IJK(3,1,0)=22 485 IJK(3,0,1)=23 486 IJK(2,2,0)=24 487 IJK(2,1,1)=25 488 IJK(2,0,2)=26 489 IJK(1,3,0)=27 490 IJK(1,2,1)=28 491 IJK(1,1,2)=29 492 IJK(1,0,3)=30 493 IJK(0,4,0)=31 494 IJK(0,3,1)=32 495 IJK(0,2,2)=33 496 IJK(0,1,3)=34 497 IJK(0,0,4)=35 498 IJK(5,0,0)=36 499 IJK(4,1,0)=37 500 IJK(4,0,1)=38 501 IJK(3,2,0)=39 502 IJK(3,1,1)=40 503 IJK(3,0,2)=41 504 IJK(2,3,0)=42 505 IJK(2,2,1)=43 506 IJK(2,1,2)=44 507 IJK(2,0,3)=45 508 IJK(1,4,0)=46 509 IJK(1,3,1)=47 510 IJK(1,2,2)=48 511 IJK(1,1,3)=49 512 IJK(1,0,4)=50 513 IJK(0,5,0)=51 514 IJK(0,4,1)=52 515 IJK(0,3,2)=53 516 IJK(0,2,3)=54 517 IJK(0,1,4)=55 518 IJK(0,0,5)=56 519c 520 do 00100 m=1,NAB 521 X=PQ(m,1) 522 Y=PQ(m,2) 523 Z=PQ(m,3) 524 X_2=X*X 525 Y_2=Y*Y 526 Z_2=Z*Z 527 X_3=X_2*X 528 Y_3=Y_2*Y 529 Z_3=Z_2*Z 530 X_4=X_3*X 531 Y_4=Y_3*Y 532 Z_4=Z_3*Z 533 X_5=X_4*X 534 Y_5=Y_4*Y 535 Z_5=Z_4*Z 536c 537 R11=R(m,1,1) 538 R21=R(m,2,1) 539 R31=R(m,3,1) 540 R41=R(m,4,1) 541 R51=R(m,5,1) 542c 543 R(M,0,2)=X*R11 544 R(M,0,3)=Y*R11 545 R(M,0,4)=Z*R11 546 R(M,0,5)=X_2*R21+R11 547 R(M,0,6)=X*Y*R21 548 R(M,0,7)=X*Z*R21 549 R(M,0,8)=Y_2*R21+R11 550 R(M,0,9)=Y*Z*R21 551 R(M,0,10)=Z_2*R21+R11 552 R(M,0,11)=X_3*R31+3.0d00*X*R21 553 R(M,0,12)=X_2*Y*R31+Y*R21 554 R(M,0,13)=X_2*Z*R31+Z*R21 555 R(M,0,14)=X*Y_2*R31+X*R21 556 R(M,0,15)=X*Y*Z*R31 557 R(M,0,16)=X*Z_2*R31+X*R21 558 R(M,0,17)=Y_3*R31+3.0d00*Y*R21 559 R(M,0,18)=Y_2*Z*R31+Z*R21 560 R(M,0,19)=Y*Z_2*R31+Y*R21 561 R(M,0,20)=Z_3*R31+3.0d00*Z*R21 562 R(M,0,21)=X_4*R41+6.0d00*X_2*R31+3.0d00*R21 563 R(M,0,22)=X_3*Y*R41+3.0d00*X*Y*R31 564 R(M,0,23)=X_3*Z*R41+3.0d00*X*Z*R31 565 R(M,0,24)=X_2*Y_2*R41+(X_2+Y_2)*R31+R21 566 R(M,0,25)=X_2*Y*Z*R41+Y*Z*R31 567 R(M,0,26)=X_2*Z_2*R41+(X_2+Z_2)*R31+R21 568 R(M,0,27)=X*Y_3*R41+3.0d00*X*Y*R31 569 R(M,0,28)=X*Y_2*Z*R41+X*Z*R31 570 R(M,0,29)=X*Y*Z_2*R41+X*Y*R31 571 R(M,0,30)=X*Z_3*R41+3.0d00*X*Z*R31 572 R(M,0,31)=Y_4*R41+6.0d00*Y_2*R31+3.0d00*R21 573 R(M,0,32)=Y_3*Z*R41+3.0d00*Y*Z*R31 574 R(M,0,33)=Y_2*Z_2*R41+(Y_2+Z_2)*R31+R21 575 R(M,0,34)=Y*Z_3*R41+3.0d00*Y*Z*R31 576 R(M,0,35)=Z_4*R41+6.0d00*Z_2*R31+3.0d00*R21 577 R(M,0,36)=X_5*R51+10.0d00* X_3*R41+15.0d00*X*R31 578 R(M,0,37)=X_4*Y*R51+6.0d00*X_2*Y*R41+3.0d00*Y*R31 579 R(M,0,38)=X_4*Z*R51+6.0d00*X_2*Z*R41+3.0d00*Z*R31 580 R(M,0,39)=X_3*Y_2*R51+(X*(X_2+Y_2)+2.0d00*X*Y_2)*R41 581 & +3.0d00*X*R31 582 R(M,0,40)=X_3*Y*Z*R51+3.0d00*X*Y*Z*R41 583 R(M,0,41)=X_3*Z_2*R51+(X*(X_2+Z_2)+2.0d00*X*Z_2)*R41 584 & +3.0d00*X*R31 585 R(M,0,42)=X_2*Y_3*R51+(3.0d00*X_2*Y+Y_3)*R41 586 & +3.0d00*Y*R31 587 R(M,0,43)=X_2*Y_2*Z*R51+(X_2*Z+Y_2*Z)*R41 588 & +Z*R31 589 R(M,0,44)=X_2*Y*Z_2*R51+(X_2*Y+Y*Z_2)*R41 590 & +Y*R31 591 R(M,0,45)=X_2*Z_3*R51+(3.0d00*X_2*Z+Z_3)*R41 592 & +3.0d00*Z*R31 593 R(M,0,46)=X*Y_4*R51+6.0d00*X*Y_2*R41+3.0d00*X*R31 594 R(M,0,47)=X*Y_3*Z*R51+3.0d00*X*Y*Z*R41 595 R(M,0,48)=X*Y_2*Z_2*R51+X*(Y_2+Z_2)*R41+ 596 & X*R31 597 R(M,0,49)=X*Y*Z_3*R51+3.0d00*X*Y*Z*R41 598 R(M,0,50)=X*Z_4*R51+6.0d00*X*Z_2*R41+3.0d00*X*R31 599 R(M,0,51)=Y_5*R51+10.0d00* Y_3*R41+15.0d00*Y*R31 600 R(M,0,52)=Y_4*Z*R51+6.0d00*Y_2*Z*R41+3.0d00*Z*R31 601 R(M,0,53)=Y_3*Z_2*R51+(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R41 602 & +3.0d00*Y*R31 603 R(M,0,54)=Y_2*Z_3*R51+(3.0d00*Y_2*Z+Z_3)*R41 604 & +3.0d00*Z*R31 605 R(M,0,55)=Y*Z_4*R51+6.0d00*Y*Z_2*R41+3.0d00*Y*R31 606 R(M,0,56)=Z_5*R51+10.0d00* Z_3*R41+15.0d00*Z*R31 60700100 continue 608c 609 end 610 subroutine hfmkr_6(R,PQ,IJK,NAB,Lr,Lr3) 611 implicit none 612c $Id$ 613c 614 integer NAB, Lr, Lr3 615 integer IJK(0:Lr,0:Lr,0:Lr) 616 double precision PQ(NAB,3) 617 double precision R(NAB,0:Lr,Lr3) 618c 619 integer m 620 double precision x,y,z 621 double precision x_2, y_2, z_2 622 double precision x_3, y_3, z_3 623 double precision x_4, y_4, z_4 624 double precision x_5, y_5, z_5 625 double precision x_6, y_6, z_6 626 double precision R11,R21,R31,R41,R51,R61 627c 628 IJK(0,0,0)=1 629 IJK(1,0,0)=2 630 IJK(0,1,0)=3 631 IJK(0,0,1)=4 632 IJK(2,0,0)=5 633 IJK(1,1,0)=6 634 IJK(1,0,1)=7 635 IJK(0,2,0)=8 636 IJK(0,1,1)=9 637 IJK(0,0,2)=10 638 IJK(3,0,0)=11 639 IJK(2,1,0)=12 640 IJK(2,0,1)=13 641 IJK(1,2,0)=14 642 IJK(1,1,1)=15 643 IJK(1,0,2)=16 644 IJK(0,3,0)=17 645 IJK(0,2,1)=18 646 IJK(0,1,2)=19 647 IJK(0,0,3)=20 648 IJK(4,0,0)=21 649 IJK(3,1,0)=22 650 IJK(3,0,1)=23 651 IJK(2,2,0)=24 652 IJK(2,1,1)=25 653 IJK(2,0,2)=26 654 IJK(1,3,0)=27 655 IJK(1,2,1)=28 656 IJK(1,1,2)=29 657 IJK(1,0,3)=30 658 IJK(0,4,0)=31 659 IJK(0,3,1)=32 660 IJK(0,2,2)=33 661 IJK(0,1,3)=34 662 IJK(0,0,4)=35 663 IJK(5,0,0)=36 664 IJK(4,1,0)=37 665 IJK(4,0,1)=38 666 IJK(3,2,0)=39 667 IJK(3,1,1)=40 668 IJK(3,0,2)=41 669 IJK(2,3,0)=42 670 IJK(2,2,1)=43 671 IJK(2,1,2)=44 672 IJK(2,0,3)=45 673 IJK(1,4,0)=46 674 IJK(1,3,1)=47 675 IJK(1,2,2)=48 676 IJK(1,1,3)=49 677 IJK(1,0,4)=50 678 IJK(0,5,0)=51 679 IJK(0,4,1)=52 680 IJK(0,3,2)=53 681 IJK(0,2,3)=54 682 IJK(0,1,4)=55 683 IJK(0,0,5)=56 684 IJK(6,0,0)=57 685 IJK(5,1,0)=58 686 IJK(5,0,1)=59 687 IJK(4,2,0)=60 688 IJK(4,1,1)=61 689 IJK(4,0,2)=62 690 IJK(3,3,0)=63 691 IJK(3,2,1)=64 692 IJK(3,1,2)=65 693 IJK(3,0,3)=66 694 IJK(2,4,0)=67 695 IJK(2,3,1)=68 696 IJK(2,2,2)=69 697 IJK(2,1,3)=70 698 IJK(2,0,4)=71 699 IJK(1,5,0)=72 700 IJK(1,4,1)=73 701 IJK(1,3,2)=74 702 IJK(1,2,3)=75 703 IJK(1,1,4)=76 704 IJK(1,0,5)=77 705 IJK(0,6,0)=78 706 IJK(0,5,1)=79 707 IJK(0,4,2)=80 708 IJK(0,3,3)=81 709 IJK(0,2,4)=82 710 IJK(0,1,5)=83 711 IJK(0,0,6)=84 712c 713 do 00100 m=1,NAB 714 X=PQ(m,1) 715 Y=PQ(m,2) 716 Z=PQ(m,3) 717 R11=R(m,1,1) 718 R(M,0,2)=X*R11 719 R(M,0,3)=Y*R11 720 R(M,0,4)=Z*R11 721 X_2=X*X 722 Y_2=Y*Y 723 Z_2=Z*Z 724 R21=R(m,2,1) 725 R(M,0,5)=X_2*R21+R11 726 R(M,0,6)=X*Y*R21 727 R(M,0,7)=X*Z*R21 728 R(M,0,8)=Y_2*R21+R11 729 R(M,0,9)=Y*Z*R21 730 R(M,0,10)=Z_2*R21+R11 731 X_3=X_2*X 732 Y_3=Y_2*Y 733 Z_3=Z_2*Z 734 R31=R(m,3,1) 735 R(M,0,11)=X_3*R31+3.0d00*X*R21 736 R(M,0,12)=X_2*Y*R31+Y*R21 737 R(M,0,13)=X_2*Z*R31+Z*R21 738 R(M,0,14)=X*Y_2*R31+X*R21 739 R(M,0,15)=X*Y*Z*R31 740 R(M,0,16)=X*Z_2*R31+X*R21 741 R(M,0,17)=Y_3*R31+3.0d00*Y*R21 742 R(M,0,18)=Y_2*Z*R31+Z*R21 743 R(M,0,19)=Y*Z_2*R31+Y*R21 744 R(M,0,20)=Z_3*R31+3.0d00*Z*R21 745 X_4=X_3*X 746 Y_4=Y_3*Y 747 Z_4=Z_3*Z 748 R41=R(m,4,1) 749 R(M,0,21)=X_4*R41+6.0d00*X_2*R31+3.0d00*R21 750 R(M,0,22)=X_3*Y*R41+3.0d00*X*Y*R31 751 R(M,0,23)=X_3*Z*R41+3.0d00*X*Z*R31 752 R(M,0,24)=X_2*Y_2*R41+(X_2+Y_2)*R31+R21 753 R(M,0,25)=X_2*Y*Z*R41+Y*Z*R31 754 R(M,0,26)=X_2*Z_2*R41+(X_2+Z_2)*R31+R21 755 R(M,0,27)=X*Y_3*R41+3.0d00*X*Y*R31 756 R(M,0,28)=X*Y_2*Z*R41+X*Z*R31 757 R(M,0,29)=X*Y*Z_2*R41+X*Y*R31 758 R(M,0,30)=X*Z_3*R41+3.0d00*X*Z*R31 759 R(M,0,31)=Y_4*R41+6.0d00*Y_2*R31+3.0d00*R21 760 R(M,0,32)=Y_3*Z*R41+3.0d00*Y*Z*R31 761 R(M,0,33)=Y_2*Z_2*R41+(Y_2+Z_2)*R31+R21 762 R(M,0,34)=Y*Z_3*R41+3.0d00*Y*Z*R31 763 R(M,0,35)=Z_4*R41+6.0d00*Z_2*R31+3.0d00*R21 764 X_5=X_4*X 765 Y_5=Y_4*Y 766 Z_5=Z_4*Z 767 R51=R(m,5,1) 768 R(M,0,36)=X_5*R51+10.0d00* X_3*R41+15.0d00*X*R31 769 R(M,0,37)=X_4*Y*R51+6.0d00*X_2*Y*R41+3.0d00*Y*R31 770 R(M,0,38)=X_4*Z*R51+6.0d00*X_2*Z*R41+3.0d00*Z*R31 771 R(M,0,39)=X_3*Y_2*R51+(X*(X_2+Y_2)+2.0d00*X*Y_2)*R41 772 & +3.0d00*X*R31 773 R(M,0,40)=X_3*Y*Z*R51+3.0d00*X*Y*Z*R41 774 R(M,0,41)=X_3*Z_2*R51+(X*(X_2+Z_2)+2.0d00*X*Z_2)*R41 775 & +3.0d00*X*R31 776 R(M,0,42)=X_2*Y_3*R51+(3.0d00*X_2*Y+Y_3)*R41 777 & +3.0d00*Y*R31 778 R(M,0,43)=X_2*Y_2*Z*R51+(X_2*Z+Y_2*Z)*R41 779 & +Z*R31 780 R(M,0,44)=X_2*Y*Z_2*R51+(X_2*Y+Y*Z_2)*R41 781 & +Y*R31 782 R(M,0,45)=X_2*Z_3*R51+(3.0d00*X_2*Z+Z_3)*R41 783 & +3.0d00*Z*R31 784 R(M,0,46)=X*Y_4*R51+6.0d00*X*Y_2*R41+3.0d00*X*R31 785 R(M,0,47)=X*Y_3*Z*R51+3.0d00*X*Y*Z*R41 786 R(M,0,48)=X*Y_2*Z_2*R51+X*(Y_2+Z_2)*R41+ 787 & X*R31 788 R(M,0,49)=X*Y*Z_3*R51+3.0d00*X*Y*Z*R41 789 R(M,0,50)=X*Z_4*R51+6.0d00*X*Z_2*R41+3.0d00*X*R31 790 R(M,0,51)=Y_5*R51+10.0d00* Y_3*R41+15.0d00*Y*R31 791 R(M,0,52)=Y_4*Z*R51+6.0d00*Y_2*Z*R41+3.0d00*Z*R31 792 R(M,0,53)=Y_3*Z_2*R51+(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R41 793 & +3.0d00*Y*R31 794 R(M,0,54)=Y_2*Z_3*R51+(3.0d00*Y_2*Z+Z_3)*R41 795 & +3.0d00*Z*R31 796 R(M,0,55)=Y*Z_4*R51+6.0d00*Y*Z_2*R41+3.0d00*Y*R31 797 R(M,0,56)=Z_5*R51+10.0d00* Z_3*R41+15.0d00*Z*R31 798 X_6=X_5 * X 799 Y_6=Y_5 * Y 800 Z_6=Z_5 * Z 801 R61=R(m,6,1) 802 R(M,0,57)=X_6*R61+15.0d00*X_4*R51+45.0d00*X_2*R41 803 & +15.0d00*R31 804 R(M,0,58)=X_5*Y*R61+10.0d00* X_3*Y*R51+15.0d00*X*Y*R41 805 R(M,0,59)=X_5*Z*R61+10.0d00* X_3*Z*R51+15.0d00*X*Z*R41 806 R(M,0,60)=X_4*Y_2*R61 807 & +(3.0d00*X_2*Y_2+X*(X*(X_2+Y_2)+2.0d00*X*Y_2))*R51 808 & +(6.0d00*X_2+3.0d00*Y_2)*R41+3.0d00*R31 809 R(M,0,61)=X_4*Y*Z*R61+6.0d00*X_2*Y*Z*R51+3.0d00*Y*Z*R41 810 R(M,0,62)=X_4*Z_2*R61 811 & +(3.0d00*X_2*Z_2+X*(X*(X_2+Z_2)+2.0d00*X*Z_2))*R51 812 & +(6.0d00*X_2+3.0d00*Z_2)*R41+3.0d00*R31 813 R(M,0,63)=X_3*Y_3*R61 814 & +(X*(3.0d00*X_2*Y+Y_3)+2.0d00*X*Y_3)*R51+9.0d00*X*Y*R41 815 R(M,0,64)=X_3*Y_2*Z*R61 816 & +(X*(X_2*Z+Y_2*Z)+2.0d00*X*Y_2*Z)*R51+3.0d00*X*Z*R41 817 R(M,0,65)=X_3*Y*Z_2*R61 818 & +(X*(X_2*Y+Y*Z_2)+2.0d00*X*Y*Z_2)*R51+3.0d00*X*Y*R41 819 R(M,0,66)=X_3*Z_3*R61 820 & +(X*(3.0d00*X_2*Z+Z_3)+2.0d00*X*Z_3)*R51+9.0d00*X*Z*R41 821 R(M,0,67)=X_2*Y_4*R61+(Y_4+6.0d00*X_2*Y_2)*R51 822 & +(3.0d00*X_2+6.0d00*Y_2)*R41+3.0d00*R31 823 R(M,0,68)=X_2*Y_3*Z*R61+(3.0d00*X_2*Y*Z+Y_3*Z)*R51 824 & +3.0d00*Y*Z*R41 825 R(M,0,69)=X_2*Y_2*Z_2*R61 826 & +(Y_2*Z_2+X_2*(Y_2+Z_2))*R51 827 & +(X_2+Y_2+Z_2)*R41+R31 828 R(M,0,70)=X_2*Y*Z_3*R61+(3.0d00*X_2*Y*Z+Y*Z_3)*R51 829 & +3.0d00*Y*Z*R41 830 R(M,0,71)=X_2*Z_4*R61+(Z_4+6.0d00*X_2*Z_2)*R51 831 & +(3.0d00*X_2+6.0d00*Z_2)*R41+3.0d00*R31 832 R(M,0,72)=X*Y_5*R61+10.0d00* X*Y_3*R51+15.0d00*X*Y*R41 833 R(M,0,73)=X*Y_4*Z*R61+6.0d00*X*Y_2*Z*R51+3.0d00*X*Z*R41 834 R(M,0,74)=X*Y_3*Z_2*R61 835 & +X*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R51+3.0d00*X*Y*R41 836 R(M,0,75)=X*Y_2*Z_3*R61+X*(3.0d00*Y_2*Z+Z_3)*R51 837 & +3.0d00*X*Z*R41 838 R(M,0,76)=X*Y*Z_4*R61+6.0d00*X*Y*Z_2*R51+3.0d00*X*Y*R41 839 R(M,0,77)=X*Z_5*R61+10.0d00* X*Z_3*R51+15.0d00*X*Z*R41 840 R(M,0,78)=Y_6*R61+15.0d00*Y_4*R51+45.0d00*Y_2*R41 841 & +15.0d00*R31 842 R(M,0,79)=Y_5*Z*R61+10.0d00* Y_3*Z*R51+15.0d00*Y*Z*R41 843 R(M,0,80)=Y_4*Z_2*R61 844 & +(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))*R51 845 & +(6.0d00*Y_2+3.0d00*Z_2)*R41+3.0d00*R31 846 R(M,0,81)=Y_3*Z_3*R61 847 & +(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)*R51+9.0d00*Y*Z*R41 848 R(M,0,82)=Y_2*Z_4*R61+(Z_4+6.0d00*Y_2*Z_2)*R51 849 & +(3.0d00*Y_2+6.0d00*Z_2)*R41+3.0d00*R31 850 R(M,0,83)=Y*Z_5*R61+10.0d00* Y*Z_3*R51+15.0d00*Y*Z*R41 851 R(M,0,84)=Z_6*R61+15.0d00*Z_4*R51+45.0d00*Z_2*R41 852 & +15.0d00*R31 85300100 continue 854c 855 end 856 subroutine hfmkr_7(R,PQ,IJK,NAB,Lr,Lr3) 857 implicit none 858c $Id$ 859c 860 integer NAB, Lr, Lr3 861 integer IJK(0:Lr,0:Lr,0:Lr) 862 double precision PQ(NAB,3) 863 double precision R(NAB,0:Lr,Lr3) 864c 865 integer m 866 double precision x,y,z 867 double precision x_2, y_2, z_2 868 double precision x_3, y_3, z_3 869 double precision x_4, y_4, z_4 870 double precision x_5, y_5, z_5 871 double precision x_6, y_6, z_6 872 double precision x_7, y_7, z_7 873 double precision R11,R21,R31,R41,R51,R61,R71 874c 875 IJK(0,0,0)=1 876 IJK(1,0,0)=2 877 IJK(0,1,0)=3 878 IJK(0,0,1)=4 879 IJK(2,0,0)=5 880 IJK(1,1,0)=6 881 IJK(1,0,1)=7 882 IJK(0,2,0)=8 883 IJK(0,1,1)=9 884 IJK(0,0,2)=10 885 IJK(3,0,0)=11 886 IJK(2,1,0)=12 887 IJK(2,0,1)=13 888 IJK(1,2,0)=14 889 IJK(1,1,1)=15 890 IJK(1,0,2)=16 891 IJK(0,3,0)=17 892 IJK(0,2,1)=18 893 IJK(0,1,2)=19 894 IJK(0,0,3)=20 895 IJK(4,0,0)=21 896 IJK(3,1,0)=22 897 IJK(3,0,1)=23 898 IJK(2,2,0)=24 899 IJK(2,1,1)=25 900 IJK(2,0,2)=26 901 IJK(1,3,0)=27 902 IJK(1,2,1)=28 903 IJK(1,1,2)=29 904 IJK(1,0,3)=30 905 IJK(0,4,0)=31 906 IJK(0,3,1)=32 907 IJK(0,2,2)=33 908 IJK(0,1,3)=34 909 IJK(0,0,4)=35 910 IJK(5,0,0)=36 911 IJK(4,1,0)=37 912 IJK(4,0,1)=38 913 IJK(3,2,0)=39 914 IJK(3,1,1)=40 915 IJK(3,0,2)=41 916 IJK(2,3,0)=42 917 IJK(2,2,1)=43 918 IJK(2,1,2)=44 919 IJK(2,0,3)=45 920 IJK(1,4,0)=46 921 IJK(1,3,1)=47 922 IJK(1,2,2)=48 923 IJK(1,1,3)=49 924 IJK(1,0,4)=50 925 IJK(0,5,0)=51 926 IJK(0,4,1)=52 927 IJK(0,3,2)=53 928 IJK(0,2,3)=54 929 IJK(0,1,4)=55 930 IJK(0,0,5)=56 931 IJK(6,0,0)=57 932 IJK(5,1,0)=58 933 IJK(5,0,1)=59 934 IJK(4,2,0)=60 935 IJK(4,1,1)=61 936 IJK(4,0,2)=62 937 IJK(3,3,0)=63 938 IJK(3,2,1)=64 939 IJK(3,1,2)=65 940 IJK(3,0,3)=66 941 IJK(2,4,0)=67 942 IJK(2,3,1)=68 943 IJK(2,2,2)=69 944 IJK(2,1,3)=70 945 IJK(2,0,4)=71 946 IJK(1,5,0)=72 947 IJK(1,4,1)=73 948 IJK(1,3,2)=74 949 IJK(1,2,3)=75 950 IJK(1,1,4)=76 951 IJK(1,0,5)=77 952 IJK(0,6,0)=78 953 IJK(0,5,1)=79 954 IJK(0,4,2)=80 955 IJK(0,3,3)=81 956 IJK(0,2,4)=82 957 IJK(0,1,5)=83 958 IJK(0,0,6)=84 959 IJK(7,0,0)=85 960 IJK(6,1,0)=86 961 IJK(6,0,1)=87 962 IJK(5,2,0)=88 963 IJK(5,1,1)=89 964 IJK(5,0,2)=90 965 IJK(4,3,0)=91 966 IJK(4,2,1)=92 967 IJK(4,1,2)=93 968 IJK(4,0,3)=94 969 IJK(3,4,0)=95 970 IJK(3,3,1)=96 971 IJK(3,2,2)=97 972 IJK(3,1,3)=98 973 IJK(3,0,4)=99 974 IJK(2,5,0)=100 975 IJK(2,4,1)=101 976 IJK(2,3,2)=102 977 IJK(2,2,3)=103 978 IJK(2,1,4)=104 979 IJK(2,0,5)=105 980 IJK(1,6,0)=106 981 IJK(1,5,1)=107 982 IJK(1,4,2)=108 983 IJK(1,3,3)=109 984 IJK(1,2,4)=110 985 IJK(1,1,5)=111 986 IJK(1,0,6)=112 987 IJK(0,7,0)=113 988 IJK(0,6,1)=114 989 IJK(0,5,2)=115 990 IJK(0,4,3)=116 991 IJK(0,3,4)=117 992 IJK(0,2,5)=118 993 IJK(0,1,6)=119 994 IJK(0,0,7)=120 995c 996 do 00100 m=1,NAB 997 X=PQ(m,1) 998 Y=PQ(m,2) 999 Z=PQ(m,3) 1000 R11=R(m,1,1) 1001 R(M,0,2)=X*R11 1002 R(M,0,3)=Y*R11 1003 R(M,0,4)=Z*R11 1004 X_2=X*X 1005 Y_2=Y*Y 1006 Z_2=Z*Z 1007 R21=R(m,2,1) 1008 R(M,0,5)=X_2*R21+R11 1009 R(M,0,6)=X*Y*R21 1010 R(M,0,7)=X*Z*R21 1011 R(M,0,8)=Y_2*R21+R11 1012 R(M,0,9)=Y*Z*R21 1013 R(M,0,10)=Z_2*R21+R11 1014 X_3=X_2*X 1015 Y_3=Y_2*Y 1016 Z_3=Z_2*Z 1017 R31=R(m,3,1) 1018 R(M,0,11)=X_3*R31+3.0d00*X*R21 1019 R(M,0,12)=X_2*Y*R31+Y*R21 1020 R(M,0,13)=X_2*Z*R31+Z*R21 1021 R(M,0,14)=X*Y_2*R31+X*R21 1022 R(M,0,15)=X*Y*Z*R31 1023 R(M,0,16)=X*Z_2*R31+X*R21 1024 R(M,0,17)=Y_3*R31+3.0d00*Y*R21 1025 R(M,0,18)=Y_2*Z*R31+Z*R21 1026 R(M,0,19)=Y*Z_2*R31+Y*R21 1027 R(M,0,20)=Z_3*R31+3.0d00*Z*R21 1028 X_4=X_3*X 1029 Y_4=Y_3*Y 1030 Z_4=Z_3*Z 1031 R41=R(m,4,1) 1032 R(M,0,21)=X_4*R41+6.0d00*X_2*R31+3.0d00*R21 1033 R(M,0,22)=X_3*Y*R41+3.0d00*X*Y*R31 1034 R(M,0,23)=X_3*Z*R41+3.0d00*X*Z*R31 1035 R(M,0,24)=X_2*Y_2*R41+(X_2+Y_2)*R31+R21 1036 R(M,0,25)=X_2*Y*Z*R41+Y*Z*R31 1037 R(M,0,26)=X_2*Z_2*R41+(X_2+Z_2)*R31+R21 1038 R(M,0,27)=X*Y_3*R41+3.0d00*X*Y*R31 1039 R(M,0,28)=X*Y_2*Z*R41+X*Z*R31 1040 R(M,0,29)=X*Y*Z_2*R41+X*Y*R31 1041 R(M,0,30)=X*Z_3*R41+3.0d00*X*Z*R31 1042 R(M,0,31)=Y_4*R41+6.0d00*Y_2*R31+3.0d00*R21 1043 R(M,0,32)=Y_3*Z*R41+3.0d00*Y*Z*R31 1044 R(M,0,33)=Y_2*Z_2*R41+(Y_2+Z_2)*R31+R21 1045 R(M,0,34)=Y*Z_3*R41+3.0d00*Y*Z*R31 1046 R(M,0,35)=Z_4*R41+6.0d00*Z_2*R31+3.0d00*R21 1047 X_5=X_4*X 1048 Y_5=Y_4*Y 1049 Z_5=Z_4*Z 1050 R51=R(m,5,1) 1051 R(M,0,36)=X_5*R51+10.0d00* X_3*R41+15.0d00*X*R31 1052 R(M,0,37)=X_4*Y*R51+6.0d00*X_2*Y*R41+3.0d00*Y*R31 1053 R(M,0,38)=X_4*Z*R51+6.0d00*X_2*Z*R41+3.0d00*Z*R31 1054 R(M,0,39)=X_3*Y_2*R51+(X*(X_2+Y_2)+2.0d00*X*Y_2)*R41 1055 & +3.0d00*X*R31 1056 R(M,0,40)=X_3*Y*Z*R51+3.0d00*X*Y*Z*R41 1057 R(M,0,41)=X_3*Z_2*R51+(X*(X_2+Z_2)+2.0d00*X*Z_2)*R41 1058 & +3.0d00*X*R31 1059 R(M,0,42)=X_2*Y_3*R51+(3.0d00*X_2*Y+Y_3)*R41 1060 & +3.0d00*Y*R31 1061 R(M,0,43)=X_2*Y_2*Z*R51+(X_2*Z+Y_2*Z)*R41 1062 & +Z*R31 1063 R(M,0,44)=X_2*Y*Z_2*R51+(X_2*Y+Y*Z_2)*R41 1064 & +Y*R31 1065 R(M,0,45)=X_2*Z_3*R51+(3.0d00*X_2*Z+Z_3)*R41 1066 & +3.0d00*Z*R31 1067 R(M,0,46)=X*Y_4*R51+6.0d00*X*Y_2*R41+3.0d00*X*R31 1068 R(M,0,47)=X*Y_3*Z*R51+3.0d00*X*Y*Z*R41 1069 R(M,0,48)=X*Y_2*Z_2*R51+X*(Y_2+Z_2)*R41+ 1070 & X*R31 1071 R(M,0,49)=X*Y*Z_3*R51+3.0d00*X*Y*Z*R41 1072 R(M,0,50)=X*Z_4*R51+6.0d00*X*Z_2*R41+3.0d00*X*R31 1073 R(M,0,51)=Y_5*R51+10.0d00* Y_3*R41+15.0d00*Y*R31 1074 R(M,0,52)=Y_4*Z*R51+6.0d00*Y_2*Z*R41+3.0d00*Z*R31 1075 R(M,0,53)=Y_3*Z_2*R51+(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R41 1076 & +3.0d00*Y*R31 1077 R(M,0,54)=Y_2*Z_3*R51+(3.0d00*Y_2*Z+Z_3)*R41 1078 & +3.0d00*Z*R31 1079 R(M,0,55)=Y*Z_4*R51+6.0d00*Y*Z_2*R41+3.0d00*Y*R31 1080 R(M,0,56)=Z_5*R51+10.0d00* Z_3*R41+15.0d00*Z*R31 1081 X_6=X_5 * X 1082 Y_6=Y_5 * Y 1083 Z_6=Z_5 * Z 1084 R61=R(m,6,1) 1085 R(M,0,57)=X_6*R61+15.0d00*X_4*R51+45.0d00*X_2*R41 1086 & +15.0d00*R31 1087 R(M,0,58)=X_5*Y*R61+10.0d00* X_3*Y*R51+15.0d00*X*Y*R41 1088 R(M,0,59)=X_5*Z*R61+10.0d00* X_3*Z*R51+15.0d00*X*Z*R41 1089 R(M,0,60)=X_4*Y_2*R61 1090 & +(3.0d00*X_2*Y_2+X*(X*(X_2+Y_2)+2.0d00*X*Y_2))*R51 1091 & +(6.0d00*X_2+3.0d00*Y_2)*R41+3.0d00*R31 1092 R(M,0,61)=X_4*Y*Z*R61+6.0d00*X_2*Y*Z*R51+3.0d00*Y*Z*R41 1093 R(M,0,62)=X_4*Z_2*R61 1094 & +(3.0d00*X_2*Z_2+X*(X*(X_2+Z_2)+2.0d00*X*Z_2))*R51 1095 & +(6.0d00*X_2+3.0d00*Z_2)*R41+3.0d00*R31 1096 R(M,0,63)=X_3*Y_3*R61 1097 & +(X*(3.0d00*X_2*Y+Y_3)+2.0d00*X*Y_3)*R51+9.0d00*X*Y*R41 1098 R(M,0,64)=X_3*Y_2*Z*R61 1099 & +(X*(X_2*Z+Y_2*Z)+2.0d00*X*Y_2*Z)*R51+3.0d00*X*Z*R41 1100 R(M,0,65)=X_3*Y*Z_2*R61 1101 & +(X*(X_2*Y+Y*Z_2)+2.0d00*X*Y*Z_2)*R51+3.0d00*X*Y*R41 1102 R(M,0,66)=X_3*Z_3*R61 1103 & +(X*(3.0d00*X_2*Z+Z_3)+2.0d00*X*Z_3)*R51+9.0d00*X*Z*R41 1104 R(M,0,67)=X_2*Y_4*R61+(Y_4+6.0d00*X_2*Y_2)*R51 1105 & +(3.0d00*X_2+6.0d00*Y_2)*R41+3.0d00*R31 1106 R(M,0,68)=X_2*Y_3*Z*R61+(3.0d00*X_2*Y*Z+Y_3*Z)*R51 1107 & +3.0d00*Y*Z*R41 1108 R(M,0,69)=X_2*Y_2*Z_2*R61 1109 & +(Y_2*Z_2+X_2*(Y_2+Z_2))*R51 1110 & +(X_2+Y_2+Z_2)*R41+R31 1111 R(M,0,70)=X_2*Y*Z_3*R61+(3.0d00*X_2*Y*Z+Y*Z_3)*R51 1112 & +3.0d00*Y*Z*R41 1113 R(M,0,71)=X_2*Z_4*R61+(Z_4+6.0d00*X_2*Z_2)*R51 1114 & +(3.0d00*X_2+6.0d00*Z_2)*R41+3.0d00*R31 1115 R(M,0,72)=X*Y_5*R61+10.0d00* X*Y_3*R51+15.0d00*X*Y*R41 1116 R(M,0,73)=X*Y_4*Z*R61+6.0d00*X*Y_2*Z*R51+3.0d00*X*Z*R41 1117 R(M,0,74)=X*Y_3*Z_2*R61 1118 & +X*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R51+3.0d00*X*Y*R41 1119 R(M,0,75)=X*Y_2*Z_3*R61+X*(3.0d00*Y_2*Z+Z_3)*R51 1120 & +3.0d00*X*Z*R41 1121 R(M,0,76)=X*Y*Z_4*R61+6.0d00*X*Y*Z_2*R51+3.0d00*X*Y*R41 1122 R(M,0,77)=X*Z_5*R61+10.0d00* X*Z_3*R51+15.0d00*X*Z*R41 1123 R(M,0,78)=Y_6*R61+15.0d00*Y_4*R51+45.0d00*Y_2*R41 1124 & +15.0d00*R31 1125 R(M,0,79)=Y_5*Z*R61+10.0d00* Y_3*Z*R51+15.0d00*Y*Z*R41 1126 R(M,0,80)=Y_4*Z_2*R61 1127 & +(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))*R51 1128 & +(6.0d00*Y_2+3.0d00*Z_2)*R41+3.0d00*R31 1129 R(M,0,81)=Y_3*Z_3*R61 1130 & +(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)*R51+9.0d00*Y*Z*R41 1131 R(M,0,82)=Y_2*Z_4*R61+(Z_4+6.0d00*Y_2*Z_2)*R51 1132 & +(3.0d00*Y_2+6.0d00*Z_2)*R41+3.0d00*R31 1133 R(M,0,83)=Y*Z_5*R61+10.0d00* Y*Z_3*R51+15.0d00*Y*Z*R41 1134 R(M,0,84)=Z_6*R61+15.0d00*Z_4*R51+45.0d00*Z_2*R41 1135 & +15.0d00*R31 1136 X_7=X_6 * X 1137 Y_7=Y_6 * Y 1138 Z_7=Z_6 * Z 1139 R71=R(m,7,1) 1140 R(M,0,85)=X_7*R71+21.0d00*X_5*R61+105.0d00*X_3*R51 1141 & +105.0d00*X*R41 1142 R(M,0,86)=X_6*Y*R71+15.0d00*X_4*Y*R61+45.0d00*X_2*Y*R51 1143 & +15.0d00*Y*R41 1144 R(M,0,87)=X_6*Z*R71+15.0d00*X_4*Z*R61+45.0d00*X_2*Z*R51 1145 & +15.0d00*Z*R41 1146 R(M,0,88)=X_5*Y_2*R71 1147 & +(X*(3.0d00*X_2*Y_2+X*(X*(X_2+Y_2)+ 1148 & 2.0d00*X*Y_2))+4.0d00*X_3*Y_2)*R61 1149 & +(X*(6.0d00*X_2+3.0d00*Y_2)+ 1150 & 4.0d00*X*(X_2+Y_2)+8.0d00*X*Y_2)*R51 1151 & +15.0d00*X*R41 1152 R(M,0,89)=X_5*Y*Z*R71+10.0d00* X_3*Y*Z*R61 1153 & +15.0d00*X*Y*Z*R51 1154 R(M,0,90)=X_5*Z_2*R71 1155 & +(X*(3.0d00*X_2*Z_2+X*(X*(X_2+Z_2)+ 1156 & 2.0d00*X*Z_2))+4.0d00*X_3*Z_2)*R61 1157 & +(X*(6.0d00*X_2+3.0d00*Z_2)+ 1158 & 4.0d00*X*(X_2+Z_2)+8.0d00*X*Z_2)*R51 1159 & +15.0d00*X*R41 1160 R(M,0,91)=X_4*Y_3*R71 1161 & +(3.0d00*X_2*Y_3+X*(X*(3.0d00*X_2*Y+Y_3)+2.0d00*X*Y_3))*R61 1162 & +(18.0d00*X_2*Y+3.0d00*Y_3)*R51+9.0d00*Y*R41 1163 R(M,0,92)=X_4*Y_2*Z*R71 1164 & +(3.0d00*X_2*Y_2*Z+X*(X*(X_2*Z+Y_2*Z)+2.0d00*X*Y_2*Z))*R61 1165 & +(6.0d00*X_2*Z+3.0d00*Y_2*Z)*R51+3.0d00*Z*R41 1166 R(M,0,93)=X_4*Y*Z_2*R71 1167 & +(X*(X*(X_2*Y+Y*Z_2)+2.0d00*X*Y*Z_2)+3.0d00*X_2*Y*Z_2)*R61 1168 & +(6.0d00*X_2*Y+3.0d00*Y*Z_2)*R51+3.0d00*Y*R41 1169 R(M,0,94)=X_4*Z_3*R71 1170 & +(3.0d00*X_2*Z_3+X*(X*(3.0d00*X_2*Z+Z_3)+2.0d00*X*Z_3))*R61 1171 & +(18.0d00*X_2*Z+3.0d00*Z_3)*R51+9.0d00*Z*R41 1172 R(M,0,95)=X_3*Y_4*R71 1173 & +(X*(Y_4+6.0d00*X_2*Y_2)+2.0d00*X*Y_4)*R61 1174 & +(X*(3.0d00*X_2+6.0d00*Y_2)+12.0d00*X*Y_2)*R51+9.0d00*X*R41 1175 R(M,0,96)=X_3*Y_3*Z*R71 1176 & +(X*(3.0d00*X_2*Y*Z+Y_3*Z)+2.0d00*X*Y_3*Z)*R61 1177 & +9.0d00*X*Y*Z*R51 1178 R(M,0,97)=X_3*Y_2*Z_2*R71 1179 & +(X*(Y_2*Z_2+X_2*(Y_2+Z_2))+2.0d00*X*Y_2*Z_2)*R61 1180 & +(X*(X_2+Y_2+Z_2)+2.0d00*X*(Y_2+Z_2))*R51+ 1181 & 3.0d00*X*R41 1182 R(M,0,98)=X_3*Y*Z_3*R71 1183 & +(X*(3.0d00*X_2*Y*Z+Y*Z_3)+2.0d00*X*Y*Z_3)*R61 1184 & +9.0d00*X*Y*Z*R51 1185 R(M,0,99)=X_3*Z_4*R71 1186 & +(X*(Z_4+6.0d00*X_2*Z_2)+2.0d00*X*Z_4)*R61 1187 & +(X*(3.0d00*X_2+6.0d00*Z_2)+12.0d00*X*Z_2)*R51+9.0d00*X*R41 1188 R(M,0,100)=X_2*Y_5*R71+(Y_5+10.0d00* X_2*Y_3)*R61 1189 & +(10.0d00*Y_3+15.0d00*X_2*Y)*R51+15.0d00*Y*R41 1190 R(M,0,101)=X_2*Y_4*Z*R71+ 1191 & (Y_4*Z+6.0d00*X_2*Y_2*Z)*R61 1192 & +(3.0d00*X_2*Z+6.0d00*Y_2*Z)*R51+3.0d00*Z*R41 1193 R(M,0,102)=X_2*Y_3*Z_2*R71 1194 & +(X_2*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+Y_3*Z_2)*R61 1195 & +(3.0d00*X_2*Y+Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R51+3.0d00*Y*R41 1196 R(M,0,103)=X_2*Y_2*Z_3*R71 1197 & +(Y_2*Z_3+X_2*(3.0d00*Y_2*Z+Z_3))*R61 1198 & +(3.0d00*X_2*Z+3.0d00*Y_2*Z+Z_3)*R51+3.0d00*Z*R41 1199 R(M,0,104)=X_2*Y*Z_4*R71+ 1200 & (6.0d00*X_2*Y*Z_2+Y*Z_4)*R61 1201 & +(3.0d00*X_2*Y+6.0d00*Y*Z_2)*R51+3.0d00*Y*R41 1202 R(M,0,105)=X_2*Z_5*R71+(Z_5+10.0d00* X_2*Z_3)*R61 1203 & +(10.0d00*Z_3+15.0d00*X_2*Z)*R51+15.0d00*Z*R41 1204 R(M,0,106)=X*Y_6*R71+15.0d00*X*Y_4*R61+ 1205 & 45.0d00*X*Y_2*R51+15.0d00*X*R41 1206 R(M,0,107)=X*Y_5*Z*R71+10.0d00* X*Y_3*Z*R61 1207 & +15.0d00*X*Y*Z*R51 1208 R(M,0,108)=X*Y_4*Z_2*R71 1209 & +X*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))*R61 1210 & +X*(6.0d00*Y_2+3.0d00*Z_2)*R51+3.0d00*X*R41 1211 R(M,0,109)=X*Y_3*Z_3*R71 1212 & +X*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)*R61+9.0d00*X*Y*Z*R51 1213 R(M,0,110)=X*Y_2*Z_4*R71+X*(Z_4+6.0d00*Y_2*Z_2)*R61 1214 & +X*(3.0d00*Y_2+6.0d00*Z_2)*R51+3.0d00*X*R41 1215 R(M,0,111)=X*Y*Z_5*R71+10.0d00* X*Y*Z_3*R61 1216 & +15.0d00*X*Y*Z*R51 1217 R(M,0,112)=X*Z_6*R71+15.0d00*X*Z_4*R61+ 1218 & 45.0d00*X*Z_2*R51 1219 & +15.0d00*X*R41 1220 R(M,0,113)=Y_7*R71+21.0d00*Y_5*R61+105.0d00*Y_3*R51 1221 & +105.0d00*Y*R41 1222 R(M,0,114)=Y_6*Z*R71+15.0d00*Y_4*Z*R61+ 1223 & 45.0d00*Y_2*Z*R51 1224 & +15.0d00*Z*R41 1225 R(M,0,115)=Y_5*Z_2*R71 1226 & +(Y*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 1227 & 2.0d00*Y*Z_2))+4.0d00*Y_3*Z_2)*R61 1228 & +(Y*(6.0d00*Y_2+3.0d00*Z_2)+4.0d00*Y*(Y_2+Z_2)+ 1229 & 8.0d00*Y*Z_2)*R51 1230 & +15.0d00*Y*R41 1231 R(M,0,116)=Y_4*Z_3*R71 1232 & +(3.0d00*Y_2*Z_3+Y*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3))*R61 1233 & +(18.0d00*Y_2*Z+3.0d00*Z_3)*R51+9.0d00*Z*R41 1234 R(M,0,117)=Y_3*Z_4*R71 1235 & +(Y*(Z_4+6.0d00*Y_2*Z_2)+2.0d00*Y*Z_4)*R61 1236 & +(Y*(3.0d00*Y_2+6.0d00*Z_2)+12.0d00*Y*Z_2)*R51+9.0d00*Y*R41 1237 R(M,0,118)=Y_2*Z_5*R71+(Z_5+10.0d00* Y_2*Z_3)*R61 1238 & +(10.0d00*Z_3+15.0d00*Y_2*Z)*R51+15.0d00*Z*R41 1239 R(M,0,119)=Y*Z_6*R71+15.0d00*Y*Z_4*R61+ 1240 & 45.0d00*Y*Z_2*R51 1241 & +15.0d00*Y*R41 1242 R(M,0,120)=Z_7*R71+21.0d00*Z_5*R61+105.0d00*Z_3*R51 1243 & +105.0d00*Z*R41 124400100 continue 1245c 1246 end 1247 subroutine hfmkr_8(R,PQ,IJK,NAB,Lr,Lr3) 1248 implicit none 1249c $Id$ 1250c 1251 integer NAB, Lr, Lr3 1252 integer IJK(0:Lr,0:Lr,0:Lr) 1253 double precision PQ(NAB,3) 1254 double precision R(NAB,0:Lr,Lr3) 1255c 1256 integer m 1257 double precision x,y,z 1258 double precision x_2, y_2, z_2 1259 double precision x_3, y_3, z_3 1260 double precision x_4, y_4, z_4 1261 double precision x_5, y_5, z_5 1262 double precision x_6, y_6, z_6 1263 double precision x_7, y_7, z_7 1264 double precision x_8, y_8, z_8 1265 double precision R11,R21,R31,R41,R51,R61,R71,R81 1266c 1267 IJK(0,0,0)=1 1268 IJK(1,0,0)=2 1269 IJK(0,1,0)=3 1270 IJK(0,0,1)=4 1271 IJK(2,0,0)=5 1272 IJK(1,1,0)=6 1273 IJK(1,0,1)=7 1274 IJK(0,2,0)=8 1275 IJK(0,1,1)=9 1276 IJK(0,0,2)=10 1277 IJK(3,0,0)=11 1278 IJK(2,1,0)=12 1279 IJK(2,0,1)=13 1280 IJK(1,2,0)=14 1281 IJK(1,1,1)=15 1282 IJK(1,0,2)=16 1283 IJK(0,3,0)=17 1284 IJK(0,2,1)=18 1285 IJK(0,1,2)=19 1286 IJK(0,0,3)=20 1287 IJK(4,0,0)=21 1288 IJK(3,1,0)=22 1289 IJK(3,0,1)=23 1290 IJK(2,2,0)=24 1291 IJK(2,1,1)=25 1292 IJK(2,0,2)=26 1293 IJK(1,3,0)=27 1294 IJK(1,2,1)=28 1295 IJK(1,1,2)=29 1296 IJK(1,0,3)=30 1297 IJK(0,4,0)=31 1298 IJK(0,3,1)=32 1299 IJK(0,2,2)=33 1300 IJK(0,1,3)=34 1301 IJK(0,0,4)=35 1302 IJK(5,0,0)=36 1303 IJK(4,1,0)=37 1304 IJK(4,0,1)=38 1305 IJK(3,2,0)=39 1306 IJK(3,1,1)=40 1307 IJK(3,0,2)=41 1308 IJK(2,3,0)=42 1309 IJK(2,2,1)=43 1310 IJK(2,1,2)=44 1311 IJK(2,0,3)=45 1312 IJK(1,4,0)=46 1313 IJK(1,3,1)=47 1314 IJK(1,2,2)=48 1315 IJK(1,1,3)=49 1316 IJK(1,0,4)=50 1317 IJK(0,5,0)=51 1318 IJK(0,4,1)=52 1319 IJK(0,3,2)=53 1320 IJK(0,2,3)=54 1321 IJK(0,1,4)=55 1322 IJK(0,0,5)=56 1323 IJK(6,0,0)=57 1324 IJK(5,1,0)=58 1325 IJK(5,0,1)=59 1326 IJK(4,2,0)=60 1327 IJK(4,1,1)=61 1328 IJK(4,0,2)=62 1329 IJK(3,3,0)=63 1330 IJK(3,2,1)=64 1331 IJK(3,1,2)=65 1332 IJK(3,0,3)=66 1333 IJK(2,4,0)=67 1334 IJK(2,3,1)=68 1335 IJK(2,2,2)=69 1336 IJK(2,1,3)=70 1337 IJK(2,0,4)=71 1338 IJK(1,5,0)=72 1339 IJK(1,4,1)=73 1340 IJK(1,3,2)=74 1341 IJK(1,2,3)=75 1342 IJK(1,1,4)=76 1343 IJK(1,0,5)=77 1344 IJK(0,6,0)=78 1345 IJK(0,5,1)=79 1346 IJK(0,4,2)=80 1347 IJK(0,3,3)=81 1348 IJK(0,2,4)=82 1349 IJK(0,1,5)=83 1350 IJK(0,0,6)=84 1351 IJK(7,0,0)=85 1352 IJK(6,1,0)=86 1353 IJK(6,0,1)=87 1354 IJK(5,2,0)=88 1355 IJK(5,1,1)=89 1356 IJK(5,0,2)=90 1357 IJK(4,3,0)=91 1358 IJK(4,2,1)=92 1359 IJK(4,1,2)=93 1360 IJK(4,0,3)=94 1361 IJK(3,4,0)=95 1362 IJK(3,3,1)=96 1363 IJK(3,2,2)=97 1364 IJK(3,1,3)=98 1365 IJK(3,0,4)=99 1366 IJK(2,5,0)=100 1367 IJK(2,4,1)=101 1368 IJK(2,3,2)=102 1369 IJK(2,2,3)=103 1370 IJK(2,1,4)=104 1371 IJK(2,0,5)=105 1372 IJK(1,6,0)=106 1373 IJK(1,5,1)=107 1374 IJK(1,4,2)=108 1375 IJK(1,3,3)=109 1376 IJK(1,2,4)=110 1377 IJK(1,1,5)=111 1378 IJK(1,0,6)=112 1379 IJK(0,7,0)=113 1380 IJK(0,6,1)=114 1381 IJK(0,5,2)=115 1382 IJK(0,4,3)=116 1383 IJK(0,3,4)=117 1384 IJK(0,2,5)=118 1385 IJK(0,1,6)=119 1386 IJK(0,0,7)=120 1387 IJK(8,0,0)=121 1388 IJK(7,1,0)=122 1389 IJK(7,0,1)=123 1390 IJK(6,2,0)=124 1391 IJK(6,1,1)=125 1392 IJK(6,0,2)=126 1393 IJK(5,3,0)=127 1394 IJK(5,2,1)=128 1395 IJK(5,1,2)=129 1396 IJK(5,0,3)=130 1397 IJK(4,4,0)=131 1398 IJK(4,3,1)=132 1399 IJK(4,2,2)=133 1400 IJK(4,1,3)=134 1401 IJK(4,0,4)=135 1402 IJK(3,5,0)=136 1403 IJK(3,4,1)=137 1404 IJK(3,3,2)=138 1405 IJK(3,2,3)=139 1406 IJK(3,1,4)=140 1407 IJK(3,0,5)=141 1408 IJK(2,6,0)=142 1409 IJK(2,5,1)=143 1410 IJK(2,4,2)=144 1411 IJK(2,3,3)=145 1412 IJK(2,2,4)=146 1413 IJK(2,1,5)=147 1414 IJK(2,0,6)=148 1415 IJK(1,7,0)=149 1416 IJK(1,6,1)=150 1417 IJK(1,5,2)=151 1418 IJK(1,4,3)=152 1419 IJK(1,3,4)=153 1420 IJK(1,2,5)=154 1421 IJK(1,1,6)=155 1422 IJK(1,0,7)=156 1423 IJK(0,8,0)=157 1424 IJK(0,7,1)=158 1425 IJK(0,6,2)=159 1426 IJK(0,5,3)=160 1427 IJK(0,4,4)=161 1428 IJK(0,3,5)=162 1429 IJK(0,2,6)=163 1430 IJK(0,1,7)=164 1431 IJK(0,0,8)=165 1432c 1433 do 00100 m=1,NAB 1434 X=PQ(m,1) 1435 Y=PQ(m,2) 1436 Z=PQ(m,3) 1437 R11=R(m,1,1) 1438 R(M,0,2)=X*R11 1439 R(M,0,3)=Y*R11 1440 R(M,0,4)=Z*R11 1441 X_2=X*X 1442 Y_2=Y*Y 1443 Z_2=Z*Z 1444 R21=R(m,2,1) 1445 R(M,0,5)=X_2*R21+R11 1446 R(M,0,6)=X*Y*R21 1447 R(M,0,7)=X*Z*R21 1448 R(M,0,8)=Y_2*R21+R11 1449 R(M,0,9)=Y*Z*R21 1450 R(M,0,10)=Z_2*R21+R11 1451 X_3=X_2*X 1452 Y_3=Y_2*Y 1453 Z_3=Z_2*Z 1454 R31=R(m,3,1) 1455 R(M,0,11)=X_3*R31+3.0d00*X*R21 1456 R(M,0,12)=X_2*Y*R31+Y*R21 1457 R(M,0,13)=X_2*Z*R31+Z*R21 1458 R(M,0,14)=X*Y_2*R31+X*R21 1459 R(M,0,15)=X*Y*Z*R31 1460 R(M,0,16)=X*Z_2*R31+X*R21 1461 R(M,0,17)=Y_3*R31+3.0d00*Y*R21 1462 R(M,0,18)=Y_2*Z*R31+Z*R21 1463 R(M,0,19)=Y*Z_2*R31+Y*R21 1464 R(M,0,20)=Z_3*R31+3.0d00*Z*R21 1465 X_4=X_3*X 1466 Y_4=Y_3*Y 1467 Z_4=Z_3*Z 1468 R41=R(m,4,1) 1469 R(M,0,21)=X_4*R41+6.0d00*X_2*R31+3.0d00*R21 1470 R(M,0,22)=X_3*Y*R41+3.0d00*X*Y*R31 1471 R(M,0,23)=X_3*Z*R41+3.0d00*X*Z*R31 1472 R(M,0,24)=X_2*Y_2*R41+(X_2+Y_2)*R31+R21 1473 R(M,0,25)=X_2*Y*Z*R41+Y*Z*R31 1474 R(M,0,26)=X_2*Z_2*R41+(X_2+Z_2)*R31+R21 1475 R(M,0,27)=X*Y_3*R41+3.0d00*X*Y*R31 1476 R(M,0,28)=X*Y_2*Z*R41+X*Z*R31 1477 R(M,0,29)=X*Y*Z_2*R41+X*Y*R31 1478 R(M,0,30)=X*Z_3*R41+3.0d00*X*Z*R31 1479 R(M,0,31)=Y_4*R41+6.0d00*Y_2*R31+3.0d00*R21 1480 R(M,0,32)=Y_3*Z*R41+3.0d00*Y*Z*R31 1481 R(M,0,33)=Y_2*Z_2*R41+(Y_2+Z_2)*R31+R21 1482 R(M,0,34)=Y*Z_3*R41+3.0d00*Y*Z*R31 1483 R(M,0,35)=Z_4*R41+6.0d00*Z_2*R31+3.0d00*R21 1484 X_5=X_4*X 1485 Y_5=Y_4*Y 1486 Z_5=Z_4*Z 1487 R51=R(m,5,1) 1488 R(M,0,36)=X_5*R51+10.0d00* X_3*R41+15.0d00*X*R31 1489 R(M,0,37)=X_4*Y*R51+6.0d00*X_2*Y*R41+3.0d00*Y*R31 1490 R(M,0,38)=X_4*Z*R51+6.0d00*X_2*Z*R41+3.0d00*Z*R31 1491 R(M,0,39)=X_3*Y_2*R51+(X*(X_2+Y_2)+2.0d00*X*Y_2)*R41 1492 & +3.0d00*X*R31 1493 R(M,0,40)=X_3*Y*Z*R51+3.0d00*X*Y*Z*R41 1494 R(M,0,41)=X_3*Z_2*R51+(X*(X_2+Z_2)+2.0d00*X*Z_2)*R41 1495 & +3.0d00*X*R31 1496 R(M,0,42)=X_2*Y_3*R51+(3.0d00*X_2*Y+Y_3)*R41 1497 & +3.0d00*Y*R31 1498 R(M,0,43)=X_2*Y_2*Z*R51+(X_2*Z+Y_2*Z)*R41 1499 & +Z*R31 1500 R(M,0,44)=X_2*Y*Z_2*R51+(X_2*Y+Y*Z_2)*R41 1501 & +Y*R31 1502 R(M,0,45)=X_2*Z_3*R51+(3.0d00*X_2*Z+Z_3)*R41 1503 & +3.0d00*Z*R31 1504 R(M,0,46)=X*Y_4*R51+6.0d00*X*Y_2*R41+3.0d00*X*R31 1505 R(M,0,47)=X*Y_3*Z*R51+3.0d00*X*Y*Z*R41 1506 R(M,0,48)=X*Y_2*Z_2*R51+X*(Y_2+Z_2)*R41+ 1507 & X*R31 1508 R(M,0,49)=X*Y*Z_3*R51+3.0d00*X*Y*Z*R41 1509 R(M,0,50)=X*Z_4*R51+6.0d00*X*Z_2*R41+3.0d00*X*R31 1510 R(M,0,51)=Y_5*R51+10.0d00* Y_3*R41+15.0d00*Y*R31 1511 R(M,0,52)=Y_4*Z*R51+6.0d00*Y_2*Z*R41+3.0d00*Z*R31 1512 R(M,0,53)=Y_3*Z_2*R51+(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R41 1513 & +3.0d00*Y*R31 1514 R(M,0,54)=Y_2*Z_3*R51+(3.0d00*Y_2*Z+Z_3)*R41 1515 & +3.0d00*Z*R31 1516 R(M,0,55)=Y*Z_4*R51+6.0d00*Y*Z_2*R41+3.0d00*Y*R31 1517 R(M,0,56)=Z_5*R51+10.0d00* Z_3*R41+15.0d00*Z*R31 1518 X_6=X_5 * X 1519 Y_6=Y_5 * Y 1520 Z_6=Z_5 * Z 1521 R61=R(m,6,1) 1522 R(M,0,57)=X_6*R61+15.0d00*X_4*R51+45.0d00*X_2*R41 1523 & +15.0d00*R31 1524 R(M,0,58)=X_5*Y*R61+10.0d00* X_3*Y*R51+15.0d00*X*Y*R41 1525 R(M,0,59)=X_5*Z*R61+10.0d00* X_3*Z*R51+15.0d00*X*Z*R41 1526 R(M,0,60)=X_4*Y_2*R61 1527 & +(3.0d00*X_2*Y_2+X*(X*(X_2+Y_2)+2.0d00*X*Y_2))*R51 1528 & +(6.0d00*X_2+3.0d00*Y_2)*R41+3.0d00*R31 1529 R(M,0,61)=X_4*Y*Z*R61+6.0d00*X_2*Y*Z*R51+3.0d00*Y*Z*R41 1530 R(M,0,62)=X_4*Z_2*R61 1531 & +(3.0d00*X_2*Z_2+X*(X*(X_2+Z_2)+2.0d00*X*Z_2))*R51 1532 & +(6.0d00*X_2+3.0d00*Z_2)*R41+3.0d00*R31 1533 R(M,0,63)=X_3*Y_3*R61 1534 & +(X*(3.0d00*X_2*Y+Y_3)+2.0d00*X*Y_3)*R51+9.0d00*X*Y*R41 1535 R(M,0,64)=X_3*Y_2*Z*R61 1536 & +(X*(X_2*Z+Y_2*Z)+2.0d00*X*Y_2*Z)*R51+3.0d00*X*Z*R41 1537 R(M,0,65)=X_3*Y*Z_2*R61 1538 & +(X*(X_2*Y+Y*Z_2)+2.0d00*X*Y*Z_2)*R51+3.0d00*X*Y*R41 1539 R(M,0,66)=X_3*Z_3*R61 1540 & +(X*(3.0d00*X_2*Z+Z_3)+2.0d00*X*Z_3)*R51+9.0d00*X*Z*R41 1541 R(M,0,67)=X_2*Y_4*R61+(Y_4+6.0d00*X_2*Y_2)*R51 1542 & +(3.0d00*X_2+6.0d00*Y_2)*R41+3.0d00*R31 1543 R(M,0,68)=X_2*Y_3*Z*R61+(3.0d00*X_2*Y*Z+Y_3*Z)*R51 1544 & +3.0d00*Y*Z*R41 1545 R(M,0,69)=X_2*Y_2*Z_2*R61 1546 & +(Y_2*Z_2+X_2*(Y_2+Z_2))*R51 1547 & +(X_2+Y_2+Z_2)*R41+R31 1548 R(M,0,70)=X_2*Y*Z_3*R61+(3.0d00*X_2*Y*Z+Y*Z_3)*R51 1549 & +3.0d00*Y*Z*R41 1550 R(M,0,71)=X_2*Z_4*R61+(Z_4+6.0d00*X_2*Z_2)*R51 1551 & +(3.0d00*X_2+6.0d00*Z_2)*R41+3.0d00*R31 1552 R(M,0,72)=X*Y_5*R61+10.0d00* X*Y_3*R51+15.0d00*X*Y*R41 1553 R(M,0,73)=X*Y_4*Z*R61+6.0d00*X*Y_2*Z*R51+3.0d00*X*Z*R41 1554 R(M,0,74)=X*Y_3*Z_2*R61 1555 & +X*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R51+3.0d00*X*Y*R41 1556 R(M,0,75)=X*Y_2*Z_3*R61+X*(3.0d00*Y_2*Z+Z_3)*R51 1557 & +3.0d00*X*Z*R41 1558 R(M,0,76)=X*Y*Z_4*R61+6.0d00*X*Y*Z_2*R51+3.0d00*X*Y*R41 1559 R(M,0,77)=X*Z_5*R61+10.0d00* X*Z_3*R51+15.0d00*X*Z*R41 1560 R(M,0,78)=Y_6*R61+15.0d00*Y_4*R51+45.0d00*Y_2*R41 1561 & +15.0d00*R31 1562 R(M,0,79)=Y_5*Z*R61+10.0d00* Y_3*Z*R51+15.0d00*Y*Z*R41 1563 R(M,0,80)=Y_4*Z_2*R61 1564 & +(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))*R51 1565 & +(6.0d00*Y_2+3.0d00*Z_2)*R41+3.0d00*R31 1566 R(M,0,81)=Y_3*Z_3*R61 1567 & +(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)*R51+9.0d00*Y*Z*R41 1568 R(M,0,82)=Y_2*Z_4*R61+(Z_4+6.0d00*Y_2*Z_2)*R51 1569 & +(3.0d00*Y_2+6.0d00*Z_2)*R41+3.0d00*R31 1570 R(M,0,83)=Y*Z_5*R61+10.0d00* Y*Z_3*R51+15.0d00*Y*Z*R41 1571 R(M,0,84)=Z_6*R61+15.0d00*Z_4*R51+45.0d00*Z_2*R41 1572 & +15.0d00*R31 1573 X_7=X_6 * X 1574 Y_7=Y_6 * Y 1575 Z_7=Z_6 * Z 1576 R71=R(m,7,1) 1577 R(M,0,85)=X_7*R71+21.0d00*X_5*R61+105.0d00*X_3*R51 1578 & +105.0d00*X*R41 1579 R(M,0,86)=X_6*Y*R71+15.0d00*X_4*Y*R61+45.0d00*X_2*Y*R51 1580 & +15.0d00*Y*R41 1581 R(M,0,87)=X_6*Z*R71+15.0d00*X_4*Z*R61+45.0d00*X_2*Z*R51 1582 & +15.0d00*Z*R41 1583 R(M,0,88)=X_5*Y_2*R71 1584 & +(X*(3.0d00*X_2*Y_2+X*(X*(X_2+Y_2)+ 1585 & 2.0d00*X*Y_2))+4.0d00*X_3*Y_2)*R61 1586 & +(X*(6.0d00*X_2+3.0d00*Y_2)+ 1587 & 4.0d00*X*(X_2+Y_2)+8.0d00*X*Y_2)*R51 1588 & +15.0d00*X*R41 1589 R(M,0,89)=X_5*Y*Z*R71+10.0d00* X_3*Y*Z*R61 1590 & +15.0d00*X*Y*Z*R51 1591 R(M,0,90)=X_5*Z_2*R71 1592 & +(X*(3.0d00*X_2*Z_2+X*(X*(X_2+Z_2)+ 1593 & 2.0d00*X*Z_2))+4.0d00*X_3*Z_2)*R61 1594 & +(X*(6.0d00*X_2+3.0d00*Z_2)+ 1595 & 4.0d00*X*(X_2+Z_2)+8.0d00*X*Z_2)*R51 1596 & +15.0d00*X*R41 1597 R(M,0,91)=X_4*Y_3*R71 1598 & +(3.0d00*X_2*Y_3+X*(X*(3.0d00*X_2*Y+Y_3)+2.0d00*X*Y_3))*R61 1599 & +(18.0d00*X_2*Y+3.0d00*Y_3)*R51+9.0d00*Y*R41 1600 R(M,0,92)=X_4*Y_2*Z*R71 1601 & +(3.0d00*X_2*Y_2*Z+X*(X*(X_2*Z+Y_2*Z)+2.0d00*X*Y_2*Z))*R61 1602 & +(6.0d00*X_2*Z+3.0d00*Y_2*Z)*R51+3.0d00*Z*R41 1603 R(M,0,93)=X_4*Y*Z_2*R71 1604 & +(X*(X*(X_2*Y+Y*Z_2)+2.0d00*X*Y*Z_2)+3.0d00*X_2*Y*Z_2)*R61 1605 & +(6.0d00*X_2*Y+3.0d00*Y*Z_2)*R51+3.0d00*Y*R41 1606 R(M,0,94)=X_4*Z_3*R71 1607 & +(3.0d00*X_2*Z_3+X*(X*(3.0d00*X_2*Z+Z_3)+2.0d00*X*Z_3))*R61 1608 & +(18.0d00*X_2*Z+3.0d00*Z_3)*R51+9.0d00*Z*R41 1609 R(M,0,95)=X_3*Y_4*R71 1610 & +(X*(Y_4+6.0d00*X_2*Y_2)+2.0d00*X*Y_4)*R61 1611 & +(X*(3.0d00*X_2+6.0d00*Y_2)+12.0d00*X*Y_2)*R51+9.0d00*X*R41 1612 R(M,0,96)=X_3*Y_3*Z*R71 1613 & +(X*(3.0d00*X_2*Y*Z+Y_3*Z)+2.0d00*X*Y_3*Z)*R61 1614 & +9.0d00*X*Y*Z*R51 1615 R(M,0,97)=X_3*Y_2*Z_2*R71 1616 & +(X*(Y_2*Z_2+X_2*(Y_2+Z_2))+2.0d00*X*Y_2*Z_2)*R61 1617 & +(X*(X_2+Y_2+Z_2)+2.0d00*X*(Y_2+Z_2))*R51+ 1618 & 3.0d00*X*R41 1619 R(M,0,98)=X_3*Y*Z_3*R71 1620 & +(X*(3.0d00*X_2*Y*Z+Y*Z_3)+2.0d00*X*Y*Z_3)*R61 1621 & +9.0d00*X*Y*Z*R51 1622 R(M,0,99)=X_3*Z_4*R71 1623 & +(X*(Z_4+6.0d00*X_2*Z_2)+2.0d00*X*Z_4)*R61 1624 & +(X*(3.0d00*X_2+6.0d00*Z_2)+12.0d00*X*Z_2)*R51+9.0d00*X*R41 1625 R(M,0,100)=X_2*Y_5*R71+(Y_5+10.0d00* X_2*Y_3)*R61 1626 & +(10.0d00*Y_3+15.0d00*X_2*Y)*R51+15.0d00*Y*R41 1627 R(M,0,101)=X_2*Y_4*Z*R71+ 1628 & (Y_4*Z+6.0d00*X_2*Y_2*Z)*R61 1629 & +(3.0d00*X_2*Z+6.0d00*Y_2*Z)*R51+3.0d00*Z*R41 1630 R(M,0,102)=X_2*Y_3*Z_2*R71 1631 & +(X_2*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+Y_3*Z_2)*R61 1632 & +(3.0d00*X_2*Y+Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R51+3.0d00*Y*R41 1633 R(M,0,103)=X_2*Y_2*Z_3*R71 1634 & +(Y_2*Z_3+X_2*(3.0d00*Y_2*Z+Z_3))*R61 1635 & +(3.0d00*X_2*Z+3.0d00*Y_2*Z+Z_3)*R51+3.0d00*Z*R41 1636 R(M,0,104)=X_2*Y*Z_4*R71+ 1637 & (6.0d00*X_2*Y*Z_2+Y*Z_4)*R61 1638 & +(3.0d00*X_2*Y+6.0d00*Y*Z_2)*R51+3.0d00*Y*R41 1639 R(M,0,105)=X_2*Z_5*R71+(Z_5+10.0d00* X_2*Z_3)*R61 1640 & +(10.0d00*Z_3+15.0d00*X_2*Z)*R51+15.0d00*Z*R41 1641 R(M,0,106)=X*Y_6*R71+15.0d00*X*Y_4*R61+ 1642 & 45.0d00*X*Y_2*R51+15.0d00*X*R41 1643 R(M,0,107)=X*Y_5*Z*R71+10.0d00* X*Y_3*Z*R61 1644 & +15.0d00*X*Y*Z*R51 1645 R(M,0,108)=X*Y_4*Z_2*R71 1646 & +X*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))*R61 1647 & +X*(6.0d00*Y_2+3.0d00*Z_2)*R51+3.0d00*X*R41 1648 R(M,0,109)=X*Y_3*Z_3*R71 1649 & +X*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)*R61+9.0d00*X*Y*Z*R51 1650 R(M,0,110)=X*Y_2*Z_4*R71+X*(Z_4+6.0d00*Y_2*Z_2)*R61 1651 & +X*(3.0d00*Y_2+6.0d00*Z_2)*R51+3.0d00*X*R41 1652 R(M,0,111)=X*Y*Z_5*R71+10.0d00* X*Y*Z_3*R61 1653 & +15.0d00*X*Y*Z*R51 1654 R(M,0,112)=X*Z_6*R71+15.0d00*X*Z_4*R61+ 1655 & 45.0d00*X*Z_2*R51 1656 & +15.0d00*X*R41 1657 R(M,0,113)=Y_7*R71+21.0d00*Y_5*R61+105.0d00*Y_3*R51 1658 & +105.0d00*Y*R41 1659 R(M,0,114)=Y_6*Z*R71+15.0d00*Y_4*Z*R61+ 1660 & 45.0d00*Y_2*Z*R51 1661 & +15.0d00*Z*R41 1662 R(M,0,115)=Y_5*Z_2*R71 1663 & +(Y*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 1664 & 2.0d00*Y*Z_2))+4.0d00*Y_3*Z_2)*R61 1665 & +(Y*(6.0d00*Y_2+3.0d00*Z_2)+4.0d00*Y*(Y_2+Z_2)+ 1666 & 8.0d00*Y*Z_2)*R51 1667 & +15.0d00*Y*R41 1668 R(M,0,116)=Y_4*Z_3*R71 1669 & +(3.0d00*Y_2*Z_3+Y*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3))*R61 1670 & +(18.0d00*Y_2*Z+3.0d00*Z_3)*R51+9.0d00*Z*R41 1671 R(M,0,117)=Y_3*Z_4*R71 1672 & +(Y*(Z_4+6.0d00*Y_2*Z_2)+2.0d00*Y*Z_4)*R61 1673 & +(Y*(3.0d00*Y_2+6.0d00*Z_2)+12.0d00*Y*Z_2)*R51+9.0d00*Y*R41 1674 R(M,0,118)=Y_2*Z_5*R71+(Z_5+10.0d00* Y_2*Z_3)*R61 1675 & +(10.0d00*Z_3+15.0d00*Y_2*Z)*R51+15.0d00*Z*R41 1676 R(M,0,119)=Y*Z_6*R71+15.0d00*Y*Z_4*R61+ 1677 & 45.0d00*Y*Z_2*R51 1678 & +15.0d00*Y*R41 1679 R(M,0,120)=Z_7*R71+21.0d00*Z_5*R61+105.0d00*Z_3*R51 1680 & +105.0d00*Z*R41 1681 X_8=X_7 * X 1682 Y_8=Y_7 * Y 1683 Z_8=Z_7 * Z 1684 R81=R(m,8,1) 1685 R(M,0,121)=X_8*R81+28.0d00*X_6*R71+210.0d00*X_4*R61 1686 & +420.0d00*X_2*R51+105.0d00*R41 1687 R(M,0,122)=X_7*Y*R81+21.0d00*X_5*Y*R71 1688 & +105.0d00*X_3*Y*R61+105.0d00*X*Y*R51 1689 R(M,0,123)=X_7*Z*R81+21.0d00*X_5*Z*R71 1690 & +105.0d00*X_3*Z*R61+105.0d00*X*Z*R51 1691 R(M,0,124)=X_6*Y_2*R81 1692 & +(5.0d00*X_4*Y_2+X*(X*(3.0d00*X_2*Y_2+ 1693 & X*(X*(X_2+Y_2)+2.0d00*X*Y_2))+4.0d00*X_3*Y_2))*R71 1694 & +(X*(X*(6.0d00*X_2+3.0d00*Y_2)+4.0d00*X*(X_2+Y_2)+ 1695 & 8.0d00*X*Y_2)+ 1696 & 5.0d00*X*(X*(X_2+Y_2)+2.0d00*X*Y_2)+15.0d00*X_2*Y_2)*R61 1697 & +(45.0d00*X_2+15.0d00*Y_2)*R51+15.0d00*R41 1698 R(M,0,125)=X_6*Y*Z*R81+15.0d00*X_4*Y*Z*R71 1699 & +45.0d00*X_2*Y*Z*R61+15.0d00*Y*Z*R51 1700 R(M,0,126)=X_6*Z_2*R81 1701 & +(5.0d00*X_4*Z_2+X*(X*(3.0d00*X_2*Z_2+X*(X*(X_2+Z_2)+ 1702 & 2.0d00*X*Z_2))+4.0d00*X_3*Z_2))*R71 1703 & +(X*(X*(6.0d00*X_2+3.0d00*Z_2)+4.0d00*X*(X_2+Z_2)+ 1704 & 8.0d00*X*Z_2)+ 1705 & 5.0d00*X*(X*(X_2+Z_2)+2.0d00*X*Z_2)+15.0d00*X_2*Z_2)*R61 1706 & +(45.0d00*X_2+15.0d00*Z_2)*R51+15.0d00*R41 1707 R(M,0,127)=X_5*Y_3*R81 1708 & +(X*(3.0d00*X_2*Y_3+X*(X*(3.0d00*X_2*Y+Y_3)+2.0d00*X*Y_3))+ 1709 & 4.0d00*X_3*Y_3)*R71 1710 & +(X*(18.0d00*X_2*Y+3.0d00*Y_3)+4.0d00*X*(3.0d00*X_2*Y+Y_3)+ 1711 & 8.0d00*X*Y_3)*R61+45.0d00*X*Y*R51 1712 R(M,0,128)=X_5*Y_2*Z*R81 1713 & +(X*(3.0d00*X_2*Y_2*Z+X*(X*(X_2*Z+Y_2*Z)+2.0d00*X*Y_2*Z))+ 1714 & 4.0d00*X_3*Y_2*Z)*R71 1715 & +(X*(6.0d00*X_2*Z+3.0d00*Y_2*Z)+4.0d00*X*(X_2*Z+Y_2*Z)+ 1716 & 8.0d00*X*Y_2*Z)*R61+15.0d00*X*Z*R51 1717 R(M,0,129)=X_5*Y*Z_2*R81 1718 & +(X*(X*(X*(X_2*Y+Y*Z_2)+2.0d00*X*Y*Z_2)+3.0d00*X_2*Y*Z_2)+ 1719 & 4.0d00*X_3*Y*Z_2)*R71 1720 & +(X*(6.0d00*X_2*Y+3.0d00*Y*Z_2)+4.0d00*X*(X_2*Y+Y*Z_2)+ 1721 & 8.0d00*X*Y*Z_2)*R61+15.0d00*X*Y*R51 1722 R(M,0,130)=X_5*Z_3*R81 1723 & +(X*(3.0d00*X_2*Z_3+X*(X*(3.0d00*X_2*Z+Z_3)+2.0d00*X*Z_3))+ 1724 & 4.0d00*X_3*Z_3)*R71 1725 & +(X*(18.0d00*X_2*Z+3.0d00*Z_3)+4.0d00*X*(3.0d00*X_2*Z+Z_3)+ 1726 & 8.0d00*X*Z_3)*R61+45.0d00*X*Z*R51 1727 R(M,0,131)=X_4*Y_4*R81 1728 & +(3.0d00*X_2*Y_4+X*(X*(Y_4+6.0d00*X_2*Y_2)+ 1729 & 2.0d00*X*Y_4))*R71 1730 & +(X*(X*(3.0d00*X_2+6.0d00*Y_2)+12.0d00*X*Y_2)+ 1731 & 18.0d00*X_2*Y_2+3.0d00*Y_4)* 1732 & R61+(18.0d00*X_2+18.0d00*Y_2)*R51+9.0d00*R41 1733 R(M,0,132)=X_4*Y_3*Z*R81 1734 & +(3.0d00*X_2*Y_3*Z+X*(X*(3.0d00*X_2*Y*Z+Y_3*Z)+ 1735 & 2.0d00*X*Y_3*Z))* 1736 & R71+(18.0d00*X_2*Y*Z+3.0d00*Y_3*Z)*R61+9.0d00*Y*Z*R51 1737 R(M,0,133)=X_4*Y_2*Z_2*R81 1738 & +(3.0d00*X_2*Y_2*Z_2+X*(X*(Y_2*Z_2+X_2*(Y_2+Z_2))+ 1739 & 2.0d00*X*Y_2*Z_2))*R71 1740 & +(X*(X*(X_2+Y_2+Z_2)+2.0d00*X*(Y_2+Z_2))+ 1741 & 3.0d00*X_2*(Y_2+Z_2)+3.0d00*Y_2*Z_2)*R61 1742 & +(6.0d00*X_2+3.0d00*Y_2+3.0d00*Z_2)*R51+3.0d00*R41 1743 R(M,0,134)=X_4*Y*Z_3*R81 1744 & +(X*(X*(3.0d00*X_2*Y*Z+Y*Z_3)+2.0d00*X*Y*Z_3)+ 1745 & 3.0d00*X_2*Y*Z_3)* 1746 & R71+(18.0d00*X_2*Y*Z+3.0d00*Y*Z_3)*R61+9.0d00*Y*Z*R51 1747 R(M,0,135)=X_4*Z_4*R81 1748 & +(3.0d00*X_2*Z_4+X*(X*(Z_4+6.0d00*X_2*Z_2)+ 1749 & 2.0d00*X*Z_4))*R71 1750 & +(X*(X*(3.0d00*X_2+6.0d00*Z_2)+12.0d00*X*Z_2)+ 1751 & 18.0d00*X_2*Z_2+3.0d00*Z_4)* 1752 & R61+(18.0d00*X_2+18.0d00*Z_2)*R51+9.0d00*R41 1753 R(M,0,136)=X_3*Y_5*R81 1754 & +(X*(Y_5+10.0d00* X_2*Y_3)+2.0d00*X*Y_5)*R71 1755 & +(X*(10.0d00*Y_3+15.0d00*X_2*Y)+20.0d00*X*Y_3)*R61+ 1756 & 45.0d00*X*Y*R51 1757 R(M,0,137)=X_3*Y_4*Z*R81 1758 & +(X*(Y_4*Z+6.0d00*X_2*Y_2*Z)+2.0d00*X*Y_4*Z)*R71 1759 & +(X*(3.0d00*X_2*Z+6.0d00*Y_2*Z)+12.0d00*X*Y_2*Z)*R61+ 1760 & 9.0d00*X*Z*R51 1761 R(M,0,138)=X_3*Y_3*Z_2*R81 1762 & +(X*(X_2*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+Y_3*Z_2)+ 1763 & 2.0d00*X*Y_3*Z_2)*R71 1764 & +(X*(3.0d00*X_2*Y+Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+ 1765 & 2.0d00*X*(Y*(Y_2+Z_2)+ 1766 & 2.0d00*Y*Z_2))*R61+9.0d00*X*Y*R51 1767 R(M,0,139)=X_3*Y_2*Z_3*R81 1768 & +(X*(Y_2*Z_3+X_2*(3.0d00*Y_2*Z+Z_3))+2.0d00*X*Y_2*Z_3)* 1769 & R71+(X*(3.0d00*X_2*Z+3.0d00*Y_2*Z+Z_3)+ 1770 & 2.0d00*X*(3.0d00*Y_2*Z+Z_3))* 1771 & R61+9.0d00*X*Z*R51 1772 R(M,0,140)=X_3*Y*Z_4*R81 1773 & +(X*(6.0d00*X_2*Y*Z_2+Y*Z_4)+2.0d00*X*Y*Z_4)*R71 1774 & +(X*(3.0d00*X_2*Y+6.0d00*Y*Z_2)+12.0d00*X*Y*Z_2)*R61+ 1775 & 9.0d00*X*Y*R51 1776 R(M,0,141)=X_3*Z_5*R81 1777 & +(X*(Z_5+10.0d00* X_2*Z_3)+2.0d00*X*Z_5)*R71 1778 & +(X*(10.0d00*Z_3+15.0d00*X_2*Z)+20.0d00*X*Z_3)*R61+ 1779 & 45.0d00*X*Z*R51 1780 R(M,0,142)=X_2*Y_6*R81+(Y_6+15.0d00*X_2*Y_4)*R71 1781 & +(45.0d00*X_2*Y_2+15.0d00*Y_4)*R61+ 1782 & (15.0d00*X_2+45.0d00*Y_2)*R51+15.0d00*R41 1783 R(M,0,143)=X_2*Y_5*Z*R81+(Y_5*Z+10.0d00* X_2*Y_3*Z)* 1784 & R71+(10.0d00*Y_3*Z+15.0d00*X_2*Y*Z)*R61+15.0d00*Y*Z*R51 1785 R(M,0,144)=X_2*Y_4*Z_2*R81 1786 & +(Y_4*Z_2+X_2*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 1787 & 2.0d00*Y*Z_2)))*R71+ 1788 & (X_2*(6.0d00*Y_2+3.0d00*Z_2)+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+ 1789 & 3.0d00*Y_2*Z_2)*R61+(3.0d00*X_2+6.0d00*Y_2+3.0d00*Z_2)*R51 1790 & +3.0d00*R41 1791 R(M,0,145)=X_2*Y_3*Z_3*R81 1792 & +(X_2*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)+Y_3*Z_3)*R71 1793 & +(9.0d00*X_2*Y*Z+Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)*R61 1794 & +9.0d00*Y*Z*R51 1795 R(M,0,146)=X_2*Y_2*Z_4*R81 1796 & +(Y_2*Z_4+X_2*(Z_4+6.0d00*Y_2*Z_2))*R71 1797 & +(X_2*(3.0d00*Y_2+6.0d00*Z_2)+6.0d00*Y_2*Z_2+Z_4)*R61 1798 & +(3.0d00*X_2+3.0d00*Y_2+6.0d00*Z_2)*R51+3.0d00*R41 1799 R(M,0,147)=X_2*Y*Z_5*R81+(10.0d00*X_2*Y*Z_3+Y*Z_5)* 1800 & R71+(15.0d00*X_2*Y*Z+10.0d00* Y*Z_3)*R61+15.0d00*Y*Z*R51 1801 R(M,0,148)=X_2*Z_6*R81+(Z_6+15.0d00*X_2*Z_4)*R71 1802 & +(45.0d00*X_2*Z_2+15.0d00*Z_4)*R61+(15.0d00*X_2+ 1803 & 45.0d00*Z_2)* 1804 & R51+15.0d00*R41 1805 R(M,0,149)=X*Y_7*R81+21.0d00*X*Y_5*R71 1806 & +105.0d00*X*Y_3*R61+105.0d00*X*Y*R51 1807 R(M,0,150)=X*Y_6*Z*R81+15.0d00*X*Y_4*Z*R71 1808 & +45.0d00*X*Y_2*Z*R61+15.0d00*X*Z*R51 1809 R(M,0,151)=X*Y_5*Z_2*R81 1810 & +X*(Y*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))+ 1811 & 4.0d00*Y_3*Z_2)*R71 1812 & +X*(Y*(6.0d00*Y_2+3.0d00*Z_2)+4.0d00*Y*(Y_2+Z_2)+ 1813 & 8.0d00*Y*Z_2)*R61 1814 & +15.0d00*X*Y*R51 1815 R(M,0,152)=X*Y_4*Z_3*R81 1816 & +X*(3.0d00*Y_2*Z_3+Y*(Y*(3.0d00*Y_2*Z+Z_3)+ 1817 & 2.0d00*Y*Z_3))*R71 1818 & +X*(18.0d00*Y_2*Z+3.0d00*Z_3)*R61+9.0d00*X*Z*R51 1819 R(M,0,153)=X*Y_3*Z_4*R81 1820 & +X*(Y*(Z_4+6.0d00*Y_2*Z_2)+2.0d00*Y*Z_4)*R71 1821 & +X*(Y*(3.0d00*Y_2+6.0d00*Z_2)+12.0d00*Y*Z_2)*R61+ 1822 & 9.0d00*X*Y*R51 1823 R(M,0,154)=X*Y_2*Z_5*R81+X*(Z_5+10.0d00* Y_2*Z_3)*R71 1824 & +X*(10.0d00*Z_3+15.0d00*Y_2*Z)*R61+15.0d00*X*Z*R51 1825 R(M,0,155)=X*Y*Z_6*R81+15.0d00*X*Y*Z_4*R71 1826 & +45.0d00*X*Y*Z_2*R61+15.0d00*X*Y*R51 1827 R(M,0,156)=X*Z_7*R81+21.0d00*X*Z_5*R71 1828 & +105.0d00*X*Z_3*R61+105.0d00*X*Z*R51 1829 R(M,0,157)=Y_8*R81+28.0d00*Y_6*R71+210.0d00*Y_4*R61 1830 & +420.0d00*Y_2*R51+105.0d00*R41 1831 R(M,0,158)=Y_7*Z*R81+21.0d00*Y_5*Z*R71 1832 & +105.0d00*Y_3*Z*R61+105.0d00*Y*Z*R51 1833 R(M,0,159)=Y_6*Z_2*R81 1834 & +(5.0d00*Y_4*Z_2+Y*(Y*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 1835 & 2.0d00*Y*Z_2))+4.0d00*Y_3*Z_2))*R71 1836 & +(Y*(Y*(6.0d00*Y_2+3.0d00*Z_2)+4.0d00*Y*(Y_2+Z_2)+ 1837 & 8.0d00*Y*Z_2)+ 1838 & 5.0d00*Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+15.0d00*Y_2*Z_2)*R61 1839 & +(45.0d00*Y_2+15.0d00*Z_2)*R51+15.0d00*R41 1840 R(M,0,160)=Y_5*Z_3*R81 1841 & +(Y*(3.0d00*Y_2*Z_3+Y*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3))+ 1842 & 4.0d00*Y_3*Z_3)*R71 1843 & +(Y*(18.0d00*Y_2*Z+3.0d00*Z_3)+4.0d00*Y*(3.0d00*Y_2*Z+Z_3)+ 1844 & 8.0d00*Y*Z_3)* 1845 & R61+45.0d00*Y*Z*R51 1846 R(M,0,161)=Y_4*Z_4*R81 1847 & +(3.0d00*Y_2*Z_4+Y*(Y*(Z_4+6.0d00*Y_2*Z_2)+ 1848 & 2.0d00*Y*Z_4))*R71 1849 & +(Y*(Y*(3.0d00*Y_2+6.0d00*Z_2)+12.0d00*Y*Z_2)+ 1850 & 18.0d00*Y_2*Z_2+3.0d00*Z_4)* 1851 & R61+(18.0d00*Y_2+18.0d00*Z_2)*R51+9.0d00*R41 1852 R(M,0,162)=Y_3*Z_5*R81 1853 & +(Y*(Z_5+10.0d00* Y_2*Z_3)+2.0d00*Y*Z_5)*R71 1854 & +(Y*(10.0d00*Z_3+15.0d00*Y_2*Z)+20.0d00*Y*Z_3)*R61+ 1855 & 45.0d00*Y*Z*R51 1856 R(M,0,163)=Y_2*Z_6*R81+(Z_6+15.0d00*Y_2*Z_4)*R71 1857 & +(45.0d00*Y_2*Z_2+15.0d00*Z_4)*R61+ 1858 & (15.0d00*Y_2+45.0d00*Z_2)*R51 1859 & +15.0d00*R41 1860 R(M,0,164)=Y*Z_7*R81+21.0d00*Y*Z_5*R71 1861 & +105.0d00*Y*Z_3*R61+105.0d00*Y*Z*R51 1862 R(M,0,165)=Z_8*R81+28.0d00*Z_6*R71+210.0d00*Z_4*R61 1863 & +420.0d00*Z_2*R51+105.0d00*R41 186400100 continue 1865c 1866 end 1867 subroutine hfmkr_9(R,PQ,IJK,NAB,Lr,Lr3) 1868 implicit none 1869c $Id$ 1870c 1871 integer NAB, Lr, Lr3 1872 integer IJK(0:Lr,0:Lr,0:Lr) 1873 double precision PQ(NAB,3) 1874 double precision R(NAB,0:Lr,Lr3) 1875c 1876 integer m 1877 double precision x,y,z 1878 double precision x_2, y_2, z_2 1879 double precision x_3, y_3, z_3 1880 double precision x_4, y_4, z_4 1881 double precision x_5, y_5, z_5 1882 double precision x_6, y_6, z_6 1883 double precision x_7, y_7, z_7 1884 double precision x_8, y_8, z_8 1885 double precision R11,R21,R31,R41,R51,R61,R71,R81,R91 1886c 1887 IJK(0,0,0)=1 1888 IJK(1,0,0)=2 1889 IJK(0,1,0)=3 1890 IJK(0,0,1)=4 1891 IJK(2,0,0)=5 1892 IJK(1,1,0)=6 1893 IJK(1,0,1)=7 1894 IJK(0,2,0)=8 1895 IJK(0,1,1)=9 1896 IJK(0,0,2)=10 1897 IJK(3,0,0)=11 1898 IJK(2,1,0)=12 1899 IJK(2,0,1)=13 1900 IJK(1,2,0)=14 1901 IJK(1,1,1)=15 1902 IJK(1,0,2)=16 1903 IJK(0,3,0)=17 1904 IJK(0,2,1)=18 1905 IJK(0,1,2)=19 1906 IJK(0,0,3)=20 1907 IJK(4,0,0)=21 1908 IJK(3,1,0)=22 1909 IJK(3,0,1)=23 1910 IJK(2,2,0)=24 1911 IJK(2,1,1)=25 1912 IJK(2,0,2)=26 1913 IJK(1,3,0)=27 1914 IJK(1,2,1)=28 1915 IJK(1,1,2)=29 1916 IJK(1,0,3)=30 1917 IJK(0,4,0)=31 1918 IJK(0,3,1)=32 1919 IJK(0,2,2)=33 1920 IJK(0,1,3)=34 1921 IJK(0,0,4)=35 1922 IJK(5,0,0)=36 1923 IJK(4,1,0)=37 1924 IJK(4,0,1)=38 1925 IJK(3,2,0)=39 1926 IJK(3,1,1)=40 1927 IJK(3,0,2)=41 1928 IJK(2,3,0)=42 1929 IJK(2,2,1)=43 1930 IJK(2,1,2)=44 1931 IJK(2,0,3)=45 1932 IJK(1,4,0)=46 1933 IJK(1,3,1)=47 1934 IJK(1,2,2)=48 1935 IJK(1,1,3)=49 1936 IJK(1,0,4)=50 1937 IJK(0,5,0)=51 1938 IJK(0,4,1)=52 1939 IJK(0,3,2)=53 1940 IJK(0,2,3)=54 1941 IJK(0,1,4)=55 1942 IJK(0,0,5)=56 1943 IJK(6,0,0)=57 1944 IJK(5,1,0)=58 1945 IJK(5,0,1)=59 1946 IJK(4,2,0)=60 1947 IJK(4,1,1)=61 1948 IJK(4,0,2)=62 1949 IJK(3,3,0)=63 1950 IJK(3,2,1)=64 1951 IJK(3,1,2)=65 1952 IJK(3,0,3)=66 1953 IJK(2,4,0)=67 1954 IJK(2,3,1)=68 1955 IJK(2,2,2)=69 1956 IJK(2,1,3)=70 1957 IJK(2,0,4)=71 1958 IJK(1,5,0)=72 1959 IJK(1,4,1)=73 1960 IJK(1,3,2)=74 1961 IJK(1,2,3)=75 1962 IJK(1,1,4)=76 1963 IJK(1,0,5)=77 1964 IJK(0,6,0)=78 1965 IJK(0,5,1)=79 1966 IJK(0,4,2)=80 1967 IJK(0,3,3)=81 1968 IJK(0,2,4)=82 1969 IJK(0,1,5)=83 1970 IJK(0,0,6)=84 1971 IJK(7,0,0)=85 1972 IJK(6,1,0)=86 1973 IJK(6,0,1)=87 1974 IJK(5,2,0)=88 1975 IJK(5,1,1)=89 1976 IJK(5,0,2)=90 1977 IJK(4,3,0)=91 1978 IJK(4,2,1)=92 1979 IJK(4,1,2)=93 1980 IJK(4,0,3)=94 1981 IJK(3,4,0)=95 1982 IJK(3,3,1)=96 1983 IJK(3,2,2)=97 1984 IJK(3,1,3)=98 1985 IJK(3,0,4)=99 1986 IJK(2,5,0)=100 1987 IJK(2,4,1)=101 1988 IJK(2,3,2)=102 1989 IJK(2,2,3)=103 1990 IJK(2,1,4)=104 1991 IJK(2,0,5)=105 1992 IJK(1,6,0)=106 1993 IJK(1,5,1)=107 1994 IJK(1,4,2)=108 1995 IJK(1,3,3)=109 1996 IJK(1,2,4)=110 1997 IJK(1,1,5)=111 1998 IJK(1,0,6)=112 1999 IJK(0,7,0)=113 2000 IJK(0,6,1)=114 2001 IJK(0,5,2)=115 2002 IJK(0,4,3)=116 2003 IJK(0,3,4)=117 2004 IJK(0,2,5)=118 2005 IJK(0,1,6)=119 2006 IJK(0,0,7)=120 2007 IJK(8,0,0)=121 2008 IJK(7,1,0)=122 2009 IJK(7,0,1)=123 2010 IJK(6,2,0)=124 2011 IJK(6,1,1)=125 2012 IJK(6,0,2)=126 2013 IJK(5,3,0)=127 2014 IJK(5,2,1)=128 2015 IJK(5,1,2)=129 2016 IJK(5,0,3)=130 2017 IJK(4,4,0)=131 2018 IJK(4,3,1)=132 2019 IJK(4,2,2)=133 2020 IJK(4,1,3)=134 2021 IJK(4,0,4)=135 2022 IJK(3,5,0)=136 2023 IJK(3,4,1)=137 2024 IJK(3,3,2)=138 2025 IJK(3,2,3)=139 2026 IJK(3,1,4)=140 2027 IJK(3,0,5)=141 2028 IJK(2,6,0)=142 2029 IJK(2,5,1)=143 2030 IJK(2,4,2)=144 2031 IJK(2,3,3)=145 2032 IJK(2,2,4)=146 2033 IJK(2,1,5)=147 2034 IJK(2,0,6)=148 2035 IJK(1,7,0)=149 2036 IJK(1,6,1)=150 2037 IJK(1,5,2)=151 2038 IJK(1,4,3)=152 2039 IJK(1,3,4)=153 2040 IJK(1,2,5)=154 2041 IJK(1,1,6)=155 2042 IJK(1,0,7)=156 2043 IJK(0,8,0)=157 2044 IJK(0,7,1)=158 2045 IJK(0,6,2)=159 2046 IJK(0,5,3)=160 2047 IJK(0,4,4)=161 2048 IJK(0,3,5)=162 2049 IJK(0,2,6)=163 2050 IJK(0,1,7)=164 2051 IJK(0,0,8)=165 2052 IJK(9,0,0)=166 2053 IJK(8,1,0)=167 2054 IJK(8,0,1)=168 2055 IJK(7,2,0)=169 2056 IJK(7,1,1)=170 2057 IJK(7,0,2)=171 2058 IJK(6,3,0)=172 2059 IJK(6,2,1)=173 2060 IJK(6,1,2)=174 2061 IJK(6,0,3)=175 2062 IJK(5,4,0)=176 2063 IJK(5,3,1)=177 2064 IJK(5,2,2)=178 2065 IJK(5,1,3)=179 2066 IJK(5,0,4)=180 2067 IJK(4,5,0)=181 2068 IJK(4,4,1)=182 2069 IJK(4,3,2)=183 2070 IJK(4,2,3)=184 2071 IJK(4,1,4)=185 2072 IJK(4,0,5)=186 2073 IJK(3,6,0)=187 2074 IJK(3,5,1)=188 2075 IJK(3,4,2)=189 2076 IJK(3,3,3)=190 2077 IJK(3,2,4)=191 2078 IJK(3,1,5)=192 2079 IJK(3,0,6)=193 2080 IJK(2,7,0)=194 2081 IJK(2,6,1)=195 2082 IJK(2,5,2)=196 2083 IJK(2,4,3)=197 2084 IJK(2,3,4)=198 2085 IJK(2,2,5)=199 2086 IJK(2,1,6)=200 2087 IJK(2,0,7)=201 2088 IJK(1,8,0)=202 2089 IJK(1,7,1)=203 2090 IJK(1,6,2)=204 2091 IJK(1,5,3)=205 2092 IJK(1,4,4)=206 2093 IJK(1,3,5)=207 2094 IJK(1,2,6)=208 2095 IJK(1,1,7)=209 2096 IJK(1,0,8)=210 2097 IJK(0,9,0)=211 2098 IJK(0,8,1)=212 2099 IJK(0,7,2)=213 2100 IJK(0,6,3)=214 2101 IJK(0,5,4)=215 2102 IJK(0,4,5)=216 2103 IJK(0,3,6)=217 2104 IJK(0,2,7)=218 2105 IJK(0,1,8)=219 2106 IJK(0,0,9)=220 2107c 2108 do 00100 m=1,NAB 2109 X=PQ(m,1) 2110 Y=PQ(m,2) 2111 Z=PQ(m,3) 2112 R11=R(m,1,1) 2113 R(M,0,2)=X*R11 2114 R(M,0,3)=Y*R11 2115 R(M,0,4)=Z*R11 2116 X_2=X*X 2117 Y_2=Y*Y 2118 Z_2=Z*Z 2119 R21=R(m,2,1) 2120 R(M,0,5)=X_2*R21+R11 2121 R(M,0,6)=X*Y*R21 2122 R(M,0,7)=X*Z*R21 2123 R(M,0,8)=Y_2*R21+R11 2124 R(M,0,9)=Y*Z*R21 2125 R(M,0,10)=Z_2*R21+R11 2126 X_3=X_2*X 2127 Y_3=Y_2*Y 2128 Z_3=Z_2*Z 2129 R31=R(m,3,1) 2130 R(M,0,11)=X_3*R31+3.0d00*X*R21 2131 R(M,0,12)=X_2*Y*R31+Y*R21 2132 R(M,0,13)=X_2*Z*R31+Z*R21 2133 R(M,0,14)=X*Y_2*R31+X*R21 2134 R(M,0,15)=X*Y*Z*R31 2135 R(M,0,16)=X*Z_2*R31+X*R21 2136 R(M,0,17)=Y_3*R31+3.0d00*Y*R21 2137 R(M,0,18)=Y_2*Z*R31+Z*R21 2138 R(M,0,19)=Y*Z_2*R31+Y*R21 2139 R(M,0,20)=Z_3*R31+3.0d00*Z*R21 2140 X_4=X_3*X 2141 Y_4=Y_3*Y 2142 Z_4=Z_3*Z 2143 R41=R(m,4,1) 2144 R(M,0,21)=X_4*R41+6.0d00*X_2*R31+3.0d00*R21 2145 R(M,0,22)=X_3*Y*R41+3.0d00*X*Y*R31 2146 R(M,0,23)=X_3*Z*R41+3.0d00*X*Z*R31 2147 R(M,0,24)=X_2*Y_2*R41+(X_2+Y_2)*R31+R21 2148 R(M,0,25)=X_2*Y*Z*R41+Y*Z*R31 2149 R(M,0,26)=X_2*Z_2*R41+(X_2+Z_2)*R31+R21 2150 R(M,0,27)=X*Y_3*R41+3.0d00*X*Y*R31 2151 R(M,0,28)=X*Y_2*Z*R41+X*Z*R31 2152 R(M,0,29)=X*Y*Z_2*R41+X*Y*R31 2153 R(M,0,30)=X*Z_3*R41+3.0d00*X*Z*R31 2154 R(M,0,31)=Y_4*R41+6.0d00*Y_2*R31+3.0d00*R21 2155 R(M,0,32)=Y_3*Z*R41+3.0d00*Y*Z*R31 2156 R(M,0,33)=Y_2*Z_2*R41+(Y_2+Z_2)*R31+R21 2157 R(M,0,34)=Y*Z_3*R41+3.0d00*Y*Z*R31 2158 R(M,0,35)=Z_4*R41+6.0d00*Z_2*R31+3.0d00*R21 2159 X_5=X_4*X 2160 Y_5=Y_4*Y 2161 Z_5=Z_4*Z 2162 R51=R(m,5,1) 2163 R(M,0,36)=X_5*R51+10.0d00* X_3*R41+15.0d00*X*R31 2164 R(M,0,37)=X_4*Y*R51+6.0d00*X_2*Y*R41+3.0d00*Y*R31 2165 R(M,0,38)=X_4*Z*R51+6.0d00*X_2*Z*R41+3.0d00*Z*R31 2166 R(M,0,39)=X_3*Y_2*R51+(X*(X_2+Y_2)+2.0d00*X*Y_2)*R41 2167 & +3.0d00*X*R31 2168 R(M,0,40)=X_3*Y*Z*R51+3.0d00*X*Y*Z*R41 2169 R(M,0,41)=X_3*Z_2*R51+(X*(X_2+Z_2)+2.0d00*X*Z_2)*R41 2170 & +3.0d00*X*R31 2171 R(M,0,42)=X_2*Y_3*R51+(3.0d00*X_2*Y+Y_3)*R41 2172 & +3.0d00*Y*R31 2173 R(M,0,43)=X_2*Y_2*Z*R51+(X_2*Z+Y_2*Z)*R41 2174 & +Z*R31 2175 R(M,0,44)=X_2*Y*Z_2*R51+(X_2*Y+Y*Z_2)*R41 2176 & +Y*R31 2177 R(M,0,45)=X_2*Z_3*R51+(3.0d00*X_2*Z+Z_3)*R41 2178 & +3.0d00*Z*R31 2179 R(M,0,46)=X*Y_4*R51+6.0d00*X*Y_2*R41+3.0d00*X*R31 2180 R(M,0,47)=X*Y_3*Z*R51+3.0d00*X*Y*Z*R41 2181 R(M,0,48)=X*Y_2*Z_2*R51+X*(Y_2+Z_2)*R41+ 2182 & X*R31 2183 R(M,0,49)=X*Y*Z_3*R51+3.0d00*X*Y*Z*R41 2184 R(M,0,50)=X*Z_4*R51+6.0d00*X*Z_2*R41+3.0d00*X*R31 2185 R(M,0,51)=Y_5*R51+10.0d00* Y_3*R41+15.0d00*Y*R31 2186 R(M,0,52)=Y_4*Z*R51+6.0d00*Y_2*Z*R41+3.0d00*Z*R31 2187 R(M,0,53)=Y_3*Z_2*R51+(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R41 2188 & +3.0d00*Y*R31 2189 R(M,0,54)=Y_2*Z_3*R51+(3.0d00*Y_2*Z+Z_3)*R41 2190 & +3.0d00*Z*R31 2191 R(M,0,55)=Y*Z_4*R51+6.0d00*Y*Z_2*R41+3.0d00*Y*R31 2192 R(M,0,56)=Z_5*R51+10.0d00* Z_3*R41+15.0d00*Z*R31 2193 X_6=X_5 * X 2194 Y_6=Y_5 * Y 2195 Z_6=Z_5 * Z 2196 R61=R(m,6,1) 2197 R(M,0,57)=X_6*R61+15.0d00*X_4*R51+45.0d00*X_2*R41 2198 & +15.0d00*R31 2199 R(M,0,58)=X_5*Y*R61+10.0d00* X_3*Y*R51+15.0d00*X*Y*R41 2200 R(M,0,59)=X_5*Z*R61+10.0d00* X_3*Z*R51+15.0d00*X*Z*R41 2201 R(M,0,60)=X_4*Y_2*R61 2202 & +(3.0d00*X_2*Y_2+X*(X*(X_2+Y_2)+2.0d00*X*Y_2))*R51 2203 & +(6.0d00*X_2+3.0d00*Y_2)*R41+3.0d00*R31 2204 R(M,0,61)=X_4*Y*Z*R61+6.0d00*X_2*Y*Z*R51+3.0d00*Y*Z*R41 2205 R(M,0,62)=X_4*Z_2*R61 2206 & +(3.0d00*X_2*Z_2+X*(X*(X_2+Z_2)+2.0d00*X*Z_2))*R51 2207 & +(6.0d00*X_2+3.0d00*Z_2)*R41+3.0d00*R31 2208 R(M,0,63)=X_3*Y_3*R61 2209 & +(X*(3.0d00*X_2*Y+Y_3)+2.0d00*X*Y_3)*R51+9.0d00*X*Y*R41 2210 R(M,0,64)=X_3*Y_2*Z*R61 2211 & +(X*(X_2*Z+Y_2*Z)+2.0d00*X*Y_2*Z)*R51+3.0d00*X*Z*R41 2212 R(M,0,65)=X_3*Y*Z_2*R61 2213 & +(X*(X_2*Y+Y*Z_2)+2.0d00*X*Y*Z_2)*R51+3.0d00*X*Y*R41 2214 R(M,0,66)=X_3*Z_3*R61 2215 & +(X*(3.0d00*X_2*Z+Z_3)+2.0d00*X*Z_3)*R51+9.0d00*X*Z*R41 2216 R(M,0,67)=X_2*Y_4*R61+(Y_4+6.0d00*X_2*Y_2)*R51 2217 & +(3.0d00*X_2+6.0d00*Y_2)*R41+3.0d00*R31 2218 R(M,0,68)=X_2*Y_3*Z*R61+(3.0d00*X_2*Y*Z+Y_3*Z)*R51 2219 & +3.0d00*Y*Z*R41 2220 R(M,0,69)=X_2*Y_2*Z_2*R61 2221 & +(Y_2*Z_2+X_2*(Y_2+Z_2))*R51 2222 & +(X_2+Y_2+Z_2)*R41+R31 2223 R(M,0,70)=X_2*Y*Z_3*R61+(3.0d00*X_2*Y*Z+Y*Z_3)*R51 2224 & +3.0d00*Y*Z*R41 2225 R(M,0,71)=X_2*Z_4*R61+(Z_4+6.0d00*X_2*Z_2)*R51 2226 & +(3.0d00*X_2+6.0d00*Z_2)*R41+3.0d00*R31 2227 R(M,0,72)=X*Y_5*R61+10.0d00* X*Y_3*R51+15.0d00*X*Y*R41 2228 R(M,0,73)=X*Y_4*Z*R61+6.0d00*X*Y_2*Z*R51+3.0d00*X*Z*R41 2229 R(M,0,74)=X*Y_3*Z_2*R61 2230 & +X*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R51+3.0d00*X*Y*R41 2231 R(M,0,75)=X*Y_2*Z_3*R61+X*(3.0d00*Y_2*Z+Z_3)*R51 2232 & +3.0d00*X*Z*R41 2233 R(M,0,76)=X*Y*Z_4*R61+6.0d00*X*Y*Z_2*R51+3.0d00*X*Y*R41 2234 R(M,0,77)=X*Z_5*R61+10.0d00* X*Z_3*R51+15.0d00*X*Z*R41 2235 R(M,0,78)=Y_6*R61+15.0d00*Y_4*R51+45.0d00*Y_2*R41 2236 & +15.0d00*R31 2237 R(M,0,79)=Y_5*Z*R61+10.0d00* Y_3*Z*R51+15.0d00*Y*Z*R41 2238 R(M,0,80)=Y_4*Z_2*R61 2239 & +(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))*R51 2240 & +(6.0d00*Y_2+3.0d00*Z_2)*R41+3.0d00*R31 2241 R(M,0,81)=Y_3*Z_3*R61 2242 & +(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)*R51+9.0d00*Y*Z*R41 2243 R(M,0,82)=Y_2*Z_4*R61+(Z_4+6.0d00*Y_2*Z_2)*R51 2244 & +(3.0d00*Y_2+6.0d00*Z_2)*R41+3.0d00*R31 2245 R(M,0,83)=Y*Z_5*R61+10.0d00* Y*Z_3*R51+15.0d00*Y*Z*R41 2246 R(M,0,84)=Z_6*R61+15.0d00*Z_4*R51+45.0d00*Z_2*R41 2247 & +15.0d00*R31 2248 X_7=X_6 * X 2249 Y_7=Y_6 * Y 2250 Z_7=Z_6 * Z 2251 R71=R(m,7,1) 2252 R(M,0,85)=X_7*R71+21.0d00*X_5*R61+105.0d00*X_3*R51 2253 & +105.0d00*X*R41 2254 R(M,0,86)=X_6*Y*R71+15.0d00*X_4*Y*R61+45.0d00*X_2*Y*R51 2255 & +15.0d00*Y*R41 2256 R(M,0,87)=X_6*Z*R71+15.0d00*X_4*Z*R61+45.0d00*X_2*Z*R51 2257 & +15.0d00*Z*R41 2258 R(M,0,88)=X_5*Y_2*R71 2259 & +(X*(3.0d00*X_2*Y_2+X*(X*(X_2+Y_2)+ 2260 & 2.0d00*X*Y_2))+4.0d00*X_3*Y_2)*R61 2261 & +(X*(6.0d00*X_2+3.0d00*Y_2)+ 2262 & 4.0d00*X*(X_2+Y_2)+8.0d00*X*Y_2)*R51 2263 & +15.0d00*X*R41 2264 R(M,0,89)=X_5*Y*Z*R71+10.0d00* X_3*Y*Z*R61 2265 & +15.0d00*X*Y*Z*R51 2266 R(M,0,90)=X_5*Z_2*R71 2267 & +(X*(3.0d00*X_2*Z_2+X*(X*(X_2+Z_2)+ 2268 & 2.0d00*X*Z_2))+4.0d00*X_3*Z_2)*R61 2269 & +(X*(6.0d00*X_2+3.0d00*Z_2)+ 2270 & 4.0d00*X*(X_2+Z_2)+8.0d00*X*Z_2)*R51 2271 & +15.0d00*X*R41 2272 R(M,0,91)=X_4*Y_3*R71 2273 & +(3.0d00*X_2*Y_3+X*(X*(3.0d00*X_2*Y+Y_3)+2.0d00*X*Y_3))*R61 2274 & +(18.0d00*X_2*Y+3.0d00*Y_3)*R51+9.0d00*Y*R41 2275 R(M,0,92)=X_4*Y_2*Z*R71 2276 & +(3.0d00*X_2*Y_2*Z+X*(X*(X_2*Z+Y_2*Z)+2.0d00*X*Y_2*Z))*R61 2277 & +(6.0d00*X_2*Z+3.0d00*Y_2*Z)*R51+3.0d00*Z*R41 2278 R(M,0,93)=X_4*Y*Z_2*R71 2279 & +(X*(X*(X_2*Y+Y*Z_2)+2.0d00*X*Y*Z_2)+3.0d00*X_2*Y*Z_2)*R61 2280 & +(6.0d00*X_2*Y+3.0d00*Y*Z_2)*R51+3.0d00*Y*R41 2281 R(M,0,94)=X_4*Z_3*R71 2282 & +(3.0d00*X_2*Z_3+X*(X*(3.0d00*X_2*Z+Z_3)+2.0d00*X*Z_3))*R61 2283 & +(18.0d00*X_2*Z+3.0d00*Z_3)*R51+9.0d00*Z*R41 2284 R(M,0,95)=X_3*Y_4*R71 2285 & +(X*(Y_4+6.0d00*X_2*Y_2)+2.0d00*X*Y_4)*R61 2286 & +(X*(3.0d00*X_2+6.0d00*Y_2)+12.0d00*X*Y_2)*R51+9.0d00*X*R41 2287 R(M,0,96)=X_3*Y_3*Z*R71 2288 & +(X*(3.0d00*X_2*Y*Z+Y_3*Z)+2.0d00*X*Y_3*Z)*R61 2289 & +9.0d00*X*Y*Z*R51 2290 R(M,0,97)=X_3*Y_2*Z_2*R71 2291 & +(X*(Y_2*Z_2+X_2*(Y_2+Z_2))+2.0d00*X*Y_2*Z_2)*R61 2292 & +(X*(X_2+Y_2+Z_2)+2.0d00*X*(Y_2+Z_2))*R51+ 2293 & 3.0d00*X*R41 2294 R(M,0,98)=X_3*Y*Z_3*R71 2295 & +(X*(3.0d00*X_2*Y*Z+Y*Z_3)+2.0d00*X*Y*Z_3)*R61 2296 & +9.0d00*X*Y*Z*R51 2297 R(M,0,99)=X_3*Z_4*R71 2298 & +(X*(Z_4+6.0d00*X_2*Z_2)+2.0d00*X*Z_4)*R61 2299 & +(X*(3.0d00*X_2+6.0d00*Z_2)+12.0d00*X*Z_2)*R51+9.0d00*X*R41 2300 R(M,0,100)=X_2*Y_5*R71+(Y_5+10.0d00* X_2*Y_3)*R61 2301 & +(10.0d00*Y_3+15.0d00*X_2*Y)*R51+15.0d00*Y*R41 2302 R(M,0,101)=X_2*Y_4*Z*R71+ 2303 & (Y_4*Z+6.0d00*X_2*Y_2*Z)*R61 2304 & +(3.0d00*X_2*Z+6.0d00*Y_2*Z)*R51+3.0d00*Z*R41 2305 R(M,0,102)=X_2*Y_3*Z_2*R71 2306 & +(X_2*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+Y_3*Z_2)*R61 2307 & +(3.0d00*X_2*Y+Y*(Y_2+Z_2)+2.0d00*Y*Z_2)*R51+3.0d00*Y*R41 2308 R(M,0,103)=X_2*Y_2*Z_3*R71 2309 & +(Y_2*Z_3+X_2*(3.0d00*Y_2*Z+Z_3))*R61 2310 & +(3.0d00*X_2*Z+3.0d00*Y_2*Z+Z_3)*R51+3.0d00*Z*R41 2311 R(M,0,104)=X_2*Y*Z_4*R71+ 2312 & (6.0d00*X_2*Y*Z_2+Y*Z_4)*R61 2313 & +(3.0d00*X_2*Y+6.0d00*Y*Z_2)*R51+3.0d00*Y*R41 2314 R(M,0,105)=X_2*Z_5*R71+(Z_5+10.0d00* X_2*Z_3)*R61 2315 & +(10.0d00*Z_3+15.0d00*X_2*Z)*R51+15.0d00*Z*R41 2316 R(M,0,106)=X*Y_6*R71+15.0d00*X*Y_4*R61+ 2317 & 45.0d00*X*Y_2*R51+15.0d00*X*R41 2318 R(M,0,107)=X*Y_5*Z*R71+10.0d00* X*Y_3*Z*R61 2319 & +15.0d00*X*Y*Z*R51 2320 R(M,0,108)=X*Y_4*Z_2*R71 2321 & +X*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))*R61 2322 & +X*(6.0d00*Y_2+3.0d00*Z_2)*R51+3.0d00*X*R41 2323 R(M,0,109)=X*Y_3*Z_3*R71 2324 & +X*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)*R61+9.0d00*X*Y*Z*R51 2325 R(M,0,110)=X*Y_2*Z_4*R71+X*(Z_4+6.0d00*Y_2*Z_2)*R61 2326 & +X*(3.0d00*Y_2+6.0d00*Z_2)*R51+3.0d00*X*R41 2327 R(M,0,111)=X*Y*Z_5*R71+10.0d00* X*Y*Z_3*R61 2328 & +15.0d00*X*Y*Z*R51 2329 R(M,0,112)=X*Z_6*R71+15.0d00*X*Z_4*R61+ 2330 & 45.0d00*X*Z_2*R51 2331 & +15.0d00*X*R41 2332 R(M,0,113)=Y_7*R71+21.0d00*Y_5*R61+105.0d00*Y_3*R51 2333 & +105.0d00*Y*R41 2334 R(M,0,114)=Y_6*Z*R71+15.0d00*Y_4*Z*R61+ 2335 & 45.0d00*Y_2*Z*R51 2336 & +15.0d00*Z*R41 2337 R(M,0,115)=Y_5*Z_2*R71 2338 & +(Y*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 2339 & 2.0d00*Y*Z_2))+4.0d00*Y_3*Z_2)*R61 2340 & +(Y*(6.0d00*Y_2+3.0d00*Z_2)+4.0d00*Y*(Y_2+Z_2)+ 2341 & 8.0d00*Y*Z_2)*R51 2342 & +15.0d00*Y*R41 2343 R(M,0,116)=Y_4*Z_3*R71 2344 & +(3.0d00*Y_2*Z_3+Y*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3))*R61 2345 & +(18.0d00*Y_2*Z+3.0d00*Z_3)*R51+9.0d00*Z*R41 2346 R(M,0,117)=Y_3*Z_4*R71 2347 & +(Y*(Z_4+6.0d00*Y_2*Z_2)+2.0d00*Y*Z_4)*R61 2348 & +(Y*(3.0d00*Y_2+6.0d00*Z_2)+12.0d00*Y*Z_2)*R51+9.0d00*Y*R41 2349 R(M,0,118)=Y_2*Z_5*R71+(Z_5+10.0d00* Y_2*Z_3)*R61 2350 & +(10.0d00*Z_3+15.0d00*Y_2*Z)*R51+15.0d00*Z*R41 2351 R(M,0,119)=Y*Z_6*R71+15.0d00*Y*Z_4*R61+ 2352 & 45.0d00*Y*Z_2*R51 2353 & +15.0d00*Y*R41 2354 R(M,0,120)=Z_7*R71+21.0d00*Z_5*R61+105.0d00*Z_3*R51 2355 & +105.0d00*Z*R41 2356 X_8=X_7 * X 2357 Y_8=Y_7 * Y 2358 Z_8=Z_7 * Z 2359 R81=R(m,8,1) 2360 R(M,0,121)=X_8*R81+28.0d00*X_6*R71+210.0d00*X_4*R61 2361 & +420.0d00*X_2*R51+105.0d00*R41 2362 R(M,0,122)=X_7*Y*R81+21.0d00*X_5*Y*R71 2363 & +105.0d00*X_3*Y*R61+105.0d00*X*Y*R51 2364 R(M,0,123)=X_7*Z*R81+21.0d00*X_5*Z*R71 2365 & +105.0d00*X_3*Z*R61+105.0d00*X*Z*R51 2366 R(M,0,124)=X_6*Y_2*R81 2367 & +(5.0d00*X_4*Y_2+X*(X*(3.0d00*X_2*Y_2+ 2368 & X*(X*(X_2+Y_2)+2.0d00*X*Y_2))+4.0d00*X_3*Y_2))*R71 2369 & +(X*(X*(6.0d00*X_2+3.0d00*Y_2)+4.0d00*X*(X_2+Y_2)+ 2370 & 8.0d00*X*Y_2)+ 2371 & 5.0d00*X*(X*(X_2+Y_2)+2.0d00*X*Y_2)+15.0d00*X_2*Y_2)*R61 2372 & +(45.0d00*X_2+15.0d00*Y_2)*R51+15.0d00*R41 2373 R(M,0,125)=X_6*Y*Z*R81+15.0d00*X_4*Y*Z*R71 2374 & +45.0d00*X_2*Y*Z*R61+15.0d00*Y*Z*R51 2375 R(M,0,126)=X_6*Z_2*R81 2376 & +(5.0d00*X_4*Z_2+X*(X*(3.0d00*X_2*Z_2+X*(X*(X_2+Z_2)+ 2377 & 2.0d00*X*Z_2))+4.0d00*X_3*Z_2))*R71 2378 & +(X*(X*(6.0d00*X_2+3.0d00*Z_2)+4.0d00*X*(X_2+Z_2)+ 2379 & 8.0d00*X*Z_2)+ 2380 & 5.0d00*X*(X*(X_2+Z_2)+2.0d00*X*Z_2)+15.0d00*X_2*Z_2)*R61 2381 & +(45.0d00*X_2+15.0d00*Z_2)*R51+15.0d00*R41 2382 R(M,0,127)=X_5*Y_3*R81 2383 & +(X*(3.0d00*X_2*Y_3+X*(X*(3.0d00*X_2*Y+Y_3)+2.0d00*X*Y_3))+ 2384 & 4.0d00*X_3*Y_3)*R71 2385 & +(X*(18.0d00*X_2*Y+3.0d00*Y_3)+4.0d00*X*(3.0d00*X_2*Y+Y_3)+ 2386 & 8.0d00*X*Y_3)*R61+45.0d00*X*Y*R51 2387 R(M,0,128)=X_5*Y_2*Z*R81 2388 & +(X*(3.0d00*X_2*Y_2*Z+X*(X*(X_2*Z+Y_2*Z)+2.0d00*X*Y_2*Z))+ 2389 & 4.0d00*X_3*Y_2*Z)*R71 2390 & +(X*(6.0d00*X_2*Z+3.0d00*Y_2*Z)+4.0d00*X*(X_2*Z+Y_2*Z)+ 2391 & 8.0d00*X*Y_2*Z)*R61+15.0d00*X*Z*R51 2392 R(M,0,129)=X_5*Y*Z_2*R81 2393 & +(X*(X*(X*(X_2*Y+Y*Z_2)+2.0d00*X*Y*Z_2)+3.0d00*X_2*Y*Z_2)+ 2394 & 4.0d00*X_3*Y*Z_2)*R71 2395 & +(X*(6.0d00*X_2*Y+3.0d00*Y*Z_2)+4.0d00*X*(X_2*Y+Y*Z_2)+ 2396 & 8.0d00*X*Y*Z_2)*R61+15.0d00*X*Y*R51 2397 R(M,0,130)=X_5*Z_3*R81 2398 & +(X*(3.0d00*X_2*Z_3+X*(X*(3.0d00*X_2*Z+Z_3)+2.0d00*X*Z_3))+ 2399 & 4.0d00*X_3*Z_3)*R71 2400 & +(X*(18.0d00*X_2*Z+3.0d00*Z_3)+4.0d00*X*(3.0d00*X_2*Z+Z_3)+ 2401 & 8.0d00*X*Z_3)*R61+45.0d00*X*Z*R51 2402 R(M,0,131)=X_4*Y_4*R81 2403 & +(3.0d00*X_2*Y_4+X*(X*(Y_4+6.0d00*X_2*Y_2)+ 2404 & 2.0d00*X*Y_4))*R71 2405 & +(X*(X*(3.0d00*X_2+6.0d00*Y_2)+12.0d00*X*Y_2)+ 2406 & 18.0d00*X_2*Y_2+3.0d00*Y_4)* 2407 & R61+(18.0d00*X_2+18.0d00*Y_2)*R51+9.0d00*R41 2408 R(M,0,132)=X_4*Y_3*Z*R81 2409 & +(3.0d00*X_2*Y_3*Z+X*(X*(3.0d00*X_2*Y*Z+Y_3*Z)+ 2410 & 2.0d00*X*Y_3*Z))* 2411 & R71+(18.0d00*X_2*Y*Z+3.0d00*Y_3*Z)*R61+9.0d00*Y*Z*R51 2412 R(M,0,133)=X_4*Y_2*Z_2*R81 2413 & +(3.0d00*X_2*Y_2*Z_2+X*(X*(Y_2*Z_2+X_2*(Y_2+Z_2))+ 2414 & 2.0d00*X*Y_2*Z_2))*R71 2415 & +(X*(X*(X_2+Y_2+Z_2)+2.0d00*X*(Y_2+Z_2))+ 2416 & 3.0d00*X_2*(Y_2+Z_2)+3.0d00*Y_2*Z_2)*R61 2417 & +(6.0d00*X_2+3.0d00*Y_2+3.0d00*Z_2)*R51+3.0d00*R41 2418 R(M,0,134)=X_4*Y*Z_3*R81 2419 & +(X*(X*(3.0d00*X_2*Y*Z+Y*Z_3)+2.0d00*X*Y*Z_3)+ 2420 & 3.0d00*X_2*Y*Z_3)* 2421 & R71+(18.0d00*X_2*Y*Z+3.0d00*Y*Z_3)*R61+9.0d00*Y*Z*R51 2422 R(M,0,135)=X_4*Z_4*R81 2423 & +(3.0d00*X_2*Z_4+X*(X*(Z_4+6.0d00*X_2*Z_2)+ 2424 & 2.0d00*X*Z_4))*R71 2425 & +(X*(X*(3.0d00*X_2+6.0d00*Z_2)+12.0d00*X*Z_2)+ 2426 & 18.0d00*X_2*Z_2+3.0d00*Z_4)* 2427 & R61+(18.0d00*X_2+18.0d00*Z_2)*R51+9.0d00*R41 2428 R(M,0,136)=X_3*Y_5*R81 2429 & +(X*(Y_5+10.0d00* X_2*Y_3)+2.0d00*X*Y_5)*R71 2430 & +(X*(10.0d00*Y_3+15.0d00*X_2*Y)+20.0d00*X*Y_3)*R61+ 2431 & 45.0d00*X*Y*R51 2432 R(M,0,137)=X_3*Y_4*Z*R81 2433 & +(X*(Y_4*Z+6.0d00*X_2*Y_2*Z)+2.0d00*X*Y_4*Z)*R71 2434 & +(X*(3.0d00*X_2*Z+6.0d00*Y_2*Z)+12.0d00*X*Y_2*Z)*R61+ 2435 & 9.0d00*X*Z*R51 2436 R(M,0,138)=X_3*Y_3*Z_2*R81 2437 & +(X*(X_2*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+Y_3*Z_2)+ 2438 & 2.0d00*X*Y_3*Z_2)*R71 2439 & +(X*(3.0d00*X_2*Y+Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+ 2440 & 2.0d00*X*(Y*(Y_2+Z_2)+ 2441 & 2.0d00*Y*Z_2))*R61+9.0d00*X*Y*R51 2442 R(M,0,139)=X_3*Y_2*Z_3*R81 2443 & +(X*(Y_2*Z_3+X_2*(3.0d00*Y_2*Z+Z_3))+2.0d00*X*Y_2*Z_3)* 2444 & R71+(X*(3.0d00*X_2*Z+3.0d00*Y_2*Z+Z_3)+ 2445 & 2.0d00*X*(3.0d00*Y_2*Z+Z_3))* 2446 & R61+9.0d00*X*Z*R51 2447 R(M,0,140)=X_3*Y*Z_4*R81 2448 & +(X*(6.0d00*X_2*Y*Z_2+Y*Z_4)+2.0d00*X*Y*Z_4)*R71 2449 & +(X*(3.0d00*X_2*Y+6.0d00*Y*Z_2)+12.0d00*X*Y*Z_2)*R61+ 2450 & 9.0d00*X*Y*R51 2451 R(M,0,141)=X_3*Z_5*R81 2452 & +(X*(Z_5+10.0d00* X_2*Z_3)+2.0d00*X*Z_5)*R71 2453 & +(X*(10.0d00*Z_3+15.0d00*X_2*Z)+20.0d00*X*Z_3)*R61+ 2454 & 45.0d00*X*Z*R51 2455 R(M,0,142)=X_2*Y_6*R81+(Y_6+15.0d00*X_2*Y_4)*R71 2456 & +(45.0d00*X_2*Y_2+15.0d00*Y_4)*R61+ 2457 & (15.0d00*X_2+45.0d00*Y_2)*R51+15.0d00*R41 2458 R(M,0,143)=X_2*Y_5*Z*R81+(Y_5*Z+10.0d00* X_2*Y_3*Z)* 2459 & R71+(10.0d00*Y_3*Z+15.0d00*X_2*Y*Z)*R61+15.0d00*Y*Z*R51 2460 R(M,0,144)=X_2*Y_4*Z_2*R81 2461 & +(Y_4*Z_2+X_2*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 2462 & 2.0d00*Y*Z_2)))*R71+ 2463 & (X_2*(6.0d00*Y_2+3.0d00*Z_2)+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+ 2464 & 3.0d00*Y_2*Z_2)*R61+(3.0d00*X_2+6.0d00*Y_2+3.0d00*Z_2)*R51 2465 & +3.0d00*R41 2466 R(M,0,145)=X_2*Y_3*Z_3*R81 2467 & +(X_2*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)+Y_3*Z_3)*R71 2468 & +(9.0d00*X_2*Y*Z+Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)*R61 2469 & +9.0d00*Y*Z*R51 2470 R(M,0,146)=X_2*Y_2*Z_4*R81 2471 & +(Y_2*Z_4+X_2*(Z_4+6.0d00*Y_2*Z_2))*R71 2472 & +(X_2*(3.0d00*Y_2+6.0d00*Z_2)+6.0d00*Y_2*Z_2+Z_4)*R61 2473 & +(3.0d00*X_2+3.0d00*Y_2+6.0d00*Z_2)*R51+3.0d00*R41 2474 R(M,0,147)=X_2*Y*Z_5*R81+(10.0d00*X_2*Y*Z_3+Y*Z_5)* 2475 & R71+(15.0d00*X_2*Y*Z+10.0d00* Y*Z_3)*R61+15.0d00*Y*Z*R51 2476 R(M,0,148)=X_2*Z_6*R81+(Z_6+15.0d00*X_2*Z_4)*R71 2477 & +(45.0d00*X_2*Z_2+15.0d00*Z_4)*R61+(15.0d00*X_2+ 2478 & 45.0d00*Z_2)* 2479 & R51+15.0d00*R41 2480 R(M,0,149)=X*Y_7*R81+21.0d00*X*Y_5*R71 2481 & +105.0d00*X*Y_3*R61+105.0d00*X*Y*R51 2482 R(M,0,150)=X*Y_6*Z*R81+15.0d00*X*Y_4*Z*R71 2483 & +45.0d00*X*Y_2*Z*R61+15.0d00*X*Z*R51 2484 R(M,0,151)=X*Y_5*Z_2*R81 2485 & +X*(Y*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))+ 2486 & 4.0d00*Y_3*Z_2)*R71 2487 & +X*(Y*(6.0d00*Y_2+3.0d00*Z_2)+4.0d00*Y*(Y_2+Z_2)+ 2488 & 8.0d00*Y*Z_2)*R61 2489 & +15.0d00*X*Y*R51 2490 R(M,0,152)=X*Y_4*Z_3*R81 2491 & +X*(3.0d00*Y_2*Z_3+Y*(Y*(3.0d00*Y_2*Z+Z_3)+ 2492 & 2.0d00*Y*Z_3))*R71 2493 & +X*(18.0d00*Y_2*Z+3.0d00*Z_3)*R61+9.0d00*X*Z*R51 2494 R(M,0,153)=X*Y_3*Z_4*R81 2495 & +X*(Y*(Z_4+6.0d00*Y_2*Z_2)+2.0d00*Y*Z_4)*R71 2496 & +X*(Y*(3.0d00*Y_2+6.0d00*Z_2)+12.0d00*Y*Z_2)*R61+ 2497 & 9.0d00*X*Y*R51 2498 R(M,0,154)=X*Y_2*Z_5*R81+X*(Z_5+10.0d00* Y_2*Z_3)*R71 2499 & +X*(10.0d00*Z_3+15.0d00*Y_2*Z)*R61+15.0d00*X*Z*R51 2500 R(M,0,155)=X*Y*Z_6*R81+15.0d00*X*Y*Z_4*R71 2501 & +45.0d00*X*Y*Z_2*R61+15.0d00*X*Y*R51 2502 R(M,0,156)=X*Z_7*R81+21.0d00*X*Z_5*R71 2503 & +105.0d00*X*Z_3*R61+105.0d00*X*Z*R51 2504 R(M,0,157)=Y_8*R81+28.0d00*Y_6*R71+210.0d00*Y_4*R61 2505 & +420.0d00*Y_2*R51+105.0d00*R41 2506 R(M,0,158)=Y_7*Z*R81+21.0d00*Y_5*Z*R71 2507 & +105.0d00*Y_3*Z*R61+105.0d00*Y*Z*R51 2508 R(M,0,159)=Y_6*Z_2*R81 2509 & +(5.0d00*Y_4*Z_2+Y*(Y*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 2510 & 2.0d00*Y*Z_2))+4.0d00*Y_3*Z_2))*R71 2511 & +(Y*(Y*(6.0d00*Y_2+3.0d00*Z_2)+4.0d00*Y*(Y_2+Z_2)+ 2512 & 8.0d00*Y*Z_2)+ 2513 & 5.0d00*Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+15.0d00*Y_2*Z_2)*R61 2514 & +(45.0d00*Y_2+15.0d00*Z_2)*R51+15.0d00*R41 2515 R(M,0,160)=Y_5*Z_3*R81 2516 & +(Y*(3.0d00*Y_2*Z_3+Y*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3))+ 2517 & 4.0d00*Y_3*Z_3)*R71 2518 & +(Y*(18.0d00*Y_2*Z+3.0d00*Z_3)+4.0d00*Y*(3.0d00*Y_2*Z+Z_3)+ 2519 & 8.0d00*Y*Z_3)* 2520 & R61+45.0d00*Y*Z*R51 2521 R(M,0,161)=Y_4*Z_4*R81 2522 & +(3.0d00*Y_2*Z_4+Y*(Y*(Z_4+6.0d00*Y_2*Z_2)+ 2523 & 2.0d00*Y*Z_4))*R71 2524 & +(Y*(Y*(3.0d00*Y_2+6.0d00*Z_2)+12.0d00*Y*Z_2)+ 2525 & 18.0d00*Y_2*Z_2+3.0d00*Z_4)* 2526 & R61+(18.0d00*Y_2+18.0d00*Z_2)*R51+9.0d00*R41 2527 R(M,0,162)=Y_3*Z_5*R81 2528 & +(Y*(Z_5+10.0d00* Y_2*Z_3)+2.0d00*Y*Z_5)*R71 2529 & +(Y*(10.0d00*Z_3+15.0d00*Y_2*Z)+20.0d00*Y*Z_3)*R61+ 2530 & 45.0d00*Y*Z*R51 2531 R(M,0,163)=Y_2*Z_6*R81+(Z_6+15.0d00*Y_2*Z_4)*R71 2532 & +(45.0d00*Y_2*Z_2+15.0d00*Z_4)*R61+ 2533 & (15.0d00*Y_2+45.0d00*Z_2)*R51 2534 & +15.0d00*R41 2535 R(M,0,164)=Y*Z_7*R81+21.0d00*Y*Z_5*R71 2536 & +105.0d00*Y*Z_3*R61+105.0d00*Y*Z*R51 2537 R(M,0,165)=Z_8*R81+28.0d00*Z_6*R71+210.0d00*Z_4*R61 2538 & +420.0d00*Z_2*R51+105.0d00*R41 2539 R91=R(m,9,1) 2540 R(M,0,166)=X_5*X_4*R91+36.0d00*X_7*R81+378.0d00*X_5*R71 2541 & +1260.0d00*X_3*R61+945.0d00*X*R51 2542 R(M,0,167)=X_8*Y*R91+28.0d00*X_6*Y*R81 2543 & +210.0d00*X_4*Y*R71+420.0d00*X_2*Y*R61+105.0d00*Y*R51 2544 R(M,0,168)=X_8*Z*R91+28.0d00*X_6*Z*R81 2545 & +210.0d00*X_4*Z*R71+420.0d00*X_2*Z*R61+105.0d00*Z*R51 2546 R(M,0,169)=X_7*Y_2*R91 2547 & +(6.0d00*X_5*Y_2+X*(5.0d00*X_4*Y_2+X*(X*(3.0d00* 2548 & X_2*Y_2+X*(X*(X_2+ 2549 & Y_2)+2.0d00*X*Y_2))+4.0d00*X_3*Y_2)))*R81 2550 & +(6.0d00*X*(3.0d00*X_2*Y_2+X*(X*(X_2+Y_2)+2.0d00*X*Y_2))+ 2551 & X*(X*(X*(6.0d00*X_2+3.0d00*Y_2)+4.0d00*X*(X_2+Y_2)+ 2552 & 8.0d00*X*Y_2)+ 2553 & 5.0d00*X*(X*(X_2+Y_2)+2.0d00*X*Y_2)+15.0d00*X_2*Y_2)+ 2554 & 24.0d00*X_3*Y_2)* 2555 & R71+(X*(45.0d00*X_2+15.0d00*Y_2)+24.0d00*X*(X_2+Y_2)+ 2556 & 48.0d00*X*Y_2+ 2557 & 6.0d00*X*(6.0d00*X_2+3.0d00*Y_2))*R61+105.0d00*X*R51 2558 R(M,0,170)=X_7*Y*Z*R91+21.0d00*X_5*Y*Z*R81 2559 & +105.0d00*X_3*Y*Z*R71+105.0d00*X*Y*Z*R61 2560 R(M,0,171)=X_7*Z_2*R91 2561 & +(6.0d00*X_5*Z_2+X*(5.0d00*X_4*Z_2+X*(X*(3.0d00*X_2*Z_2+ 2562 & X*(X*(X_2+ 2563 & Z_2)+2.0d00*X*Z_2))+4.0d00*X_3*Z_2)))*R81+(6.0d00*X* 2564 & (3.0d00*X_2*Z_2 2565 & +X*(X*(X_2+Z_2)+2.0d00*X*Z_2))+X*(X*(X*(6.0d00*X_2+ 2566 & 3.0d00*Z_2)+4.0d00*X* 2567 & (X_2+Z_2)+8.0d00*X*Z_2)+5.0d00*X*(X*(X_2+Z_2)+ 2568 & 2.0d00*X*Z_2)+ 2569 & 15.0d00*X_2*Z_2)+24.0d00*X_3*Z_2)*R71+(X*(45.0d00*X_2+ 2570 & 15.0d00*Z_2)+ 2571 & 24.0d00*X*(X_2+Z_2)+48.0d00*X*Z_2+6.0d00*X*(6.0d00*X_2+ 2572 & 3.0d00*Z_2))*R61 2573 & +105.0d00*X*R51 2574 R(M,0,172)=X_6*Y_3*R91 2575 & +(5.0d00*X_4*Y_3+X*(X*(3.0d00*X_2*Y_3+ 2576 & X*(X*(3.0d00*X_2*Y+Y_3)+ 2577 & 2.0d00*X*Y_3))+4.0d00*X_3*Y_3))*R81+(X*(X*(18.0d00*X_2*Y+ 2578 & 3.0d00*Y_3)+4.0d00*X*(3.0d00*X_2*Y+Y_3)+8.0d00*X*Y_3)+ 2579 & 5.0d00*X*(X*(3.0d00*X_2*Y+ 2580 & Y_3)+2.0d00*X*Y_3)+15.0d00*X_2*Y_3)*R71+(135.0d00*X_2*Y+ 2581 & 15.0d00*Y_3)*R61+45.0d00*Y*R51 2582 R(M,0,173)=X_6*Y_2*Z*R91 2583 & +(5.0d00*X_4*Y_2*Z+X*(X*(3.0d00*X_2*Y_2*Z+ 2584 & X*(X*(X_2*Z+Y_2*Z)+ 2585 & 2.0d00*X*Y_2*Z))+4.0d00*X_3*Y_2*Z))*R81+(X*(X* 2586 & (6.0d00*X_2*Z+ 2587 & 3.0d00*Y_2*Z)+4.0d00*X*(X_2*Z+Y_2*Z)+ 2588 & 8.0d00*X*Y_2*Z)+5.0d00*X*(X*(X_2*Z+ 2589 & Y_2*Z)+2.0d00*X*Y_2*Z)+15.0d00*X_2*Y_2*Z)*R71 2590 & +(45.0d00*X_2*Z+15.0d00*Y_2*Z)*R61+15.0d00*Z*R51 2591 R(M,0,174)=X_6*Y*Z_2*R91 2592 & +(5.0d00*X_4*Y*Z_2+X*(X*(X*(X*(X_2*Y+Y*Z_2)+ 2593 & 2.0d00*X*Y*Z_2)+ 2594 & 3.0d00*X_2*Y*Z_2)+4.0d00*X_3*Y*Z_2))*R81+ 2595 & (5.0d00*X*(X*(X_2*Y+ 2596 & Y*Z_2)+2.0d00*X*Y*Z_2)+X*(X*(6.0d00*X_2*Y+ 2597 & 3.0d00*Y*Z_2)+4.0d00*X*(X_2*Y+ 2598 & Y*Z_2)+8.0d00*X*Y*Z_2)+15.0d00*X_2*Y*Z_2)*R71+ 2599 & (45.0d00*X_2*Y+ 2600 & 15.0d00*Y*Z_2)*R61+15.0d00*Y*R51 2601 R(M,0,175)=X_6*Z_3*R91 2602 & +(5.0d00*X_4*Z_3+X*(X*(3.0d00*X_2*Z_3+X*(X* 2603 & (3.0d00*X_2*Z+Z_3)+ 2604 & 2.0d00*X*Z_3))+4.0d00*X_3*Z_3))*R81+(X*(X*(18.0d00*X_2*Z+ 2605 & 3.0d00*Z_3)+4.0d00*X*(3.0d00*X_2*Z+Z_3)+8.0d00*X*Z_3)+ 2606 & 5.0d00*X*(X*(3.0d00*X_2*Z+ 2607 & Z_3)+2.0d00*X*Z_3)+15.0d00*X_2*Z_3)*R71+(135.0d00*X_2*Z+ 2608 & 15.0d00*Z_3)*R61+45.0d00*Z*R51 2609 R(M,0,176)=X_5*Y_4*R91 2610 & +(X*(3.0d00*X_2*Y_4+X*(X*(Y_4+6.0d00*X_2*Y_2)+ 2611 & 2.0d00*X*Y_4))+ 2612 & 4.0d00*X_3*Y_4)*R81+(X*(X*(X*(3.0d00*X_2+ 2613 & 6.0d00*Y_2)+12.0d00*X*Y_2)+ 2614 & 18.0d00*X_2*Y_2+3.0d00*Y_4)+4.0d00*X*(Y_4+ 2615 & 6.0d00*X_2*Y_2)+8.0d00*X*Y_4)* 2616 & R71+(X*(18.0d00*X_2+18.0d00*Y_2)+ 2617 & 4.0d00*X*(3.0d00*X_2+6.0d00*Y_2)+ 2618 & 48.0d00*X*Y_2)*R61+45.0d00*X*R51 2619 R(M,0,177)=X_5*Y_3*Z*R91 2620 & +(X*(3.0d00*X_2*Y_3*Z+X*(X*(3.0d00*X_2*Y*Z+Y_3*Z)+ 2621 & 2.0d00*X*Y_3*Z))+ 2622 & 4.0d00*X_3*Y_3*Z)*R81+(X*(18.0d00*X_2*Y*Z+3.0d00*Y_3*Z)+ 2623 & 4.0d00*X*(3.0d00*X_2*Y*Z+Y_3*Z)+8.0d00*X*Y_3*Z)*R71+ 2624 & 45.0d00*X*Y*Z*R61 2625 R(M,0,178)=X_5*Y_2*Z_2*R91 2626 & +(X*(3.0d00*X_2*Y_2*Z_2+X*(X*(Y_2*Z_2+X_2*(Y_2+Z_2))+ 2627 & 2.0d00*X*Y_2*Z_2))+4.0d00*X_3*Y_2*Z_2)*R81+(X*(X*(X*(X_2+ 2628 & Y_2+Z_2)+2.0d00*X*(Y_2+Z_2))+3.0d00*X_2*(Y_2+Z_2)+ 2629 & 3.0d00*Y_2* 2630 & Z_2)+4.0d00*X*(Y_2*Z_2+X_2*(Y_2+Z_2))+8.0d00*X*Y_2*Z_2)* 2631 & R71+(X*(6.0d00*X_2+3.0d00*Y_2+3.0d00*Z_2)+ 2632 & 4.0d00*X*(X_2+Y_2+Z_2)+ 2633 & 8.0d00*X*(Y_2+Z_2))*R61+15.0d00*X*R51 2634 R(M,0,179)=X_5*Y*Z_3*R91 2635 & +(X*(X*(X*(3.0d00*X_2*Y*Z+Y*Z_3)+2.0d00*X*Y*Z_3)+ 2636 & 3.0d00*X_2*Y*Z_3)+ 2637 & 4.0d00*X_3*Y*Z_3)*R81+(X*(18.0d00*X_2*Y*Z+ 2638 & 3.0d00*Y*Z_3)+4.0d00*X* 2639 & (3.0d00*X_2*Y*Z+Y*Z_3)+8.0d00*X*Y*Z_3)*R71+ 2640 & 45.0d00*X*Y*Z*R61 2641 R(M,0,180)=X_5*Z_4*R91 2642 & +(X*(3.0d00*X_2*Z_4+X*(X*(Z_4+6.0d00*X_2*Z_2)+ 2643 & 2.0d00*X*Z_4))+ 2644 & 4.0d00*X_3*Z_4)*R81+(X*(X*(X*(3.0d00*X_2+6.0d00*Z_2)+ 2645 & 12.0d00*X*Z_2)+ 2646 & 18.0d00*X_2*Z_2+3.0d00*Z_4)+4.0d00*X*(Z_4+6.0d00*X_2*Z_2)+ 2647 & 8.0d00*X*Z_4)* 2648 & R71+(X*(18.0d00*X_2+18.0d00*Z_2)+4.0d00*X*(3.0d00*X_2+ 2649 & 6.0d00*Z_2)+48.0d00*X*Z_2)*R61+45.0d00*X*R51 2650 R(M,0,181)=X_4*Y_5*R91 2651 & +(3.0d00*X_2*Y_5+X*(X*(Y_5+10.0d00* X_2*Y_3)+ 2652 & 2.0d00*X*Y_5))*R81 2653 & +(X*(X*(10.0d00*Y_3+15.0d00*X_2*Y)+20.0d00*X*Y_3)+ 2654 & 30.0d00*X_2*Y_3+ 2655 & 3.0d00*Y_5)*R71+(90.0d00*X_2*Y+30.0d00*Y_3)*R61+ 2656 & 45.0d00*Y*R51 2657 R(M,0,182)=X_4*Y_4*Z*R91+(3.0d00*X_2*Y_4*Z+X*(X*(Y_4*Z+ 2658 & 6.0d00*X_2*Y_2*Z)+2.0d00*X*Y_4*Z))*R81+(X*(X*( 2659 & 3.0d00*X_2*Z+ 2660 & 6.0d00*Y_2*Z)+12.0d00*X*Y_2*Z)+18.0d00*X_2*Y_2*Z+ 2661 & 3.0d00*Y_4*Z)*R71 2662 & +(18.0d00*X_2*Z+18.0d00*Y_2*Z)*R61+9.0d00*Z*R51 2663 R(M,0,183)=X_4*Y_3*Z_2*R91 2664 & +(3.0d00*X_2*Y_3*Z_2+X*(X*(X_2*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+ 2665 & Y_3*Z_2)+2.0d00*X*Y_3*Z_2))*R81+(3.0d00*X_2*(Y*(Y_2+ 2666 & Z_2)+2.0d00*Y*Z_2)+X*(X*(3.0d00*X_2*Y+Y*(Y_2+Z_2)+ 2667 & 2.0d00*Y*Z_2)+ 2668 & 2.0d00*X*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))+3.0d00*Y_3*Z_2)*R71+ 2669 & (18.0d00*X_2*Y+3.0d00*Y*(Y_2+Z_2)+6.0d00*Y*Z_2)*R61+ 2670 & 9.0d00*Y*R51 2671 R(M,0,184)=X_4*Y_2*Z_3*R91 2672 & +(3.0d00*X_2*Y_2*Z_3+X*(X*(Y_2*Z_3+X_2*(3.0d00*Y_2*Z+Z_3))+ 2673 & 2.0d00*X*Y_2*Z_3))*R81+(X*(X*(3.0d00*X_2*Z+ 2674 & 3.0d00*Y_2*Z+Z_3)+ 2675 & 2.0d00*X*(3.0d00*Y_2*Z+Z_3))+ 2676 & 3.0d00*X_2*(3.0d00*Y_2*Z+Z_3)+3.0d00*Y_2*Z_3)* 2677 & R71+(9.0d00*Y_2*Z+3.0d00*Z_3+18.0d00*X_2*Z)*R61+ 2678 & 9.0d00*Z*R51 2679 R(M,0,185)=X_4*Y*Z_4*R91 2680 & +(X*(X*(6.0d00*X_2*Y*Z_2+Y*Z_4)+2.0d00*X*Y*Z_4)+ 2681 & 3.0d00*X_2*Y*Z_4)* 2682 & R81+(X*(X*(3.0d00*X_2*Y+6.0d00*Y*Z_2)+12.0d00*X*Y*Z_2)+ 2683 & 18.0d00*X_2*Y*Z_2+3.0d00*Y*Z_4)*R71 2684 & +(18.0d00*X_2*Y+18.0d00*Y*Z_2)*R61+9.0d00*Y*R51 2685 R(M,0,186)=X_4*Z_5*R91 2686 & +(3.0d00*X_2*Z_5+X*(X*(Z_5+10.0d00* X_2*Z_3)+ 2687 & 2.0d00*X*Z_5))*R81 2688 & +(X*(X*(10.0d00*Z_3+15.0d00*X_2*Z)+20.0d00*X*Z_3)+ 2689 & 30.0d00*X_2*Z_3+ 2690 & 3.0d00*Z_5)*R71+(90.0d00*X_2*Z+30.0d00*Z_3)*R61+ 2691 & 45.0d00*Z*R51 2692 R(M,0,187)=X_3*Y_6*R91 2693 & +(X*(Y_6+15.0d00*X_2*Y_4)+2.0d00*X*Y_6)*R81 2694 & +(X*(45.0d00*X_2*Y_2+15.0d00*Y_4)+30.0d00*X*Y_4)*R71 2695 & +(X*(15.0d00*X_2+45.0d00*Y_2)+90.0d00*X*Y_2)*R61+ 2696 & 45.0d00*X*R51 2697 R(M,0,188)=X_3*Y_5*Z*R91 2698 & +(X*(Y_5*Z+10.0d00* X_2*Y_3*Z)+2.0d00*X*Y_5*Z)*R81 2699 & +(X*(10.0d00*Y_3*Z+15.0d00*X_2*Y*Z)+20.0d00*X*Y_3*Z)*R71 2700 & +45.0d00*X*Y*Z*R61 2701 R(M,0,189)=X_3*Y_4*Z_2*R91 2702 & +(X*(Y_4*Z_2+X_2*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 2703 & 2.0d00*Y*Z_2)))+2.0d00*X*Y_4*Z_2)*R81 2704 & +(X*(X_2*(6.0d00*Y_2+ 2705 & 3.0d00*Z_2)+Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+ 2706 & 3.0d00*Y_2*Z_2)+2.0d00*X*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 2707 & 2.0d00*Y*Z_2)))*R71+(X*(3.0d00*X_2+6.0d00*Y_2+3.0d00*Z_2)+ 2708 & 2.0d00*X*(6.0d00*Y_2+3.0d00*Z_2))*R61+9.0d00*X*R51 2709 R(M,0,190)=X_3*Y_3*Z_3*R91 2710 & +(X*(X_2*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)+Y_3*Z_3)+ 2711 & 2.0d00*X*Y_3*Z_3)*R81 2712 & +(X*(9.0d00*X_2*Y*Z+Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)+ 2713 & 2.0d00*X*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3))*R71+ 2714 & 27.0d00*X*Y*Z*R61 2715 R(M,0,191)=X_3*Y_2*Z_4*R91 2716 & +(X*(Y_2*Z_4+X_2*(Z_4+6.0d00*Y_2*Z_2))+2.0d00*X*Y_2*Z_4)* 2717 & R81+(X*(X_2*(3.0d00*Y_2+6.0d00*Z_2)+6.0d00*Y_2*Z_2+Z_4)+ 2718 & 2.0d00*X*(Z_4+6.0d00*Y_2*Z_2))*R71 2719 & +(X*(3.0d00*X_2+3.0d00*Y_2+6.0d00*Z_2)+ 2720 & 2.0d00*X*(3.0d00*Y_2+6.0d00*Z_2))*R61 2721 & +9.0d00*X*R51 2722 R(M,0,192)=X_3*Y*Z_5*R91 2723 & +(X*(10.0d00*X_2*Y*Z_3+Y*Z_5)+2.0d00*X*Y*Z_5)*R81 2724 & +(X*(15.0d00*X_2*Y*Z+10.0d00* Y*Z_3)+20.0d00*X*Y*Z_3)*R71 2725 & +45.0d00*X*Y*Z*R61 2726 R(M,0,193)=X_3*Z_6*R91 2727 & +(X*(Z_6+15.0d00*X_2*Z_4)+2.0d00*X*Z_6)*R81 2728 & +(X*(45.0d00*X_2*Z_2+15.0d00*Z_4)+30.0d00*X*Z_4)*R71 2729 & +(X*(15.0d00*X_2+45.0d00*Z_2)+90.0d00*X*Z_2)*R61+ 2730 & 45.0d00*X*R51 2731 R(M,0,194)=X_2*Y_7*R91+(Y_7+21.0d00*X_2*Y_5)*R81 2732 & +(105.0d00*X_2*Y_3+21.0d00*Y_5)*R71 2733 & +(105.0d00*X_2*Y+105.0d00*Y_3)*R61+105.0d00*Y*R51 2734 R(M,0,195)=X_2*Y_6*Z*R91+(Y_6*Z+15.0d00*X_2*Y_4*Z)* 2735 & R81+(45.0d00*X_2*Y_2*Z+15.0d00*Y_4*Z)*R71 2736 & +(15.0d00*X_2*Z+45.0d00*Y_2*Z)*R61+15.0d00*Z*R51 2737 R(M,0,196)=X_2*Y_5*Z_2*R91 2738 & +(Y_5*Z_2+X_2*(Y*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 2739 & 2.0d00*Y*Z_2))+4.0d00*Y_3*Z_2))*R81+(Y*(3.0d00*Y_2*Z_2+ 2740 & Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))+X_2*(Y*(6.0d00*Y_2+ 2741 & 3.0d00*Z_2)+ 2742 & 4.0d00*Y*(Y_2+Z_2)+8.0d00*Y*Z_2)+4.0d00*Y_3*Z_2)*R71 2743 & +(15.0d00*X_2*Y+4.0d00*Y*(Y_2+Z_2)+8.0d00*Y*Z_2+ 2744 & Y*(6.0d00*Y_2+3.0d00*Z_2))* 2745 & R61+15.0d00*Y*R51 2746 R(M,0,197)=X_2*Y_4*Z_3*R91 2747 & +(Y_4*Z_3+X_2*(3.0d00*Y_2*Z_3+Y*(Y*(3.0d00*Y_2*Z+Z_3)+ 2748 & 2.0d00*Y*Z_3)))*R81+(X_2*(18.0d00*Y_2*Z+3.0d00*Z_3)+ 2749 & Y*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)+3.0d00*Y_2*Z_3)*R71 2750 & +(18.0d00*Y_2*Z+3.0d00*Z_3+9.0d00*X_2*Z)*R61+9.0d00*Z*R51 2751 R(M,0,198)=X_2*Y_3*Z_4*R91 2752 & +(X_2*(Y*(Z_4+6.0d00*Y_2*Z_2)+2.0d00*Y*Z_4)+Y_3*Z_4)*R81 2753 & +(X_2*(Y*(3.0d00*Y_2+6.0d00*Z_2)+ 2754 & 12.0d00*Y*Z_2)+Y*(Z_4+6.0d00*Y_2*Z_2)+ 2755 & 2.0d00*Y*Z_4)*R71 2756 & +(9.0d00*X_2*Y+Y*(3.0d00*Y_2+6.0d00*Z_2)+12.0d00*Y*Z_2)*R61 2757 & +9.0d00*Y*R51 2758 R(M,0,199)=X_2*Y_2*Z_5*R91 2759 & +(Y_2*Z_5+X_2*(Z_5+10.0d00* Y_2*Z_3))*R81 2760 & +(X_2*(10.0d00*Z_3+15.0d00*Y_2*Z)+10.0d00* Y_2*Z_3+Z_5)*R71 2761 & +(15.0d00*Y_2*Z+10.0d00* Z_3+15.0d00*X_2*Z)*R61+ 2762 & 15.0d00*Z*R51 2763 R(M,0,200)=X_2*Y*Z_6*R91+ 2764 & (15.0d00*X_2*Y*Z_4+Y*Z_6)*R81 2765 & +(45.0d00*X_2*Y*Z_2+15.0d00*Y*Z_4)*R71 2766 & +(15.0d00*X_2*Y+45.0d00*Y*Z_2)*R61+15.0d00*Y*R51 2767 R(M,0,201)=X_2*Z_7*R91+(Z_7+21.0d00*X_2*Z_5)*R81 2768 & +(105.0d00*X_2*Z_3+21.0d00*Z_5)*R71 2769 & +(105.0d00*X_2*Z+105.0d00*Z_3)*R61+105.0d00*Z*R51 2770 R(M,0,202)=X*Y_8*R91+28.0d00*X*Y_6*R81 2771 & +210.0d00*X*Y_4*R71+420.0d00*X*Y_2*R61+105.0d00*X*R51 2772 R(M,0,203)=X*Y_7*Z*R91+21.0d00*X*Y_5*Z*R81 2773 & +105.0d00*X*Y_3*Z*R71+105.0d00*X*Y*Z*R61 2774 R(M,0,204)=X*Y_6*Z_2*R91 2775 & +X*(5.0d00*Y_4*Z_2+Y*(Y*(3.0d00*Y_2*Z_2+Y*(Y*(Y_2+Z_2)+ 2776 & 2.0d00*Y*Z_2))+4.0d00*Y_3*Z_2))*R81 2777 & +X*(Y*(Y*(6.0d00*Y_2+3.0d00*Z_2)+4.0d00*Y*(Y_2+Z_2)+ 2778 & 8.0d00*Y*Z_2)+ 2779 & 5.0d00*Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2)+15.0d00*Y_2*Z_2)*R71 2780 & +X*(45.0d00*Y_2+15.0d00*Z_2)*R61+15.0d00*X*R51 2781 R(M,0,205)=X*Y_5*Z_3*R91 2782 & +X*(Y*(3.0d00*Y_2*Z_3+Y*(Y*(3.0d00*Y_2*Z+Z_3)+ 2783 & 2.0d00*Y*Z_3))+4.0d00*Y_3*Z_3)*R81 2784 & +X*(Y*(18.0d00*Y_2*Z+3.0d00*Z_3)+ 2785 & 4.0d00*Y*(3.0d00*Y_2*Z+Z_3)+8.0d00*Y*Z_3)*R71+ 2786 & 45.0d00*X*Y*Z*R61 2787 R(M,0,206)=X*Y_4*Z_4*R91 2788 & +X*(3.0d00*Y_2*Z_4+Y*(Y*(Z_4+6.0d00*Y_2*Z_2)+ 2789 & 2.0d00*Y*Z_4))* 2790 & R81+X*(Y*(Y*(3.0d00*Y_2+6.0d00*Z_2)+12.0d00*Y*Z_2)+ 2791 & 18.0d00*Y_2*Z_2+ 2792 & 3.0d00*Z_4)*R71+X*(18.0d00*Y_2+18.0d00*Z_2)*R61+ 2793 & 9.0d00*X*R51 2794 R(M,0,207)=X*Y_3*Z_5*R91 2795 & +X*(Y*(Z_5+10.0d00* Y_2*Z_3)+2.0d00*Y*Z_5)*R81 2796 & +X*(Y*(10.0d00*Z_3+15.0d00*Y_2*Z)+20.0d00*Y*Z_3)*R71 2797 & +45.0d00*X*Y*Z*R61 2798 R(M,0,208)=X*Y_2*Z_6*R91+X*(Z_6+15.0d00*Y_2*Z_4)*R81 2799 & +X*(45.0d00*Y_2*Z_2+15.0d00*Z_4)*R71 2800 & +X*(15.0d00*Y_2+45.0d00*Z_2)*R61+15.0d00*X*R51 2801 R(M,0,209)=X*Y*Z_7*R91+21.0d00*X*Y*Z_5*R81 2802 & +105.0d00*X*Y*Z_3*R71+105.0d00*X*Y*Z*R61 2803 R(M,0,210)=X*Z_8*R91+28.0d00*X*Z_6*R81 2804 & +210.0d00*X*Z_4*R71+420.0d00*X*Z_2*R61+105.0d00*X*R51 2805 R(M,0,211)=Y_5*Y_4*R91+36.0d00*Y_7*R81+378.0d00*Y_5*R71 2806 & +1260.0d00*Y_3*R61+945.0d00*Y*R51 2807 R(M,0,212)=Y_8*Z*R91+28.0d00*Y_6*Z*R81 2808 & +210.0d00*Y_4*Z*R71+420.0d00*Y_2*Z*R61+105.0d00*Z*R51 2809 R(M,0,213)=Y_7*Z_2*R91 2810 & +(6.0d00*Y_5*Z_2+Y*(5.0d00*Y_4*Z_2+Y*(Y*(3.0d00*Y_2*Z_2+ 2811 & Y*(Y*(Y_2+ 2812 & Z_2)+2.0d00*Y*Z_2))+4.0d00*Y_3*Z_2)))*R81+(6.0d00*Y* 2813 & (3.0d00*Y_2*Z_2 2814 & +Y*(Y*(Y_2+Z_2)+2.0d00*Y*Z_2))+Y*(Y*(Y*(6.0d00*Y_2+ 2815 & 3.0d00*Z_2)+4.0d00*Y* 2816 & (Y_2+Z_2)+8.0d00*Y*Z_2)+5.0d00*Y*(Y*(Y_2+Z_2)+ 2817 & 2.0d00*Y*Z_2)+15.0d00* 2818 & Y_2*Z_2)+24.0d00*Y_3*Z_2)*R71 2819 & +(Y*(45.0d00*Y_2+15.0d00*Z_2)+24.0d00*Y*(Y_2+Z_2)+ 2820 & 48.0d00*Y*Z_2+ 2821 & 6.0d00*Y*(6.0d00*Y_2+3.0d00*Z_2))*R61+105.0d00*Y*R51 2822 R(M,0,214)=Y_6*Z_3*R91 2823 & +(5.0d00*Y_4*Z_3+Y*(Y*(3.0d00*Y_2*Z_3+Y*(Y* 2824 & (3.0d00*Y_2*Z+Z_3)+ 2825 & 2.0d00*Y*Z_3))+4.0d00*Y_3*Z_3))*R81 2826 & +(Y*(Y*(18.0d00*Y_2*Z+3.0d00*Z_3)+4.0d00*Y*(3.0d00*Y_2*Z+ 2827 & Z_3)+8.0d00*Y*Z_3)+ 2828 & 5.0d00*Y*(Y*(3.0d00*Y_2*Z+Z_3)+2.0d00*Y*Z_3)+ 2829 & 15.0d00*Y_2*Z_3)*R71 2830 & +(135.0d00*Y_2*Z+15.0d00*Z_3)*R61+45.0d00*Z*R51 2831 R(M,0,215)=Y_5*Z_4*R91 2832 & +(Y*(3.0d00*Y_2*Z_4+Y*(Y*(Z_4+6.0d00*Y_2*Z_2)+ 2833 & 2.0d00*Y*Z_4))+ 2834 & 4.0d00*Y_3*Z_4)*R81+(Y*(Y*(Y*(3.0d00*Y_2+6.0d00*Z_2)+ 2835 & 12.0d00*Y*Z_2)+18.0d00*Y_2*Z_2+3.0d00*Z_4)+4.0d00*Y*(Z_4+ 2836 & 6.0d00*Y_2*Z_2)+ 2837 & 8.0d00*Y*Z_4)*R71 2838 & +(Y*(18.0d00*Y_2+18.0d00*Z_2)+4.0d00*Y*(3.0d00*Y_2+ 2839 & 6.0d00*Z_2)+48.0d00*Y*Z_2)* 2840 & R61+45.0d00*Y*R51 2841 R(M,0,216)=Y_4*Z_5*R91 2842 & +(3.0d00*Y_2*Z_5+Y*(Y*(Z_5+10.0d00* Y_2*Z_3)+ 2843 & 2.0d00*Y*Z_5))*R81 2844 & +(Y*(Y*(10.0d00*Z_3+15.0d00*Y_2*Z)+20.0d00*Y*Z_3)+ 2845 & 30.0d00*Y_2*Z_3+3.0d00*Z_5) 2846 & *R71+(90.0d00*Y_2*Z+30.0d00*Z_3)*R61+45.0d00*Z*R51 2847 R(M,0,217)=Y_3*Z_6*R91 2848 & +(Y*(Z_6+15.0d00*Y_2*Z_4)+2.0d00*Y*Z_6)*R81 2849 & +(Y*(45.0d00*Y_2*Z_2+15.0d00*Z_4)+30.0d00*Y*Z_4)*R71 2850 & +(Y*(15.0d00*Y_2+45.0d00*Z_2)+90.0d00*Y*Z_2)*R61+ 2851 & 45.0d00*Y*R51 2852 R(M,0,218)=Y_2*Z_7*R91+(Z_7+21.0d00*Y_2*Z_5)*R81 2853 & +(105.0d00*Y_2*Z_3+21.0d00*Z_5)*R71 2854 & +(105.0d00*Y_2*Z+105.0d00*Z_3)*R61+105.0d00*Z*R51 2855 R(M,0,219)=Y*Z_8*R91+28.0d00*Y*Z_6*R81 2856 & +210.0d00*Y*Z_4*R71+420.0d00*Y*Z_2*R61+105.0d00*Y*R51 2857 R(M,0,220)=Z_5*Z_4*R91+36.0d00*Z_7*R81+378.0d00*Z_5*R71 2858 & +1260.0d00*Z_3*R61+945.0d00*Z*R51 285900100 continue 2860c 2861 end 2862#endif 2863