function [y]=rk4(f,x,y0)
% Beräkna en approximation av lösningen till begynnelsevärdes-
% problemet y'=f(x,y),  y(x(1))=y0  med hjälp av Runge-Kuttas 
% metod och steglängden h=x(n+1)-x(n).
%
N=length(x); h=x(2)-x(1);
y(1)=y0;
for n=1:N-1,
  k1=h*feval(f,x(n),y(n));
  k2=h*feval(f,x(n)+h/2,y(n)+k1/2);
  k3=h*feval(f,x(n)+h/2,y(n)+k2/2);
  k4=h*feval(f,x(n)+h,y(n)+k3);
  y(n+1)=y(n)+(1/6)*(k1+2*k2+2*k3+k4);
end

