import math


def yield_strength(impactor_density_kgpm3):
    """
    Yield strength equation for breakup altitude calculation. Only valid for density range 1000 to 8000.
    :param impactor_density_kgpm3: Impactor density in kg/m^3
    :returns: Yield Strength in Pascals.
    :Reference: EarthImpactEffect.pdf, Equation 10
    """

    return 10 ** (2.107 + (0.0624 * math.sqrt(impactor_density_kgpm3)))