clear all
aa=load('../Données/Soi.txt');
SOI=reshape(aa(:,2:13)',1,1692);
dm=[1,31,28,31,30,31,30,31,31,30,31,30];
sdm=(cumsum(dm)+15)/365;
for j=1:12
   t(j:12:12*size(aa,1))=aa(:,1)'+sdm(j);
end
t=t(SOI>-5);
SOI=SOI(SOI>-5);
soi2=sign(SOI).*SOI.^2;
p=sqrt(var(soi2))/(length(soi2)-1)

