Iteration 1
augmented
  1 -2 1
  3 -1 3

 solution =  1 -0
 p    =  2 1
 f(p) =  0 4

Iteration 2
augmented
  1 -2 -0
  12 -1 -4

 solution =  -0.347826 -0.173913
 p    =  1.65217 0.826087
 f(p) =  -0.0302457 0.683817

Iteration 3
augmented
  1 -1.65217 0.0302457
  8.18904 -1 -0.683817

 solution =  -0.0925824 -0.0743433
 p    =  1.55959 0.751744
 f(p) =  -0.00552693 0.0416912

Iteration 4
augmented
  1 -1.50349 0.00552693
  7.29698 -1 -0.0416912

 solution =  -0.00684081 -0.00822604
 p    =  1.55275 0.743518
 f(p) =  -6.76677e-05 0.000218631

Iteration 5
augmented
  1 -1.48704 6.76677e-05
  7.2331 -1 -0.000218631

 solution =  -4.02608e-05 -7.25797e-05
 p    =  1.55271 0.743445
 f(p) =  -5.26781e-09 7.55066e-09

Iteration 6
augmented
  1 -1.48689 5.26781e-09
  7.23273 -1 -7.55066e-09

 solution =  -1.69103e-09 -4.68013e-09
 p    =  1.55271 0.743445
 f(p) =  -2.02041e-16 -8.23126e-16