In this talk we formulate and prove results for the exponential matrix representing the dynamics of the Fermi-Bose model in an undepleted bosonic field approximation. A recent application of this model is molecular dimmers dissociating into its atomic compounds. The problem is solved in D spatial dimensions by dividing the system matrix into blocks with generalizations of Hankel matrices, here refered to as D-block Hankel-matrices. The method is practically useful for treating large systems, i.e. dense computational grids or higher spatial dimensions, either on a single standard computer or a cluster. In particular the results can be used for studies of threedimensional physical systems of arbitrary geometry. We illustrate the generality of our approach by giving numerical results for the dynamics of Glauber type atomic pair correlation functions for a non-isotropic three-dimensional harmonically trapped molecular Bose-Einstein condensate.