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