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
|
fp t58,
t56,
t64,
t54,
t55,
t63,
t60,
t57,
t50;
t58 = x*x;
t56 = t58*t58;
t64 = RATIONAL(2.0,3.0)*t58+RATIONAL(-1.0,6.0)*t56;
t54 = RATIONAL(1.0,24.0);
t55 = RATIONAL(-1.0,24.0);
t63 = t54*t56+t55*t58;
t60 = t58*x;
t57 = t58*t60;
t50 = t55*t60;
coeff_I_m2 = RATIONAL(1.0,20.0)*x+t50+RATIONAL(-1.0,120.0)*t57+t63;
coeff_I_m1 = RATIONAL(-1.0,2.0)*x+t50+t54*t57+t64;
coeff_I_0 = RATIONAL(1.0,1.0)+RATIONAL(-1.0,3.0)*x+RATIONAL(-5.0,4.0)*t58
+RATIONAL(5.0,12.0)*t60+RATIONAL(1.0,4.0)*t56+RATIONAL(-1.0,12.0)*t57;
coeff_I_p1 = x+RATIONAL(-7.0,12.0)*t60+RATIONAL(1.0,12.0)*t57+t64;
coeff_I_p2 = RATIONAL(-1.0,4.0)*x+RATIONAL(7.0,24.0)*t60+t55*t57+t63;
coeff_I_p3 = RATIONAL(1.0,30.0)*x+RATIONAL(1.0,120.0)*t57+t50;
|