Complete description of Lagrangian subproblems

Refactoring: change variable names.
Heat storage variables are now 's', while the Lagrangian
variable 'w^q' are now 'u'.

Complete description of Lagrangian subproblems,

Add description of Subgradient Algorithm for the Lagrangian Dual.