def machine_shop():
    p = MixedIntegerLinearProgram()
    x = p.new_variable(real=True, nonnegative=True)

    p.set_objective(100 * x[1] + 150 * x[2])

    p.add_constraint(8000 * x[1] + 4000 * x[2] <= 40000)
    p.add_constraint(15 * x[1] + 30 * x[2] <= 200)

    p.show()
    
    print(round(p.solve(), 2))