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

[Stavros M. <mac...@gm...>]

L'enjeu est plutôt d'assurer que les matheux soient conscients que Maxima
reste un outil pratique pour tout type de calcul.

Honnêtement, je suis un peu perplexe que OEIS estime cette série digne
d'inclusion.... Mais là c'est un autre débat.

         -s

On Tue, Nov 15, 2022 at 12:48 PM Thomas Baruchel <bar...@gm...> wrote:

> 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-discuss mailing list
> Max...@li...
> https://lists.sourceforge.net/lists/listinfo/maxima-discuss
>