import math


def solar_equation_of_center(julian_century, solar_geometric_mean_anomaly):
    """Returns the Solar Equation of Center with Julian Century, julian_century and 
    Solar Geometric Mean Anomaly, solar_geometric_mean_anomaly."""
    solar_equation_of_center = (
        math.sin(math.radians(solar_geometric_mean_anomaly))
        * (1.914602 - julian_century * (0.004817 + 0.000014 * julian_century))
        + math.sin(math.radians(2 * solar_geometric_mean_anomaly))
        * (0.019993 - 0.000101 * julian_century)
        + math.sin(math.radians(3 * solar_geometric_mean_anomaly)) * 0.000289
    )
    return solar_equation_of_center