global myrandomvalue

function myrandom(min, max)
  myrandomvalue = ((myrandomvalue*123+125)^2)  mod (256^4)
  return (myrandomvalue mod (max-min)) + min
end function

function seedmyrandom(val)
 myrandomvalue = val
end function

Function RndDouble!()
   Const TWO27! = 134217728.0      '2 ^ 27
   Const TWO29! = 536870912.0      '2 ^ 29

   rnd_state=RND_A*(rnd_state Mod RND_Q)-RND_R*(rnd_state/RND_Q)
   If rnd_state<0 rnd_state=rnd_state+RND_M
   Local r_hi! = rnd_state & $1ffffffc

   rnd_state=RND_A*(rnd_state Mod RND_Q)-RND_R*(rnd_state/RND_Q)
   If rnd_state<0 rnd_state=rnd_state+RND_M
   Local r_lo! = rnd_state & $1ffffff8

   Return (r_hi + r_lo/TWO27)/TWO29
End Function

AppTitle "Zufälle gibt's..."
Graphics 640, 480, 16, 2
Dim zahlen(50)

SeedRnd MilliSecs()

Repeat
SeedRnd MilliSecs()
Print "Wähle ein Modi:"
Print "1 - Blitzbasic Zufall"
Print "2 - eigener Zufall"
Print "3 - Statistischer vergleich beider Funktionen mit großen Zahlen"
wahl = Input(":|> ")
Cls
If wahl = 0 Then End
If wahl = 1 Then Gosub bbzufall
If wahl = 2 Then Gosub meinzufall
If wahl = 3 Then Gosub bignum

Forever

End

.bbzufall
Print "100 Zahlen zwischen 2 und 10:"
For i = 1 To 100
zahl = Rand(2,10)
If zahl = 2 Then zahlen(2) = zahlen(2) + 1
If zahl = 3 Then zahlen(3) = zahlen(3) + 1
If zahl = 4 Then zahlen(4) = zahlen(4) + 1
If zahl = 5 Then zahlen(5) = zahlen(5) + 1
If zahl = 6 Then zahlen(6) = zahlen(6) + 1
If zahl = 7 Then zahlen(7) = zahlen(7) + 1
If zahl = 8 Then zahlen(8) = zahlen(8) + 1
If zahl = 9 Then zahlen(9) = zahlen(9) + 1
If zahl = 10 Then zahlen(10) = zahlen(10) + 1
Next
For i = 2 To 10
Print "Zahl "+i+" "+zahlen(i)+" mal."
zahlen(i) = 0
Next
WaitKey()
Cls
Return

.meinzufall
Print "370 Zahlen zwischen 0 und 36:"
For i = 1 To 370
zahl = randomite(0,36)
For j = 0 To 36
If zahl = j Then zahlen(j) = zahlen(j) + 1
Next
Next
For i = 0 To 36 Step 2
Print "Zahl "+i+" "+zahlen(i)+" mal.          Zahl "+(i+1)+" "+zahlen(i+1)+" mal."
Gsammt = Gsammt + zahlen(i) + zahlen(i+1)
zahlen(i) = 0
zahlen(i+1) = 0
Next
Print "Summe: "+Gsammt
Print "(Taste)"
WaitKey()
Gsammt = 0
Cls
Return

.bignum
Print "5000 Zahlen zwischen 5 und 50:"
Print "Blitzbasic"
For i = 1 To 5000
zahl = Rand(5,50)
For j = 5 To 50
If zahl = j Then zahlen(j) = zahlen(j) + 1
Next
Next
For i = 5 To 50 Step 2
Print "Zahl "+i+" "+zahlen(i)+" mal.          Zahl "+(i+1)+" "+zahlen(i+1)+" mal."
Gsammt = Gsammt + zahlen(i) + zahlen(i+1)
zahlen(i) = 0
zahlen(i+1) = 0
Next
Print "Summe: "+Gsammt
Print "(Taste)"
WaitKey()
Gsammt = 0
Cls

Print "5000 Zahlen zwischen 5 und 50:"
Print "Eigene funktion"
For i = 1 To 5000
zahl = randomite(5,50)
For j = 5 To 50
If zahl = j Then zahlen(j) = zahlen(j) + 1
Next
Next
For i = 5 To 50 Step 2
Print "Zahl "+i+" "+zahlen(i)+" mal.          Zahl "+(i+1)+" "+zahlen(i+1)+" mal."
Gsammt = Gsammt + zahlen(i) + zahlen(i+1)
zahlen(i) = 0
zahlen(i+1) = 0
Next
Print "Summe: "+Gsammt
Print "(Taste)"
WaitKey()
Gsammt = 0
Cls
Return

Function randomite(start=0,ende=1)
bis = (ende+1) - start
a = 7^5
c = 1020
zahl = ((a * RndSeed() + c) Mod bis) + (start)
SeedRnd zahl
Return zahl
End Function