Properties of the hybrid of block-pulse functions and Lagrange
polynomials based on the Legendre-Gauss-type points are
investigated and utilized to define the composite interpolation
operator as an extension of the well-known Legendre interpolation
operator. The uniqueness and interpolating properties are
discussed and the corresponding differentiation matrix is also
introduced. The applicability and effectiveness of the method are
illustrated by two numerical experiments.