In this paper a numerical method for the computation of the acoustic field surrounding a set of vibrating bodies and coupled, forced shells (or shields) is introduced. The method is derived through the reformulation of the Helmholtz equation, which governs the acoustic field, as an integral equation termed the boundary and shell integral equation. Collocation is applied to this equation and a method termed the boundary and shell element method (BSEM) results. Modal modelling of the shell is then included to complete the model of the full acoustic-structure system. The method is implemented for general three-dimensional problems in a Fortran subroutine ABSEMGEN. The subroutine is applied to test problems and results are demonstrated.