def equilibrium_temperature(T_star, R_star, albedo, semi_major_axis):
    """
    Calculate the equilibrium temperature of a planet.
    """
    return T_star * (1 - albedo)**.25 * (R_star/(2*semi_major_axis))**.5