How to “Estimate the model several times”

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