n,theta=var('n,theta')
A = Matrix([[cos(n*theta),sin(n*theta)],[-sin(n*theta),cos(n*theta)]])
print A
print
print
print 'Determinant =',A.determinant()

SageDays/Kreyszig-6.6-6 (last edited 2010-08-12 20:26:04 by bhanukiran)