getc.m

来自「污染物连续衰变浓度变化 基于matlab的编程」· M 代码 · 共 27 行

function Con=GetC(num,aa,k)
global Con
% [v,D,num,c0,k]=Readfile;
% a0=Geta0(num,c0,k);
% t=200;
Con=zeros(num,1);
product=1;
sum=0;
Con(1)=aa(1);

for i=2:num
    for j=1:i-1
        for l=j:i-1
            product=product*k(l)/(k(l)-k(i));
        end
        sum=sum+product*Con(j);
        product=1;
    end
    Con(i)=aa(i)-sum;
    sum=0;
end