Numerical Integration of the Equations of Motion

Hi, I am wondering how to write a code to represent the discrete equation of motion. The input consists of stiffness and force matrix and the output will be a matrix of displacement. Any suggestion regarding this matter, Your cooperation is highly appreciated