function y = AB3(a,b,h,y0,f)

  x = a:h:b;
  n = length(x);
  m = length(y0);
  
  y = zeros(m,n);
  y(:,1) = y0;
  
  % zunächst zwei Kutta-Schritte (3. Ordnung) durchführen:
  for i=1:2
      
      k1 = f(x(i),y(:,i));
      k2 = f(x(i)+1/2*h, y(:,i) + 1/2*h*k1);
      k3 = f(x(i)+h, y(:,i) + h*(-k1 + 2*k2));
      y(:,i+1) = y(:,i) + h*(1/6*k1+2/3*k2+1/6*k3);
      
  end
  
  % Startwerte stören später nicht vergessen:
  % ...
  
  % dann Adams-Bashforth 3-Schritt-Verfahren:
  for i=3:n-1
      
      y(:,i+1) = y(:,i) + h/12*(23*f(x(i),y(:,i)) - 16*f(x(i-1),y(:,i-1)) + 5*f(x(i-2),y(:,i-2)));
      
  end  
  
end