为什么索引超出数组边界

f=@(a,b,c) a.sin(a)+b.^2+3c;
N=500;
d=3;
ger=500;
w=0.9;
limit=[-2,2;-2,2;-2,2];
vlimit=[-1.5,1.5;-1.5,1.5;-1.5,1.5];
c_2 = 0.5;
c_3 = 0.5;
tic;
for i=1:d
x(:,1)=limit(i,1)+(limit(i,2)-limit(i,1))*rand(N,1);
end

v=rand(N,d);
xm=x;
ym=zeros(1,d);
fxm=zeros(N,1);
fym=-inf;
plot3(xm(:,1),xm(:,2),f(xm(:,1),xm(:,2)),"mo");
title("种群初始状态分布图");
hold on;