A=matrix([[3,3,4],[6,8,7],[-3,5,-9]])
L=matrix([[1,0,0],[2,1,0],[-1,4,1]])
U=matrix([[3,3,4],[0,2,-1],[0,0,-1]])
a=L*U
print det(A)
detL=det(L)
print detL
detU=det(U)
print detU
l=L**(-1) 
u=U**(-1)
detlu=det(l*u)
print detlu
detula=det(l*u*A)
print detula

Strang-4.2 27-U (last edited 2010-12-18 06:23:47 by shwetaradha)