getf("GG1.sci");
a=grand(1,20000,"exp",1);          // exponentielle parametre 1
d=grand(1,20000,"exp",1/1.001);        // exponentielle parametre 2
xset("window",0); xbasc();
moy=GG1(a,d)
xset("window",1); xbasc();
d=grand(1,1000,"exp",2);          // exponentielle parametre 0.5
moy=GG1(a,d)
