function [x,y] = odj_euler(f,a,b,y0,n)
% implementacija Eulerove metode za rješavanje 
% inicijalnog problema ODJ
% Matematika 2

% ULAZNI PARAMETRI:   - f   => funkcija
%                     - a   => donji rub intervala
%                     - b   => gornji rub intervala
%                     - y0  => početni uvjet
%                     - n   => broj čvorova

% IZLAZNI PARAMETRI:  - x   => čvorovi na apscisi
%                     - y   => aproksimacije rješenja u čvorovima

h = (b-a)/(n-1);		% duljina podintervala

x(1) = a;		% početni čvor je a
y(1) = y0;		% početna vrijednost u donjem rubu

% koraci Eulerove metode
for i = 1:n-1
    x(i+1) = a + h*i;
    y(i+1) = y(i) + h*f(x(i),y(i));
end