623.sagews

#an  ill-conditioned problem
#taken from Tiomthy Sauer's Numerical Analysis page 94
m2 = matrix([[1,1],[1.0001,1]])
#show that this has a large condition number
print condition(m2)

40004.0001

b = vector([2,2.0001])
x = m2.solve_right(b)
#answer should be (1,1)
print x

(1.00000000000222, 0.999999999997779)