Durante más de 5 décadas, los cálculos perturbativos en Física de Altas Energías se han realizado "a la Feynman". La enorme cantidad de datos que se obtienen en el Gran Colisionador de Hadrones (LHC, CERN), demanda predicciones teóricas cada vez más precisas.

En los años treinta los trabajos de los astrónomos Zwiky y Oort sugirieron por primera vez una posible discrepancia entre el potencial gravitacional observado y el esperado dada la masa en estrellas y gas.

Este trabajo presenta, en concreto, una breve descripción del formalismo post-Newtoniano, la descripción del método numérico utilizado, el análisis y validación de las soluciones numéricas obtenidas del código creado para simular sistemas de N cuerpos incluyendo todos los términos de las ecuaciones de Einstein-Infeld-Hoffmann.