diff --git a/fdahpca.m b/fdahpca.m
index b64d0f9..355070d 100644
--- a/fdahpca.m
+++ b/fdahpca.m
@@ -110,7 +110,7 @@
             obj.vec = obj.vec.';
             vm = mean(obj.vec,2);
             for jj = 1:no
-                for ii = 1:size(gam,2)
+                for ii = 1:size(gam,1)
                     c(ii,jj) = (obj.vec(:,ii)-vm).'*obj.U(:,jj);
                 end
             end