test_solve.py  a x = y 
a=
[[ 0.5         0.33333333  0.25        0.2       ]
 [ 0.33333333  0.25        0.2         0.16666667]
 [ 0.25        0.2         0.16666667  0.14285714]
 [ 0.2         0.16666667  0.14285714  0.125     ]]
y=
[1 2 3 4]
solve x=
[ -380.  3240. -7140.  4480.]
 
check y = y1 = a*x =
[1 2 3 4]
 
check y2 = a*x =
[[ -190.          1080.         -1785.           896.        ]
 [ -126.66666667   810.         -1428.           746.66666667]
 [  -95.           648.         -1190.           640.        ]
 [  -76.           540.         -1020.           560.        ]]
check for i in range(4): sum(y2[i])
1.0
2.0
3.0
4.0