Abstract
In this paper a matrix separation of the variables method for solving initial-boundary value problems for coupled systems of second order partial differential equations is presented. Finite computable approximate solutions and upper error bounds for them in terms of data are given.