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