In this paper,we investigate the numerical differentiation of higher order.Based on orthogonal polynomial theory,integral operator method is utilized.Using the proposed method we can estimate any order derivatives of approximately specified functions.The new method also generalized the result of the Lanczos' s.