Bonjour,
je fais appel à vous car j'ai un problème avec le calcul d'une intégrale à savoir : intégrale de (-1) à +inf u^((1/x)-1)e^(-u)du.
Cela ressemble à la fonction gamma(1/x) mis à part que la borne inférieure de l'intégrale est -1.
La macro calcule l'intégrale pour un x compris entre 0 et 1.
Le calcul fonctionne très bien si je choisis x=0.5. En revanche pour x=0.4, voici l'erreur que me renvoie SAS :
J'aurais aimé savoir si l'erreur peut venir du code ou plutôt du domaine de définition de l'intégrale ?
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 NOTE: IML Ready NOTE: Module FONCTION defined. WARNING: Invalid argument resulted in missing value result. operation : ** at line 28 column 1 operands : t, _TEM1002 t 1 row 1 col (numeric) -0.264241 _TEM1002 1 row 1 col (numeric) 2.3333333 statement : ASSIGN at line 28 column 1 traceback : module FONCTION at line 28 column 1 ERROR: Invalid Operation in QUAD. operation : QUAD at line 28 column 1 operands : *LIT1006, a *LIT1006 1 row 1 col (character, size 8) fonction a 1 row 2 cols (numeric) -1 P statement : CALL at line 28 column 1 WARNING: All data set variables are unvalued. No APPEND done.
Merci beaucoup !
Partager