The transfer matrix method based is extensively used and well-validated for predicting the transmission loss of multilayer structures. However, this method leads to poor results at low frequencies due to the infinite extent assumption it is based on. This paper presents an efficient implementation of a Rayleigh-integral based method to account for the finite size effects. The accuracy of the method is illustrated by various examples.