A = [ 12667.67 100 ; 100 1 ]
 b = [ 24700 ; 220 ]
 linsolve(A,b)