vertrauen ist wichtig

%% Vertrauensintervall
clc; clear all;

% Variablen
%x=[0.23,0.23,0.23,0.30,0.32,0.32,0.34,0.34,0.34 , 0.43,0.43,0.43,0.43,0.45,0.45,0.45,0.45,0.45 , 0.54,0.54,0.54,0.54,0.54,0.54,0.56,0.56,0.56, 0.62,0.65,0.65,0.65,0.67,0.67,0.68,0.76,0.76 ];
%N = length(x);
N = 18;
% Irrtumswahrscheinlichkeit
a = 0.1;
% Vetrauensintervall (0.95 ändern)
%a = 1 - 0.95;
% Vetrauensintervall länge
vl = 0.5;


% Mittelwert & Standardabweichung
%x_mean = mean(x)
%x_std = std(x) % =s
x_mean = 90;
x_std = 1.2;

if x_std > 3
    fprintf("Näherung durch Normalverteilung zulässig: %.1f > %.0f\n", x_std, 3)

    % Normalverteilung
    disp ' '
    disp '------------------------'
    disp 'Normalverteilung Z-Raum'
    disp '------------------------'
    z_oben = norminv(1-a/2, 0, 1)
    z_unten = x_mean - (z_oben * x_std) / sqrt(N)
    z_oben = x_mean + (z_oben * x_std) / sqrt(N)

else
    fprintf("Näherung durch Normalverteilung NICHT zulässig: %.1f < %.0f\n", x_std, 3)

    % Student-t
    disp ' '
    disp '------------------------'
    disp 'Student-t'
    disp '------------------------'
    % Annahmebereich bestimmen
    t_oben = tinv(1-a/2, N-1) % zweiseitiger Test
    %t_oben = tinv(1-a, N-1) % einseitiger Test
    t_unten = -t_oben;

    % Vertrauensintervall
    x_unten = x_mean - (t_oben * x_std) / sqrt(N)
    x_oben = x_mean + (t_oben * x_std) / sqrt(N)
end

laenge = 4* x_std^2 / (0.5)^2