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
|
fp t25,
t28,
t20,
t27,
t18,
t34,
t33,
t23,
t22,
t17,
t32,
t15,
t31,
t30,
t29,
t26,
t21,
t19,
t16,
t14;
t25 = RATIONAL(1.0,6.0);
t28 = y*y;
t20 = t25*t28;
t27 = x*x;
t18 = t25*t27;
t34 = t20+RATIONAL(-1.0,9.0)+t18;
t33 = x*y;
t23 = RATIONAL(2.0,9.0);
t22 = RATIONAL(-1.0,3.0);
t17 = t22*t27;
t32 = t20+t23+t17;
t15 = t22*t28;
t31 = t15+t23+t18;
t30 = RATIONAL(-1.0,4.0)*t33+t34;
t29 = RATIONAL(1.0,4.0)*t33+t34;
t26 = RATIONAL(-1.0,6.0);
t21 = t26*y;
t19 = t26*x;
t16 = t25*x;
t14 = t25*y;
coeff_I_m1_m1 = t21+t19+t29;
coeff_I_0_m1 = t21+t32;
coeff_I_p1_m1 = t21+t16+t30;
coeff_I_m1_0 = t19+t31;
coeff_I_0_0 = t15+t17+RATIONAL(5.0,9.0);
coeff_I_p1_0 = t16+t31;
coeff_I_m1_p1 = t14+t19+t30;
coeff_I_0_p1 = t14+t32;
coeff_I_p1_p1 = t14+t16+t29;
|