% Trekker realisasjoner fra normalfordelingen, og plotter histogrammer med
% pdf-kurver for sammenlikning.


N = 1000; % Sett antall realisasjoner som skal trekkes
X = random('normal',0,1,N,1); % Trekk realisasjoner

% Sett antall realisasjoner som skal tas med i figuren
n1 = 100; n2 = N;
% Bestem antall underintervaller, rund ned til nærmeste heltall
nbins1 = sqrt(n1); nbins2 = floor(sqrt(n2));

% Lag figur
figure(1)
subplot(1,2,1) % Venstre halvdel av plottet
histfit(X(1:n1),nbins1,'normal')
title([num2str(n1) ' realisasjoner'])
xlabel('x')
ylabel('Frekvens')

subplot(1,2,2) % Høyre halvdel av plottet
histfit(X(1:n2),nbins2,'normal') 
title([num2str(n2) ' realisasjoner'])
xlabel('x')
ylabel('Frekvens')
