Published Paper Details:

In recent years there has been an increasing demand for solutions of the general three-body problem in various astrophysical situations. For example, binaries and their interactions with single stars play a major role in the evolution of star clusters. Triple stellar systems are another obvious astrophysical three-body problem. Many other astrophysical bodies, ranging from compact bodies to galaxies, occur in triple systems. At present time, three body astrophysics is primarily motivated by the need to understand the role of binaries in the evolution of stellar systems because mathematically ideal solution for the three-body problem is of no use i.e., Sundman's series solution which has extremely slow convergence. In this paper we have implemented the numerical integration method in python for two body problem for understanding the implementation and then extended it to higher body problem (three and four body). Here we have used Odeint package of Scipy library in Python to solve the equation of motion for bodies under gravitation. By observations of the results of computer simulations, it is found that this problem got some patterns when observed over a very long interval of time. Particularly in problems of astrophysical importance, one may almost always identify a binary and a third body. A binary can be treated as a single entity with certain internal properties (component masses, energy and angular momentum, as well as orientation in space) and this system and third body makes themselves effective binary.