partDéc() qui renvoie 1

[iequa]

Je ne comprends pas bien, mettre K^(1/3) dans une variable ? C'est ce que j'ai fais dès le début (voir le code dans le premier message), sans succès

[Bisam]

C'est malheureusement un "bug" impossible à supprimer.
En effet, cela vient directement de la façon de stocker des réels.

Ici, après le calcul (approché) de 64^(1/3), c'est un nombre très proche de 4, mais légèrement inférieur, qui est obtenu. (En fait, ce sera le cas pour toutes les valeurs suivantes : (n^3)^(1/3) renvoie toujours une valeur très légèrement plus petite que n...)
Sa partie fractionnaire vaut donc 0.9999999999999.. avec plein de 9, mais la calculatrice n'affiche pas tous les chiffres qu'elle a en mémoire, donc elle affiche une valeur proche, à savoir 1.

C'est même plus qu'un problème d'affichage puisque l'on ne peut pas faire des calculs avec ces nombres impossibles à afficher :

Code: Select all

1-2^(-40)-1

renvoie 0
Pourtant

Code: Select all

fPart(1-2^-40)

renvoie 1 et prouve que le résultat en mémoire est différent du nombre 1.
et

Code: Select all

1-2^(-40)-1/2-1/2

ne renvoie pas 0... mais bien 2^-40 !

On peut en déduire que les calculettes z80 effectuent les calculs sur 40 bits de mantisse.

Enfin, pour résoudre le problème, il suffit de ne pas tester l'égalité à 0 ! On peut le contourner ainsi :

Code: Select all

Prompt K
K^(1/3)->R
If abs(partDéc(R)-partDéc(R)^2)<E-13
Then
Disp R
End

On teste en fait si la partie décimale est très proche de 0 ou bien très proche de 1...

[Adriweb]

Jolie analyse