On an operational matrix method based on generalized Bernoulli polynomials of level m