I have the simplified version of what I want to do:

restart: with(LinearAlgebra):
d:=4:

product(MatrixExponential(I*Matrix(d,{(l,l)=1})*A[l,l]),l=1..d):

 

But I get the following error:

Error, (in Matrix) integer indices required for Matrix

Help appreciated.