function raeuberbeute

h = 0.05;
tfin = 10;
T = 0:h:tfin;

a = 3;
b = 1.5;
c = 0.8;
d = 1.5;

q0 = [2 1];

[T,Q] = ode45(@(t,q) rhs(t,q,a,b,c,d), T, q0);

plot(T,Q)
xlabel('Zeit')
ylabel('u(t), v(t)')
title('Raeuber Beute Modellierung')
legend('Beute','Raeuber','location','northeast')
set(gcf,'paperpositionmode','auto')
set(gca,'xlim',[0 tfin],'ylim',[0 4.5],'fontsize',15)
print('-depsc','-loose','rb1');

function dq = rhs(t,q,a,b,c,d)
dq = q;
dq(1) = q(1)*(a-b*q(2));
dq(2) = q(2)*(c*q(1)-d);