clc
close all
clear all

%------------------------------------------------- 
global alpha beta gamma delta;
    % PARAMETERZUWEISUNG
    alpha = 5;
    beta = 0.1;
    gamma = 3;
    delta = 0.05;
%-------------------------------------------------

% Startwerte
x = [10; 100];

% Festlegung d. Zeit und Lösung der Differentialgleichung
time = [0 5];
[t,y] = ode45(@lotka_volterra, time, x);

% Plot
figure, plot( t,y(:,1), t,y(:,2),'--' );
title('Raeuber-Beute Model nach Lotka-Volterra', 'fontsize', 16, 'fontweight', 'b');
xlabel({'Zeit'; 't [Jahre]'},'fontsize',14,'fontweight','l','color','black');
ylabel({'Population'; 'Räuber / Beute'},'fontsize',14,'fontweight','l','color','black');
legend('Beute', 'Räuber');

% AUSTEHENDE AUFGABE:
% - Simulation mit versch. Modellparameter