Se presenta un modelo no lineal con precisión robusta para el análisis del margen de estabilidad transitoria en sistemas eléctricos de potencia multimáquinas. Para evaluar las condiciones críticas de estabilidad dinámica, la técnica propuesta combina los métodos de interpolación de Lagrange con el de diferencias divididas de Newton, el resultado de la combinación es un método de interpolación de alto orden. El modelo usado para aplicar la técnica corresponde al sistema interconectado de 5-áreas de la red mexicana con 45 generadores. Como caso particular, en el contenido del trabajo se muestra la aplicación del procedimiento de alto orden para identificar el modo de baja frecuencia producido por una contingencia crítica. Ejemplos numéricos ilustran el método y demuestran la habilidad para extraer el comportamiento modal.