version()
'SageMath version 9.2, Release Date: 2020-10-24'
We get some data from the kerrgeodesic_gw
SageMath package:
from kerrgeodesic_gw.astro_data import *
1 radian in microarcseconds:
rad_muas = n(180/pi * 3600 * 1e6)
rad_muas
2.06264806247096e11
1/rad_muas
4.84813681109536e-12
SgrA_mass_sol
4100000.0
SgrA_distance_pc
8120.0
s = SgrA_mass_m / SgrA_distance_m
s
2.416366617764889e-11
s*rad_muas
4.98411392235226
M87_mass_sol
6200000000.0
M87_distance_pc
16700000.0
s = M87_mass_m / M87_distance_m
s
1.776684109898577e-11
s*rad_muas
3.66467403690525
CygX1_mass_sol = 15
CygX1_mass_m = CygX1_mass_sol * solar_mass_m
CygX1_mass_sol
15
CygX1_distance_pc = 1.86e3
CygX1_distance_m = CygX1_distance_pc * pc
CygX1_distance_pc
1860.00000000000
s = CygX1_mass_m / CygX1_distance_m
s
3.85934243435305e-16
s*rad_muas
0.0000796046519463029
s*rad_muas*1e5
7.96046519463029
M31_mass_sol = 1.5e8
M31_mass_m = M31_mass_sol * solar_mass_m
M31_mass_sol
1.50000000000000e8
M31_distance_pc = 7.6e5
M31_distance_m = M31_distance_pc * pc
M31_distance_pc
760000.000000000
s = M31_mass_m / M31_distance_m
s
9.44523279986404e-12
s*rad_muas
1.94821911342268