Dear Maple Support Team and Community,

    I would like to report what seems to be a bug in Maple's MatrixExponential function in the LinearAlgebra package.

    Consider the following 8×8 skew-Hermitian matrix B:
 

B := Matrix(8, 8, [[0, I, 0, -I, 0, I, 0, -I], 
                   [-I, 0, I, 0, -I, 0, I, 0], 
                   [0, -I, 0, I, 0, -I, 0, I], 
                   [I, 0, -I, 0, I, 0, -I, 0], 
                   [0, I, 0, -I, 0, I, 0, -I], 
                   [-I, 0, I, 0, -I, 0, I, 0], 
                   [0, -I, 0, I, 0, -I, 0, I], 
                   [I, 0, -I, 0, I, 0, -I, 0]]);

    Now compute the matrix exponential exp(I*t*B) using MatrixExponential

H := LinearAlgebra:-MatrixExponential(I * t * B);

However, this result is inconsistent with the matrix exponential computed using MATLAB and SageMath, both of which give the same result, different from Maple's.

I would appreciate it if someone from the development team could look into this.