clear all
load('../Données/OceanTide1000.mat')
p1=p1(~isnan(p1));
nb=10;
[k,x]=hist(p1,nb);
obs=k;
att=0*x+length(p1)/length(k);
U=sum((obs-att).^2./att)
tabl=[x;k; att];
chisqp(U,nb-1)