Code : Tout sélectionner
usepackage("siunitx","locale=DE");
usepackage("fourier","upright");
locale("FR_fr");
import gsl;
// import math;
import geometry;
import graph;
///////////////////////////////////////////////
void FnIn(picture pic=currentpicture, int n, real p, real ua){
real pk,spk;
real dp=ua*sqrt(p*(1-p))/sqrt(n);
for(int k=0; k<=n; ++k){
//pk=choose(n,k)*p^k*(1-p)^(n-k);
pk=pdf_binomial(k,p,n);
if((k/n>=p-dp)&&(k/n<=p+dp)) spk+=pk;
draw(pic,rotate(90)*scale(.2)*Label(format((n>140)?"$\hphantom{\num{%.10f}}$":"$\num{%.10f}$",pk),EndPoint),
(k/n,0)--(k/n,pk),
blue+squarecap);
}
xaxis(pic,xmin=0,xmax=1,RightTicks(scale(.6)*Label(),step=1),above=true);
// label(pic,format("$%.3f$",p+dp),(p+dp,-.5));
distance(pic,"$I_{"+string(n)+"}$",(p-dp,0),(p+dp,0),7mm);
distance(pic,scale(.3)*Label("Intervalle de fluctuation asymptotique au seuil de 95\% :",align=realEpsilon*E),(0,0),(.01,0),7mm,nullpen);
distance(pic,scale(.5)*Label(format("$\approx\num{%.6f}$",p-dp),align=W),(p-dp-.1,0),(p-dp+.1,0),7mm,nullpen);
distance(pic,scale(.5)*Label(format("$\approx\num{%.6f}$",p+dp),align=E),(p+dp-.1,0),(p+dp+.1,0),7mm,nullpen);
distance(pic,format("$P\left(F_{"+string(n)+"}\in I_{"+string(n)+"}\right)\approx%.7f$",spk),(p-.1,0),(p+.1,0),1.2cm,nullpen);
label(pic,"$F_{"+string(n)+"}=\frac{X_{"+string(n)+"}}{"+string(n)+"}$
avec
$X_{"+string(n)+"}\hookrightarrow \mathcal{B}("+string(n)+";\num{"+string(p)+"})$
",truepoint(pic,NW),SE);
label(pic,scale(.5)*"GM",truepoint(pic,SW),NE,lightgray);
}
///////////////////////////////////////////////
real p=.7;
int N[]={30,40,50,60,70,80,90,100,110,120,130,140,150,200,250,300};
picture pict;
size(pict,8cm,4.5cm,IgnoreAspect);
for(int n:N){
FnIn(pict,n,p,1.96);
add(pict.fit());
erase(pict);
newpage();
}
Je n'y ai plus mis le nez dedans depuis les messages que tu as pu lire... et donc, il faudra que je regarde le code avant... car je ne diffuse pas un fichier où je mets mon nom, s'il y a des problèmes ou s'il y a des bouts de code de test qui trainent dans un coin. Mais dans l'instant, je n'ai vraiment pas le temps de m'y re-plonger.