% Lager kvantilplott og Q-Q plott av n realisasjoner trukket fra standard normalfordelingen
clc, clear all, close all

n = 24;                         % Antall realisasjoner som skal trekkes
x = random('normal',0,1,1,24);  % Trekk realisasjoner fra standard normalfordelingen
x_sorted = sort(x);             % Sorter tallene
p = zeros(1,n);                 % Opprett vektor med plass til plotteposisjonene             
for m = 1:n                     
    p(m) = m/(n+1);             % Beregn plotteposisjonene
end

% KVANTILPLOTT
figure(1), hold on              % Åpne figurvindu, sett på hold
plot(x_sorted,p,'kx','LineWidth',2)   % Plott punktene som svarte kryss
xlim([-2.6 2.6])                      % Juster kanten av figurutsnittet
hold on
t = linspace(-2.6,2.6,100);        
F = cdf('normal',t,0,1);              % Evaluer cdf-en til standard normalfordelingen på samme intervall
plot(t,F,'r-')                           % Plott cdf-en som en rød linje
legend('(x_{(m)},m/25)','\Phi(x)','Location','NorthWest')   % Lag tegnforklaring
title('Kvantilplott','FontSize',14)           % Lag tittel, bruk skriftstørelse 14
hold off

% KVANTIL-KVANTIL (Q-Q) PLOTT
figure(2), hold on              % Åpne nytt figurvindu
inv_x_sorted = icdf('normal',p,0,1);    % Evaluer den inverse cdf-en til standard normalfordelingen i plotteposisjonene
plot(x_sorted,inv_x_sorted,'kx','LineWidth',2)  % Plott punktene
plot(t,t,'r-')                                  % Plott linja y = x
xlim([-2.6 2.6]); ylim([-2.6 2.6])              % Juster utsnitt
legend('(x_{(m)},\Phi^{-1}(m/25))','x','Location','NorthWest') 
title('Q-Q plott','FontSize',14)
hold off
