The present study considers the solutions of hyperbolic partial differential equations. For this, an approximate method based on Bernoulli polynomials is developed. This method transforms the equation into the matrix equation and the unknown of this equation is a Bernoulli coefficients matrix. To demostrate the validity and applicability of the method, an error analysis developed based on residual function. Also examples are presented to illustrate the accuracy of the method.