wait, I think I might be able to solve my own suggestion.
I'm working on a cubic interpolation function.
function cubicInterpolate1D(%pa,%pb,%pc,%pd,%x)
{
%ax3 = (((-1/2) * %pa) + ((3/2) * %pb) + ((-3/2) * %pc) + ((1/2) * %pd)) * mPow(%x, 3) //generate ax3
%bx2 = (%pa + ((-5/2) * %pb) + (2 * %pc) + ((-1/2) * %pd)) * mPow(%x, 2); //generate bx2
%cx = ((-1/2) * %pa + (1/2) * %pc) * %x) //generate cx
return %ax3 + %bx2 + %cx + %pa
}idk if this generates correct results.