 function y = agyu(dy0)
    global h L
    [tv,yv] = euler2(h,0,L,0,dy0);
    plot(tv,yv,'o-', 'LineWidth',2);
    y = yv(L/h+1);	    
    return;
%% Az explicit Euler módszer    
    function [tv,yv] = euler2(h,t0,tmax,y0,dy0)   

    g = 9.8;	
    y = y0;	    
    dy = dy0;	 
    tv = [t0];
    yv = [y0];
    v=1;   % a konstans sebesség
    for t = t0:h:tmax
   %% Euler     
	y = y+dy*h; 
	dy = dy+(-g*h)/(v*v);
  %%  
	tv = [tv; t+h];
	yv = [yv; y];
    end
    return;