clear all, clc;
global m g;
m=75;
g=9.81;

delta=1;
kf=m*g/2500;
ks=m*g/25;

H=[2000; 0]
allaH=H;
for t=0:delta:30-delta
    H=H+delta*fEXS711(t,H,kf);
    allaH=[allaH H];
end

delta=0.01*delta
while H(1)>0
    H=H+delta*fEXS711(t,H,ks);
    allaH=[allaH H];
end 
allaT=[0:30 30+delta:delta:+30+(length(allaH)-31)*delta];
plot(allaT,allaH);