This paper presents a new method for finding the Block Pulse series coefficients, deriving the Block Pulse integration operational matrices and generalizing the integration operational matrices which are necessary for the control fields using the Block Pulse functions. In order to apply the Block Pulse function technique to the problems of state estimation or parameter identification more efficiently, it is necessary to find the more exact value of the Block Pulse series coefficients and integral operational matrices. This paper presents the method for improving the accuracy of the Block Pulse series coefficients and derives the related integration operational matrices and generalized integration operational matrix by using the Lagrange second order interpolation polynomial.