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
