jacobi eigenvalue algorithm

jae is implemented as described here. the signature of the only exported method is (matrix, threshold). it would not examine the parameters for correctness, so check them yourself. threshold should be > 0. as it would internally calculate the number of iterations for given threshold and it would do ln(threshold), do not supply 0 as threshold - you can get infinite loop.