This paper shows a new approach for non-linear least squares fitting with NURBS as curves and surfaces to measured data by the Gauss-Newton method. A Trust Region algorithm is used to reach global convergence as well as variable substitution and simple bounds.