function [t,y] = euler (f, y0, a, b, h)
% K�z�ns�ges differenci�l egyenletek kezdeti �rt�k probl�m�j�nak megold�sa
% Euler m�dszerrel
% Bemenet: f - diff. egy. rsz., y0 - kezdeti �rt�kek, 
% a,b - intervallum eleje/v�ge, h - feloszt�s nagys�ga
% Kimenet: t-f�ggetlen v�ltoz� �rt�kei, y - f�gg� v�ltoz� �s deriv�ltjai
n = round((b - a)/h);
t(1) = a;
y(1) = y0;
for i = 1 : n
    y(i + 1) = y(i) + h*f(t(i), y(i));
    t(i + 1) = t(i) + h;
end