Maximum likelihood estimation of factor analysis models - fitting function