Re: [Maxima-lang-fr] [Maxima-discuss] D'après vous, est-ce-correct comme code Maxima pour donner un exemple concernant oeis.org/A173272 ?

[Thomas B. <bar...@gm...>]

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