def monomer_from_anisotropy(a, Am, Ad, b):
    """
    Calculate monomer fraction from anisotropy, monomer and dimer anisotropy
    and brightness relation.
    """
    return (b * a - Ad*b) / ((Am - b * Ad) - (1-b) * a)