First method:
Convert the output to ordinary loading matrices:

> [A1,B1,C1] = fac2let(Factors1,DimX,Fac);

and

> [A2,B2,C2] = fac2let(Factors2,DimX,Fac);

Then verify that A1 equals A2, for example using plot (plot([A1 A2])) or monitoring the correlation coefficients (corrcoef([A1 A2])). Similar can be done for the second and third mode loadings.

Second method:
Use the multiple cosine to monitor the similarity of all loading matrices simultaneously. This measure is more or less the same, but it applies to all loading matrices:

MultipleCosine = ncosine(Factors1,Factors2,DimX);