Operational matrices of integration are used to convert the given partial differential equations with initial and boundary conditions into system of algebraic equations.Then we can easily find the solution by solving simultaneous algebraic equations.