function [x,y] = euler(fun,y0,a,b,n)
% [x,y] = euler(fun,y0,a,b,n)
%
% Funkce pocita diferencialni rovnici 
%      y' = f(x,y)   na <a,b>
% Eulerovou metodou
%
% fun ... nazev m-filu obsahujici funkci
% y0 ... pocatecni podminka
% a ... dolni mez intervalu
% b ... horni mez intervalu
% n ... pocet deleni
%
% x ... dikretizace <a,b>
% y ... hledana funkce na siti

h = (b-a)/n;
x = a:h:b;

y(1) = y0;
for k = 1:n
 y(k+1) = y(k) + h*feval(fun,x(k),y(k));
end

plot(x,y)

%end
