% Aufgabe 3

clc
clear all

[x,y]=ode45(@igelschnecke,[0,5],[12,4]);
plot(x,y(:,1))
hold on
grid on
plot(x,y(:,2),'r')
legend('Schnecken','Igel')
xlabel('t/s')
ylabel('N(t)')
title('Raeuber-Beute-System')
hold off

% Startwerte für Simulink
a1=5;
a2=5;
b1=0.5;
b2=0.5;

function dy=igelschnecke(t,y)
dy=zeros(2,1);
a1=5;
a2=5;
b1=0.5;
b2=0.5;
dy(1)=a1*y(1)-b1*y(1)*y(2);
dy(2)=b2*y(1)*y(2)-a2*y(2);

end