On Tue, 15 Nov 2022, Robert Dodier wrote:
Le record donné par G. C. Greubel est de 10 000 digits.
J’en donne 150 parce que MAPLE en donne 120 (comme
Mathematica)
/* OEIS A173272 Decimal expansion of the positive
solution of sqrt((2-x)(2+x)) + sqrt((3-x)(a+x)) =
sqrt((2-x)(2+x))*sqrt((3-x)(3+x)). */
/* After MAPLE and MATHEMATICA examples, Patrick
Guillemin is adding a MAXIMA example with 150 Digits */
find_root_abs:1/10^150$ find_root_rel:1/10^150$
fpprec:150$ fpprintprec:150$
x0:bf_find_root(x^8-22*x^6+163*x^4-454*x^2+385, x,
1.1b0, 1.3b0);
/* end of code */
Bonjour Patrick,
le code semble correct, mais si l'enjeu est la course aux
nombres, peut-être vaut-il mieux utiliser Pari-GP, plus
adapté, qui calcule de façon instantané plusieurs dizaines
de milliers de chiffres :
(18:46) gp > \p200
realprecision = 211 significant digits (200 digits displayed)
(18:47) gp > solve(x=1.1, 1.3, x^8-22*x^6+163*x^4-454*x^2+385)
%2 =
%1.2311857237786688299627058347697888745686490269976349243438469028632788354636825802070220761365423157787386759254111932030711775173725696924568433951120227051506561962475322075369285896297975417856324
Ici, j'ai tapé \p200 (200 chiffres) pour ne pas alourdir le message,
mais on peut sans hésiter tenter \p20000 pour passer à 20000 chiffres.
Cordialement,
--
Thomas Baruchel
_______________________________________________
Maxima-lang-fr mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/maxima-lang-fr
