Abstract
This paper focus on research of numerical methods for solving systems of differential equations that contain fractional order derivatives. At first, the definitions of fractional order derivatives were presented, based on these definitions, an approximation formula for calculating fractional order derivatives was built and programmed. Finally, the program that has been built are used to calculate vibration of car model in which viscoelastic material contain fractional order derivatives.