function solve_dgl()

% Set our preferences for ode45
% The default relative tolerance is 1e-3.
% To set our output to non-negative, we provide an array listing
%   each population that we want to constrain.  Since this example
%   has two populations, we pass the array [1 2]
options = odeset('RelTol', 1e-4, 'NonNegative', [1 2]);

% Use ode45 to solve our ODE
% Place the time points in a vector 't'
% Place the solution in a vector 'x'
[t,x] = ode45('lotka_volterra', [0 20], [10 10], options);



xLabel = 'Zeit t';
yLabel = 'Räuber/Beute';
Title = 'Räuber und Beute über Zeit';

h=figure('Visible', 'off'); % Es wird kein Plot-Fenster geöffnet
plot(t,x);
%set(p,'Color',[0.02,0.08,0.5]);
xlabel(xLabel);
ylabel(yLabel);
title(Title);
legend('Beute', 'Räuber');
print(h,'-dpng','rbt.png'); % Plot in Datei geschrieben


xLabel = 'Beute';
yLabel = 'Räuber';
Title = 'Phasenraum Beute(x), Räuber(y)';

h=figure('Visible', 'off'); % Es wird kein Plot-Fenster geöffnet
plot(x(:,1), x(:,2));
%set(p,'Color',[0.02,0.08,0.5]);
xlabel(xLabel);
ylabel(yLabel);
title(Title);
print(h,'-dpng','rbp.png'); % Plot in Datei geschrieben