The numerical integration of partial difference equations (PDEs) resulting from passive physical problems can be performed by simulating the actual system with multidimensional (MD) passive wave digital filters. The new result of this paper is a synthesis procedure for automatic generation of the algorithms, based on multidimensional wave digital filters for the numerical integration of partial differential equations describing the linear hyperbolic passive systems.