Numerical Determination of a Canonical Form of a Symplectic Matrix