def Mass_function1(asini, porb):
    """
    Returns the mass function of a binary star.

    asini: Projected semi-major axis in light-second.
    porb: Orbital period in second.
    """
    mfunc = 8015123.37129 * asini**3 / porb**2
    return mfunc