An automatic method for data smoothing and numerical differentiation, based on mollification techniques, trigonometric interpolation and the principle of generalized cross validation is presented. With data measured at a discrete set of equidistant points in a given interval, the method allows for the approximate recovery of the function and its first derivative on the entire domain.