% przed użyciem funkcji należy ją zapisać do pliku o nazwie `euler.m`

%%%%%%%%%%%%%%%%
%
% Metoda Eulera 
%
%%%%%%%%%%%%%%%%

function x = euler(f, x0, t0, T, N)
% Dla równania $x' = f(t,x), \quad t_0 < t \leq T$ z warunkiem początkowym $x(t_0) = x_0$  
% wyznacza x(T)
% metodą Eulera z krokiem $h = (T-t_0)/N$
x = x0; t = t0; h = (T-t0)/N;
for i = 1:N
	x = x + h * feval(f,t,x); % nowy x to x(t+h)
	t = t+h;
end
end
