A new operational matrix based on Bernoulli polynomials
  In this research, the Bernoulli polynomials are introduced. The properties of these polynomials are employed to construct the operational matrices of integration together with the derivative and product. These properties are then utilized to transform the differential equation to a matrix equation which corresponds to a system of algebraic equations with unknown Bernoulli coefficients. This method can be used for many problems such as differential equations, integral equations and so on. Numerical examples show the method is computationally simple and also illustrate the efficiency and accuracy of the method. 
J.A. Rad, S. Kazem, M. Shaban, K. Parand
cs.NA math.NA