function haseluchs
figure(1)
clc
clear all
close all

lt=[1850:1897,1899,1900]; %es fehlt der Wert für 1989
ht=[1850:1900]; %Zeit läuft von 1850 - 1900
%Hasenfelle in 1000 einlesen von den Jahren 1850 bis 1900
hasen=[50 80 80 90 69 80 91 70 24 12 13 32 5 150 140 105 42 20 3 4 7 6 52 40 44 95 80 55 15 8 14 13 11 38 50 130 131 90 30 20 41 43 53 49 62 79 40 13 4 1 12];
%Luchsfelle einlesen von den Jahren 1850 - 1900   
luchse=[6 5 6 12 23 31 33 25 19 10 4 5 5 18 32 62 65 38 21 9 5 4 15 19 21 40 30 18 13 11 10 9 29 50 71 74 32 20 16 6 4 11 20 34 45 35 26 17 8 10];

%Luchsfelle plotten über Zeit

plot(lt,luchse,'r.-')
hold on
%Hasenfelle plotten über Zeit
plot(ht,hasen,'b.-')

%Mittelwerte berechnen
Hm=mean(hasen)
Lm=mean(luchse)
% Die Mittelwerte sind hm = 48.7451 und lm=22.64


%geschätzte Schwingungszeit T=9

omega=(2*pi)/9;
lambda=1;
%Startwerte
H0=20;
L0=10;


[t,z]=ode45(@HALU,[1850 1900],[H0 L0],[],Hm,Lm,omega,lambda)

plot(t,z(:,1),'b')%Hasen
plot(t,z(:,2),'r')%Luchse
legend('Luchse','Hasen','Hasen*','Luchse*')
title('Zeitdiagramm')

%figure(2)
%H=[30:5:70];
%L=[30:5:70];
%[h,l]=meshgrid(H,L);
%dh=omega*lambda*h.*(1-(l/Lm));
%dl=(omega/lambda).*l*((h/Hm)-1);
%quiver(h,l,dh,dl)
%Axis..ggf

function dHLdt = HALU(t,z,Hm,Lm,omega,lambda)
H=z(1)
L=z(2)

%ODEs
DH=omega*lambda*H*(1-(L/Lm)); %Änderungen Hasen mit der Zeit
DL=(omega/lambda)*L*((H/Hm)-1); %Änderungen Luchse mit der Zeit

dHLdt=[DH;DL]