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
|
fp t23,
t24,
t22,
t20,
t19,
t18,
t17,
t16,
t15,
t14,
t13;
t23 = x*x;
t24 = x*t23;
t22 = t24*t23;
t20 = t23*t23;
t19 = RATIONAL(-1.0,24.0);
t18 = RATIONAL(1.0,12.0);
t17 = RATIONAL(-7.0,12.0);
t16 = RATIONAL(5.0,12.0);
t15 = RATIONAL(-1.0,12.0);
t14 = RATIONAL(13.0,12.0);
t13 = t19*t24;
coeff_I_m2 = t18*x+RATIONAL(-1.0,8.0)*t23+t13+RATIONAL(1.0,8.0)*t20+t19*
t22;
coeff_I_m1 = RATIONAL(-2.0,3.0)*x+t14*t23+t13+t17*t20+RATIONAL(5.0,24.0)*
t22;
coeff_I_0 = RATIONAL(-25.0,12.0)*t23+t16*t24+t14*t20+RATIONAL(1.0,1.0)+
RATIONAL(-5.0,12.0)*t22;
coeff_I_p1 = t17*t24+RATIONAL(2.0,3.0)*x+RATIONAL(3.0,2.0)*t23+t16*t22-
t20;
coeff_I_p2 = t15*x+RATIONAL(-11.0,24.0)*t23+RATIONAL(7.0,24.0)*t24+
RATIONAL(-5.0,24.0)*t22+RATIONAL(11.0,24.0)*t20;
coeff_I_p3 = t13+RATIONAL(1.0,24.0)*t22+t18*t23+t15*t20;
|
