Is it possible to get the Lagrangian multiplier for the constraints from my solution?
The Lagrange multipliers, also known as marginals, can be accessed after a solve as follows.
Say the constraint of interest is called CON(j,k) and you want to set lamda(j,k) to the value of its marginal after the solve. Then anywhere after the solve statement, just write:
lamda(j,k) = CON.M(j,k);