 A=vpa([0.003 59.14; 5.291 -6.13],4)
 b=vpa([59.17; 46.78],4)
 
 format long
 c=vpa(A(2,1)/A(1,1),4)
vpa(A(2,:)-vpa(c*vpa(A(1,:),4),4),4)
 A=vpa(A,4)
 b(2)=vpa(b(2)-vpa(c*vpa(b(1),4),4),4)
 b=vpa(b,4)
 y=vpa(b(2)/A(2,2),4)
 x=vpa((vpa(b(1)-vpa(A(1,2)*y,4),4))/A(1,1),4)