function ratten

% ratten (T8H2ab_Ute)
% Simulation des Schlaf-Wach-Verhaltens der Ratten, es wird jeweils einmal
% in der Stunde bestimmt, wie viele Ratten wach sind.

%% Datenspeicher definieren

rat=ones(1,5); % zeigt, ob eine Ratte jeweils wach ist oder schläft
zaehler=zeros(1,5); % zählt, wie häufig jede Ratte schläft
prozent=zeros(1,5); % Schlafanteil in Prozent

%% Ausgangssituation herstellen (per Zufall)

% wach = 1
% schlafen = 0

a=randperm(5);

rat(a(1))=0;
rat(a(2))=0;

zaehler(a(1))=1;
zaehler(a(2))=1;

anzahl=sum(rat);
disp([num2str(anzahl) ' Ratten sind wach.'])

stund=1;

ben=input('Wiederholen? Nein (n) oder Ja (j). ','s');

%% Simulation

while ben~='n'

    for i=1:5
        if anzahl==0 | anzahl==1
           a5=randperm(5);
            a=a5(1);
            if a==1          % 20% Wahrscheinlichkeit für wach
                rat(i)=1;
            else             % 80% Wahrscheinlichkeit für schlaf
                rat(i)=0;
                zaehler(i)=zaehler(i)+1;
            end
        else
          a5=randperm(5);
            a=a5(1);
            if a==1          % 20% Wahrscheinlichkeit für schlaf
                rat(i)=0;
                zaehler(i)=zaehler(i)+1;
            else             % 80% Wahrscheinlichkeit für wach
                rat(i)=1;
            end
        end
    end

anzahl=sum(rat);
if anzahl==1
    disp('1 Ratte ist wach.')
else
    disp([num2str(anzahl) ' Ratten sind wach.'])
end

stund=stund+1;

ben=input('Wiederholen? Nein (n) oder Ja (j). ','s');

end

%% Resultat

disp(['Insgesamt wurden ' num2str(stund) ' Stunden simuliert.'])

for i=1:5
    prozent(i)=100*zaehler(i)/stund;
    disp(['Ratte ' num2str(i) ' hat ' num2str(prozent(i)) ' Prozent der Zeit geschlafen.'])
end

disp(['Das ergibt im Durchschnitt einen Schlafanteil von ' num2str(mean(prozent)) ' Prozent.'])