A new approach with the ability to use the multiple observations based on the least square approach has been proposed for initial orbit determination. This approach considers the Earth’s Oblateness by using the developed Lagrange coefficients. The efficiency of the proposed method has been tested in two scenarios. The first scenario is to use the simulated and the second one is to utilize the r...