The coupled vibration functions for a half railway vehicle model running on a simply
sported beam bridge were derived. Based on Ruge-Kutta method, the coupled vibration functions
were solved with the ODE serial functions of MATLAB. Compared with Nemark-βmethod, the
proposed method has higher computing efficiency with the same accuracy.