This study presents an energy preserving partially implicit scheme for the simulation of wave propagation in homogeneous and heterogeneous mediums. Despite its implicit nature, the developed scheme does not require any explicit numerical or analytical inversion of the coefficient matrix. Theoretical analysis and numerical experiments are performed to validate the energy preserving properties of the fully-discrete scheme. Convergence analysis is also performed to assess the rate of convergence of the developed scheme. The efficiency and accuracy of the developed scheme are validated by numerical solutions of wave propagation in layered heterogeneous mediums. Furthermore, simulations of soliton propagation following nonlinear sine-Gordon and Klein-Gordon equations in homogeneous and heterogeneous mediums are discussed. Numerical solutions are also compared with the results available in the literature. The present method accurately resolves the physical characteristics of the chosen problems, competing well with existing multi-stage time-integration methods. Moreover, it has significantly lower computational complexity than the four-stage, fourth-order Runge-Kutta-Nyström method.