PB Debutant Simulation VBA

Function AléatN(ByVal ÉT As Double, ParamArray Moy() As Variant) As Variant
Dim Ax As Double, Ay As Double, S As Double
Ax = Rnd * 2 - 1
Do: Ay = Ax: Ax = Rnd * 2 - 1: S = Ax * Ax + Ay * Ay: Loop Until S <= 1
S = Sqr(-2 * Log(S) / S) * ÉT
If UBound(Moy) = 1 Then
      AléatN = Array(Moy(0) + Ax * S, Moy(1) + Ay * S)
Else
      AléatN = Moy(0) + Ax * S
      End If
End Function

Function DistrQsN(Rnd0à1 As Double, Moyenne As Double, ÉcartType As Double) As Double
Rem. Distribution quasi normale à part que le nombre engendré ne fuira la moyenne de plus de 4 fois l'écart type
DistrQsN = (Rnd0à1 ^ 0.18148 - (1 - Rnd0à1) ^ 0.18148) * 4 * ÉcartType + Moyenne
End Function

Function DistrATH(Rnd0à1 As Double, Moyenne As Double, ÉcartMoitié As Double) As Double
Const Inv2ATanH½ = 191646774 / 210545501
Dim RndDispersé As Double: RndDispersé = Log(Rnd0à1 / (1 - Rnd0à1)) * Inv2ATanH½ ' = ATanH(2*Rnd0à1-1) / AtanH(0.5)
DistrATH = RndDispersé * ÉcartMoitié + Moyenne            ' Transformée de Fisher :  ATanH(x)=ln((1+x)/(1-x))/2
' ÉcartMoitié : Fourchette en + et - de la moyenne contenant statitistiquement la moitié des nombres à engendrer.
End Function