rotated scores after varimax rotation

Dear all,

I'd like to verify if my approach is right. I'd like to obtain rotated subject scores after varimax rotation.

I have matrix M and perform pca:

[coeff, scores, latent, tsquared explained] = pca(M);

Then I rotate the first two PCs:

[RotatedCoeff,RotMatrix] = rotatefactors(coeff(:,1:2));

And compute rotated subject scores:

RotatedScores = scores(:,1:2)*RotMatrix;

Is this the right way to obtain rotated scores?

Thanks,

Best

Hans