diff options
Diffstat (limited to 'src/ksq_odd.x')
| -rw-r--r-- | src/ksq_odd.x | 403 |
1 files changed, 403 insertions, 0 deletions
diff --git a/src/ksq_odd.x b/src/ksq_odd.x new file mode 100644 index 0000000..3f116b9 --- /dev/null +++ b/src/ksq_odd.x @@ -0,0 +1,403 @@ + o1 = mm**2 + o2 = 5.0000000000000d1*o1 + o3 = -8.0000000000000d1 + o2 + o4 = 1/o3 + o5 = 2.00000000000000d0*qgrd(j) + o6 = cos(o5) + o7 = 5.0000000000000d0*o6 + o8 = 3.00000000000000d0 + o7 + o9 = dnnnnngtil + o10 = mm*phigrd(k) + o11 = sin(o10) + o12 = sin(qgrd(j)) + o13 = o12**2 + o14 = o13**2 + o15 = 2.50000000000000d1*o1 + o16 = -4.0000000000000d1 + o15 + o17 = 1/o16 + o18 = dnnngtil + o19 = 2.00000000000000d1*o1 + o20 = -3.2000000000000d1 + o19 + o21 = 1/o20 + o22 = dngtil + o23 = 5.0000000000000d0*o1 + o24 = -8.0000000000000d0 + o23 + o25 = 1/o24 + o26 = cos(o10) + o27 = cos(qgrd(j)) + o28 = 1.03000000000000d2*o1 + o29 = -9.2800000000000d2 + o28 + o30 = o12*o13 + o31 = 3.5000000000000d1*o1 + o32 = -8.3000000000000d1 + o31 + o33 = -1.60000000000000d1 + o1 + o34 = dnnnnnnngtil + o35 = 4.0000000000000d0*qgrd(j) + o36 = cos(o35) + o37 = dnnnngtil + o38 = 1.00000000000000d1*o1 + o39 = -1.60000000000000d1 + o38 + o40 = 1/o39 + o41 = dnngtil + o42 = dnnnnnngtil + o43 = o27**2 + o44 = gtil + o45 = o13*o14 + o46 = 1.07000000000000d2*o1 + o47 = 5.3000000000000d1*o1 + o48 = sin(o5) + o49 = byJ**2 + o50 = o24**2 + o51 = 1/o50 + o52 = o26**2 + o53 = o37**2 + o54 = o1**2 + o55 = o6**2 + o56 = o36**2 + o57 = o39**2 + o58 = 1/o57 + o59 = o41**2 + o60 = o9**2 + o61 = o16**2 + o62 = 1/o61 + o63 = o18**2 + o64 = o20**2 + o65 = 1/o64 + o66 = o22**2 + o67 = o11**2 + o68 = o1*o54 + o69 = o3**2 + o70 = 1/o69 + o71 = o42**2 + o72 = o34**2 + o73 = o14**2 + o74 = o48**2 + o75 = o43**2 + o76 = o44**2 + t1 = -1.00800000000000d2*byJ*o13*o25*o26*o37 - 1.39500000000000d1* + & byJ*o1*o13*o25*o26*o37 - 9.7440000000000d2*byJ*o13*o25*o26*o36*o37 + + & 1.08150000000000d2*byJ*o1*o13*o25*o26*o36*o37 + 1.33200000000000d + & 3*byJ*o13*o26*o40*o41 - 4.6800000000000d2*byJ*o1*o13*o26*o40*o41 + 6 + & .9720000000000d3*byJ*o13*o26*o36*o40*o41 - 2.94000000000000d3*byJ*o1 + & *o13*o26*o36*o40*o41 - 9.6000000000000d2*byJ*o14*o26*o4*o42 + 6.00 + & 00000000000d1*byJ*o1*o14*o26*o4*o42 - 1.80000000000000d2*byJ*o14*o26 + & *o44 + 7.2000000000000d2*byJ*o13*o26*o43*o44 - 2.70000000000000d3* + & byJ*o14*o26*o43*o44 + 4.5000000000000d2*byJ*o26*o44*o45 + 1.80000000 + & 000000d1*o13*o49 - 3.7296000000000d3*o13*o25*o37*o40*o41*o52 + 7 + & .9425000000000d2*o1*o13*o25*o37*o40*o41*o52 - 5.5574400000000d4* + & o13*o25*o36*o37*o40*o41*o52 + 2.21991000000000d4*o1*o13*o25*o36* + & o37*o40*o41*o52 + 2.68800000000000d3*o14*o25*o37*o4*o42*o52 + 2. + & 04000000000000d2*o1*o14*o25*o37*o4*o42*o52 + 2.59840000000000d4* + & o14*o25*o36*o37*o4*o42*o52 - 4.5080000000000d3*o1*o14*o25*o36*o3 + & 7*o4*o42*o52 - 3.5520000000000d4*o14*o4*o40*o41*o42*o52 + 1.4700 + & 0000000000d4*o1*o14*o4*o40*o41*o42*o52 - 1.85920000000000d5*o14* + & o36*o4*o40*o41*o42*o52 + 9.0020000000000d4*o1*o14*o36*o4*o40*o41 + & *o42*o52 + t2 = -1.10400000000000d5*o14*o17*o18*o21*o22*o43*o52 + 2.9160000 + & 0000000d4*o1*o14*o17*o18*o21*o22*o43*o52 + 5.0400000000000d2*o14 + & *o25*o37*o44*o52 + 6.9750000000000d1*o1*o14*o25*o37*o44*o52 + 4. + & 8720000000000d3*o14*o25*o36*o37*o44*o52 - 5.4075000000000d2*o1*o + & 14*o25*o36*o37*o44*o52 - 6.6600000000000d3*o14*o40*o41*o44*o52 + + & 2.34000000000000d3*o1*o14*o40*o41*o44*o52 - 3.4860000000000d4*o + & 14*o36*o40*o41*o44*o52 + 1.47000000000000d4*o1*o14*o36*o40*o41*o + & 44*o52 - 2.01600000000000d3*o13*o25*o37*o43*o44*o52 - 2.79000000 + & 000000d2*o1*o13*o25*o37*o43*o44*o52 + 7.5600000000000d3*o14*o25* + & o37*o43*o44*o52 + 1.04625000000000d3*o1*o14*o25*o37*o43*o44*o52 + & - 1.94880000000000d4*o13*o25*o36*o37*o43*o44*o52 + 2.16300000000 + & 000d3*o1*o13*o25*o36*o37*o43*o44*o52 + 7.3080000000000d4*o14*o25 + & *o36*o37*o43*o44*o52 - 8.1112500000000d3*o1*o14*o25*o36*o37*o43* + & o44*o52 + 2.66400000000000d4*o13*o40*o41*o43*o44*o52 - 9.3600000 + & 000000d3*o1*o13*o40*o41*o43*o44*o52 - 9.9900000000000d4*o14*o40* + & o41*o43*o44*o52 + 3.5100000000000d4*o1*o14*o40*o41*o43*o44*o52 + + & 1.39440000000000d5*o13*o36*o40*o41*o43*o44*o52 - 5.880000000000 + & 0d4*o1*o13*o36*o40*o41*o43*o44*o52 - 5.2290000000000d5*o14*o36*o + & 40*o41*o43*o44*o52 + 2.20500000000000d5*o1*o14*o36*o40*o41*o43*o + & 44*o52 - 1.92000000000000d4*o14*o4*o42*o43*o44*o52 + t3 = 1.20000000000000d3*o1*o14*o4*o42*o43*o44*o52 - 7.6800000000 + & 000d3*o17*o21*o22*o34*o43*o45*o52 + 1.44000000000000d3*o1*o17*o2 + & 1*o22*o34*o43*o45*o52 - 1.26000000000000d3*o25*o37*o44*o45*o52 - + & 1.74375000000000d2*o1*o25*o37*o44*o45*o52 - 1.21800000000000d4* + & o25*o36*o37*o44*o45*o52 + 1.35187500000000d3*o1*o25*o36*o37*o44* + & o45*o52 + 1.66500000000000d4*o40*o41*o44*o45*o52 - 5.85000000000 + & 00d3*o1*o40*o41*o44*o45*o52 + 8.7150000000000d4*o36*o40*o41*o44* + & o45*o52 - 3.6750000000000d4*o1*o36*o40*o41*o44*o45*o52 + 4.80000 + & 00000000d3*o4*o42*o44*o45*o52 - 3.00000000000000d2*o1*o4*o42*o44 + & *o45*o52 + 7.2000000000000d4*o4*o42*o43*o44*o45*o52 - 4.50000000 + & 00000d3*o1*o4*o42*o43*o44*o45*o52 + 1.41120000000000d2*o13*o51*o + & 52*o53 + 3.9060000000000d1*o1*o13*o51*o52*o53 + 2.72832000000000 + & d3*o13*o36*o51*o52*o53 + 7.4760000000000d1*o1*o13*o36*o51*o52*o5 + & 3 + 1.81350000000000d2*o13*o25*o37*o40*o41*o52*o54 - 2.667000000 + & 00000d2*o13*o25*o36*o37*o40*o41*o52*o54 - 2.32500000000000d1*o14 + & *o25*o37*o4*o42*o52*o54 + 1.80250000000000d2*o14*o25*o36*o37*o4* + & o42*o52*o54 - 7.8000000000000d2*o14*o4*o40*o41*o42*o52*o54 - 4.9 + & 000000000000d3*o14*o36*o4*o40*o41*o42*o52*o54 - 1.92000000000000 + & d3*o14*o17*o18*o21*o22*o43*o52*o54 - 6.0000000000000d1*o17*o21*o + & 22*o34*o43*o45*o52*o54 + t4 = 2.70281250000000d0*o13*o51*o52*o53*o54 - 4.1908125000000d1* + & o13*o36*o51*o52*o53*o54 - 2.74432000000000d4*o13*o25*o37*o40*o41 + & *o52*o55 + 9.5432000000000d3*o1*o13*o25*o37*o40*o41*o52*o55 + 1. + & 14688000000000d4*o14*o25*o37*o4*o42*o52*o55 - 4.2560000000000d2* + & o1*o14*o25*o37*o4*o42*o52*o55 - 1.20064000000000d5*o14*o4*o40*o4 + & 1*o42*o52*o55 + 5.2304000000000d4*o1*o14*o4*o40*o41*o42*o52*o55 + & - 1.99200000000000d6*o14*o17*o18*o21*o22*o43*o52*o55 + 2.0850000 + & 0000000d6*o1*o14*o17*o18*o21*o22*o43*o52*o55 + 1.31072000000000d + & 3*o13*o51*o52*o53*o55 + 6.6560000000000d1*o1*o13*o51*o52*o53*o55 + & + 2.60000000000000d2*o13*o25*o37*o40*o41*o52*o54*o55 - 1.820000 + & 00000000d1*o14*o25*o37*o4*o42*o52*o54*o55 - 2.80000000000000d3*o + & 14*o4*o40*o41*o42*o52*o54*o55 - 5.2500000000000d5*o14*o17*o18*o2 + & 1*o22*o43*o52*o54*o55 + 8.4500000000000d-1*o13*o51*o52*o53*o54*o + & 55 - 1.88708800000000d5*o13*o25*o37*o40*o41*o52*o56 + 1.00521050 + & 000000d5*o1*o13*o25*o37*o40*o41*o52*o56 + 1.31868800000000d4*o13 + & *o51*o52*o53*o56 - 2.92726000000000d3*o1*o13*o51*o52*o53*o56 - 8 + & .8322500000000d3*o13*o25*o37*o40*o41*o52*o54*o56 + 1.62450312500 + & 000d2*o13*o51*o52*o53*o54*o56 + 2.46420000000000d4*o13*o52*o58*o + & 59 - 1.73160000000000d4*o1*o13*o52*o58*o59 + 2.57964000000000d5* + & o13*o36*o52*o58*o59 - 1.99416000000000d5*o1*o13*o36*o52*o58*o59 + t5 = 3.04200000000000d3*o13*o52*o54*o58*o59 + 3.8220000000000d4* + & o13*o36*o52*o54*o58*o59 + 1.43648000000000d5*o13*o52*o55*o58*o59 + & - 1.07200000000000d5*o1*o13*o52*o55*o58*o59 + 2.00000000000000d + & 4*o13*o52*o54*o55*o58*o59 + 6.7512200000000d5*o13*o52*o56*o58*o5 + & 9 - 5.6938000000000d5*o1*o13*o52*o56*o58*o59 + 1.20050000000000d + & 5*o13*o52*o54*o56*o58*o59 - 3.07200000000000d2*byJ*o13*o25*o26*o37 + & *o6 - 7.8000000000000d0*byJ*o1*o13*o25*o26*o37*o6 + 3.216000000000 + & 0d3*byJ*o13*o26*o40*o41*o6 - 1.20000000000000d3*byJ*o1*o13*o26*o40*o + & 41*o6 - 1.34400000000000d3*byJ*o14*o26*o4*o42*o6 + 8.4000000000000 + & d1*byJ*o1*o14*o26*o4*o42*o6 - 2.03712000000000d4*o13*o25*o37*o40*o + & 41*o52*o6 + 5.8188000000000d3*o1*o13*o25*o37*o40*o41*o52*o6 - 1. + & 46540800000000d5*o13*o25*o36*o37*o40*o41*o52*o6 + 6.571880000000 + & 0d4*o1*o13*o25*o36*o37*o40*o41*o52*o6 + 1.19552000000000d4*o14*o + & 25*o37*o4*o42*o52*o6 - 1.84000000000000d1*o1*o14*o25*o37*o4*o42* + & o52*o6 + 3.6377600000000d4*o14*o25*o36*o37*o4*o42*o52*o6 - 6.311 + & 2000000000d3*o1*o14*o25*o36*o37*o4*o42*o52*o6 - 1.35488000000000 + & d5*o14*o4*o40*o41*o42*o52*o6 + 5.7940000000000d4*o1*o14*o4*o40*o + & 41*o42*o52*o6 - 2.60288000000000d5*o14*o36*o4*o40*o41*o42*o52*o6 + & + 1.26028000000000d5*o1*o14*o36*o4*o40*o41*o42*o52*o6 + 9.50400 + & 00000000d5*o14*o17*o18*o21*o22*o43*o52*o6 + t6 = -6.3960000000000d5*o1*o14*o17*o18*o21*o22*o43*o52*o6 + 1.53 + & 600000000000d3*o14*o25*o37*o44*o52*o6 + 3.9000000000000d1*o1*o14 + & *o25*o37*o44*o52*o6 - 1.60800000000000d4*o14*o40*o41*o44*o52*o6 + & + 6.0000000000000d3*o1*o14*o40*o41*o44*o52*o6 - 6.1440000000000d + & 3*o13*o25*o37*o43*o44*o52*o6 - 1.56000000000000d2*o1*o13*o25*o37 + & *o43*o44*o52*o6 + 2.30400000000000d4*o14*o25*o37*o43*o44*o52*o6 + & + 5.8500000000000d2*o1*o14*o25*o37*o43*o44*o52*o6 + 6.4320000000 + & 000d4*o13*o40*o41*o43*o44*o52*o6 - 2.40000000000000d4*o1*o13*o40 + & *o41*o43*o44*o52*o6 - 2.41200000000000d5*o14*o40*o41*o43*o44*o52 + & *o6 + 9.0000000000000d4*o1*o14*o40*o41*o43*o44*o52*o6 - 2.688000 + & 00000000d4*o14*o4*o42*o43*o44*o52*o6 + 1.68000000000000d3*o1*o14 + & *o4*o42*o43*o44*o52*o6 + 3.8400000000000d4*o17*o21*o22*o34*o43*o + & 45*o52*o6 - 2.64000000000000d4*o1*o17*o21*o22*o34*o43*o45*o52*o6 + & - 3.8400000000000d3*o25*o37*o44*o45*o52*o6 - 9.7500000000000d1* + & o1*o25*o37*o44*o45*o52*o6 + 4.0200000000000d4*o40*o41*o44*o45*o5 + & 2*o6 - 1.50000000000000d4*o1*o40*o41*o44*o45*o52*o6 + 6.72000000 + & 00000d3*o4*o42*o44*o45*o52*o6 - 4.2000000000000d2*o1*o4*o42*o44* + & o45*o52*o6 + 1.00800000000000d5*o4*o42*o43*o44*o45*o52*o6 - 6.30 + & 00000000000d3*o1*o4*o42*o43*o44*o45*o52*o6 + 8.6016000000000d2*o + & 13*o51*o52*o53*o6 + 1.40880000000000d2*o1*o13*o51*o52*o53*o6 + t7 = 8.3148800000000d3*o13*o36*o51*o52*o53*o6 - 7.1176000000000d + & 2*o1*o13*o36*o51*o52*o53*o6 + 5.6640000000000d2*o13*o25*o37*o40* + & o41*o52*o54*o6 - 2.96800000000000d3*o13*o25*o36*o37*o40*o41*o52* + & o54*o6 - 4.5550000000000d1*o14*o25*o37*o4*o42*o52*o54*o6 + 2.523 + & 50000000000d2*o14*o25*o36*o37*o4*o42*o52*o54*o6 - 3.092000000000 + & 00d3*o14*o4*o40*o41*o42*o52*o54*o6 - 6.8600000000000d3*o14*o36*o + & 4*o40*o41*o42*o52*o54*o6 + 6.9000000000000d4*o14*o17*o18*o21*o22 + & *o43*o52*o54*o6 + 1.50000000000000d3*o17*o21*o22*o34*o43*o45*o52 + & *o54*o6 + 3.02250000000000d0*o13*o51*o52*o53*o54*o6 - 2.34325000 + & 000000d1*o13*o36*o51*o52*o53*o54*o6 + 1.18992000000000d5*o13*o52 + & *o58*o59*o6 - 8.6208000000000d4*o1*o13*o52*o58*o59*o6 + 6.228320 + & 0000000d5*o13*o36*o52*o58*o59*o6 - 4.9504000000000d5*o1*o13*o36* + & o52*o58*o59*o6 + 1.56000000000000d4*o13*o52*o54*o58*o59*o6 + 9.8 + & 000000000000d4*o13*o36*o52*o54*o58*o59*o6 + 6.1952000000000d2*o1 + & 4*o43*o51*o52*o60 - 1.93600000000000d2*o1*o14*o43*o51*o52*o60 + + & 1.51250000000000d1*o14*o43*o51*o52*o54*o60 + 4.3059200000000d3*o + & 14*o43*o51*o52*o55*o60 - 9.5584000000000d2*o1*o14*o43*o51*o52*o5 + & 5*o60 + 5.3045000000000d1*o14*o43*o51*o52*o54*o55*o60 - 3.266560 + & 0000000d3*o14*o43*o51*o52*o6*o60 + 8.7296000000000d2*o1*o14*o43* + & o51*o52*o6*o60 - 5.6650000000000d1*o14*o43*o51*o52*o54*o6*o60 + t8 = 1.47200000000000d4*o18*o34*o43*o45*o52*o62 - 2.968000000000 + & 00d3*o1*o18*o34*o43*o45*o52*o62 + 1.28000000000000d2*o18*o34*o43 + & *o45*o52*o54*o62 - 5.3120000000000d4*o18*o34*o43*o45*o52*o6*o62 + & + 2.57200000000000d4*o1*o18*o34*o43*o45*o52*o6*o62 - 1.400000000 + & 00000d3*o18*o34*o43*o45*o52*o54*o6*o62 + 1.05800000000000d5*o14* + & o43*o52*o62*o63 - 2.94400000000000d4*o1*o14*o43*o52*o62*o63 + 2. + & 04800000000000d3*o14*o43*o52*o54*o62*o63 + 1.37780000000000d6*o1 + & 4*o43*o52*o55*o62*o63 - 1.16200000000000d6*o1*o14*o43*o52*o55*o6 + & 2*o63 + 2.45000000000000d5*o14*o43*o52*o54*o55*o62*o63 - 7.63600 + & 00000000d5*o14*o43*o52*o6*o62*o63 + 4.2824000000000d5*o1*o14*o43 + & *o52*o6*o62*o63 - 4.4800000000000d4*o14*o43*o52*o54*o6*o62*o63 + + & 2.88000000000000d4*o14*o43*o52*o65*o66 - 7.2000000000000d3*o1*o + & 14*o43*o52*o65*o66 + 4.5000000000000d2*o14*o43*o52*o54*o65*o66 + + & 7.2000000000000d5*o14*o43*o52*o55*o65*o66 - 9.0000000000000d5*o + & 1*o14*o43*o52*o55*o65*o66 + 2.81250000000000d5*o14*o43*o52*o54*o + & 55*o65*o66 - 2.88000000000000d5*o14*o43*o52*o6*o65*o66 + 2.16000 + & 000000000d5*o1*o14*o43*o52*o6*o65*o66 - 2.25000000000000d4*o14*o + & 43*o52*o54*o6*o65*o66 - 2.10600000000000d4*o1*o14*o17*o18*o21*o2 + & 2*o67 - 6.7562000000000d4*o1*o14*o17*o18*o22*o25*o67 - 7.0180000 + & 000000d4*o1*o14*o17*o18*o22*o40*o67 + 1.92000000000000d3*o1*o14* + & o17*o42*o43*o44*o67 + t9 = -3.3600000000000d3*o1*o17*o42*o43*o44*o45*o67 - 5.440000000 + & 0000d2*o1*o17*o25*o27*o30*o37*o42*o48*o67 + 1.98400000000000d3*o + & 1*o17*o27*o30*o40*o41*o42*o48*o67 - 1.02000000000000d3*o1*o12*o2 + & 5*o27*o37*o44*o48*o67 + 1.78500000000000d3*o1*o25*o27*o30*o37*o4 + & 4*o48*o67 + 3.7200000000000d3*o1*o12*o27*o40*o41*o44*o48*o67 - 6 + & .5100000000000d3*o1*o27*o30*o40*o41*o44*o48*o67 + 4.625200000000 + & 0d4*o14*o17*o18*o22*o25*o54*o67 + 6.6494000000000d4*o14*o17*o18* + & o22*o40*o54*o67 - 1.20000000000000d2*o14*o17*o42*o43*o44*o54*o67 + & + 2.10000000000000d2*o17*o42*o43*o44*o45*o54*o67 + 1.2200000000 + & 0000d2*o17*o25*o27*o30*o37*o42*o48*o54*o67 - 5.7200000000000d2*o + & 17*o27*o30*o40*o41*o42*o48*o54*o67 + 1.65000000000000d2*o12*o25* + & o27*o37*o44*o48*o54*o67 - 2.88750000000000d2*o25*o27*o30*o37*o44 + & *o48*o54*o67 - 8.4000000000000d2*o12*o27*o40*o41*o44*o48*o54*o67 + & + 1.47000000000000d3*o27*o30*o40*o41*o44*o48*o54*o67 - 5.850000 + & 0000000d4*o1*o14*o17*o18*o21*o22*o55*o67 - 1.42506000000000d5*o1 + & *o14*o17*o18*o22*o25*o55*o67 - 1.31652000000000d5*o1*o14*o17*o18 + & *o22*o40*o55*o67 + 9.2016000000000d4*o14*o17*o18*o22*o25*o54*o55 + & *o67 + 1.25082000000000d5*o14*o17*o18*o22*o40*o54*o55*o67 - 7.02 + & 00000000000d4*o1*o14*o17*o18*o21*o22*o6*o67 - 1.96332000000000d5 + & *o1*o14*o17*o18*o22*o25*o6*o67 - 1.92408000000000d5*o1*o14*o17*o + & 18*o22*o40*o6*o67 - 1.29920000000000d3*o1*o17*o25*o27*o30*o37*o4 + & 2*o48*o6*o67 + 1.33760000000000d4*o1*o17*o27*o30*o40*o41*o42*o48 + & *o6*o67 + t10 = -2.43600000000000d3*o1*o12*o25*o27*o37*o44*o48*o6*o67 + 4. + & 2630000000000d3*o1*o25*o27*o30*o37*o44*o48*o6*o67 + 2.5080000000 + & 0000d4*o1*o12*o27*o40*o41*o44*o48*o6*o67 - 4.3890000000000d4*o1* + & o27*o30*o40*o41*o44*o48*o6*o67 + 1.30692000000000d5*o14*o17*o18* + & o22*o25*o54*o6*o67 + 1.82664000000000d5*o14*o17*o18*o22*o40*o54* + & o6*o67 + 1.08400000000000d2*o17*o25*o27*o30*o37*o42*o48*o54*o6*o + & 67 - 6.5320000000000d3*o17*o27*o30*o40*o41*o42*o48*o54*o6*o67 + + & 5.1000000000000d1*o12*o25*o27*o37*o44*o48*o54*o6*o67 - 8.9250000 + & 000000d1*o25*o27*o30*o37*o44*o48*o54*o6*o67 - 1.06800000000000d4 + & *o12*o27*o40*o41*o44*o48*o54*o6*o67 + 1.86900000000000d4*o27*o30 + & *o40*o41*o44*o48*o54*o6*o67 + 4.0960000000000d3*o1*o14*o60*o62*o + & 67 - 5.1200000000000d2*o14*o54*o60*o62*o67 + 9.2160000000000d3*o + & 1*o14*o55*o60*o62*o67 - 1.15200000000000d3*o14*o54*o55*o60*o62*o + & 67 + 1.22880000000000d4*o1*o14*o6*o60*o62*o67 - 1.53600000000000 + & d3*o14*o54*o6*o60*o62*o67 + 2.63438000000000d5*o1*o14*o62*o63*o6 + & 7 - 1.49158000000000d5*o14*o54*o62*o63*o67 + 5.8030200000000d5*o + & 1*o14*o55*o62*o63*o67 - 3.2531400000000d5*o14*o54*o55*o62*o63*o6 + & 7 + 7.8147600000000d5*o1*o14*o6*o62*o63*o67 - 4.4056800000000d5* + & o14*o54*o6*o62*o63*o67 + 6.8890000000000d3*o1*o14*o51*o66*o67 - + & 5.8100000000000d3*o14*o51*o54*o66*o67 + 1.36890000000000d4*o1*o1 + & 4*o51*o55*o66*o67 + t11 = -1.05300000000000d4*o14*o51*o54*o55*o66*o67 + 1.4641000000 + & 0000d4*o1*o14*o58*o66*o67 - 1.69400000000000d4*o14*o54*o58*o66*o + & 67 + 2.52810000000000d4*o1*o14*o55*o58*o66*o67 - 2.8620000000000 + & 0d4*o14*o54*o55*o58*o66*o67 + 1.94220000000000d4*o1*o14*o51*o6*o + & 66*o67 - 1.56600000000000d4*o14*o51*o54*o6*o66*o67 + 3.847800000 + & 0000d4*o1*o14*o58*o6*o66*o67 - 4.4040000000000d4*o14*o54*o58*o6* + & o66*o67 + 8.1000000000000d3*o1*o14*o65*o66*o67 + 2.2500000000000 + & 0d4*o1*o14*o55*o65*o66*o67 + 2.70000000000000d4*o1*o14*o6*o65*o6 + & 6*o67 - 7.4900000000000d3*o14*o17*o18*o22*o25*o67*o68 - 1.498000 + & 00000000d4*o14*o17*o18*o22*o40*o67*o68 - 5.5000000000000d0*o17*o + & 25*o27*o30*o37*o42*o48*o67*o68 + 2.80000000000000d1*o17*o27*o30* + & o40*o41*o42*o48*o67*o68 - 1.43100000000000d4*o14*o17*o18*o22*o25 + & *o55*o67*o68 - 2.86200000000000d4*o14*o17*o18*o22*o40*o55*o67*o6 + & 8 - 2.07600000000000d4*o14*o17*o18*o22*o25*o6*o67*o68 - 4.152000 + & 0000000d4*o14*o17*o18*o22*o40*o6*o67*o68 - 1.70000000000000d0*o1 + & 7*o25*o27*o30*o37*o42*o48*o6*o67*o68 + 3.5600000000000d2*o17*o27 + & *o30*o40*o41*o42*o48*o6*o67*o68 + 1.60000000000000d1*o14*o60*o62 + & *o67*o68 + 3.6000000000000d1*o14*o55*o60*o62*o67*o68 + 4.8000000 + & 000000d1*o14*o6*o60*o62*o67*o68 + 2.28980000000000d4*o14*o62*o63 + & *o67*o68 + 5.0562000000000d4*o14*o55*o62*o63*o67*o68 + t12 = 6.8052000000000d4*o14*o6*o62*o63*o67*o68 + 1.2250000000000 + & 0d3*o14*o51*o66*o67*o68 + 2.02500000000000d3*o14*o51*o55*o66*o67 + & *o68 + 4.9000000000000d3*o14*o58*o66*o67*o68 + 8.1000000000000d3 + & *o14*o55*o58*o66*o67*o68 + 3.15000000000000d3*o14*o51*o6*o66*o67 + & *o68 + 1.26000000000000d4*o14*o58*o6*o66*o67*o68 + 1.42420000000 + & 000d4*o1*o14*o60*o67*o70 - 1.90400000000000d3*o14*o54*o60*o67*o7 + & 0 + 3.15540000000000d4*o1*o14*o55*o60*o67*o70 - 4.2480000000000d + & 3*o14*o54*o55*o60*o67*o70 + 4.2396000000000d4*o1*o14*o6*o60*o67* + & o70 - 5.6880000000000d3*o14*o54*o6*o60*o67*o70 + 6.4000000000000 + & d1*o14*o60*o67*o68*o70 + 1.44000000000000d2*o14*o55*o60*o67*o68* + & o70 + 1.92000000000000d2*o14*o6*o60*o67*o68*o70 + 5.120000000000 + & 0d2*o1*o43*o45*o62*o67*o71 - 6.4000000000000d1*o43*o45*o54*o62*o + & 67*o71 + 2.00000000000000d0*o43*o45*o62*o67*o68*o71 + 1.28000000 + & 000000d4*o45*o52*o70*o71 - 1.60000000000000d3*o1*o45*o52*o70*o71 + & + 5.0000000000000d1*o45*o52*o54*o70*o71 + 2.50880000000000d4*o4 + & 5*o52*o55*o70*o71 - 3.13600000000000d3*o1*o45*o52*o55*o70*o71 + + & 9.8000000000000d1*o45*o52*o54*o55*o70*o71 + 3.5840000000000d4*o4 + & 5*o52*o6*o70*o71 - 4.4800000000000d3*o1*o45*o52*o6*o70*o71 + 1.4 + & 0000000000000d2*o45*o52*o54*o6*o70*o71 + t13 = -1.20000000000000d4*o4*o42*o44*o52*o73 + 7.5000000000000d2 + & *o1*o4*o42*o44*o52*o73 - 1.68000000000000d4*o4*o42*o44*o52*o6*o7 + & 3 + 1.05000000000000d3*o1*o4*o42*o44*o52*o6*o73 + 5.120000000000 + & 0d2*o43*o52*o62*o72*o73 - 6.4000000000000d1*o1*o43*o52*o62*o72*o + & 73 + 2.00000000000000d0*o43*o52*o54*o62*o72*o73 - 1.054000000000 + & 00d3*o1*o25*o37*o40*o41*o67*o74 + 1.44500000000000d2*o1*o51*o53* + & o67*o74 + 4.0850000000000d2*o25*o37*o40*o41*o54*o67*o74 - 4.6750 + & 000000000d1*o51*o53*o54*o67*o74 - 1.69708000000000d4*o1*o25*o37* + & o40*o41*o55*o67*o74 + 8.2418000000000d2*o1*o51*o53*o55*o67*o74 + + & 7.5821000000000d3*o25*o37*o40*o41*o54*o55*o67*o74 - 3.451000000 + & 0000d1*o51*o53*o54*o55*o67*o74 + 1.92200000000000d3*o1*o58*o59*o + & 67*o74 - 8.6800000000000d2*o54*o58*o59*o67*o74 + 8.7362000000000 + & d4*o1*o55*o58*o59*o67*o74 - 7.4404000000000d4*o54*o55*o58*o59*o6 + & 7*o74 - 9.6232000000000d3*o1*o25*o37*o40*o41*o6*o67*o74 + 6.9020 + & 000000000d2*o1*o51*o53*o6*o67*o74 + 4.7966000000000d3*o25*o37*o4 + & 0*o41*o54*o6*o67*o74 - 1.26100000000000d2*o51*o53*o54*o6*o67*o74 + & + 2.59160000000000d4*o1*o58*o59*o6*o67*o74 - 1.68880000000000d4 + & *o54*o58*o59*o6*o67*o74 - 3.8500000000000d1*o25*o37*o40*o41*o67* + & o68*o74 + 3.7812500000000d0*o51*o53*o67*o68*o74 + t14 = -1.51300000000000d2*o25*o37*o40*o41*o55*o67*o68*o74 + 3.61 + & 25000000000d-1*o51*o53*o55*o67*o68*o74 + 9.8000000000000d1*o58*o + & 59*o67*o68*o74 + 1.58420000000000d4*o55*o58*o59*o67*o68*o74 - 5. + & 0140000000000d2*o25*o37*o40*o41*o6*o67*o68*o74 + 2.3375000000000 + & 0d0*o51*o53*o6*o67*o68*o74 + 2.49200000000000d3*o58*o59*o6*o67*o + & 68*o74 - 3.6000000000000d3*o14*o43*o52*o76 + 4.5000000000000d2*o + & 45*o52*o76 + 2.25000000000000d4*o43*o45*o52*o76 + 1.800000000000 + & 00d3*o1*o13*o43*o67*o76 - 6.3000000000000d3*o1*o14*o43*o67*o76 + + & 5.5125000000000d3*o1*o43*o45*o67*o76 - 2.25000000000000d3*o52*o + & 73*o76 + 2.81250000000000d3*o13*o52*o73*o76 - 3.3750000000000d4* + & o43*o52*o73*o76 + 7.2000000000000d3*o13*o52*o75*o76 - 5.40000000 + & 00000d4*o14*o52*o75*o76 + 1.01250000000000d5*o45*o52*o75*o76 - 1 + & .61920000000000d4*o14*o17*o18*o25*o43*o52*o9 + 4.7828000000000d3 + & *o1*o14*o17*o18*o25*o43*o52*o9 + 8.4480000000000d3*o14*o21*o22*o + & 25*o43*o52*o9 - 2.37600000000000d3*o1*o14*o21*o22*o25*o43*o52*o9 + & - 1.12640000000000d3*o17*o25*o34*o43*o45*o52*o9 + 2.46400000000 + & 000d2*o1*o17*o25*o34*o43*o45*o52*o9 - 3.5200000000000d2*o14*o17* + & o18*o25*o43*o52*o54*o9 + 1.65000000000000d2*o14*o21*o22*o25*o43* + & o52*o54*o9 + t15 = -1.10000000000000d1*o17*o25*o34*o43*o45*o52*o54*o9 - 1.540 + & 48000000000d5*o14*o17*o18*o25*o43*o52*o55*o9 + 8.2058000000000d4 + & *o1*o14*o17*o18*o25*o43*o52*o55*o9 + 1.11360000000000d5*o14*o21* + & o22*o25*o43*o52*o55*o9 - 8.1960000000000d4*o1*o14*o21*o22*o25*o4 + & 3*o52*o55*o9 - 7.2100000000000d3*o14*o17*o18*o25*o43*o52*o54*o55 + & *o9 + 7.7250000000000d3*o14*o21*o22*o25*o43*o52*o54*o55*o9 + 1.0 + & 1120000000000d5*o14*o17*o18*o25*o43*o52*o6*o9 - 4.4447200000000d + & 4*o1*o14*o17*o18*o25*o43*o52*o6*o9 - 6.4512000000000d4*o14*o21*o + & 22*o25*o43*o52*o6*o9 + 3.8256000000000d4*o1*o14*o21*o22*o25*o43* + & o52*o6*o9 + 2.96960000000000d3*o17*o25*o34*o43*o45*o52*o6*o9 - 5 + & .1520000000000d2*o1*o17*o25*o34*o43*o45*o52*o6*o9 + 4.5092000000 + & 000d3*o14*o17*o18*o25*o43*o52*o54*o6*o9 - 4.4340000000000d3*o14* + & o21*o22*o25*o43*o52*o54*o6*o9 + 2.06000000000000d1*o17*o25*o34*o + & 43*o45*o52*o54*o6*o9 + 1.06240000000000d4*o1*o14*o17*o22*o25*o67 + & *o9 - 7.1126000000000d4*o1*o14*o17*o18*o4*o67*o9 + 1.62000000000 + & 000d3*o1*o14*o21*o22*o4*o67*o9 + 2.87980000000000d4*o1*o14*o22*o + & 4*o40*o67*o9 - 5.1440000000000d3*o14*o17*o22*o25*o54*o67*o9 + 3. + & 01060000000000d4*o14*o17*o18*o4*o54*o67*o9 - 1.85960000000000d4* + & o14*o22*o4*o40*o54*o67*o9 + 2.24640000000000d4*o1*o14*o17*o22*o2 + & 5*o55*o67*o9 - 1.52406000000000d5*o1*o14*o17*o18*o4*o55*o67*o9 + + & 4.5000000000000d3*o1*o14*o21*o22*o4*o55*o67*o9 + 5.628600000000 + & 0d4*o1*o14*o22*o4*o40*o55*o67*o9 + t16 = -1.00440000000000d4*o14*o17*o22*o25*o54*o55*o67*o9 + 6.622 + & 2000000000d4*o14*o17*o18*o4*o54*o55*o67*o9 - 3.5676000000000d4*o + & 14*o22*o4*o40*o54*o55*o67*o9 + 3.09120000000000d4*o1*o14*o17*o22 + & *o25*o6*o67*o9 - 2.08212000000000d5*o1*o14*o17*o18*o4*o6*o67*o9 + & + 5.4000000000000d3*o1*o14*o21*o22*o4*o6*o67*o9 + 8.067600000000 + & 0d4*o1*o14*o22*o4*o40*o6*o67*o9 - 1.44120000000000d4*o14*o17*o22 + & *o25*o54*o6*o67*o9 + 8.9304000000000d4*o14*o17*o18*o4*o54*o6*o67 + & *o9 - 5.1648000000000d4*o14*o22*o4*o40*o54*o6*o67*o9 - 5.2096000 + & 000000d4*o1*o14*o18*o62*o67*o9 + 1.69520000000000d4*o14*o18*o54* + & o62*o67*o9 - 1.16928000000000d5*o1*o14*o18*o55*o62*o67*o9 + 3.78 + & 36000000000d4*o14*o18*o54*o55*o62*o67*o9 - 1.56096000000000d5*o1 + & *o14*o18*o6*o62*o67*o9 + 5.0652000000000d4*o14*o18*o54*o6*o62*o6 + & 7*o9 + 2.80000000000000d2*o14*o17*o22*o25*o67*o68*o9 - 1.7120000 + & 0000000d3*o14*o17*o18*o4*o67*o68*o9 + 1.12000000000000d3*o14*o22 + & *o4*o40*o67*o68*o9 + 5.4000000000000d2*o14*o17*o22*o25*o55*o67*o + & 68*o9 - 3.8160000000000d3*o14*o17*o18*o4*o55*o67*o68*o9 + 2.1600 + & 0000000000d3*o14*o22*o4*o40*o55*o67*o68*o9 + 7.8000000000000d2*o + & 14*o17*o22*o25*o6*o67*o68*o9 - 5.1120000000000d3*o14*o17*o18*o4* + & o6*o67*o68*o9 + 3.12000000000000d3*o14*o22*o4*o40*o6*o67*o68*o9 + & - 8.5600000000000d2*o14*o18*o62*o67*o68*o9 - 1.90800000000000d3* + & o14*o18*o55*o62*o67*o68*o9 - 2.55600000000000d3*o14*o18*o6*o62*o + & 67*o68*o9 + ksq(i,j,k) = t1 + t2 + t3 + t4 + t5 + t6 + t7 + t8 + t9 + t10 + + & t11 + t12 + t13 + t14 + t15 + t16 |