A note on high order Bernoulli numbers and polynomials using differential equations

An ordinary differential equation is constructed to determine coefficients of a recurrence formula related to the generating function of Bernoulli numbers. This construction is more complicated work than the case of Eulerian numbers and polynomials. Solving this differential equation, we derive some identities on Bernoulli numbers and polynomials of higher order.