Bill Allombert on Wed, 08 May 2013 22:14:56 +0200


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: silly long list of suggestions 

On Tue, May 07, 2013 at 11:42:50PM +0200, Christian Hoffmann wrote:
> 
> >? bestapprPade(sin(x),5)
> >%6 = (911*x^11 - 250734*x^9 + 29604960*x^7 - 1665195840*x^5 + 39178339200*x^3 - 250637587200*x)/(-8981280*x^4 - 2594592000*x^2 - 250637587200)
> >
> >Good!
> This gives :
> 
> ? bestapprPade(sin(x),5)
>   ***   at top-level: bestapprPade(sin(x),
>   ***                 ^--------------------
>   ***   not a function in function call
> 
> on my machine:
> 
>           GP/PARI CALCULATOR Version 2.5.3 (released)
>           i386 running darwin (x86-64/GMP-5.0.5 kernel) 64-bit version
>                     compiled: Mar 19 2013, gcc-4.7.1 (GCC)

On Pari 2.5, you have to use the bestappr() function 

? bestappr(Mod(truncate(sin(x)),x^16),11,4)
%10 = (911*x^11-250734*x^9+29604960*x^7-1665195840*x^5+39178339200*x^3-250637587200*x)/(-8981280*x^4-2594592000*x^2-250637587200)

We improved the interface in 2.6.

Cheers,
Bill.