Pulp the LP solver

hello_pulp.py

import pulp

problem = pulp.LpProblem('hamburg and omelet', pulp.LpMaximize)

# create variables
x = pulp.LpVariable('x', cat = 'Integer')
y = pulp.LpVariable('y', cat = 'Integer')

# maximize
problem += 400 * x + 300 * y

# subjects
problem += 60 * x + 40 * y <= 3800
problem += 20 * x + 30 * y <= 2100
problem += 20 * x + 10 * y <= 1200

print problem

status = problem.solve()
print pulp.LpStatus[status]
print "Result"
print "x", x.value()
print "y", y.value()

max_value = 400 * x.value() + 300 * y.value()

print "MAX VALUE", "400 * %d + 300 * %d = %d" % (x.value(), y.value(), max_value)