Numerically solving differential equations of motion