Two boundary element method formulations are presented for the analysis of the one-dimensional scalar wave propagation problem in multi-region domains. One of the formulations employs the time-domain fundamental solution; the other, the fundamental solution related to the static problem. The problem domain is constituted of sub-domains with different material properties and, consequently, with different wave propagation velocities. Thus, the sub-region technique, akin to the boundary element method, is used to deal with the non-homogeneous problem. At the end of the article, some examples are presented, illustrating the characteristics of each formulation and their accuracy.