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
|
fp t878,
t835,
t877,
t851,
t876,
t842,
t875,
t852,
t874,
t837,
t873,
t848,
t872,
t830,
t871,
t831,
t870,
t869,
t841,
t794,
t868,
t799,
t867,
t822,
t866,
t789,
t840,
t865,
t864,
t824,
t863,
t805,
t796,
t862,
t861,
t783,
t860,
t819,
t791,
t859,
t838,
t858,
t823,
t857,
t834,
t856,
t846,
t855,
t854,
t850,
t849,
t847,
t845,
t844,
t843,
t839,
t836,
t833,
t832,
t829,
t828,
t827,
t826,
t825,
t820,
t803,
t802,
t801,
t800,
t798,
t797,
t795,
t793,
t790,
t788,
t787,
t786,
t785,
t784,
t781;
t878 = x*y;
t835 = RATIONAL(-1.0,24.0);
t877 = t835*x;
t851 = y*y;
t876 = x*t851;
t842 = RATIONAL(2.0,3.0);
t875 = t842*y;
t852 = t851*y;
t874 = t852+y;
t837 = RATIONAL(7.0,24.0);
t873 = x*t837;
t848 = RATIONAL(-2.0,3.0);
t872 = x*t848;
t830 = RATIONAL(1.0,24.0);
t871 = x*t830;
t831 = RATIONAL(-7.0,24.0);
t870 = x*t831;
t869 = x*t852;
t841 = RATIONAL(1.0,6.0);
t794 = t837*t869;
t868 = t794+t841*t851;
t799 = t842*t869;
t867 = t799+RATIONAL(7.0,9.0)*t851;
t822 = RATIONAL(5.0,24.0);
t866 = t794+t822*t851;
t789 = RATIONAL(14.0,3.0)*t869;
t840 = RATIONAL(7.0,3.0);
t865 = t789+t840*t851;
t864 = RATIONAL(-49.0,24.0)*t869+RATIONAL(-5.0,4.0)*t851;
t824 = RATIONAL(-1.0,36.0);
t863 = t835*t869+t824*t851;
t805 = RATIONAL(-35.0,6.0)*t851;
t796 = t848*t869;
t862 = t796+RATIONAL(-5.0,9.0)*t851;
t861 = RATIONAL(-32.0,3.0)*t869+RATIONAL(-70.0,9.0)*t851;
t783 = t831*t869;
t860 = t783+RATIONAL(-5.0,12.0)*t851;
t819 = RATIONAL(-14.0,3.0);
t791 = t819*t869;
t859 = t791+t805;
t838 = RATIONAL(-7.0,18.0);
t858 = t796+t838*t851;
t823 = RATIONAL(5.0,18.0);
t857 = t799+t823*t851;
t834 = RATIONAL(-1.0,12.0);
t856 = t783+t834*t851;
t846 = RATIONAL(-5.0,3.0);
t855 = t791+t846*t851;
t854 = t789+RATIONAL(25.0,6.0)*t851;
t850 = RATIONAL(-2.0,9.0);
t849 = RATIONAL(2.0,9.0);
t847 = RATIONAL(-1.0,3.0);
t845 = RATIONAL(1.0,2.0);
t844 = RATIONAL(4.0,3.0);
t843 = RATIONAL(-4.0,3.0);
t839 = RATIONAL(-5.0,18.0);
t836 = RATIONAL(1.0,72.0);
t833 = RATIONAL(7.0,12.0);
t832 = RATIONAL(-7.0,12.0);
t829 = RATIONAL(-1.0,72.0);
t828 = RATIONAL(1.0,36.0);
t827 = RATIONAL(-5.0,24.0);
t826 = RATIONAL(1.0,12.0);
t825 = RATIONAL(7.0,18.0);
t820 = RATIONAL(10.0,3.0);
t803 = t847*t878;
t802 = t840*t878;
t801 = RATIONAL(-7.0,3.0)*t878;
t800 = RATIONAL(1.0,3.0)*t878;
t798 = x*t875;
t797 = y*t872;
t795 = y*t873;
t793 = y*t877;
t790 = RATIONAL(-16.0,3.0)*t878;
t788 = RATIONAL(16.0,3.0)*t878;
t787 = RATIONAL(32.0,3.0)*t869;
t786 = t830*t869;
t785 = y*t870;
t784 = y*t871;
t781 = RATIONAL(49.0,24.0)*t869;
coeff_dxx_m2_m2 = t786+t784+t874*t824+(RATIONAL(1.0,18.0)+t834*x)*t851;
coeff_dxx_m1_m2 = t833*t876+t785+t874*t822+t860;
coeff_dxx_0_m2 = t798+t843*t876+t874*t838+t867;
coeff_dxx_p1_m2 = t844*t876+t797+t874*t823+t862;
coeff_dxx_p2_m2 = t832*t876+t795+t874*t834+t868;
coeff_dxx_p3_m2 = t826*t876+t793+t874*t836+t863;
coeff_dxx_m2_m1 = t803+RATIONAL(5.0,8.0)*t876+t849*y+RATIONAL(7.0,36.0)*
t852+t860;
coeff_dxx_m1_m1 = t802+t846*y+RATIONAL(-35.0,24.0)*t852+t781+(RATIONAL(
25.0,8.0)+RATIONAL(-35.0,8.0)*x)*t851;
coeff_dxx_0_m1 = RATIONAL(49.0,18.0)*t852+t790+RATIONAL(28.0,9.0)*y+
RATIONAL(10.0,1.0)*t876+t859;
coeff_dxx_p1_m1 = RATIONAL(-35.0,18.0)*t852+RATIONAL(-20.0,9.0)*y+
RATIONAL(-10.0,1.0)*t876+t788+t854;
coeff_dxx_p2_m1 = t801+t833*t852+t875+RATIONAL(35.0,8.0)*t876+t864;
coeff_dxx_p3_m1 = t800+RATIONAL(-7.0,72.0)*t852+RATIONAL(-1.0,9.0)*y+
RATIONAL(-5.0,8.0)*t876+t866;
coeff_dxx_m2_0 = t847+RATIONAL(-4.0,9.0)*t852+(t845+RATIONAL(-7.0,6.0)*
t851)*x+t867;
coeff_dxx_m1_0 = t820*t852+RATIONAL(5.0,2.0)+(RATIONAL(49.0,6.0)*t851+
RATIONAL(-7.0,2.0))*x+t859;
coeff_dxx_0_0 = t819+t787+RATIONAL(98.0,9.0)*t851+RATIONAL(-56.0,9.0)*
t852+(RATIONAL(8.0,1.0)+RATIONAL(-56.0,3.0)*t851)*x;
coeff_dxx_p1_0 = t820+RATIONAL(40.0,9.0)*t852+(RATIONAL(-8.0,1.0)+
RATIONAL(56.0,3.0)*t851)*x+t861;
coeff_dxx_p2_0 = RATIONAL(-1.0,1.0)+t843*t852+(RATIONAL(-49.0,6.0)*t851+
RATIONAL(7.0,2.0))*x+t865;
coeff_dxx_p3_0 = t841+t849*t852+(RATIONAL(-1.0,2.0)+RATIONAL(7.0,6.0)*
t851)*x+t858;
coeff_dxx_m2_p1 = t800+RATIONAL(4.0,9.0)*t852+RATIONAL(5.0,6.0)*t876+t850
*y+t862;
coeff_dxx_m1_p1 = t801+RATIONAL(-10.0,3.0)*t852+RATIONAL(5.0,3.0)*y+x*
t805+t854;
coeff_dxx_0_p1 = RATIONAL(-28.0,9.0)*y+RATIONAL(40.0,3.0)*t876+t788+
RATIONAL(56.0,9.0)*t852+t861;
coeff_dxx_p1_p1 = RATIONAL(20.0,9.0)*y+t790+RATIONAL(-40.0,9.0)*t852+t787
+(RATIONAL(-40.0,3.0)*x+RATIONAL(50.0,9.0))*t851;
coeff_dxx_p2_p1 = RATIONAL(35.0,6.0)*t876+t802+t844*t852+t848*y+t855;
coeff_dxx_p3_p1 = t803+RATIONAL(-5.0,6.0)*t876+RATIONAL(1.0,9.0)*y+t850*
t852+t857;
coeff_dxx_m2_p2 = t793+RATIONAL(-1.0,4.0)*t876+RATIONAL(-7.0,36.0)*t852+
t828*y+t868;
coeff_dxx_m1_p2 = t827*y+t795+RATIONAL(7.0,4.0)*t876+RATIONAL(35.0,24.0)*
t852+t864;
coeff_dxx_0_p2 = RATIONAL(-4.0,1.0)*t876+t825*y+RATIONAL(-49.0,18.0)*t852
+t797+t865;
coeff_dxx_p1_p2 = t839*y+t798+RATIONAL(35.0,18.0)*t852+RATIONAL(4.0,1.0)*
t876+t855;
coeff_dxx_p2_p2 = t832*t852+t826*y+t785+t781+(RATIONAL(-7.0,4.0)*x+t845)*
t851;
coeff_dxx_p3_p2 = t784+RATIONAL(1.0,4.0)*t876+RATIONAL(7.0,72.0)*t852+
t829*y+t856;
coeff_dxx_m2_p3 = t851*t871+t828*t852+t863;
coeff_dxx_m1_p3 = t851*t870+t827*t852+t866;
coeff_dxx_0_p3 = t825*t852+t842*t876+t858;
coeff_dxx_p1_p3 = t839*t852+t851*t872+t857;
coeff_dxx_p2_p3 = t826*t852+t851*t873+t856;
coeff_dxx_p3_p3 = t829*t852+t786+(t877+t836)*t851;
|
