Bonjour à tous,
Le problème a pourtant l'air simple, mais j'ai des résultats surprenants...
Mon polynome est le suivant :
où les "e" signifie "x 10^" : c'est du c++...
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6 ax^2 + bx + c : a = 4 - 1e12; b = 1e12 * 0.0101; c = -1e6; delta = b*b-4*a*c; sol = (-b+sqrt(delta))/2/a;
En appliquant la formule du delta, j'obtiens pour racine :
0.00009999999999999564087843
et le polynome évalué vaut :
-0.00000000314321368932723999
Si j'applique un brave algo de dichotomie, j'obtiens pour racine :
0.00009999999999999595258655
et le polynome évalué vaut :
0.00000000000000000000000000
Une méthode de Newton converge comme la dichotomie. Apparemment, le calcul "analytique" par delta fonctionne moins bien que les méthodes d'approximation numérique : est-ce normal ? Ai-je fait une bêtise ? Toute suggestion serait la bienvenue !!
Merci !
Marc
Partager