Générateur de nombres pseudo-aléatoires

rem ============================================================================
rem          Générateur de nombres pseudo-aléatoires
rem                      Par Papydall
rem ============================================================================
rem Génération de nombres pseudo-aléatoires distribués uniformément entre 0 et 1
rem selon la formule de Lehmer et l algorithme et les paramètres de
rem Parker and Miller (1988):
rem
rem              ____________________________
rem             |                            |
rem             | x(i+1) = a * x(i) modulo c |
rem             |____________________________|
rem
rem
rem Avec :
rem    a = 7^5      = 16807
rem    c = 2^31 - 1 = 2147483647
rem Sa période est fixée par c
rem ============================================================================

rem Appel :
rem Aleat() <===  retourne un nombre pseudo-aléatoire [0..1[ dans la variable globale RAND

rem ============================================================================
dim i
Application_Title "Générateur aléatoire Par Papydall"
list 10 : top 10,50 : left 10,100 : height 10,380 : item_add 10,"RND Papydall" : item_add 10,""
list 20 : top 20,50 : left 20,300 : height 20,380 : item_add 20,"RND Panoramic": item_add 20,""
font_bold 10 : font_bold 20 : width 10,150 : width 20,150
for i = 1 to 25 : Aleat() : item_add 10,rand : item_add 20,rnd(1) : pause 1000 : next i
Verif_Aleat(10000)
terminate
end
rem ============================================================================

SUB Random()
    dim_local a,c,x,y
    a = 16807 : c = 2147483647 : x = a * seed : seed = mod(x,c)
    y = seed/c : '  normalisation par c
    rand = y
END_SUB
rem ============================================================================
' Génération de la constante d'initialisation du générateur Random
' à partir de la date courante selon Anderson (1990)
SUB GenerationSeed()
    if variable("ret") = 0 then dim ret
    dim_local t0,t1,t2,t3,t4,t5,lseed,t$,d$

    t$ = time$  : ' Récupération de l'heure actuelle
    d$ = date$  : ' Récupération de la date actuelle

    t0 = val(mid$(t$,7,8)) : ' secondes
    t1 = val(mid$(t$,4,2)) : ' minutes
    t2 = val(mid$(t$,1,2)) : ' heures

    t3 = val(mid$(d$,1,2)) : ' Jour
    t4 = val(mid$(d$,4,2)) : ' Mois
    t5 = val(mid$(d$,7,4)) : ' année

' calcul du seed d'après l'algo d'Anderson (seed varie entre 0 et 2^31-1
    lseed = t5 + 70*(t4 + 12*(t3 + 31*(t2 + 23*(t1 + 59*t0))))
' s'assurer que seed est impair
    if even(lseed) = 1 then lseed = lseed - 1
    ret =  lseed
END_SUB
rem ============================================================================
SUB Aleat()
    if variable("rand") = 0 then dim rand
    if variable("seed") = 0 then dim seed
    GenerationSeed() : seed = ret : Random()
END_SUB
rem ============================================================================
SUB Verif_Aleat(NbTirage)
    dim_local i,Serie(NbTirage),Somme, Moyenne, EcartType,tex$
    print "  ..... Veuillez patienter, je calcule ....."
' Génération d'une série aléatoire de distribution uniforme
    for i = 0 to NbTirage : Aleat() : Serie(i) = rand : display : next i
' Vérification du caractère aléatoire de la série
' Pour une distribution aléatoire uniforme entre 0 et 1, la moyenne de
' la distribution doit être de 0.5 et son écart 0,289 (1/12)
    Somme = 0
    for i = 0 to NbTirage : Somme = Somme + Serie(i) : next i
    Moyenne = Somme / NbTirage

    Somme = 0
    for i = 0 to NbTirage : Somme = Somme + (Serie(i) - Moyenne)*(Serie(i) - Moyenne) : next i
    EcartType = sqr(Somme/Nbtirage)
    
    tex$ = "Moyenne : " + str$(Moyenne) + " pour " + str$(NbTirage) + " tirages" + chr$(13)
    tex$ = tex$ + "Ecart type : " + str$(EcartType) + " pour " + str$(NbTirage)+" tirages"+ chr$(13)+chr$(13)
    tex$ = tex$ + "Pour une distribution aleatoire uniforme entre 0 et 1," +chr$(13)
    tex$ = tex$ + "pour un nombre infini de tirages" + chr$(13)
    tex$ = tex$ + "- la moyenne vaut 0.5" + chr$(13)
    tex$ = tex$ + "- l'ecart-type vaut 1/12 soit 0.289"
    message tex$

END_SUB
rem ============================================================================

 for i = 1 to 25 : Aleat() : item_add 10,int(rand*100) : item_add 20,int(rnd(100)) : pause 1000 : next i

dim i
list 10 : top 10,10 : left 10,020 : height 10,400 : item_add 10,"RandomValue,100" : item_add 10,""
dll_on "kgf.dll"
for i = 1 to 25
    item_add 10,DLL_call1("RandomValue",100)
    ' On obtient toujours la même série : 1 , 4 , 87 , 21 , 28 , 68 , 32 , 17 etc.
next i

label start
dim i, m1%
m1% = 0 - 1
list 10 : top 10,10 : left 10,020 : height 10,400
button  1 : left 1,200 : caption 1,"Start" : on_click 1,start
dll_on "kgf.dll"
end
start:
clear 10
item_add 10,"RandomValue,100" : item_add 10,""
i = DLL_call1("RandomValue",m1%)
for i = 1 to 25
    item_add 10,DLL_call1("RandomValue",100)
    ' On obtient toujours la même série : 1 , 4 , 87 , 21 , 28 , 68 , 32 , 17 etc.
next i
return

rem ============================================================================
list 10 : top 10,10 : left 10,020 : height 10,400
dll_on "kgf.dll"
 Randomize()
 Random(100)
end
rem ============================================================================
SUB Random(n%)
    dim_local i
    item_add 10,"RandomValue,"+str$(n%) : item_add 10,""
    for i = 1 to 25
        item_add 10,DLL_call1("RandomValue",n%)
    next i
END_SUB
rem ============================================================================
SUB Randomize()
    dim_local bidon
    bidon = DLL_call1("RandomValue",0-1)
END_SUB
rem ============================================================================