N体シミュレーション
その後も1979年には Efstathiou & Jones が N = 500 {\displaystyle N=500} 体による銀河回転の研究など、計算機の発達に伴ってより大規模なシミュレーションがなされるようになっていった。
より大規模なシミュレーションの要求は強く、1986年に Barnes & Hut はツリー法を導入し、同時期に PM 法も確立した。1989年にはGRAPEプロジェクトがスタートしている。
一方で積分スキームに関する研究も進められた。天体力学の分野からは1990年に吉田春夫によりシンプレクティック積分子の一般的な構成方法が示された。その翌年牧野はエルミート積分子を導入した。やがて対称型公式の有用性が認められるようになった。
2005年のミレニアム・シミュレーション(英語版)では N = 1.0 × 10 10 {\displaystyle N=1.0\times 10^{10}} の宇宙論的構造形成シミュレーションが遂行されるに至った。