def orbital_params(m1, m2, p):
    Msolar = 1.989e+33 # g
    mt = m1 + m2
    mt = Msolar * mt
    G = 6.67430e-8 # cgs
    a = ((p**2 * G * mt)/(4*np.pi**2))**(1/3) # cm
    return a