function solve_dgl_own()

time = 0;
i = 1;
tinc = .0001;
tend = 30;
x1 = 10;
x2 = 10;

x = zeros(tend/tinc, 2);

while(time < tend)

    x(i,1) = x1;
    x(i,2) = x2;
    t(i) = time;
    dx = lotka_volterra(time, [x1, x2]);
    x1 = x1 + dx(1)*tinc;
    x2 = x2 + dx(2)*tinc;
    time = time + tinc;
    i = i + 1;

end





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
