import pyomo.environ as pyo
model = pyo.ConcreteModel(name='Bubbles Brewery')

# Decision variables
model.light_beer = pyo.Var(within=pyo.NonNegativeReals)
model.dark_beer = pyo.Var(within=pyo.NonNegativeReals)

# Objective function
model.objective = pyo.Objective(expr=13*model.light_beer + 23*model.dark_beer, sense=pyo.maximize)

# Constraints
model.corn = pyo.Constraint(expr=5*model.light_beer + 15*model.dark_beer <= 480)
model.malt = pyo.Constraint(expr=4*model.light_beer + 4*model.dark_beer <= 160)
model.hops = pyo.Constraint(expr=35*model.light_beer + 20*model.dark_beer <= 1190)

# Export the model
model.write('bubbles_brewery.mps', io_options={'symbolic_solver_labels': True})

# Optionally print the model to verify
model.pprint()