@Preben Alsholm When I run the code it gives me a vector of length 16 and I know that the dimension should be 1, so thats good. There are a lot of errors in between though.

Thanks a lot!

@Preben Alsholm Yeah, I do get m:=4. And sorry, I should have mentioned this from the beginning. Matrices have complex entires, however, the eigenvalues must be real.

@Preben Alsholm when I try to use my matrices, for some reason the code stops at m:=nops(L); and doesn't give me the basis. Here is the modufied code that I am using

A := foldl(LinearAlgebra:-KroneckerProduct, X, Z, X, X);

B := foldl(LinearAlgebra:-KroneckerProduct, X, Z, Id, Z);

C := foldl(LinearAlgebra:-KroneckerProduct, Z, Id, Z, Id);

E := foldl(LinearAlgebra:-KroneckerProduct, Z, Z, Id, Z);

L := [A, B, C, E];

LE := `~`[[Eigenvectors]](L);

LE[1];

LEV := map2(op, 2, LE);

LEV[1]; seq(LEV[1][() .. (), i], i = 1 .. n);

m := nops(L);

IntersectionBasis([seq([seq(LEV[j][() .. (), i], i = 1 .. n)], j = 1 .. m)]);

for some prespeciefied X,Y,X and Id. Any idea what's going wrong and/or suggestions how to fix it?

Also my matrices are 16 by 16 and even the could would give me the answer I won't be able to see the basis. How can I address this?

@Preben Alsholm I am a biginner with maple. My matrices are very big, I cannot see the eigenvectors and I have no clue how to use the eigenvectors within the intersectbasis command. My use of the word recursive might have been incorrect. Either way, I am not able to get the answer.