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