geta0.m

污染物连续衰变浓度变化 基于matlab的编程

function a0=Geta0(num,c0,k)
global a0
a0=zeros(num,1);
product=1;
sum=0;
a0(1)=c0(1);
for i=2:num
    for j=1:i-1
        for l=j:i-1
            product=k(l)/(k(l)-k(i))*product
        end
        sum=sum+c0(j)*product;
        product=1;
    end
    a0(i)=c0(i)+sum;
    sum=0;
end