Resumen: Runge–Kutta–Nyström (RKN) methods for the numerical solution of inhomogeneous linear initial value problems with constant coefficients are considered. A general procedure to construct explicit-stage RKN methods with maximal order [fórmula], similar to the developed by the authors (Montijano et al., 2023) for the class of second order IVP under consideration, depending on the nodes [fórmula] is presented. This procedure requires only the solution of successive linear equations in the elementij, [fórmula], of the matrix of coefficients A of the RKN method and avoids the solution of non linear equations. The remarkable fact is that using as free parameters the nodes [fórmula] with a quadrature relation, the [fórmula] elements of matrix can be computed by solving successively linear systems with coefficients depending on the nodes, so that if they are non-singular we get a unique 8-stage method with maximal order [fórmula]. We obtain an optimized six-stage seventh-order RKN method in the sense that the nodes are chosen so that minimize the leading term of the local error. Finally, some numerical experiments are presented to test the behaviour of the optimized RKN method with others with Radau and Lobatto nodes. Idioma: Inglés DOI: 10.1016/j.cam.2023.115533 Año: 2024 Publicado en: Journal of Computational and Applied Mathematics 438 (2024), 115533 [19 pp.] ISSN: 0377-0427 Financiación: info:eu-repo/grantAgreement/ES/MICINN/PID2019-109045GB-C31 Tipo y forma: Artículo (Versión definitiva) Área (Departamento): Área Matemática Aplicada (Dpto. Matemática Aplicada)