This file contains Matlab code useful with the paper:

Gradient Projection Algorithms and Software for Arbitrary
Rotation Criteria in Factor Analysis.

By:

Coen A. Bernaards and Robert I. Jennrich.

Website: http://www.stat.ucla.edu/research

The Matlab function "randorth" will produce a random k by k 
orthogonal matrix. That is a random orthogonal matrix T uniformily 
distributed over the group of k by k orthogonal matricies. The code 
is very brief:

    function T=randorth(k)

    [T,R]=qr(randn(k,k));