for (i = 0; i <= 3; i++)
{
t2 = t2 * t;
y = y + p* * t2;
}
if (quad == 2 || quad == 1)
y = -y; // correct sign
for (i = 0; i <= 3; i++)
{
t2 = t2 * t;
y = y + p* * t2;
}
if (quad == 2 || quad == 1)
y = -y; // correct sign