function [xup,xdown]=confintmean(x,alpha);
%
% Soit un ensemble de données supposées distribuées comme une population
% normale. On veut estimer un intervalle de confiance pour la moyenne 
% avec un degré de confiance alpha. On suppose que l'on ne sait rien de 
% la variabilité des données
%
% Input: x, la série de données
%        alpha, le degré de confiance 
% Output: xup, la borne supérieure de l'intervalle de confiance
%         xdown, la borne inférieure de l'intervalle de confiance
%
%  [xup,xdown]=confintmean(x,alpha)
%
xb=mean(x);
n=length(x);
s=sqrt((n*sum(x.^2)-sum(x)^2)/n/(n-1));
t=tq(1-alpha/2,n-1);
xup=xb+t*s/sqrt(n);
xdown=xb-t*s/sqrt(n);