reof for eof ，at first you can assum the countribution([new_loads, rotmax] = varimax( loadings, normalize, tolerance, it_max )
[yPr,T] = rotatefactors(rpcdata,'method','varimax','Normalize','off');
you can chose one)