The natural frequency and main vibration mode of multi-freedom vibration system are the base of the research on the vibration characteristics, matrix iterative algorithm is one of the basic methods to caculate these dynamic features. It programmes matrix iterative algorithm on MATLAB. This program can quickly get the natural frequency and the main vibration mode of multi-freedom vibration systems, provide theory foundation for designers to avoid system resonance.