UPDATED: Performs Multivariate Polynomial Regression on multidimensional data. May not be able to handle ill-conditioned data, so make sure your columns are independent. The fits are limited to standard polynomial bases with minor modification options. Feel free to implement a term reduction heuristic.
The functionality is explained in hopefully sufficient detail within the m.file. Feel free to post a comment or inquiry.
No longer requires the statistical toolbox. Optional Functionality requires Symbolic Toolbox.