function v=euler_h(f,x0,u0,h,k)
%pomocna funkcija koja tabelira funkciju u(x) na intervalu [x0,x0+(k-1)*h] modifikovanom Ojlerovom metodom 
v=zeros(1,k);
v(1)=u0;
for i=1:k-1
    u_z=v(i)+h/2*f(x0+(i-1)*h,v(i)); %pomocna vrednost u*
    v(i+1)=v(i)+h*f(x0+(i-1/2)*h,u_z);% racunamo vrednost u(x0+i*h)=u_i
end