clc
clear all
close all

% Dati
x = [6.24, 6.31, 6.28, 6.30, 6.25, 6.26, 6.24, 6.29, 6.22, 6.28];
y = [6.27, 6.25, 6.33, 6.27, 6.24, 6.31, 6.28, 6.29, 6.34, 6.27];

% Sistema di ipotesi
% H0: mu_x = mu_y
% H1: mu_x diverso da mu_y

% Livello di significatività
alpha = 0.05;

% Deviazioni standard note
sigma_x = 0.05;
sigma_y = 0.05;

% Z test a due campioni

% Valore osservato per la statistica del test
z_oss = (mean(x) -mean(y))/sqrt(sigma_x^2/length(x)+sigma_y^2/length(y))

% Valore critico 
z_c = norminv(1-alpha/2)

% Poichè |z_oss| < z_c non possiamo rifiutare H0

% P-value
p_value = 2*(1-normcdf(abs(z_oss)))
% Anche il p-value conferma che non possiamo rifiutare H0 perchè risulta
% maggiore del livello di significatività scelto