This paper presents a new method for finding the Block Pulse series coefficients and deriving the Block Pulse 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 continuous-time dynamic systems more efficiently, it is necessary to find the more exact value of the Block Pulse series coefficients and drives the related integration operational matrices by using the Lagrange second order interpolation polynomial.