11 messages in org.openoffice.fr.progRe: [prog] Précision numérique
FromSent OnAttachments
xavier roguiezDec 11, 2008 8:52 am 
xavier roguiezDec 11, 2008 10:08 am 
Laurent GodardDec 11, 2008 11:10 am 
xavier roguiezDec 11, 2008 11:19 am 
Laurent GodardDec 11, 2008 11:15 pm 
Jean-Baptiste FaureDec 11, 2008 11:52 pm 
Bernard MarcellyDec 12, 2008 4:38 am 
xavier roguiezDec 12, 2008 11:10 am.ods
Jean-Baptiste FaureDec 12, 2008 12:20 pm 
xavier roguiezDec 12, 2008 12:59 pm 
xavier roguiezDec 12, 2008 1:17 pm 
Actions with this message:
Paste this link in email or IM:
Paste this link in email or IM:
Atom feed for this thread
Paste this URL into your reader:
Subject:Re: [prog] Précision numériqueActions...
From:xavier roguiez (xr@geo.hmg.inpg.fr)
Date:Dec 12, 2008 11:10:36 am
List:org.openoffice.fr.prog
Attachments:
test.ods - 23k - 1 page

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é.

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

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 -->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.);

-->//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

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.

D'avance merci xavier

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, 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.

Bernard