Discussion :

[proimane]

salut tous le monde comment écrire racine cubique dans fortron car dans mon programme j'ai rencontre quelque problème tel que:
j'ai:
dp=-23.
k0=12d-1
l=100.
f=dp/(k0*l)
g=f**(1./3.)
quand j'exécute il affiche NaN pour g mais il affiche pour f = -0.916666666666667 c'est une vrai valeur
et Merci

[Sylvain Bergeron]

Le problème vient du fait que dans a ** b, a ne peut pas être négatif.

Il est facile d'admettre que pour des exposants entiers, ça va bien. Mais pour des exposants fractionnaires, l'évaluation devient problématique. Si mathématiquement, il est facile de résoudre certains cas triviaux, par exemple (-1)^(1/3), que vaut (-1)^(1/3+Epsilon) ? Sachant que 1/3 n'est pas représentable directement par une variable réel, ton calcul est équivalent à (-1)^(1/3+-Epsilon). Il est donc normal que la réponse soit NaN (Not a number).