The associated paper presents a novel computational homogenization method enabling efficient and accurate modelling of wave propagation phenomena in locally resonant acoustic metamaterial (LRAM) panels, which was not possible for direct numerical simulations as well as available homogenization methods. The proposed method provides a considerable advantage for fast analyses and design of metamaterial panels. The method above is implemented in MATLAB and the provided MATLAB scripts include:1) Calculation of the frequency-independent dynamic effective material properties for a given LRAM unit cell design. Finite element method is employed.2) Dispersion analysis of an infinite LRAM beam, using plane wave transformation.3) Frequency domain analysis of a finite LRAM beam, based on harmonic solution assumption. Isogeometric analysis is employed.4) Transient analysis of a finite LRAM beam. Isogeometric analysis is employed. Time integration is performed using Newmark method. Reference solutions by direct numerical simulations using COMSOL are given as well.