关于for循环绘图的问题

clc
clear
n=4;
m=1000;
s=1;
for Y=-3:0.001:3;
for k=1:s
X=randn(n,m);
A=mean(X);
B=sort(X);
w=0;
for j=1:m
h=0;
for i=1:n
if A(j)>=B(i,j)
h=h+1;
end
end
u=min(B(h+1,j)-n*(A(j)-B(h+1,j)),A(j));
l=B(h,j);
if Y>l && Y<u
w=w+1;
end
end
p=w/m;
end
P=1/s*sum(p);
hold on
plot(Y,P);
axis([-3 3 0 1])
xlabel('X0')
ylabel('概率')
grid on
end