function v=runge_kutta_h(f,x0,u0,h,k)


%pomocna funkcija koja tabelira funkciju u(x) na intervalu [x0,x0+(k-1)*h]  
%runge kuta metodom 

v=zeros(1,k);
v(1)=u0;
for i=1:k-1
    k1=h*f(x0+(i-1)*h,v(i));
    k2=h*f(x0+(i-1/2)*h,v(i)+k1/2);
    k3=h*f(x0+(i-1/2)*h,v(i)+k2/2);
    k4=h*f(x0+i*h,v(i)+k3);
    v(i+1)=v(i)+(k1+2*k2+2*k3+k4)/6;% racunamo vrednost u(x0+i*h)
end