ASY-Forum

Code : Tout sélectionner

int n=100;
for(int k=0; k<=100; ++k) {
  real p_k=choose(n,k);
  write(k);
  write(p_k);
}

Code : Tout sélectionner

14
4.41869426773236e+16
15
2.53338471349989e+17
echantillonnage-1-graphique.asy: 20.3: Integer overflow

Compilation exited abnormally with code 1 at Wed May  8 09:51:10

Code : Tout sélectionner

import gsl;
real p=.7, pX;
int  n=300, nbsucces=100;
...
  
   pX=pdf_binomial(nbsucces,p,n);

Code : Tout sélectionner

import gsl;

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();
}