Often it's necessary to fit or evaluate a function with 2 dimensions. Use polyFit2D.m to fit a 2D polynomila to data and polyVal2D.m to evaluate it. Fun!
These functions are explained in detail here:
They are also hosted in a Gist on Github here:
now with polyDer2D.m