I wanted to do find the polynomial function in multivariable with set of data points. Something similar to the one done in mathematica that is the function InterpolatePolynomial. In one variable there is the function in sage called .lagrange_polynomial which works pretty good. I was wonder if there is something in multivariable ?