I’d like to know how it’s possible to obtain the lagrange multipliers from an optimal solution in a Concrete model solved with glpk? Thanks! Answer This works for me in gurobi, try it out and tell me if it works with glpk. You need to prompt getting the Lagrange multipliers / dual variable in pyomo by putting the following line