v1=[3,1]
v2=[2,1]
S1=matrix([v1,v2])
S=transpose(S1)
D=matrix([[1,0],[0,4]])
A=S*D*S^(-1)
print A
[-5 18]
[-3 10]

Strang-250-2-U (last edited 2010-12-17 05:05:08 by rupampy)