function f = tcl(i,j,n,p,nb,step,A)
N=n/step;
x=sum(rbeta([p,N,step],i,j),3);
E=step*i/(i+j);
sigma=(i*j/((i+j+1)*(i+j)^2))*step;
m=(cumsum(x')'.*(ones([1,p])'*(1./(1:N)))-E)/sqrt(sigma);
rf=pnorm(-A:2*A/nb:A,0,1);
index=cumsum(ones([1,nb])/(nb+2));
for k=1:N,
vecteur=(m(:,k))*sqrt(k);
e=quantile(vecteur,index);
plot(-A:2*A/nb:A,rf,e,index);
xlabel(sprintf('Fonction de repartition de la moyenne de
%d vas beta(%d,%d)',k*step,i,j));
if (k ~= N)
sprintf('Appuyez sur une touche pour la
suite'), pause;
end
end;
|