blob: 3f116b934083d66a1287dfdfd266d1d4f1f962da (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
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
|