function dx=ahf( t,x )
global gg alpha arhem;
dx=zeros(2,1);
dx(1)=x(2);
dx(2)=gg-alpha*abs(x(2))*x(2)-arhem;
end

