import math


def approx_atmospheric_refraction(solar_elevation_angle):
    """Returns Approximate Atmospheric Refraction in degrees with Solar Elevation 
    Angle, solar_elevation_angle."""
    if solar_elevation_angle > 85:
        approx_atmospheric_refraction = 0
    elif solar_elevation_angle > 5:
        approx_atmospheric_refraction = (
            58.1 / math.tan(math.radians(solar_elevation_angle))
            - 0.07 / pow(math.tan(math.radians(solar_elevation_angle)), 3)
            + 0.000086 / pow(math.tan(math.radians(solar_elevation_angle)), 5)
        ) / 3600
    elif solar_elevation_angle > -0.575:
        approx_atmospheric_refraction = (
            1735
            + solar_elevation_angle
            * (
                -518.2
                + solar_elevation_angle
                * (
                    103.4
                    + solar_elevation_angle * (-12.79 + solar_elevation_angle * 0.711)
                )
            )
        ) / 3600
    else:
        approx_atmospheric_refraction = (
            -20.772 / math.tan(math.radians(solar_elevation_angle))
        ) / 3600
    return approx_atmospheric_refraction