Final Ver. (I hope)
Lagrangian relax. It is slow...
Original data
runs slow