Re: Doublefactorial and Pochhammer

From	Joe Riel <joer@san.rr.com>
Newsgroups	comp.soft-sys.math.maple
Subject	Re: Doublefactorial and Pochhammer
Date	2015-11-09 15:44 -0800
Organization	A noiseless patient Spider
Message-ID	<87bnb2u4an.fsf@san.rr.com> (permalink)
References	<18156841-3104-434f-8cdc-1999bb794b97@googlegroups.com>

Show all headers | View raw

Peter Luschny <peter.luschny@gmail.com> writes:

> a := n -> doublefactorial(2*n-1)-2^n*pochhammer(1/2,n):
> seq(a(n), n=0..99);
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
>
> a(100)
> 6666308670072953744411215006735034163324489389674388736363184954745922258576896518414625915283128424390474317708176893511841954015267176587405666801912441638268530971962511539459228515625-1267650600228229401496703205376*GAMMA(201/2)/sqrt(Pi)
>
> Can somebody explain this sudden jump?
>

GAMMA has a preset limit at which it will apply
a particular simplification to an integer input.
You can extend this by modifying a global constant:

`GAMMA/magic` := 2*`GAMMA/magic`:

a(100);
             0

The constant actually is an expression sequence of four integers.

-- 
Joe Riel

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

Thread

Doublefactorial and Pochhammer Peter Luschny <peter.luschny@gmail.com> - 2015-11-05 15:00 -0800
  Re: Doublefactorial and Pochhammer jfh <john.harper@vuw.ac.nz> - 2015-11-05 16:05 -0800
    Re: Doublefactorial and Pochhammer Peter Luschny <peter.luschny@gmail.com> - 2015-11-06 00:40 -0800
  Re: Doublefactorial and Pochhammer Joe Riel <joer@san.rr.com> - 2015-11-09 15:44 -0800

csiph-web