test_psimeq_dbg.py
x_data[ 0 ]= 0.0 , y_data[ 0 ]= 0.0
x_data[ 1 ]= 0.1 , y_data[ 1 ]= 6.0
x_data[ 2 ]= 0.2 , y_data[ 2 ]= 14.1
x_data[ 3 ]= 0.3 , y_data[ 3 ]= 5.01
x_data[ 4 ]= 0.4 , y_data[ 4 ]= 4.326
x_data[ 5 ]= 0.5 , y_data[ 5 ]= 4.326
x_data[ 6 ]= 0.6 , y_data[ 6 ]= 4.326
x_data[ 7 ]= 0.7 , y_data[ 7 ]= 4.326
x_data[ 8 ]= 0.8 , y_data[ 8 ]= 4.326
x_data[ 9 ]= 0.9 , y_data[ 9 ]= 4.326
x_data[ 10 ]= 1.0 , y_data[ 10 ]= 4.326
x_data[ 11 ]= 1.1 , y_data[ 11 ]= 4.326
x_data[ 12 ]= 1.2 , y_data[ 12 ]= 4.326
x_data[ 13 ]= 1.3 , y_data[ 13 ]= 4.326
x_data[ 14 ]= 1.4 , y_data[ 14 ]= 4.326
x_data[ 15 ]= 1.5 , y_data[ 15 ]= 4.326
x_data[ 16 ]= 1.6 , y_data[ 16 ]= 4.326
x_data[ 17 ]= 1.7 , y_data[ 17 ]= 4.326
x_data[ 18 ]= 1.8 , y_data[ 18 ]= 4.326
x_data[ 19 ]= 1.9 , y_data[ 19 ]= 0.0
A=
[[ 0.  0.  0.  0.]
 [ 0.  0.  0.  0.]
 [ 0.  0.  0.  0.]
 [ 0.  0.  0.  0.]]
C=
[ 0.  0.  0.  0.]
fit data to  3  degree polynomial
A[ 0 ][ 0 ]= 20.0
A[ 0 ][ 1 ]= 19.0
A[ 0 ][ 2 ]= 24.7
A[ 0 ][ 3 ]= 36.1
Y[ 0 ]= 90.0
A[ 1 ][ 0 ]= 19.0
A[ 1 ][ 1 ]= 24.7
A[ 1 ][ 2 ]= 36.1
A[ 1 ][ 3 ]= 56.2666
Y[ 1 ]= 76.302
A[ 2 ][ 0 ]= 24.7
A[ 2 ][ 1 ]= 36.1
A[ 2 ][ 2 ]= 56.2666
A[ 2 ][ 3 ]= 91.333
Y[ 2 ]= 91.7046
A[ 3 ][ 0 ]= 36.1
A[ 3 ][ 1 ]= 56.2666
A[ 3 ][ 2 ]= 91.333
A[ 3 ][ 3 ]= 152.45581
Y[ 3 ]= 126.5949
C[ 0 ]= 5.07760587239
C[ 1 ]= 1.45271951094
C[ 2 ]= -2.0730744555
C[ 3 ]= 0.333826894602
xmin= 0.0 , xmax= 1.9 , ymin= 0.0 , ymax= 14.1
max= 8.81210258849  at i= 2 , xbad= 0.2 , ybad= 14.1
C[ 0 ]= 5.07760587239
C[ 1 ]= 1.45271951094
C[ 2 ]= -2.0730744555
C[ 3 ]= 0.333826894602
max_err= 8.81210258849
avg_err= 1.48478125305
rms_err= 2.45605363066
A=
[[ 0.  0.  0.  0.  0.]
 [ 0.  0.  0.  0.  0.]
 [ 0.  0.  0.  0.  0.]
 [ 0.  0.  0.  0.  0.]
 [ 0.  0.  0.  0.  0.]]
C=
[ 0.  0.  0.  0.  0.]
fit data to  4  degree polynomial
A[ 0 ][ 0 ]= 20.0
A[ 0 ][ 1 ]= 19.0
A[ 0 ][ 2 ]= 24.7
A[ 0 ][ 3 ]= 36.1
A[ 0 ][ 4 ]= 56.2666
Y[ 0 ]= 90.0
A[ 1 ][ 0 ]= 19.0
A[ 1 ][ 1 ]= 24.7
A[ 1 ][ 2 ]= 36.1
A[ 1 ][ 3 ]= 56.2666
A[ 1 ][ 4 ]= 91.333
Y[ 1 ]= 76.302
A[ 2 ][ 0 ]= 24.7
A[ 2 ][ 1 ]= 36.1
A[ 2 ][ 2 ]= 56.2666
A[ 2 ][ 3 ]= 91.333
A[ 2 ][ 4 ]= 152.45581
Y[ 2 ]= 91.7046
A[ 3 ][ 0 ]= 36.1
A[ 3 ][ 1 ]= 56.2666
A[ 3 ][ 2 ]= 91.333
A[ 3 ][ 3 ]= 152.45581
A[ 3 ][ 4 ]= 259.72867
Y[ 3 ]= 126.5949
A[ 4 ][ 0 ]= 56.2666
A[ 4 ][ 1 ]= 91.333
A[ 4 ][ 2 ]= 152.45581
A[ 4 ][ 3 ]= 259.72867
A[ 4 ][ 4 ]= 449.40730666
Y[ 4 ]= 187.0537932
C[ 0 ]= 2.63290909089
C[ 1 ]= 34.9075161531
C[ 2 ]= -86.0911301232
C[ 3 ]= 70.239372117
C[ 4 ]= -18.3961961112
xmin= 0.0 , xmax= 1.9 , ymin= 0.0 , ymax= 14.1
max= 7.39675182025  at i= 2 , xbad= 0.2 , ybad= 14.1
C[ 0 ]= 2.63290909089
C[ 1 ]= 34.9075161531
C[ 2 ]= -86.0911301232
C[ 3 ]= 70.239372117
C[ 4 ]= -18.3961961112
max_err= 7.39675182025
avg_err= 1.3871928077
rms_err= 2.05223962573