Revert previous change, and describe subproblems

Revert previous commit: it actually made things more complicated.
Change name I^{cog} to I^{inc} (not all cogeneration units are
meant to fulfil the incentive constraints).
Correct wrong delta lagrangian term.
Add description of lagrangian function organizing the variables
by d.
Add description of lagrangian subproblems.