We report the results of numerical experiments performed to study the
orbital circularization time of binary galaxies. We find that the time
scale is quite long (larger than the Hubble time), confirming earlier
calculations. The results depend on the initial conditions. From our
simulations we obtained a fitting formula for the circularization time
as a function of the initial orbital parameters like the pericentric
distance, mass ratio, and eccentricity.