A17= matrix([[3.2,1.6,0,0],[1.6,-0.8,2.4,0],[0,2.4,-4.8,3.6],[0,0,3.6,2.4]]) 
b17= vector([-0.8,16.0,-39,10.2]) 
x17= A17.solve_right(b17) 
print x17 

SageDays/Kreyszig-18.1-17 (last edited 2010-08-11 13:18:27 by SWATI AHUJA)