diff --git a/matlab/k_order_pert.m b/matlab/k_order_pert.m
index fe0973964e84307923b047134b364a4a7866dba7..0fb2355ecf4c44cf4f0c9a8651fa1745a7f21568 100644
--- a/matlab/k_order_pert.m
+++ b/matlab/k_order_pert.m
@@ -159,7 +159,7 @@ for i=1:n1
             xx = x(:,m);
             y(:,(i-1)*n1*n2+(j-1)*n2+k) = xx;
             if j ~= i
-                y(:,(j-1)*n1*n2+(i-1)*n2+i) = xx;
+                y(:,(j-1)*n1*n2+(i-1)*n2+k) = xx;
             end
         end
     end