To estimate the model several times without orthogonality constraints do the following:

T1=[];
for i=1:5
  tic
  f1=parafac(X,5);
  t1=toc;
  T1=[T1 t1];
end

To do the same using orthogonality, e.g., in the emission mode do

T2=[];
for i=1:5
  tic
  f2=parafac(X,5,[],[0 1 0]);
  t2=toc;
  T2=[T2 t2];
end