Dimensional interpolation

Suppose you have a list of points symmetrically distributed around x=0 at integers...

For this method you need (n-1)/2 variables so in this case for 5 points we need 2, let's call them x and y....
Consider the following polynomial in 2 variables with 9 a(n) coefficients:

I generate the following table based on the points to be interpolated, x and y can each have values [-1,0,1], and the sum x+y is related to the range of the function. So x=-1 y=-1 x+y = -2, from the list of points the value of the function at -2 is 3 so A is three there. 

So I can solve for the 9 unknowns a1..a9 with 9 equations made from the table above:

So we can vary x between -1..1 and y between -1..1 with every possible pairwise combination to generate a plot:

Note that it is not single valued, because there are a lot of ways to add to 0 with x equaling a number between -1 and 1 and y also a number from -1 to 1. But it still does nicely interpolate the points other than that.
I haven't tried but I imagine wider ranges with more variables would work as well. I also suppose you could narrow the thickness of the band by ranging one variable over a smaller interval like y from -.5...5 and x from -1.5..1.5