A=([[0,2],[3,4]])
Z=matrix(2,2)
var('x,y')
L=matrix([[1,0],[x,1]])
U=matrix([[1,y],[0,1]])
Z=L*U
for i in range(2):
    for j in range (2):
        Z[i][j]==A[i][j]
f = function('f', x)
f = L
print f(3)
g=function('g',y)
g=U
print g(2)
R=f(3)*g(2)
print R
bool(R==A)

Strang-1.5-2 (last edited 2010-12-18 10:58:18 by sai kumar)