We consider the numerical treatment of systems of second order differential-algebraic equations (DAEs). The classical approach of transforming a second order system to first order by introducing new variables can lead to difficulties such as an increase in the index or the loss of structure. We show how we can compute an equivalent strangeness-free second order system using the derivative array approach and we present Runge-Kutta methods for the direct numerical solution of second order DAEs.We discuss the numerical solution of the linear second order differential-algebraic initial value problemwhere M, C, K ∈ C(I, C m×n ) and f (t) ∈ C(I, C m ) are sufficiently smooth together with given initial values x 0 ,ẋ 0 ∈ C n . Usually, in the classical theory of ordinary differential equations, high order systems are turned into first order systems by introducing new variables for the derivatives. However, it has been discussed in several publications, see e.g. [1,4], that for differential-algebraic equations (DAEs) this classical approach has to be performed with great care, since it may lead to a number of mathematical difficulties. In [3,5] several examples show that the classical approach of introducing new variables may lead to higher smoothness requirements for the inhomogeneity, corresponding to an increase in the index of the DAE, or may result in a loss of structures in the system. Further, it is shown in [1, 2, 4, 6] that numerical methods can fail. Therefore, the direct discretization and numerical solution of the second order differential-algebraic system can yield better numerical results. The theoretical analysis of linear second order differential-algebraic equations is studied in [3,5], where condensed forms for matrix triple and corresponding invariants under equivalence transformations are derived. One of the main results is the equivalence to a so-called strangeness-free system, that separates the given second order DAE into uncoupled subsystems of second order differential, first order differential and algebraic equations together with some consistency conditions for the right-hand side. The algebraic approach described in [3,5] gives the theoretical analysis of second order DAEs including decisions about existence and uniqueness of solutions and consistency of initial values, but it cannot be used for the development of numerical methods as it is not feasible for numerical computations. Here, we introduce a numerical computable way for the determination of the characteristic invariants of a given DAE, i.e., the number of second order and the number of first order differential equations, the number of algebraic equations and the number of consistency conditions in the strangeness-free system, as well as for the extraction of an equivalent strangeness-free second order differential-algebraic system. The basic idea is to differentiate the DAE (1) a number of times and to put the original DAE together with the derivatives up to a sufficiently high orderμ into an enlarged system. In this way we obtain t...