‎‎در مطالعه حاضر، با استفاده از روش عددی در محیط ‎$ mathematica $‎، مسئله چند جسمی در مکانیک سماوی محاسبه شده است. سپس کاربردهای این حل را در پیش بینی برخی پدیده های نجومی امروزین (ماه و خورشید گرفتگی) و نیز مدل سازی فرایندهای کهن در سامانه های سیاره ای ( نظیر تحول مدار سیارات و سیارکها، پایداری حرکت سیارکهای تروژان ) آزموده ایم. بدین منظور نخست به حل مسئله دو جسمی و سه جسمی در دستگاه مرکز جرم در دو بعد پرداخته و در گامهای بعدی مسئله را در سه بعد و برای تعداد ذرات بیشتر حل نموده ایم. با تغییر پارامترهای مرتبط با نسبت جرم ذرات و فواصل و سرعتهای اولیه آنها در معادلات بی بعد شده حرکت تحت نیروی گرانش نیوتنی‎‎‎ و حل عددی به پیش بینی ماه گرفتگی و خورشید گرفتگی و در ادامه به موضوعاتی نظیر سرنوشت سیارکها در سامانه خورشیدی، پاکسازی سامانه خورشیدی از خرده سیاره ها و پدیده نجومی قمرهای هم مدار پرداخته ایم.