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
|
fp t620,
t649,
t618,
t648,
t619,
t647,
t646,
t598,
t639,
t535,
t645,
t593,
t644,
t610,
t640,
t553,
t583,
t643,
t642,
t641,
t638,
t637,
t636,
t587,
t566,
t546,
t635,
t634,
t584,
t633,
t602,
t614,
t631,
t556,
t632,
t605,
t545,
t630,
t582,
t530,
t629,
t628,
t611,
t627,
t604,
t568,
t626,
t591,
t570,
t586,
t625,
t606,
t578,
t624,
t613,
t623,
t622,
t617,
t616,
t615,
t612,
t609,
t608,
t607,
t603,
t601,
t600,
t599,
t597,
t596,
t595,
t594,
t592,
t590,
t589,
t585,
t579,
t577,
t576,
t575,
t573,
t572,
t571,
t569,
t567,
t565,
t564,
t563,
t561,
t560,
t559,
t555,
t554,
t551,
t548,
t547,
t544,
t543,
t541,
t537,
t536,
t531,
t529,
t528,
t524,
t521;
t620 = x*x;
t649 = y*t620;
t618 = y*y;
t648 = x*t618;
t619 = t618*y;
t647 = x*t619;
t646 = t619+y;
t598 = RATIONAL(-1.0,12.0);
t639 = RATIONAL(-7.0,48.0)*t620;
t535 = t619*t639;
t645 = t535+t598*t648;
t593 = RATIONAL(1.0,48.0);
t644 = t593*t620;
t610 = RATIONAL(1.0,3.0);
t640 = t610*t620;
t553 = t619*t640;
t583 = RATIONAL(5.0,18.0);
t643 = t553+t583*t648;
t642 = t620*t619;
t641 = t618*t620;
t638 = RATIONAL(-16.0,3.0)*t642+RATIONAL(-70.0,9.0)*t648;
t637 = RATIONAL(7.0,9.0)*t648+t553;
t636 = RATIONAL(-1.0,48.0)*t620;
t587 = RATIONAL(7.0,18.0);
t566 = t587*t619;
t546 = RATIONAL(-7.0,3.0)*t642;
t635 = t546+RATIONAL(-35.0,6.0)*t648;
t634 = RATIONAL(7.0,48.0)*t620;
t584 = RATIONAL(-1.0,36.0);
t633 = t619*t636+t584*t648;
t602 = RATIONAL(-7.0,18.0);
t614 = RATIONAL(-1.0,3.0);
t631 = t614*t620;
t556 = t619*t631;
t632 = t556+t602*t648;
t605 = RATIONAL(7.0,3.0);
t545 = t605*t642;
t630 = t605*t648+t545;
t582 = RATIONAL(5.0,24.0);
t530 = t619*t634;
t629 = t582*t648+t530;
t628 = x*t646;
t611 = RATIONAL(-4.0,9.0);
t627 = RATIONAL(25.0,6.0)*t648+t545+t611*y;
t604 = RATIONAL(1.0,18.0);
t568 = t604*y;
t626 = t556+RATIONAL(-5.0,9.0)*t648+t568;
t591 = RATIONAL(-1.0,18.0);
t570 = t591*y;
t586 = RATIONAL(-5.0,12.0);
t625 = t570+t535+t586*t648;
t606 = RATIONAL(1.0,6.0);
t578 = RATIONAL(-1.0,144.0);
t624 = t578*y+t530+t606*t648;
t613 = RATIONAL(-5.0,3.0);
t623 = t613*t648+t546+t570;
t622 = RATIONAL(-49.0,48.0)*t642+RATIONAL(-5.0,4.0)*t648+t568;
t617 = RATIONAL(-2.0,9.0);
t616 = RATIONAL(2.0,9.0);
t615 = RATIONAL(-2.0,3.0);
t612 = RATIONAL(-1.0,9.0);
t609 = RATIONAL(2.0,3.0);
t608 = RATIONAL(4.0,9.0);
t607 = RATIONAL(1.0,9.0);
t603 = RATIONAL(-5.0,18.0);
t601 = RATIONAL(1.0,72.0);
t600 = RATIONAL(7.0,36.0);
t599 = RATIONAL(-1.0,24.0);
t597 = RATIONAL(7.0,12.0);
t596 = RATIONAL(-7.0,12.0);
t595 = RATIONAL(1.0,24.0);
t594 = RATIONAL(-1.0,72.0);
t592 = RATIONAL(1.0,36.0);
t590 = RATIONAL(-5.0,24.0);
t589 = RATIONAL(1.0,12.0);
t585 = RATIONAL(-7.0,36.0);
t579 = RATIONAL(10.0,3.0);
t577 = RATIONAL(1.0,144.0);
t576 = t612*t619;
t575 = t608*y;
t573 = t607*t619;
t572 = RATIONAL(-8.0,9.0)*t619;
t571 = RATIONAL(8.0,9.0)*t619;
t569 = t591*t619;
t567 = t604*t619;
t565 = t602*t619;
t564 = t577*t619;
t563 = t577*y;
t561 = t578*t619;
t560 = RATIONAL(-7.0,144.0)*t619;
t559 = RATIONAL(7.0,144.0)*t619;
t555 = y*t631;
t554 = y*t640;
t551 = RATIONAL(-1.0,6.0)*t649;
t548 = RATIONAL(7.0,6.0)*t649;
t547 = RATIONAL(-8.0,3.0)*t649;
t544 = RATIONAL(8.0,3.0)*t649;
t543 = RATIONAL(-7.0,6.0)*t649;
t541 = t606*t649;
t537 = y*t636;
t536 = y*t639;
t531 = y*t634;
t529 = t593*t642;
t528 = y*t644;
t524 = RATIONAL(16.0,3.0)*t642;
t521 = RATIONAL(49.0,48.0)*t642;
coeff_dx_m2_m2 = t529+t564+t528+t563+(t599*t620+t594)*t618+(t604*t618+
t646*t584)*x;
coeff_dx_m1_m2 = t569+t536+(t607+RATIONAL(7.0,24.0)*t620)*t618+t582*t628+
t625;
coeff_dx_0_m2 = t615*t641+t554+t602*t628+t637;
coeff_dx_p1_m2 = t555+t567+(t612+t609*t620)*t618+t583*t628+t626;
coeff_dx_p2_m2 = t561+t531+(t601+RATIONAL(-7.0,24.0)*t620)*t618+t598*t628
+t624;
coeff_dx_p3_m2 = t537+t595*t641+t601*t628+t633;
coeff_dx_m2_m1 = t551+t560+(RATIONAL(5.0,48.0)+RATIONAL(5.0,16.0)*t620)*
t618+(t600*t619+t616*y)*x+t625;
coeff_dx_m1_m1 = t548+t575+t566+t521+(RATIONAL(-35.0,16.0)*t620+RATIONAL(
-5.0,6.0))*t618+(t613*y+RATIONAL(25.0,8.0)*t618+RATIONAL(-35.0,24.0)*t619)*x;
coeff_dx_0_m1 = t547+RATIONAL(5.0,1.0)*t641+(RATIONAL(28.0,9.0)*y+
RATIONAL(49.0,18.0)*t619)*x+t635;
coeff_dx_p1_m1 = t565+t544+(RATIONAL(-5.0,1.0)*t620+RATIONAL(5.0,6.0))*
t618+(RATIONAL(-35.0,18.0)*t619+RATIONAL(-20.0,9.0)*y)*x+t627;
coeff_dx_p2_m1 = t543+t559+(RATIONAL(35.0,16.0)*t620+RATIONAL(-5.0,48.0))
*t618+(t597*t619+t609*y)*x+t622;
coeff_dx_p3_m1 = t541+RATIONAL(-5.0,16.0)*t641+(RATIONAL(-7.0,72.0)*t619+
t612*y)*x+t629;
coeff_dx_m2_0 = t589+t573+RATIONAL(1.0,4.0)*t620+(t585+t596*t620)*t618+(
t614+t611*t619)*x+t637;
coeff_dx_m1_0 = RATIONAL(-7.0,4.0)*t620+t615+t572+(RATIONAL(14.0,9.0)+
RATIONAL(49.0,12.0)*t620)*t618+(t579*t619+RATIONAL(5.0,2.0))*x+t635;
coeff_dx_0_0 = t524+(RATIONAL(-28.0,3.0)*t618+RATIONAL(4.0,1.0))*t620+(
RATIONAL(-14.0,3.0)+RATIONAL(98.0,9.0)*t618+RATIONAL(-56.0,9.0)*t619)*x;
coeff_dx_p1_0 = RATIONAL(-4.0,1.0)*t620+t571+t609+(RATIONAL(28.0,3.0)*
t620+RATIONAL(-14.0,9.0))*t618+(t579+RATIONAL(40.0,9.0)*t619)*x+t638;
coeff_dx_p2_0 = t598+RATIONAL(-4.0,3.0)*t647-x+t576+RATIONAL(7.0,4.0)*
t620+(RATIONAL(-49.0,12.0)*t620+t600)*t618+t630;
coeff_dx_p3_0 = (t597*t618+RATIONAL(-1.0,4.0))*t620+(t616*t619+t606)*x+
t632;
coeff_dx_m2_p1 = t576+t541+(RATIONAL(5.0,12.0)*t620+RATIONAL(5.0,36.0))*
t618+(t617*y+t608*t619)*x+t626;
coeff_dx_m1_p1 = t571+t543+(RATIONAL(-10.0,9.0)+RATIONAL(-35.0,12.0)*t620
)*t618+(RATIONAL(-10.0,3.0)*t619+RATIONAL(5.0,3.0)*y)*x+t627;
coeff_dx_0_p1 = t544+RATIONAL(20.0,3.0)*t641+(RATIONAL(56.0,9.0)*t619+
RATIONAL(-28.0,9.0)*y)*x+t638;
coeff_dx_p1_p1 = t575+t547+t572+t524+(RATIONAL(-20.0,3.0)*t620+RATIONAL(
10.0,9.0))*t618+(RATIONAL(50.0,9.0)*t618+RATIONAL(20.0,9.0)*y+RATIONAL(-40.0,
9.0)*t619)*x;
coeff_dx_p2_p1 = t548+t573+(RATIONAL(-5.0,36.0)+RATIONAL(35.0,12.0)*t620)
*t618+(RATIONAL(4.0,3.0)*t619+t615*y)*x+t623;
coeff_dx_p3_p1 = t551+t586*t641+(t607*y+t617*t619)*x+t643;
coeff_dx_m2_p2 = t537+t559+(RATIONAL(-1.0,8.0)*t620+t599)*t618+(t592*y+
t585*t619)*x+t624;
coeff_dx_m1_p2 = t565+t531+(RATIONAL(7.0,8.0)*t620+t610)*t618+(RATIONAL(
35.0,24.0)*t619+t590*y)*x+t622;
coeff_dx_0_p2 = t555+RATIONAL(-2.0,1.0)*t641+(RATIONAL(-49.0,18.0)*t619+
t587*y)*x+t630;
coeff_dx_p1_p2 = t554+t566+(RATIONAL(2.0,1.0)*t620+t614)*t618+(t603*y+
RATIONAL(35.0,18.0)*t619)*x+t623;
coeff_dx_p2_p2 = t563+t521+t536+t560+(RATIONAL(-7.0,8.0)*t620+t595)*t618+
(t596*t619+t589*y+RATIONAL(1.0,2.0)*t618)*x;
coeff_dx_p3_p2 = t528+RATIONAL(1.0,8.0)*t641+(t594*y+RATIONAL(7.0,72.0)*
t619)*x+t645;
coeff_dx_m2_p3 = t592*t647+t561+(t644+t577)*t618+t633;
coeff_dx_m1_p3 = t590*t647+t567+(t639+t591)*t618+t629;
coeff_dx_0_p3 = x*t566+t618*t640+t632;
coeff_dx_p1_p3 = t603*t647+t569+(t604+t631)*t618+t643;
coeff_dx_p2_p3 = t564+t589*t647+(t634+t578)*t618+t645;
coeff_dx_p3_p3 = t529+t618*t636+(t594*t619+t601*t618)*x;
|