How to calculate determinant of matrices with symbolic entries fast?

I am using the det function in MATLAB directly to calculate determinant of matrices with symbolic entries but it is very slow and doesn't yield result in a considerable time. Are there any other options/functions to do this task fast? For example, symengine or MUPAD? I don't know much about these.

syms z1 z2 z3 z4

A=[0        1       0       0;
   0        0       1       0;
   0        0       0       1;
   -29.17   -56   -36.7   -10.1];

B1=[-1.55   1       0       0;
    -1    -0.3      0       0;
    0       0       0.5     0;
    -0.7    0       -0.34   -2.6];

B2=[0       0       0       0;
    1      1.5      4       0;
    0       0       0       0;
  -0.33     0       0      -1.1];
B3=[0       0       0       0;
    0       0       0       0;
    0       0       0       0;
  -0.08    -0.7     0       -1];

B4=zeros(4,4);B4(4,3)=-3;

n=4;

det(kron(A,eye(n))+kron(eye(n),A)+kron(B1,eye(n))*z1+kron(B2,eye(n))*z2+kron(B3,eye(n))*z3+kron(B4,eye(n))*z4+kron(eye(n),B1)*z1^(-1)+kron(eye(n),B2)*z2^(-1)+kron(eye(n),B3)*z3^(-1)+kron(eye(n),B4)*z4^(-1))