We propose a variant of the multiconfiguration time-dependent Hartree (MCTDH) method within the framework of Hermite-distributed approximating functional (HDAF) method. The discretized Hamiltonian is a highly banded Toeplitz matrix which significantly reduces computational cost in terms of both storage and number of operations. The method proposed is employed to carry out the study of tunnelling dynamics of two coupled double well oscillators. We have calculated the orthogonality time \tau , which is a measure of the time interval for an initial state to evolve into its orthogonal state. It is observed that the coupling has a significant effect on \tau .