A new method for shape modification of NURBS surfaces is proposed in this paper. Explicit formulae for computing new control points are derived by using constrained opimization method. Examples are also given to compare the results of Piegl's method with those of the new method.