Fitting data in space by surfaces in parametric representation with polynomial components

We consider fitting measured data points in space in the total least squares sense by surfaces in parametric representation with polynomial components in two variables. A well-known descent algorithm is suitably modified. Numerical examples are given.