In this paper the vertical vibration response of a 9 degrees-of-freedom (DOF) 4-axle freight wagon, modeled as a multi-body system is investigated. The Lagrange's equation is used to obtain the equations of motion, and the corresponding system matrices are then assembled. A combined modal analysis and numerical integration of the convolution integral is adopted, which is achieved by programming a MATLAB script. Modal analysis is used to uncouple the differential equations and thereafter, the convolution integrals are numerically integrated using piecewise constants as interpolating functions. The Complex-Exponential Method (CEM) is then used to extract damping ratio, mode shapes and the system's natural frequencies from test data, this information is then used to compute the numerical solution. The developed models are then validated with measured data.