A bicubic model for local smoothing of surfaces is constructed on the base of pivot points. Such an approach allows reducing the dimension of matrix of normal equations more than twice. The model enables to increase essentially the speed and stability of calculations. The algorithms, constructed by the aid of the offered model, can be used both in applications and the development of global methods for smoothing and approximation of surfaces.