Calcul d'intégrale IML

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 :
Code:

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.
J'aurais aimé savoir si l'erreur peut venir du code ou plutôt du domaine de définition de l'intégrale ?

Merci beaucoup !