function result = mf(X,model,mode,options);
%MF Matrix factorization for simplifying PCA and Tucker models.
%  Takes a PCA or Tucker model and rotates one mode to optimal simplicity.
%  The input is the data used to fit the original model, the model 
%  structure and the mode to be rotated. The remainder of the model is 
%  counterrotated so that the fit of the model is the same. In PCA 
%  the scores or loadings are counterrotated whereas in Tucker, the core 
%  is counterrotated.
%  The rotation is done by orthomax where default is to set the parameter 
%  gamma to one (corresponds to varimax) and where setting it to
%  zero would correspond to quartimax.
%  Optional structure input (options) can contain any of the fields:
%
%    display   [ {'off'} | 'on']   Governs screen display.
%    gamma     between 0 and 1. One is default.
%
%I/O: result = mf(X,model,mode,options);
%
%See also: PCA, TUCKER

% Copyright