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;