clc
clear all
CA=zeros(20+1,1);
B=zeros(20+1,1);
CA(1)=1/2;
CA(2)=1/2;
l=2;
for D=0.1:0.1:0.5
    CL = []; % in each D loop, CL should be initialized to empty matrix.
    for t=0.1:0.1:1
        for i=2:20
            B(i)=-CA(i)*D;
        end
        CA=B*5
        CL(l)=CA(20)
        l=l+1
    end
    a=numel(CL)
end

CL = []; % in each D loop, CL should be initialized to empty matrix.

l = 2;

3楼: Originally posted by dnp at 2014-07-03 23:17:52
上面那个程序还有问题，l也要重新再定义。
即在

CL = []; % in each D loop, CL should be initialized to empty matrix.

后面加上

l = 2;