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)