Computation of gyroscopic systems and symplectic eigensolutions of skew-symmetric matrices