说明：  if derros(n)<derros(n-1) studyspace(1,n)=amlf*studyspace(1,n-1); end if derros(n)>=kkk*derros(n-1) studyspace(1,n)=betat*studyspace(1,n-1); end for j=1:HN wincrease(1,j)=-studyspace(1,n)*wdel(1,j)+mc*(w(1,j,n-1)-w(1,j,n-2)); w(1,j,n)=w(1,j,n-1)+wincrease(1,j); for i=1:IN vincrease(j,i)=-studyspace(1,n)*vdel(j,i)+mc*(v(j,i,n-1)-v(j,i,n-2)); v(j,i,n)=v(j,i,n-1)+vincrease(j,i); end aincrease(1,j)=-studyspace(1,n)*adel(1,j)+mc*(a(1,j,n-1)-a(1,j,n-2)); a(1,j,n)=a(1,j,n-1)+aincrease(1,j); bincrease(1,j)=-studyspace(1,n)*bdel(1,j)+mc*(b(1,j,n-1)-b(1,j,n-2)); b(1,j,n)=b(1,j,n-1)+bincrease(1,j); end end