a=[[0.00000000000001, -50.0, -3.0, -90.0], [-85.0, 2.0, -25.0, -6.0], [79.0, 5.0, 30.0, -1.0]] ?> gj(a) gj(a) => [0.0, 1.8, 0.0] >> gj_pivot(a) gj_pivot(a) => [0.99314840191104, 1.97870449242621, -2.97840820710344](元々はgj(a)がNANになる例だったので改作しました)
Err(pi1(n)) = 0.410545842036582/√n Err(pi2(n)) = 0.199265567320772/√nあたりになるはずである.したがって,「pi2の誤差が小さい」でも良いし, 数回繰り返した程度では有意な差が得られないはずなので,「それほど違い がない」でも良い.