How Can I slove this problem in sage?
asked 2015-01-17 09:50:35 -0500
This post is a wiki. Anyone with karma >750 is welcome to improve it.
Sorry for my bad English.
Write a procedure that accepts as input an orthogonal matrix A |- O(2m), with m >= 1, such that A^2= -I(2m) [subscript 2m]
and give as output an orthonormal basis of R^2m type:
u1,Au1,u2,Au2, ... , um,Aum`
Test the procedure on matrices of type A = tMJM [t apex] where M |- O(2M) is a random matrix and J |- O(2m) is the block matrix 2x2 all equal to
J = [column 1] ( 0 1) [column 2] (-1 0)