def get_mu(X, Y, Z):
    """ Find the mean molecular weight

    Positional Arguments:
        X -> Hydrogen Abundence
        Y -> Helium Abundence
        Z -> Metallicity
    Returns:
        The mean molecular weight
    """

    return (1/(2*X+(3/4)*Y+(1/2)*Z))