%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                            Michael Pokojovy                             % 
%                                                                         %
% Nachkommastellen von Pi: http://www.eveandersson.com/pi/digits/100000   %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

fid = fopen('pi.txt');
PI = '';
tline = fgetl(fid);
while ischar(tline)
    PI = [PI tline];
    tline = fgetl(fid);
end
fclose(fid);

N = length(PI); % Anzahl der Nachkommastellen
M = 100;        % Anzahl der Einzelbilder
delay = 0.05;

figure(1);

set(gcf, 'PaperUnits', 'centimeters');
xSize = 21; ySize = 10;
xLeft = (21 - xSize)/2; yTop = (30 - ySize)/2;
set(gcf,'PaperPosition', [xLeft yTop xSize ySize]);
set(gcf,'Position',[0 0 xSize*50 ySize*50]);

for i = 1:M
    clf;
    hold on;
    freq = zeros(1, 10);    
    
    l = N*i/M;
    
    a = floor(a);
    b = floor(b);
    
    for k = 1:10
        freq(k) = length(find(PI(1:l) == num2str(k-1)))/l;
    end
    
    t = title(['Histogramm f\"ur die ersten ', num2str(l), ' Nachkommastellen von $\pi$'], 'interpreter', 'latex');
    set(t, 'FontSize', 16);
    
    bar(0:1:9, freq, 1, 'g');
    axis([-0.5 9.5 0 0.12]);
    xlabel('Ziffer');
    ylabel('relative Häufigkeit');
    
    pause(delay);
end