global gg alpha arhem;
x0=[0.;0];
p=100000;mu=29;tt=300;ru=8314;r=0.25;
m=0.01;gg=9.814;ro=p*mu/(ru*tt);
alpha=0.47*pi*r*r*ro/(2*m);arhem=ro*gg*4*pi*r*r*r/m;
options = odeset('RelTol',1e-8,'AbsTol',1e-8,'MaxStep',0.1);
[t,x]=ode45(@ahf,[0,20],x0,options);
hold on
plot(t,x(:,1))

%plot(t,x(:,1),'r');
%a=load('outh.dat');
