#2
A2= matrix([[4,-1,0],[-1,4,-1],[0,-1,4]])
b2= vector([21,-45,33])
x2= A2.solve_right(b2)
print x2

SageDays/Kreyszig-18.3-2 (last edited 2010-08-12 12:37:25 by SWATI AHUJA)