A=matrix([[1,1,1],[1,2,2],[1,2,3]])
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([u == 1,v == 1,x == 1,u*x + 1 == 2,v*x + w == 2,y == 1,u*y + z == 2,v*y + w*z + 1 == 3],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-3 (last edited 2010-12-18 11:28:09 by sai kumar)