damfail.m

2005MCMA题程序

function dZ=damfail(t,Z)
hour=60*60;
As=2*10^8;
r=0.001;
%r=0.007;
a=1.7*hour;
Z0=52.9;
H0=53;
b=1000;
%b=2000;
q=As*r/b;
c=H0-Z0-q;
dZ=[];
dZ(1)=-r*a*sqrt(q+c*exp(-(Z0-Z(1))/q));