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