Discussion :

[PtitSuisse]

Bonjour, je cherche à calculer une intégrale double sur [0;1]².
J'ai simplifié le truc en faisant une double somme à partir de nombre aléatoires x et y générés au hasard entre 0 et 1.

Mon code ne marche pas.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
Sub Digamma()
Dim b As Long, a As Long, i As Integer, j As Integer, nu As Long, x As Long, y As Long
b = 0
a = 0
nu = Range("B7").Value
For j = 1 To 10000
x = Rnd(1)
b = b + a
    For i = 1 To 10000
    y = Rnd(1)
    a = a - (1 - x) * Log(x * y) * (x * y) ^ (nu - 1) / (1 - x * y)
    Next i
Next j
Range("AU33").Value = b
End Sub

Je vous remercie.

Ptit Suisse

[marsupilami34]

Bonjour

Qd je parcours ton code, dans la première boucle, tu as x=1 et y=1
Donc l'expression suivante

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
(1 - x * y)

vaut 0

La formule ci desous fais donc une div par 0

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
    a = a - (1 - x) * Log(x * y) * (x * y) ^ (nu - 1) / (1 - x * y)

[PtitSuisse]

c'est posisble
mais du coup si j'écris rnd(1), j'e veux un nombre compris entre 0 et 1 réels.
Si j'ajoute par exemple la condition x différent de y comme ca je suis sure que le dénominateur ne vaut pas 0, ca irait ?

Merci