test_lsfit2.adb nx= 15, step= 1.42857142857143E-01 ny= 15, step= 1.42857142857143E-01 a^ 0 b^ 0 a^ 1 b^ 0 a^ 0 b^ 1 C( 1)= 1.0000E+00 x^0 y^0 C( 2)= 2.0000E+00 x^1 y^0 C( 3)= 3.0000E+00 x^0 y^1 power 1 test, maxerr= 4.88498130835069E-15 a^ 0 b^ 0 a^ 1 b^ 0 a^ 0 b^ 1 a^ 2 b^ 0 a^ 1 b^ 1 a^ 0 b^ 2 C( 1)= 1.0000E+00 x^0 y^0 C( 2)= 2.0000E+00 x^1 y^0 C( 3)= 3.0000E+00 x^0 y^1 C( 4)= 4.0000E+00 x^2 y^0 C( 5)= 5.0000E+00 x^1 y^1 C( 6)= 6.0000E+00 x^0 y^2 power 2 test, maxerr= 1.42108547152020E-14 a^ 0 b^ 0 a^ 1 b^ 0 a^ 0 b^ 1 a^ 2 b^ 0 a^ 1 b^ 1 a^ 0 b^ 2 a^ 3 b^ 0 a^ 2 b^ 1 a^ 1 b^ 2 a^ 0 b^ 3 C( 1)= 1.0000E+00 x^0 y^0 C( 2)= 2.0000E+00 x^1 y^0 C( 3)= 3.0000E+00 x^0 y^1 C( 4)= 4.0000E+00 x^2 y^0 C( 5)= 5.0000E+00 x^1 y^1 C( 6)= 6.0000E+00 x^0 y^2 C( 7)= 7.0000E+00 x^3 y^0 C( 8)= 8.0000E+00 x^2 y^1 C( 9)= 9.0000E+00 x^1 y^2 C( 10)= 1.0000E+01 x^0 y^3 power 3 test, maxerr= 3.55271367880050E-14 a^ 0 b^ 0 a^ 1 b^ 0 a^ 0 b^ 1 a^ 2 b^ 0 a^ 1 b^ 1 a^ 0 b^ 2 a^ 3 b^ 0 a^ 2 b^ 1 a^ 1 b^ 2 a^ 0 b^ 3 a^ 4 b^ 0 a^ 3 b^ 1 a^ 2 b^ 2 a^ 1 b^ 3 a^ 0 b^ 4 C( 1)= 1.0000E+00 x^0 y^0 C( 2)= 2.0000E+00 x^1 y^0 C( 3)= 3.0000E+00 x^0 y^1 C( 4)= 4.0000E+00 x^2 y^0 C( 5)= 5.0000E+00 x^1 y^1 C( 6)= 6.0000E+00 x^0 y^2 C( 7)= 7.0000E+00 x^3 y^0 C( 8)= 8.0000E+00 x^2 y^1 C( 9)= 9.0000E+00 x^1 y^2 C( 10)= 1.0000E+01 x^0 y^3 C( 11)= 1.1000E+01 x^4 y^0 C( 12)= 1.2000E+01 x^3 y^1 C( 13)= 1.3000E+01 x^2 y^2 C( 14)= 1.4000E+01 x^1 y^3 C( 15)= 1.5000E+01 x^0 y^4 power 4 test, maxerr= 1.45661260830821E-13 a^ 0 b^ 0 a^ 1 b^ 0 a^ 0 b^ 1 a^ 2 b^ 0 a^ 1 b^ 1 a^ 0 b^ 2 a^ 3 b^ 0 a^ 2 b^ 1 a^ 1 b^ 2 a^ 0 b^ 3 a^ 4 b^ 0 a^ 3 b^ 1 a^ 2 b^ 2 a^ 1 b^ 3 a^ 0 b^ 4 a^ 5 b^ 0 a^ 4 b^ 1 a^ 3 b^ 2 a^ 2 b^ 3 a^ 1 b^ 4 a^ 0 b^ 5 C( 1)= 1.0000E+00 x^0 y^0 C( 2)= 2.0000E+00 x^1 y^0 C( 3)= 3.0000E+00 x^0 y^1 C( 4)= 4.0000E+00 x^2 y^0 C( 5)= 5.0000E+00 x^1 y^1 C( 6)= 6.0000E+00 x^0 y^2 C( 7)= 7.0000E+00 x^3 y^0 C( 8)= 8.0000E+00 x^2 y^1 C( 9)= 9.0000E+00 x^1 y^2 C( 10)= 1.0000E+01 x^0 y^3 C( 11)= 1.1000E+01 x^4 y^0 C( 12)= 1.2000E+01 x^3 y^1 C( 13)= 1.3000E+01 x^2 y^2 C( 14)= 1.4000E+01 x^1 y^3 C( 15)= 1.5000E+01 x^0 y^4 C( 16)= 1.6000E+01 x^5 y^0 C( 17)= 1.7000E+01 x^4 y^1 C( 18)= 1.8000E+01 x^3 y^2 C( 19)= 1.9000E+01 x^2 y^3 C( 20)= 2.0000E+01 x^1 y^4 C( 21)= 2.1000E+01 x^0 y^5 power 5 test, maxerr= 5.64881474929280E-13 a^ 0 b^ 0 a^ 1 b^ 0 a^ 0 b^ 1 a^ 2 b^ 0 a^ 1 b^ 1 a^ 0 b^ 2 a^ 3 b^ 0 a^ 2 b^ 1 a^ 1 b^ 2 a^ 0 b^ 3 a^ 4 b^ 0 a^ 3 b^ 1 a^ 2 b^ 2 a^ 1 b^ 3 a^ 0 b^ 4 a^ 5 b^ 0 a^ 4 b^ 1 a^ 3 b^ 2 a^ 2 b^ 3 a^ 1 b^ 4 a^ 0 b^ 5 a^ 6 b^ 0 a^ 5 b^ 1 a^ 4 b^ 2 a^ 3 b^ 3 a^ 2 b^ 4 a^ 1 b^ 5 a^ 0 b^ 6 C( 1)= 1.0000E+00 x^0 y^0 C( 2)= 2.0000E+00 x^1 y^0 C( 3)= 3.0000E+00 x^0 y^1 C( 4)= 4.0000E+00 x^2 y^0 C( 5)= 5.0000E+00 x^1 y^1 C( 6)= 6.0000E+00 x^0 y^2 C( 7)= 7.0000E+00 x^3 y^0 C( 8)= 8.0000E+00 x^2 y^1 C( 9)= 9.0000E+00 x^1 y^2 C( 10)= 1.0000E+01 x^0 y^3 C( 11)= 1.1000E+01 x^4 y^0 C( 12)= 1.2000E+01 x^3 y^1 C( 13)= 1.3000E+01 x^2 y^2 C( 14)= 1.4000E+01 x^1 y^3 C( 15)= 1.5000E+01 x^0 y^4 C( 16)= 1.6000E+01 x^5 y^0 C( 17)= 1.7000E+01 x^4 y^1 C( 18)= 1.8000E+01 x^3 y^2 C( 19)= 1.9000E+01 x^2 y^3 C( 20)= 2.0000E+01 x^1 y^4 C( 21)= 2.1000E+01 x^0 y^5 C( 22)= 2.2000E+01 x^6 y^0 C( 23)= 2.3000E+01 x^5 y^1 C( 24)= 2.4000E+01 x^4 y^2 C( 25)= 2.5000E+01 x^3 y^3 C( 26)= 2.6000E+01 x^2 y^4 C( 27)= 2.7000E+01 x^1 y^5 C( 28)= 2.8000E+01 x^0 y^6 power 6 test, maxerr= 2.33058017329313E-12 just boundary to data, power 4 a^ 0 b^ 0 a^ 1 b^ 0 a^ 0 b^ 1 a^ 2 b^ 0 a^ 1 b^ 1 a^ 0 b^ 2 a^ 3 b^ 0 a^ 2 b^ 1 a^ 1 b^ 2 a^ 0 b^ 3 a^ 4 b^ 0 a^ 3 b^ 1 a^ 2 b^ 2 a^ 1 b^ 3 a^ 0 b^ 4 C( 1)= -1.2000E+01 x^0 y^0 C( 2)= 2.0000E+00 x^1 y^0 C( 3)= 3.0000E+00 x^0 y^1 C( 4)= 1.7000E+01 x^2 y^0 C( 5)= 5.0000E+00 x^1 y^1 C( 6)= 1.9000E+01 x^0 y^2 C( 7)= 7.0000E+00 x^3 y^0 C( 8)= 8.0000E+00 x^2 y^1 C( 9)= 9.0000E+00 x^1 y^2 C( 10)= 1.0000E+01 x^0 y^3 C( 11)= 1.1000E+01 x^4 y^0 C( 12)= 1.2000E+01 x^3 y^1 C( 13)= 0.0000E+00 x^2 y^2 C( 14)= 1.4000E+01 x^1 y^3 C( 15)= 1.5000E+01 x^0 y^4 just boundary test4, maxerr= 4.40536496171262E-13 all internal test4, maxerr= 1.30000000000000E+01 check Bpoly power 4, dimension 2 a^ 0 b^ 0 a^ 1 b^ 0 a^ 0 b^ 1 a^ 2 b^ 0 a^ 1 b^ 1 a^ 0 b^ 2 a^ 3 b^ 0 a^ 2 b^ 1 a^ 1 b^ 2 a^ 0 b^ 3 a^ 4 b^ 0 a^ 3 b^ 1 a^ 2 b^ 2 a^ 1 b^ 3 a^ 0 b^ 4 C( 1)= -1.2000E+01 x^0 y^0 C( 2)= 2.0000E+00 x^1 y^0 C( 3)= 3.0000E+00 x^0 y^1 C( 4)= 1.7000E+01 x^2 y^0 C( 5)= 5.0000E+00 x^1 y^1 C( 6)= 1.9000E+01 x^0 y^2 C( 7)= 7.0000E+00 x^3 y^0 C( 8)= 8.0000E+00 x^2 y^1 C( 9)= 9.0000E+00 x^1 y^2 C( 10)= 1.0000E+01 x^0 y^3 C( 11)= 1.1000E+01 x^4 y^0 C( 12)= 1.2000E+01 x^3 y^1 C( 13)= 0.0000E+00 x^2 y^2 C( 14)= 1.4000E+01 x^1 y^3 C( 15)= 1.5000E+01 x^0 y^4 Bpoly fit, maxerr= 4.42109890363973E-13 internal Bpoly fit, maxerr= 1.30000000000000E+01 test_lsfit2.adb finished