def mean_obliquity_ecliptic(julian_century):
    """Returns the Mean Obliquity of Ecliptic in Degrees with Julian Century, 
    julian_century."""
    mean_obliquity_of_ecliptic = (
        23
        + (
            26
            + (
                (
                    21.448
                    - julian_century
                    * (46.815 + julian_century * (0.00059 - julian_century * 0.001813))
                )
            )
            / 60
        )
        / 60
    )
    return mean_obliquity_of_ecliptic