恩格尔曲线

def solve(____):
    a, b = utility_func
    engel = []
    for income in income_values:
        x_opt = a * income / ((____) * px)
        engel.append({"____": income, "____": round(____)})
    return {"____": engel}