m=0.;
t0=[0 400];
x0=[0 1];
n=20;
[t,y]=ode45(@(t,y)gyak4(t,y,m), t0,x0);
plot(y(:,1), y(:,2));
 
%körrõl indít kp körül, hova megy ,idõ fordítva


