hdavid16

using JuMP, HiGHS, Polyhedra, CDDLib, Plots

# solve LP
m = Model(HiGHS.Optimizer)
x, y = @variable(m, v[i=1:2], lower_bound = 0)
@constraint(m, c[i=1:4], [6,1,-1,0][i]*x + [4,2,1,1][i]*y ≤ [24,6,1,2][i])
@objective(m, Max, x+y)
optimize!(m)
xopt, yopt = value.(v)