theta=var('theta')
A=matrix([[1,0,0],[0,cos(theta),sin(theta)],[0,-sin(theta),cos(theta)]])
Ain = A.inverse()
y1, y2,y3 = var('y1,y2,y3')
y = matrix([[y1],[y2],[y3]])
Ain*y

SageDays/Kreyszig-6.3-30 (last edited 2010-08-12 21:26:09 by bhanukiran)