A=matrix([[1,0,1],[2,2,2],[3,4,5]])
Z=matrix(3,3)
var('y,z,u,v,w')
L=matrix([[1,0,0],[x,1,0],[y,z,1]])
U=matrix([[1,u,v],[0,1,w],[0,0,1]])
Z=L*U
for i in range(3):
    for j in range (3):
        Z[i][j]==A[i][j]        
solve([1== 1,u == 0,v == 1,x == 2,u*x + 1 == 2,v*x + w == 2,y == 3,u*y + z == 4,v*y + w*z + 1 == 5],x,y,z,u,v,w)
f=function('f',x,y,z)
f=L
print f(1,1,1)
g=function('g',u,v,w)
g=U
print g(1,1,1)
P=f(1,1,1)*g(1,1,1)
print P
d=bool(P==A)
print d

*sadhana,student,snist *sai kumar,student,mriet

Strang-1.5-24 (last edited 2010-12-18 12:08:20 by sai kumar)