%% Eksempel 1
    % Parametere
    k = 3;
    p = 0.7;
    
    % Regn ut punktsannsynligheter
    xs = 0:20;
    P = pdf('nbin', xs, k, p);
    
    % Plott punktsannsynligheter
    figure(1);
    x = 0:20;
    bar(x+3, P, 0.2);
    xlabel('Antall syke pasienter', 'FontSize', 26);
    ylabel('Sannsynlighet', 'FontSize', 26)
    title('Eks 1: Punktsannsynlighet', 'FontSize', 30);
    set(gca, 'FontSize', 26);
    xlim([-0.5, 20.5]);
    print(1, '-depsc2', 'Eks11Punkt.eps');
    
    % Plott kumulativ fordeling
    xs = linspace(-10, 21, 1000);
    F = cdf('nbin', xs, k, p);
    figure(2);
    plot(xs+3, F, 'LineWidth', 4);
    set(gca, 'FontSize', 26, 'xlim', [-1, 21], 'ylim', [-0.2, 1.2]);
    xlabel('Antall syke pasienter', 'FontSize', 26);
    ylabel('Kumulativ sannsynlighet', 'FontSize', 26);
    title('Eks 1: Kumulativ fordeling', 'FontSize', 30);
    print(2, '-depsc2', 'Eks11Kum.eps');
    
%% Eksempel 2: Forskjellige parametere i binomisk fordeling
    % Iterer gjennom forskellige verdier av n
    Ks = [3, 6, 9];
    p = 0.7;
    h = figure(3);
    set(h, 'pos', [10 10 650 1000], 'PaperPositionMode', 'auto')
    
    % Iterer gjennom forskjellige verdier av k
    for i=1:3
        subplot(3, 2, i);
        k = Ks(i);
        xs = 0:25;
        pF = pdf('nbin', xs, k, p);
        bar(xs+k, pF, 0.3)
        xlim([-0.5, 25+0.5]);
        title(strcat('k = ', num2str(k), ', 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.3, 0.4, 0.5];
    for i=1:3
        subplot(3, 2, i+3);
        p = P(i);
        k = 3;
        xs = 0:25;
        pF = pdf('nbin', xs, k, p);
        bar(xs+k, pF, 0.3)
        xlim([-0.5, 25+0.5]);
        title(strcat('k = ', num2str(k), ', 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', 'Eks12Punkt.eps')
 
     % Iterer gjennom forskellige verdier av k
    Ks = [3, 6, 9];
    p = 0.7;
    h = figure(4);
    set(h, 'pos', [10 10 650 1000], 'PaperPositionMode', 'auto')
    for i=1:3
        subplot(3, 2, i);
        k = Ks(i);
        xs = linspace(-10, 25, 1000);
        pF = cdf('nbin', xs, k, p);
        plot(xs+k, pF, 'LineWidth', 4)
        xlim([-0.5, 25+0.5]);
        title(strcat('k = ', num2str(k), ', 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
    P = [0.3, 0.4, 0.5];
    for i=1:3
        subplot(3, 2, i+3);
        p = P(i);
        k = 3;
        xs = linspace(-10, 25, 1000);
        pF = cdf('nbin', xs, k, p);
        plot(xs+k, pF, 'LineWidth', 4)
        xlim([-0.5, 25+0.5]);
        title(strcat('k = ', num2str(k), ', 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', 'Eks12Kum.eps')
    
%% Eksempel 4: Geometrisk fordeling
    % Parametere
    p = 0.7;
    
    % Regn ut punktsannsynligheter
    xs = 0:20;
    P = pdf('geo', xs, p);
    
    % Plott punktsannsynligheter
    figure(5);
    bar(xs+1, P, 0.2);
    xlabel('Antall pasienter', 'FontSize', 26);
    ylabel('Sannsynlighet', 'FontSize', 26)
    title('Eks 4: Punktsannsynlighet', 'FontSize', 30);
    set(gca, 'FontSize', 26);
    xlim([-1, 21]);
    ylim([0, 1]);
    print(5, '-depsc2', 'Eks14Punkt.eps');
    
    % Plott kumulativ fordeling
    xs = linspace(-10, 25, 1000);
    F = cdf('geo', xs, p);
    figure(6);
    stairs(xs+1, F, 'LineWidth', 4);
    set(gca, 'FontSize', 26, 'xlim', [-1, 21], 'ylim', [-0.2, 1.2]);
    xlabel('Antall placebo', 'FontSize', 26);
    ylabel('Kumulativ sannsynlighet', 'FontSize', 26);
    title('Eks 4: Kumulativ fordeling', 'FontSize', 30);
    print(6, '-depsc2', 'Eks14Kum.eps');
    
%% Eksempel 5: Forskjellige parametere i geometrisk fordeling
    % Iterer gjennom forskellige verdier av p
    Ps = [0.2, 0.3, 0.4, 0.5];
    
    h = figure(7);
    set(h, 'pos', [10 10 600 800], 'PaperPositionMode', 'auto')
    for i=1:4
        subplot(2, 2, i);
        xs = 0:20;
        p = Ps(i);
        pF = pdf('geo', xs, p);
        bar(xs+1, pF, 0.3)
        xlim([-0.5, 20+0.5]);
        title(strcat('p = ', num2str(p)), 'FontSize', 20);
        xlabel('Antall', 'FontSize', 16);
        ylabel('Sannsynlighet', 'FontSize', 16);
        set(gca, 'FontSize', 20);
        ylim([0, 1.1])
    end
     
    print(7, '-depsc2', 'Eks15Punkt.eps')
 
    % Iterer gjennom forskellige verdier av p
    h = figure(8);
    set(h, 'pos', [10 10 600 800], 'PaperPositionMode', 'auto')
    for i=1:4
        subplot(2, 2, i);
        p = Ps(i);
        xs = linspace(-1, 20+1, 1000);
        Cf = cdf('geo', xs, p);
        plot(xs+1, Cf, 'LineWidth', 4);
        xlim([-0.5, 20+0.5]);
        title(strcat('p = ', num2str(p)), 'FontSize', 20);
        xlabel('Antall', 'FontSize', 16);
        ylabel('Sannsynlighet', 'FontSize', 16);
        set(gca, 'FontSize', 20);
        ylim([-0.2, 1.2])
    end
    
    print(8, '-depsc2', 'Eks15Kum.eps')
    
%% Eksempel 7: F?rste premie i Lotto
    % Parametere
    p = 1.86e-7;
    
    % Regn ut Kumulativ fordeling
    xs = linspace(1, 10*5.38e6, 1000);
    P1 = cdf('geo', xs, p);
    
    % Plott kumulativ fordeling
    figure(9);
    plot((xs+1)/(52*1000), P1, 'LineWidth', 4);
    xlabel('Antall millennium', 'FontSize', 26);
    ylabel('Kumulativ sannsynlighet', 'FontSize', 26)
    title('Eks 7: Kumulativ fordeling', 'FontSize', 30);
    set(gca, 'FontSize', 26);
    xlim([-1, max(xs/(52*1000))]);
    print(9, '-depsc2', 'Eks18Kum.eps');
    
%% Eksempel 9: SMS-eksempel
    mu = 10;
    
    % Punktsannsynlighet
    figure(10);
    xs = 0:30;
    pF = pdf('poiss', xs, mu);
    bar(xs, pF, 0.3)
    xlim([-0.5, max(xs)+0.5]);
    title('Eksempel 9: Punktsannsynlighet', 'FontSize', 20);
    xlabel('Antall', 'FontSize', 16);
    ylabel('Sannsynlighet', 'FontSize', 16);
    set(gca, 'FontSize', 20);
 %   ylim([0, 1.1])
     
    print(10, '-depsc2', 'Eks19Punkt.eps')
 
    % kumulativ fordeling
    h = figure(11);
    xs = linspace(-1, 30+1, 1000);
    Cf = cdf('poiss', xs, mu);
    plot(xs, Cf, 'LineWidth', 4);
    xlim([-0.5, 30+0.5]);
    title('Eksempel 9: Kumulativ fordeling', 'FontSize', 20);
    xlabel('Antall', 'FontSize', 16);
    ylabel('Sannsynlighet', 'FontSize', 16);
    set(gca, 'FontSize', 20);
    ylim([-0.2, 1.2])
    
    print(11, '-depsc2', 'Eks19Kum.eps')
    
%% Eksempel 11: Forskjellige parametere i Poissonfordeling
    % Iterer gjennom forskellige verdier av mu
    muS = [0.5, 1, 2, 4];
    
    h = figure(12);
    set(h, 'pos', [10 10 600 800], 'PaperPositionMode', 'auto')
    for i=1:4
        subplot(2, 2, i);
        xs = 0:20;
        mu = muS(i);
        pF = pdf('poiss', xs, mu);
        bar(xs, pF, 0.3)
        xlim([-0.5, 20+0.5]);
        title(strcat('mu = ', num2str(mu)), 'FontSize', 20);
        xlabel('Antall', 'FontSize', 16);
        ylabel('Sannsynlighet', 'FontSize', 16);
        set(gca, 'FontSize', 20);
        ylim([0, 1.1])
    end
     
    print(12, '-depsc2', 'Eks21Punkt.eps')
 
    % Iterer gjennom forskellige verdier av mu
    h = figure(13);
    set(h, 'pos', [10 10 600 800], 'PaperPositionMode', 'auto')
    for i=1:4
        subplot(2, 2, i);
        mu = muS(i);
        xs = linspace(-1, 20+1, 1000);
        Cf = cdf('poiss', xs, mu);
        plot(xs, Cf, 'LineWidth', 4);
        xlim([-0.5, 20+0.5]);
        title(strcat('mu = ', num2str(mu)), 'FontSize', 20);
        xlabel('Antall', 'FontSize', 16);
        ylabel('Sannsynlighet', 'FontSize', 16);
        set(gca, 'FontSize', 20);
        ylim([-0.2, 1.2])
    end
    
    print(13, '-depsc2', 'Eks21Kum.eps')
    