Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]

Groups > comp.soft-sys.math.maple > #1199

Re: Convert to factorial, bug?

Show all headers | View raw

On 01.11.2015 19:40, Peter Luschny wrote:
> A:=(n,k)->(n!/k!)*binomial(k,n-k)*2^(2*k-n);
> convert(A(n,k), factorial);
> B:=(n,k)->n!*2^(2*k-n)/((n-k)!*(2*k-n)!);
>
> A(1,0); # gives 0
> B(1,0); # Error, div by zero
>
> Peter
>

I guess that Maple uses 1 diveded by GAMMA and not the function 1/GAMMA,
which is holomorphic and zero in non-negative integers.

B(n,0) = 2^(-n)/(-n)! = 2^(-n)/GAMMA(-n+1)

To work around that weakness here one can use the limit to get the result:

'Limit(B(n,0), n=1)';  ``= %;  #convert(%, GAMMA);
value(%); # gives 0

Back to comp.soft-sys.math.maple | Previous | Next — Previous in thread | Next in thread | Find similar

Thread

Convert to factorial, bug? Peter Luschny <peter.luschny@gmail.com> - 2015-11-01 10:40 -0800
  Re: Convert to factorial, bug? Axel Vogt <&noreply@axelvogt.de> - 2015-11-02 08:44 +0100
    Re: Convert to factorial, bug? Peter Luschny <peter.luschny@gmail.com> - 2015-11-02 00:15 -0800

csiph-web