Polynomial Surface Fitting