杨小玲对教育事业:请问高手,我的MATLAB程序怎么运行不了啊

% Pareto distribution generate methods
r=1.2;
a=rand(1,4000);
%end
for i=1:4000;
p(i)=(a(i)^(-1/r));
%sum=sum+p(i);
if p(i)>100
p(i)=100;
end
end

c=p*20;
d=0;
ii=0;
b=rand(1,1)
if b>0.5
ii=1;
end
if ii==1
for i=1:4000;
if mod(i,2)==0;
d0=zeros(1,fix(c(i)));
else
d0=ones(1,fix(c(i)));
end
d=[d,d0];
end
else
for i=1:4000;
if mod(i,2)==0;
d0=ones(1,fix(c(i)));
else
d0=zeros(1,fix(c(i)));
end
d=[d,d0];
end
end
d=d(2:max(size(d)));
for j=1: 50
pareto_0
suma(j,1:20000)=d(1:20000);
j
end
self0=sum(suma(:,1:20000));
for k=1:20000/20
self1(k)=sum(self0((k-1)*20+1:k*20));
end
return;