% djur, kaniner och rävar
  clear, clf
  t=0; y=[300 150]; T=[]; Y=[]; tslut=9; 
  h=0.1;                                 % egen rungekutta
  while t<tslut-h/2
    f1=fkr(y);        f2=fkr(y+h*f1/2);
    f3=fkr(y+h*f2/2); f4=fkr(y+h*f3);
    y=y+h/6*(f1+2*f2+2*f3+f4); t=t+h;  
    T=[T; t]; Y=[Y; y]; 
  end
  [T,Y]=ode45('fdjur',[0 tslut],[300 150]');       % alternativet ode45
  kanin=Y(:,1); raev=Y(:,2);
  subplot(2,2,1), plot(T,kanin, T,raev,'--')
  title('Kaniner och rävar som funk.av tiden')
  subplot(2,2,2), plot(300,150,'x', kanin,raev), axis equal
  xlabel('Antal kaniner'), ylabel('Antal rävar')



