%% Eksempel 2
    % Parametere
    n = 4;
    p = 0.7;
    
    % Regn ut punktsannsynligheter
    P = zeros(n+1, 1);
    for i = 1:(n+1)
        P(i) = nchoosek(n, i-1)*p^(i-1)*(1-p)^(n-i+1);
    end
    
    % Plott punktsannsynligheter
    figure(1);
    x = 0:n;
    bar(x, P, 0.2);
    xlabel('Antall syke pasienter', 'FontSize', 26);
    ylabel('Sannsynlighet', 'FontSize', 26)
    title('Eks 2: Punktsannsynlighet', 'FontSize', 30);
    set(gca, 'FontSize', 26);
    xlim([-1, 5]);
    print(1, '-depsc2', 'Eks2Punkt.eps');
    
    % Plott kumulativ fordeling
    xExtended = [-1, x, 10];
    F = cumsum(P);
    Fextended = [0; F; 1];
    figure(2);
    stairs(xExtended, Fextended, 'LineWidth', 4);
    set(gca, 'FontSize', 26, 'xlim', [-1, 5], 'ylim', [-0.2, 1.2]);
    xlabel('Antall syke pasienter', 'FontSize', 26);
    ylabel('Kumulativ sannsynlighet', 'FontSize', 26);
    title('Eks 2: Kumulativ fordeling', 'FontSize', 30);
    print(2, '-depsc2', 'Eks2Kum.eps');
    
%% Eksempel 3: Forskjellige parametere i binomisk fordeling
    % Iterer gjennom forskellige verdier av n
    N = [4, 10, 40];
    p = 0.7;
    h = figure(3);
    set(h, 'pos', [10 10 650 1000], 'PaperPositionMode', 'auto')
    for i=1:3
        subplot(3, 2, i);
        n = N(i);
        xs = 0:n;
        pF = pdf('binomial', xs, n, p);
        bar(xs, pF, 0.3)
        xlim([-0.5, n+0.5]);
        title(strcat('n = ', num2str(n), ', p = ', num2str(p)), 'FontSize', 20);
        xlabel('Antall', 'FontSize', 16);
        ylabel('Sannsynlighet', 'FontSize', 16);
        set(gca, 'FontSize', 20);
        ylim([0, 0.6])
    end
    
    % Iterer gjennom forskjellige verdier av p
    n = 10;
    P = [0.1, 0.5, 0.9];
    for i=1:3
        subplot(3, 2, i+3);
        p = P(i);
        xs = 0:n;
        pF = pdf('binomial', xs, n, p);
        bar(xs, pF, 0.3)
        xlim([-0.5, n+0.5]);
        title(strcat('n = ', num2str(n), ', p = ', num2str(p)), 'FontSize', 20);
        xlabel('Antall', 'FontSize', 16);
        ylabel('Sannsynlighet', 'FontSize', 16);
        set(gca, 'FontSize', 20);
        ylim([0, 0.6])
    end
   
    print(3, '-depsc2', 'Eks3Punkt.eps')
 
    % Iterer gjennom forskellige verdier av n
    N = [4, 10, 40];
    p = 0.7;
    h = figure(4);
    set(h, 'pos', [10 10 650 1000], 'PaperPositionMode', 'auto')
    for i=1:3
        subplot(3, 2, i);
        n = N(i);
        xs = linspace(-1, n+1, 1000);
        Cf = cdf('binomial', xs, n, p);
        plot(xs, Cf, 'LineWidth', 4);
        xlim([-0.5, n+0.5]);
        title(strcat('n = ', num2str(n), ', p = ', num2str(p)), 'FontSize', 20);
        xlabel('Antall', 'FontSize', 16);
        ylabel('Sannsynlighet', 'FontSize', 16);
        set(gca, 'FontSize', 20);
        ylim([-0.2, 1.2])
    end
    
    % Iterer gjennom forskjellige verdier av p
    n = 10;
    P = [0.1, 0.5, 0.9];
    for i=1:3
        subplot(3, 2, i+3);
        p = P(i);
        xs = linspace(-1, n+1, 1000);
        Cf = cdf('binomial', xs, n, p);
        plot(xs, Cf, 'LineWidth', 4);
        xlim([-0.5, n+0.5]);
        title(strcat('n = ', num2str(n), ', p = ', num2str(p)), 'FontSize', 20);
        xlabel('Antall', 'FontSize', 16);
        ylabel('Sannsynlighet', 'FontSize', 16);
        set(gca, 'FontSize', 20);
        ylim([-0.2, 1.2])
    end
   
    print(4, '-depsc2', 'Eks3Kum.eps')
    
%% Eksempel 7
    % Parametere
    N = 20;
    n = 4;
    k = 6;
    
    % Regn ut punktsannsynligheter
    xs = 0:n;
    P = pdf('hyge', xs, N, n, k);
    
    % Plott punktsannsynligheter
    figure(5);
    bar(xs, P, 0.2);
    xlabel('Antall placebo', 'FontSize', 26);
    ylabel('Sannsynlighet', 'FontSize', 26)
    title('Eks 7: Punktsannsynlighet', 'FontSize', 30);
    set(gca, 'FontSize', 26);
    xlim([-1, 5]);
    print(5, '-depsc2', 'Eks7Punkt.eps');
    
    % Plott kumulativ fordeling
    xs = linspace(-1, 5, 1000);
    F = cdf('hyge', xs, N, n, k);
    figure(6);
    stairs(xs, F, 'LineWidth', 4);
    set(gca, 'FontSize', 26, 'xlim', [-1, 5], 'ylim', [-0.2, 1.2]);
    xlabel('Antall placebo', 'FontSize', 26);
    ylabel('Kumulativ sannsynlighet', 'FontSize', 26);
    title('Eks 7: Kumulativ fordeling', 'FontSize', 30);
    print(6, '-depsc2', 'Eks7Kum.eps');
    
%% Eksempel 8: Forskjellige parametere i hyper-geometrisk fordeling
    % Iterer gjennom forskellige verdier av n
    N = 20;
    nAll = [4, 10, 15, 20];
    k = 6;
    
    h = figure(7);
    set(h, 'pos', [10 10 600 800], 'PaperPositionMode', 'auto')
    for i=1:4
        subplot(2, 2, i);
        n = nAll(i);
        xs = 0:k;
        pF = pdf('hyge', xs, N, n, k);
        bar(xs, pF, 0.3)
        xlim([-0.5, k+0.5]);
        title(strcat('N = ', num2str(N), ', n = ', num2str(n), ', k = ', num2str(k)), 'FontSize', 20);
        xlabel('Antall', 'FontSize', 16);
        ylabel('Sannsynlighet', 'FontSize', 16);
        set(gca, 'FontSize', 20);
        ylim([0, 1.1])
    end
     
    print(7, '-depsc2', 'Eks8Punkt.eps')
 
    % Iterer gjennom forskellige verdier av n
    N = 20;
    nAll = [4, 10, 15, 20];
    k = 6;
    
    h = figure(8);
    set(h, 'pos', [10 10 600 800], 'PaperPositionMode', 'auto')
    for i=1:4
        subplot(2, 2, i);
        n = nAll(i);
        xs = linspace(-1, k+1, 1000);
        Cf = cdf('hyge', xs, N, n, k);
        plot(xs, Cf, 'LineWidth', 4);
        xlim([-0.5, k+0.5]);
        title(strcat('N = ', num2str(N), ', n = ', num2str(n), ', k = ', num2str(k)), 'FontSize', 20);
        xlabel('Antall', 'FontSize', 16);
        ylabel('Sannsynlighet', 'FontSize', 16);
        set(gca, 'FontSize', 20);
        ylim([-0.2, 1.2])
    end
    
    print(8, '-depsc2', 'Eks8Kum.eps')
    