function eu=euler(x0,y0,x1,h)
% Eulerova metoda pro reseni ODR
global n x y

N = round((x1 - x0)/h);
x = zeros(N,1);
y= zeros(N,2);

x(1) = x0;
y(1,:) = y0;
n = 1;
while x(n) <= x1
  y(n+1,:) = y(n,:) + h*f(x(n),y(n,:));
  x(n+1) = x(n) + h;
  n = n+1;
end;
plot(x,y(:,1));

function ff = f(x,y)
ff = [y(2),exp(y(1))];




