# Polynomial interpolation

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 ?