f(i) = 1 / ((2 * i - 2) * (2 * i - 3) * a * a);

float x, rez_one, rez_two, buf;
pt >> x;
buf = x * ((-x) * (4 * x - 3) - 2);
rez_one = 1 + buf;
rez_two = 1 - buf;
pt << rez_one;
pt << rez_two;

y = x*x;       // 1
z = 3*y+1;     // 2,3
t = (4*y+2)*x; // 4,5,6
y = z - t;     // 7
z = z + t;     //8