function y=ojlerm2(imef,h)
% h - korak integracije
% Izlazni podatak:
% y - priblizno resenje
%
% Test: y=ojler('primer1',0.5)
%close all % zatvori sve grafike korisno za desktop varijantu
run(imef) % ucitaj skript fajl, tj. pocetni problem
%% formiranje mreze
x=x0:h:xm; % formiramo ravnomernu mrezu
n=length(x); % duzina mreze
%% inicijalizacija
y=zeros(1,n); % incijalizacija aproksimacije
y(1)=y0; % dodeljujemo pocetnu vrednost
%% iteracije
for j=1:n-1
yp=y(j)+h/2*f(x(j),y(j));
y(j+1)=y(j)+h*f(x(j),yp);
 % Druga modifikovana Ojlerova metoda
end
%% grafik aproksimacije
if nargout==0 % ako smo funkciju pozvali bez izlaznih argumenata, onda neka isctra grafik
plot(x,y, '-or')
end
 % u isom grafickom prozoru iscrtavamo i grafik tacnog resenja
hold on
x=linspace(x0,xm);
plot(x,ye(x), 'b')
legend('aproskimacija', 'tacno')
hold off
