from pulp import *

## max 3 x + 5 y
## s.t.  x         <= 4
##              2y <= 12
##     3 x + 2 y <= 18
##       x         >= 0
##               y >= 0

prob = LpProblem("WynCor", LpMaximize)
# Decision variables
x = LpVariable("x", 0)
y = LpVariable("y", 0)
# Constraints
prob += x <=4
prob += 2*y <=12
prob += 3*x + 2*y <=18
# Objective
prob += 3*x + 5*y
# Solve
prob.solve()
LpStatus[prob.status]
# Results
x.varValue
y.varValue
value(prob.objective)