Karim Belabas on Sun, 07 Apr 2019 15:58:46 +0200
|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
- To: Dirk Laurie <dirk.laurie@gmail.com>
- Subject: Re: Padé approximants?
- From: Karim Belabas <Karim.Belabas@math.u-bordeaux.fr>
- Date: Sun, 7 Apr 2019 15:58:42 +0200
- Cc: pari-users@pari.math.u-bordeaux.fr
- Delivery-date: Sun, 07 Apr 2019 15:58:46 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=simple/simple; d=math.u-bordeaux.fr; s=mail; t=1554645524; bh=JpVM2TeUp1ngqJD8e+gVTozHZ28cNuOfngnFsElH7Rg=; h=Date:From:To:Cc:Subject:References:In-Reply-To:From; b=0+3kVZbmuB1Pczij/PcA1Aivr2J+bEohVxS/+Zi0VwAAE/0FSRfHSoWFaLca68qUe JtIR2qFXZSWL2y0SpH/GUXiOY9WAtfuTh6RuyVa/vVsZkTYRaw6E3LAGppcYk/a8PY egDjnb+jLtULjD6aLz+W6XD0+Pf/qhV0w7mFPv7s=
- In-reply-to: <CABcj=tmUx5MadzryURxMTw3amN0FwSV5hw1bi=EMGoe-dfT2fw@mail.gmail.com>
- Mail-followup-to: Dirk Laurie <dirk.laurie@gmail.com>, pari-users@pari.math.u-bordeaux.fr
- References: <CABcj=tmUx5MadzryURxMTw3amN0FwSV5hw1bi=EMGoe-dfT2fw@mail.gmail.com>
- User-agent: NeoMutt/20180716
* Dirk Laurie [2019-04-07 14:46]:
> I can't find a built-in function that computes Padé approximants of a
> power series.
>
> e.g. pade(1 + x + 1/2*x^2 + 1/6*x^3 + 1/24*x^4 + O(x^5)) → (x^2 + 6*x
> + 12)/(x^2 - 6*x + 12)
>
> It is of course quite easy to program for oneself in the nondegenerate
> case (epsilon algorithm) but I wonder if it may be disguised as a
> special case of a more general function.
See ??bestapprPade
Cheers,
K.B.
--
Karim Belabas, IMB (UMR 5251) Tel: (+33) (0)5 40 00 26 17
Universite de Bordeaux Fax: (+33) (0)5 40 00 21 23
351, cours de la Liberation http://www.math.u-bordeaux.fr/~kbelabas/
F-33405 Talence (France) http://pari.math.u-bordeaux.fr/ [PARI/GP]
`