Almost correct

Bugs:
1. If use flat start, will converge to different solution
2. the results are a little different from the standard results