function dx=ahf2(t,x)
global gg;
dx=zeros(4,1);
dx(1)=x(3);
dx(2)=x(4);
vv=sqrt((x(3)-vt)^2+x(4)^2);
dx(3)=0.;
dx(4)=-gg;
end

