%zad 1----------

xi = [11.0552, 10.6477, 9.63204, 9.88118, 9.37967, 8.24137, 10.73796, 9.22886, 8.82199, 9.192863];
war=0.9523 % podana w zadaniu
n=length(xi);
p1=0.95;
p2=0.99;
a1=1-p1;
a2=1-p2;
sr = mean(xi) % -> do protokołu
warianc = var(xi) % z próbki -> do protokołu
odch_std = sqrt(war);

z1=1-(a1/2);
z2=1-(a2/2);


up1 = 1.96; %odczytane z tabeli 2 dla swojego z
up2 = 2.58; %odczytane z tabeli 2 dla swojego z

dol1 = sr-up1*odch_std/sqrt(n)
gor1 = sr+up1*odch_std/sqrt(n)

dol2 = sr-up2*odch_std/sqrt(n)
gor2 = sr+up2*odch_std/sqrt(n)


%zad 2-----
xi = [11.0552, 10.6477, 9.63204, 9.88118, 9.37967, 8.24137, 10.73796, 9.22886, 8.82199, 9.192863];
n=length(xi);
sr = mean(xi) % -> do protokołu
war = var(xi) % z próbki -> do protokołu
odch_std = sqrt(war);

p1=0.95;
p2=0.99;
a1=1-p1;
a2=1-p2;

z1=1-(a1/2)
z2=1-(a2/2)

t_z1 = 2.2622; %odczytane z tabeli 3 dla v=n-1 i z1
t_z2 = 3.2498; %odczytane z tabeli 3 dla v=n-1 i z2

dol1 = sr-t_z1*odch_std/sqrt(n)
gor1 = sr+t_z1*odch_std/sqrt(n)
dol2 = sr-t_z2*odch_std/sqrt(n)
gor2 = sr+t_z2*odch_std/sqrt(n)


%zad3-----

xi = [8.591,8.35633,8.95066,9.2165,9.26367,8.85530,9.94805,9.69714,9.06728,9.46711];
n=length(xi);
srednia = mean(xi) % -> do protokołu
wariancja = var(xi) % z próbki -> do protokołu

p1=0.95;
p2=0.99;
a1=1-p1;
a2=1-p2;

z1=(a1/2)
z2=1-(a1/2)

x_z1 = 2.7004; %odczytane z tabeli 4 dla v=n-1 i z1
x_z2 = 19.0228; %odczytane z tabeli 4 dla v=n-1 i z2

k1_p1 = ((n-1)*wariancja)/x_z2
k2_p1 = ((n-1)*wariancja)/x_z1

z3=a2/2
z4=1-(a2/2)

x_z3 = 1.7349; %odczytane z tabeli 4 dla v=n-1 i z1
x_z4 = 23.5894; %odczytane z tabeli 4 dla v=n-1 i z2

k1_p2 = ((n-1)*wariancja)/x_z4
k2_p2 = ((n-1)*wariancja)/x_z3