human
I wrote the below code

    beam = beam_area.to_value(si.sr)
    nu = disp.to(si.GHz, spectral())
    def convert_Jy_to_K(x_jybm):
        factor = ((((2 * _si.k_B) * si.K) * (nu ** 2)) / (_si.c ** 2)).to_value(astrophys.Jy)
        return ((x_jybm / beam) / factor)
    def convert_K_to_Jy(x_K):
        factor = (astrophys.Jy / (((2 * _si.k_B) * (nu ** 2)) / (_si.c ** 2))).to_value(si.K)
        return ((x_K * beam) / factor)
    return [(astrophys.Jy, si.K, convert_Jy_to_K, convert_K_to_Jy)]

What's a good function header?
 ... Here's a good function header?
assistant
#Here's a good function header.

def brightness_temperature(beam_area, disp):