# Copyright (c) 2025 Jakub Szyndler # # This program is free software: you can redistribute it and/or modify # it under the terms of the GNU General Public License as published by # the Free Software Foundation, either version 3 of the License, or # (at your option) any later version. # # This program is distributed in the hope that it will be useful, # but WITHOUT ANY WARRANTY; without even the implied warranty of # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the # GNU General Public License for more details. # # You should have received a copy of the GNU General Public License # along with this program. If not, see #constants Year = 365.25 #days Day = 24*3600 AU = 1.495978707e11#m pc = 3.0857e16 #m solar_mass = 1.9891e30 #kg mass_conversion = 1.3466e27/solar_mass #kg time_conversion= 3.3356e-9#s G = 2.96e-4 c = 173 #AU/days #mercury mass 1.652e-7 m_sun #Input paraeters m_1 = 1.4*G/c**2 # masses in solar mass ( we work in geometrized units) m_2 = 1.4*G/c**2 M = m_1 + m_2 eta = m_1*m_2/M r_0 = 746600000/AU # initial separation( we operate in Astronomicals units) t0 = 35*Year*c # Time in which we want to evolve the configuration (in years*days*c) H = 1#7.751938773864*60/4/Year