I'll give you a couple of answers. I gave a huge hint in the terminology that I used for the actual method I used.
For simplicity's sake, let's assume a quadratic equation, for 3 values.
5,8,14
Think of it as 3 ordered pairs, (1,5),(2,8),(3,14)
Since they are non-linear, you can fit a parabola to the three points.
y=ax^2+bx+c
plug in 1 for x, 5 for y,
then 2 for x, 8 for y
3 for x 14 for y
and you get:
a+b+c=5
4a+2b+c=8
9a+3b+c=14
3 equations, 3 unknowns, solved very simply and quickly using linear algebra and matrices.
Alternatively, you can use the x-values of 0,1,and 2, which simply shifts the parabola over 1 unit. Thus, when x=0, c=y. At that point, you can do mental math to solve for a and b. Of course, doing it this way, and shifting it back over 1 unit, would require just as much work as using 1,2,3 for x's.
Thus, with the original 4 values, you end up with 4 equations and 4 unknowns.
If you assume whatever value you want for the 4th term, then you have 5 equations and 5 unknowns - pretty trivial stuff to solve using linear algebra, especially if you have a calculator available that handles matrices.
But, if you grab, for example, a TI83+ that handles matrices fairly well, you may as well simply enter the coordinate pairs in the lists and have the calculator do a cubic or quartic regression, since the best fit curve will, of course, pass exactly through those points (provided there is one more value than the degree of the polynomial.)