Erreur d'exécution 5

**aurama** · 20/03/2013, 05h09

Bonjour!

J'aimerais avoir un petit coup de main avec le code que je viens de créer.

J'essaie de calculer Beta(1.5,4.3,5,2) mais ça ne fonctionne pas. Ça m'indique l'erreur d'exécution 5 (Argument ou appel de procédure incorrect).

Avec le débogage, l'erreur d'exécution 5 apparaît dans la Function appelée fonction (ligne 27).

Cependant, si j'essaie de calculer Beta(1.5,4,5,2), ça fonctionne. On dirait que mon erreur vient du fait que le 4.3 n'est pas un entier.....

Si quelqu'un pourrait m'aider à trouver d'où provient mon erreur, ce serait très apprécié.

Merci d'avance et si vous avez besoin de plus d'informations n'hésitez pas à me faire signe.

Mon code est le suivant:

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
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
 
'********** Tâche # 1 : Programmation de fonctions VBA **********
 
 
Function z_i(x_i As Double) As Double
    'Fonction permettant de calculer la nouvelle variable z_i (suite au changement de variable)
 
    'Déclaration des variables
    Dim C As Double, D As Double
    C = 0 'Borne inférieure de l'intégrale
    D = 1 'Borne supérieure de l'intégrale
 
    'Déclaration d'une nouvelle variable
    Quotient = (D - C) / 2
 
    'Calcul de la nouvelle variable
    z_i = (Quotient * x_i) + ((C + D) / 2)
 
'Fin de la procédure
End Function
 
 
Function fonction(z_i, a As Double, b As Double) As Double
    'Fonction permettant de calculer le bêta
 
    'Calcul
    fonction = (z_i ^ (a - 1)) * ((1 - z_i) ^ (b - 1))
 
'Fin de la procédure
End Function
 
 
Function Beta(a As Double, b As Double, Optional n As Long = 10, Optional m As Byte = 5) As Variant
    'Fonction permettant de calculer la fonction bêta
    'Référence : http://en.wikipedia.org/wiki/Gaussian_quadrature
 
    'Déclaration des variables
    Dim w_i_1 As Double, w_i_2 As Double, w_i_3 As Double, w_i_4 As Double, w_i_5 As Double 'Poids
    Dim c_i As Double
    Dim x_i_1 As Double, x_i_2 As Double, x_i_3 As Double, x_i_4 As Double, x_i_5 As Double 'Points
    Dim i As Long
    Dim Somme As Double, Somme1 As Double, Somme2 As Double, Somme3 As Double, Somme4 As Double, Somme5 As Double
    Dim Quotient As Double, delta_x As Double
    Dim C As Double, D As Double
    C = 0 'Borne inférieure de l'intégrale
    D = 1 'Borne supérieure de l'intégrale
 
    'Message si a est négatif
    If a < 0 Then
        Message1 = MsgBox("La valeur de a est négative.", vbRetryCancel + vbCritical, "Problématique")
        If Message1 = vbCancel Then Exit Function
        If Message1 = vbRetry Then
            Message2 = MsgBox("Veuillez réessayer votre calcul avec un a positif.", vbOKOnly + vbInformation, "Problématique")
            If Message2 = vbOK Then Exit Function
        End If
    End If
 
    'Message si b est négatif
    If b < 0 Then
        Message3 = MsgBox("La valeur de b est négative.", vbRetryCancel + vbCritical, "Problématique")
        If Message3 = vbCancel Then Exit Function
        If Message3 = vbRetry Then
            Message4 = MsgBox("Veuillez réessayer votre calcul avec un b positif.", vbOKOnly + vbInformation, "Problématique")
            If Message4 = vbOK Then Exit Function
        End If
    End If
 
    'Message si a est nul
    If a = 0 Then
        Message5 = MsgBox("La valeur de a est nulle.", vbRetryCancel + vbCritical, "Problématique")
        If Message5 = vbCancel Then Exit Function
        If Message5 = vbRetry Then
            Message6 = MsgBox("Veuillez réessayer votre calcul avec un a positif.", vbOKOnly + vbInformation, "Problématique")
            If Message6 = vbOK Then Exit Function
        End If
    End If
 
    'Message si b est nul
    If b = 0 Then
        Message7 = MsgBox("La valeur de b est nulle.", vbRetryCancel + vbCritical, "Problématique")
        If Message7 = vbCancel Then Exit Function
        If Message7 = vbRetry Then
            Message8 = MsgBox("Veuillez réessayer votre calcul avec un b positif.", vbOKOnly + vbInformation, "Problématique")
            If Message8 = vbOK Then Exit Function
        End If
    End If
 
    'Longueur des n sous-intervalles de l'intervalle [C,D]
    delta_x = (D - C) / n 'n est le nombre de sous-intervalles
 
    'Déclaration d'une nouvelle variable
    Quotient = (D - C) / 2
 
    For i = 1 To n
        'Borne inférieure des n sous-intervalles
        c_i = C + (i - 1) * delta_x
 
        If m = 1 Then 'm est le nombre de points de la quadrature gaussienne
            'Déclaration des points et des poids qui sont déterminés grâce aux polynômes de Legendre
            x_i_1 = c_i + 0
            w_i_1 = 2
            'Somme du produit des poids et de f(z_i)
            Somme1 = w_i_1 * fonction(z_i(x_i_1), a, b)
 
        ElseIf m = 2 Then
            'Déclaration des points et des poids
            x_i_1 = c_i - ((3 ^ 0.5) / 3)
            w_i_1 = 1
            x_i_2 = c_i + ((3 ^ 0.5) / 3)
            w_i_2 = 1
            'Somme du produit des poids et de f(z_i)
            Somme2 = w_i_1 * fonction(z_i(x_i_1), a, b) + w_i_2 * fonction(z_i(x_i_2), a, b)
 
        ElseIf m = 3 Then
            'Déclaration des points et des poids
            x_i_1 = c_i + 0
            w_i_1 = 8 / 9
            x_i_2 = c_i - ((3 / 5) ^ 0.5)
            w_i_2 = 5 / 9
            x_i_3 = c_i + ((3 / 5) ^ 0.5)
            w_i_3 = 5 / 9
            'Somme du produit des poids et de f(z_i)
            Somme3 = w_i_1 * fonction(z_i(x_i_1), a, b) + w_i_2 * fonction(z_i(x_i_2), a, b) + w_i_3 * fonction(z_i(x_i_3), a, b)
 
        ElseIf m = 4 Then
            'Déclaration des points et des poids
            x_i_1 = c_i + (((3 - (2 * (6 / 5) ^ 0.5)) / 7) ^ 0.5)
            w_i_1 = (18 + (30 ^ 0.5)) / 36
            x_i_2 = c_i - (((3 - (2 * (6 / 5) ^ 0.5)) / 7) ^ 0.5)
            w_i_2 = (18 + (30 ^ 0.5)) / 36
            x_i_3 = c_i + (((3 + (2 * (6 / 5) ^ 0.5)) / 7) ^ 0.5)
            w_i_3 = (18 - (30 ^ 0.5)) / 36
            x_i_4 = c_i - (((3 + (2 * (6 / 5) ^ 0.5)) / 7) ^ 0.5)
            w_i_4 = (18 - (30 ^ 0.5)) / 36
            'Somme du produit des poids et de f(z_i)
            Somme4 = w_i_1 * fonction(z_i(x_i_1), a, b) + w_i_2 * fonction(z_i(x_i_2), a, b) + w_i_3 * fonction(z_i(x_i_3), a, b) + w_i_4 * fonction(z_i(x_i_4), a, b)
 
        ElseIf m = 5 Then
            'Déclaration des points et des poids
            x_i_1 = c_i + 0
            w_i_1 = 128 / 225
            x_i_2 = c_i + (((5 - (2 * (10 / 7) ^ 0.5)) ^ 0.5) / 3)
            w_i_2 = (322 + (13 * (70 ^ 0.5))) / 900
            x_i_3 = c_i - (((5 - (2 * (10 / 7) ^ 0.5)) ^ 0.5) / 3)
            w_i_3 = (322 + (13 * (70 ^ 0.5))) / 900
            x_i_4 = c_i + (((5 + (2 * (10 / 7) ^ 0.5)) ^ 0.5) / 3)
            w_i_4 = (322 - (13 * (70 ^ 0.5))) / 900
            x_i_5 = c_i - (((5 + (2 * (10 / 7) ^ 0.5)) ^ 0.5) / 3)
            w_i_5 = (322 - (13 * (70 ^ 0.5))) / 900
            'Somme du produit des poids et de f(z_i)
            Somme5 = w_i_1 * fonction(z_i(x_i_1), a, b) + w_i_2 * fonction(z_i(x_i_2), a, b) + w_i_3 * fonction(z_i(x_i_3), a, b) + w_i_4 * fonction(z_i(x_i_4), a, b) + w_i_5 * fonction(z_i(x_i_5), a, b)
 
        End If
 
        'Somme de la quadrature gaussienne
        Somme = Somme + Somme1 + Somme2 + Somme3 + Somme4 + Somme5
 
    Next i
    'Approximation de la fonction bêta avec la quadrature gaussienne
    Beta = Somme * Quotient
 
'Fin de la procédure
End Function

**ProgElecT** · 20/03/2013, 07h27

Salut, bienvenue sur DVP

Donner comme nom fonction à une Function n'est pas des plus lisible.
Lors de l'appel à cette fonction tu passes 3 paramètres, alors quelle n'en comporte que 2.

Code :

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

Function fonction(z_i, a As Double, b As Double) As Double

Différents appels, lignes 103, 112, 123 ....
Somme1 = w_i_1 * fonction(z_i(x_i_1) , a, b)
.........
Somme2 = w_i_1 * fonction(z_i(x_i_1) , a, b) + w_i_2 * fonction(z_i(x_i_2) , a, b)
..........
Somme3 = w_i_1 * fonction(z_i(x_i_1) , a, b) + w_i_2 * fonction(z_i(x_i_2) , a, b) + w_i_3 * fonction(z_i(x_i_3) , a, b)

**aurama** · 21/03/2013, 04h04

Bonjour!

Merci beaucoup ProgElecT pour ta réponse! Tu as raison, ce n'était pas très lisible, j'ai également changée cette fonction pour qu'elle ait 3 paramètres, comme tu me l'as suggéré!

Malheureusement, j'ai toujours le même problème...

Je cherche à calculer Beta(1.5,4.3,n,m) et cela avec différents n=5,10,50 et 100 et m=1,2,3,4 et 5.

Mon calcul fonctionne seulement si m=1, mais pour toutes les autres valeurs de m (peut importe le n utilisé) ça ne fonctionne pas.

Ça m'indique l'erreur d'exécution 5 (Argument ou appel de procédure incorrect).

Avec le débogage, l'erreur d'exécution 5 apparaît dans la Function appelée f (ligne 18).

Si quelqu'un pouvait m'aider à trouver d'où provient mon erreur, ce serait très apprécié.

Merci d'avance et si vous avez besoin de plus d'informations n'hésitez pas à me faire signe.

Voici mon nouveau code:

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
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
'********** Tâche # 1 : Programmation de fonctions VBA **********
 
 
Function f(k As Double, a As Double, b As Double) As Double
    'Fonction permettant de calculer le bêta
 
    'Déclaration des variables
    Dim z As Double
    Dim C As Double, D As Double
    C = 0 'Borne inférieure de l'intégrale
    D = 1 'Borne supérieure de l'intégrale
 
    'Déclaration d'une nouvelle variable
    Quotient = (D - C) / 2
 
    'Calculs
    z = ((Quotient * k) + ((C + D) / 2))
    f = (z ^ (a - 1)) * ((1 - z) ^ (b - 1))
 
'Fin de la procédure
End Function
 
 
Function Beta(a As Double, b As Double, Optional n As Long = 10, Optional m As Byte = 5) As Variant
    'Fonction permettant de calculer la fonction bêta
    'Référence : http://en.wikipedia.org/wiki/Gaussian_quadrature
 
    'Déclaration des variables
    Dim w1 As Double, w2 As Double, w3 As Double, w4 As Double, w5 As Double, _
    w6 As Double, w7 As Double, w8 As Double, w9 As Double, w10 As Double, _
    w11 As Double, w12 As Double, w13 As Double, w14 As Double, w15 As Double 'Poids
    Dim x1 As Double, x2 As Double, x3 As Double, x4 As Double, x5 As Double, _
    x6 As Double, x7 As Double, x8 As Double, x9 As Double, x10 As Double, _
    x11 As Double, x12 As Double, x13 As Double, x14 As Double, x15 As Double 'Points
    Dim j As Long
    Dim c_i As Double
    Dim Somme As Double, Somme1 As Double, Somme2 As Double, Somme3 As Double, Somme4 As Double, Somme5 As Double
    Dim Quotient As Double, delta_x As Double
    Dim C As Double, D As Double
    C = 0 'Borne inférieure de l'intégrale
    D = 1 'Borne supérieure de l'intégrale
 
    'Message si a est négatif
    If a < 0 Then
        Message1 = MsgBox("La valeur de a est négative.", vbRetryCancel + vbCritical, "Problématique")
        If Message1 = vbCancel Then Exit Function
        If Message1 = vbRetry Then
            Message2 = MsgBox("Veuillez réessayer votre calcul avec un a positif.", vbOKOnly + vbInformation, "Problématique")
            If Message2 = vbOK Then Exit Function
        End If
    End If
 
    'Message si b est négatif
    If b < 0 Then
        Message3 = MsgBox("La valeur de b est négative.", vbRetryCancel + vbCritical, "Problématique")
        If Message3 = vbCancel Then Exit Function
        If Message3 = vbRetry Then
            Message4 = MsgBox("Veuillez réessayer votre calcul avec un b positif.", vbOKOnly + vbInformation, "Problématique")
            If Message4 = vbOK Then Exit Function
        End If
    End If
 
    'Message si a est nul
    If a = 0 Then
        Message5 = MsgBox("La valeur de a est nulle.", vbRetryCancel + vbCritical, "Problématique")
        If Message5 = vbCancel Then Exit Function
        If Message5 = vbRetry Then
            Message6 = MsgBox("Veuillez réessayer votre calcul avec un a positif.", vbOKOnly + vbInformation, "Problématique")
            If Message6 = vbOK Then Exit Function
        End If
    End If
 
    'Message si b est nul
    If b = 0 Then
        Message7 = MsgBox("La valeur de b est nulle.", vbRetryCancel + vbCritical, "Problématique")
        If Message7 = vbCancel Then Exit Function
        If Message7 = vbRetry Then
            Message8 = MsgBox("Veuillez réessayer votre calcul avec un b positif.", vbOKOnly + vbInformation, "Problématique")
            If Message8 = vbOK Then Exit Function
        End If
    End If
 
    'Longueur des n sous-intervalles de l'intervalle [C,D]
    delta_x = (D - C) / n 'n est le nombre de sous-intervalles
 
    'Déclaration d'une nouvelle variable
    Quotient = (D - C) / 2
 
    For j = 1 To n
        'Borne inférieure des n sous-intervalles
        c_i = C + (j - 1) * delta_x
 
        'Déclaration des points et des poids qui sont déterminés grâce aux polynômes de Legendre
        'Référence : http://en.wikipedia.org/wiki/Gaussian_quadrature
        x1 = c_i + 0
        w1 = 2
        x2 = c_i - ((3 ^ 0.5) / 3)
        w2 = 1
        x3 = c_i + ((3 ^ 0.5) / 3)
        w3 = 1
        x4 = c_i + 0
        w4 = 8 / 9
        x5 = c_i - ((3 / 5) ^ 0.5)
        w5 = 5 / 9
        x6 = c_i + ((3 / 5) ^ 0.5)
        w6 = 5 / 9
        x7 = c_i + (((3 - (2 * (6 / 5) ^ 0.5)) / 7) ^ 0.5)
        w7 = (18 + (30 ^ 0.5)) / 36
        x8 = c_i - (((3 - (2 * (6 / 5) ^ 0.5)) / 7) ^ 0.5)
        w8 = (18 + (30 ^ 0.5)) / 36
        x9 = c_i + (((3 + (2 * (6 / 5) ^ 0.5)) / 7) ^ 0.5)
        w9 = (18 - (30 ^ 0.5)) / 36
        x10 = c_i - (((3 + (2 * (6 / 5) ^ 0.5)) / 7) ^ 0.5)
        w10 = (18 - (30 ^ 0.5)) / 36
        x11 = c_i + 0
        w11 = 128 / 225
        x12 = c_i + (((5 - (2 * (10 / 7) ^ 0.5)) ^ 0.5) / 3)
        w12 = (322 + (13 * (70 ^ 0.5))) / 900
        x13 = c_i - (((5 - (2 * (10 / 7) ^ 0.5)) ^ 0.5) / 3)
        w13 = (322 + (13 * (70 ^ 0.5))) / 900
        x14 = c_i + (((5 + (2 * (10 / 7) ^ 0.5)) ^ 0.5) / 3)
        w14 = (322 - (13 * (70 ^ 0.5))) / 900
        x15 = c_i - (((5 + (2 * (10 / 7) ^ 0.5)) ^ 0.5) / 3)
        w15 = (322 - (13 * (70 ^ 0.5))) / 900
 
        If m = 1 Then 'm est le nombre de points de la quadrature gaussienne
            'Somme du produit des poids et de f
            Somme1 = w1 * f(x1, a, b)
 
        ElseIf m = 2 Then
            'Somme du produit des poids et de f
            Somme2 = (w2 * f(x2, a, b)) + (w3 * f(x3, a, b))
 
        ElseIf m = 3 Then
            'Somme du produit des poids et de f
            Somme3 = (w4 * f(x4, a, b)) + (w5 * f(x5, a, b)) + (w6 * f(x6, a, b))
 
        ElseIf m = 4 Then
            'Somme du produit des poids et de f
            Somme4 = (w7 * f(x7, a, b)) + (w8 * f(x8, a, b)) + (w9 * f(x9, a, b)) + (w10 * f(x10, a, b))
 
        ElseIf m = 5 Then
            'Somme du produit des poids et de f
            Somme5 = (w11 * f(x11, a, b)) + (w12 * f(x12, a, b)) + (w13 * f(x13, a, b)) + (w14 * f(x14, a, b)) + (w15 * f(x15, a, b))
 
        End If
 
        Somme = Somme + Somme1 + Somme2 + Somme3 + Somme4 + Somme5
 
    Next
    'Approximation de la fonction bêta avec la quadrature gaussienne
    Beta = Somme * Quotient
 
'Fin de la procédure
End Function

**ProgElecT** · 21/03/2013, 13h00

Salut

Cause de l'erreur.

result = NumbeR^ExponenT
le signe ^, Permet d'élever un nombre à une puissance.
L'argument NumbeR ne peut être négatif que si l'argument ExponenT est un entier

Appliqué à ton code:
dans cette partie du calcul cela donne: result = ((1 - z) ^ (b - 1)) il arrive que
(1 - z) = 1 - 1.0886... = -8.8675...E-02 donc NumbeR est négatif
(b - 1) = 4.3 - 1 = 3.3 donc ExponenT n'est pas un entier

Donc, soit il faut modifier l'algorithme dans la fonction f, soit ce cas de figure ne devrait pas arrivé, donc il y a une erreur ailleurs avant l'appel à la fonction f.

Tu pourrais éventuellement te faire aider dans le forum Algorithmes

**aurama** · 24/03/2013, 02h33

Merci ProgElecT!

Grâce à tes explications, j'ai pu régler mon problème.

Erreur d'exécution 5

Discussions diverses VB

Discussions similaires

Partager

Partager