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;
dt = 0.01    % Zeitschritt fuer Euler-Verfahren in Jahren

H=H0;
L=L0;

pH=[];
pL=[];
for t=1850:dt:1900
    
    DH=(omega*lambda*H*(1-(L/Lm)))*dt; %Änderungen Hasen mit der Zeit
    DL=((omega/lambda)*L*((H/Hm)-1))*dt; %Änderungen Luchse mit der Zeit
   
    
 H=H+DH;
 L=L+DL;
 
pH=[pH H];
pL=[pL L];

 
end

 plot(pH)  
 plot(pL, 'r')