matlab衰减指数的程序是,衰减方程的数值计算程序（matlab和fortran)

衰减方程如下：

matlab程序：

C0=1;

k=0.36;

t=linspace(0,15,90);

dt=t(2)-t(1);

C=C0*exp(-k*t);%解析解

%下面求解数值解

C1=zeros(1,length(t));

C1(1)=C0;

%下面用显式求解

for i=2:length(t)

C1(i)=C1(i-1)-k*dt*C1(i-1);

end

%下面用隐式求解

C2=C1;

for i=2:length(t)

C2(i)=C2(i-1)/(1+k*dt);

end

%四阶龙格库塔

lambda=-k;

C3=C1;

for i=2:length(t)

C3(i)=C3(i-1)*(1+dt*lambda+0.5*(dt*lambda)^2+1/6*(dt*lambda)^3+1/24*(dt*lambda)^4);

end

%半隐半显

C4=C1;

for i=2:length(t)

C4(i)=C4(i-1)*(1-0.5*k*dt)/(1