utilisation fonction integration

[hermes]

Bonsoir,

j'ai éssayé de faire retrouver une borne d'une integrale à la TI89 cependant les resultats qu'elle me renvois sont aberrant notament quand je lui demande l'equivalent d'une lecture inverse de la loi normale centré reduite soit lorsque je demande le resultat de
integral·(1/(2*pi)^(1/2)*e^(-x^2/2),x,-inf,1.34)
c'est egal à 0.9909
Alors que si je lui demande la borne correspondante a 0.9909 grace à une fonction solve la calculette me renvoit 2.36 via la fonction que j'ai mis en PJ alors que la même chose dans l'ecran home me renvois le bon resultat.

Cdt

Hermes

[diot]

bizarre en effet. Le bug peuvent exister sur les calculatrices. Peut-être que c'en est pas un mais je vais avoir du mal à t'aider, car je n'utilise pas cette calculatrice.

[Bisam]

Mince, je voulais regarder et j'ai pas allumé l'ordi... J'y vais de ce pas.

(30 minutes passent...)

Après quelques essais, ma calculette ne fournit jamais de résultat... elle tourne sans s'arrêter pendant plus de 5 minutes sans fournir de réponse et ensuite, je "Break" par pitié pour ses piles.

Il serait bien plus efficace de construire une table de la fonction et de la compléter par interpolation.

Pour commencer, on peut aussi considérablement réduire le temps de calcul en faisant calculer l'intégrale de 0 à t et en enlevant 1/2 à la valeur cherchée...

Code: Select all

if c>0.5 then
exp>list(solve(1/racine(2*pi)* int(exp(-x^2/2),x,0,t)=c-0.5,t),t)
else
-exp>list(solve(1/racine(2*pi)* int(exp(-x^2/2),x,0,t)=0.5-c,t),t)
endif


ou encore en simplifiant :

Code: Select all

sign(c-0.5)*exp>list(solve(1/racine(2*pi)* int(exp(-x^2/2),x,0,t)=abs(c-0.5),t),t)

[hermes]

Bonjour,

Merci pour l'astuce pour que mon programme aille plus vlite (il aurait suffit que je comprenne un peu mon cours et j'aurais pu y pensé tous seul mais bon les proba c'est pas ma tasse de thé) cependant le problème reste le même, c'est a dire en fonction si le calcule est entré dans HOME ou dans un programme le resultat n'est pas le même. Bizzar non?

Cdt

Hermès

[Bisam]

Je n'ai pas du tout constaté ce problème en faisant mes tests.
Es-tu certain d'avoir bien tapé la même chose dans les deux cas ?
Peut-être y a-t-il une variable qui n'a pas la même valeur...