Re: [prog] Précision numérique

[xavier roguiez]

Désolé pour le premier mail , mais je me suis rendu compte que les valeurs approchées proposées en copier coller donnait effectivement le résultat affiché.

Voici une vérification avec mapple :
bw:=0.25; d:=0.917; dp:=0.067;A:=2945.*10^(-6);  Ap:=490.*10^(-6);
> Es:=200000.;  Ecm:=22000.*(28./10.)^0.3;alpha:=Es/Ecm*(1.+2.);
                             bw := .25
                             d := .917
                             dp := .067
                         A := .002945000000
                        Ap := .0004900000000
                           Es := 200000.
                         Ecm := 29961.95105
                        alpha := 20.02539818

> solve(bw*y^2/2+alpha*(A+Ap)*y-alpha*(d*A+dp*Ap),y);

                     .4415138105, -.9918117525
Le résultat est donc 0.441538105

Voici une vérification avec Scilab

2.);

-->//solve(bw*y^2/2+alpha*(A+Ap)*y-alpha*(d*A+dp*Ap),y)

-->a=bw/2.;  b=alpha*(A+Ap); c=-alpha*(d*A+dp*Ap);

-->Delta=b*b-4*a*c;

-->(-b+Delta^0.5)/(2*a)
ans  =

   0.4415138


D'avance merci
xavier

Maintenant je cherche mon erreur, je remets pas en cause OO, mais plutôt mon approche. J'ai rajouté des 28. sur mes chiffres mais rien n'y change -->bw=0.25; d=0.917; dp=0.067; A=2945.*10^(-6); Ap=490.*10^(-6); Es=200000.; Ecm=22000.*(28./10.)^0.3; alpha=Es/Ecm*(1.+ Sinon ci joint le document open office avec deux versions sur la même feuille : la version basic (macro fonction) et la version tableur, mais j'ai le même résultat. Cela vient certainement de moi, mais j'avoue atteint mes limites en terme de recherche d'erreur. C'est pourquoi je me tourne vers la liste.

Mais la résolution avec d'autres outils numériques : C, mapple ou autre donne la vraie solution de x1 qui est x1 := .4415138104

Bonjour,


   Bernard

---------------------------------------------------------------------
To unsubscribe, e-mail: prog-unsubscr...@fr.openoffice.org
For additional commands, e-mail: prog-h...@fr.openoffice.org

Vérifie plutôt tes calculs avec tes autres outils. Peut-être une erreur dans les données entrées. Avec ce que tu donnes, j'ai refait les calculs (péniblement) avec la calculette scientifique de Windows, qui est plus précise. Je retrouve les valeurs du Delta et du résultat : 0,43820369. Le calcul de la racine carrée de Delta par la fonction Basic SQR() est aussi correct: je retrouve le delta avec une multiplication par lui-même.

test.ods
Description: application/vnd.oasis.opendocument.spreadsh

---------------------------------------------------------------------
To unsubscribe, e-mail: prog-unsubscr...@fr.openoffice.org
For additional commands, e-mail: prog-h...@fr.openoffice.org