dett_m=1/12; %monthly data
alpha=0.1126;
mu=0.0287;
sigma=0.0021;
r0=0.027227045583;
N=1000;
n=12;
array_list_r=zeros(n,N);
for k=1:N
array_list_r(1,k)=r0;
end
rng(1,'twister');
for i=1:N %N roads,here 1000
for j=1:n-1 %n months,here 12
array_list_r(j+1,i)=array_list_r(j,i)+alpha*(mu-array_list_r(j,i))dett_m+sigmasqrt...
(array_list_r(j,i))*normrnd(0,sqrt(dett_m));
for array_list_r(j+l,i)<0
array_list_r(j+1,i)=0;
end
end
hold on;
for c=1:N
plot(array_list_r(1,c))
xlim([1 n]);
end
save('array_list_r.mat','array_list_r');
而且他说我for array_list_r(j+l,i)<0这个部分是错误的 ,可不可以好心人帮我解答一下!!!!