Re: Numerical calculation of a double sum (Appell's function F4)

• To: mathgroup at smc.vnet.net
• Subject: [mg75352] Re: Numerical calculation of a double sum (Appell's function F4)
• From: dh <dh at metrohm.ch>
• Date: Thu, 26 Apr 2007 03:31:41 -0400 (EDT)
• References: <f0n8j5\$s6s\$1@smc.vnet.net>

```
Hi Markus,

Are you aware that Sum can do double sums?

Daniel

Markus Huber wrote:

> Hello,

>

> I am working currently with double hypergeometric sums. Unfortunately

> the one I need (Appell's F4) is not implemented in M. So I decided to

> write my own F4.

> I discovered quickly that this a far more complicated task than

> expected, which is due to many special cases of the arguments (e.g.

> the series can be truncated, transform into another series or have

> singularities in Gamma functions; there are also cases where these

> singularities cancel). I managed to overcome all those difficulties,

> but still I wonder, if there is an easier solution to the following:

>

> F4 is a double sum. My approach is a simple While loop (ok, it's more

> complicated, because cancelling singularities can lead to intermediate

> terms that are zero and so you have t think carefully about getting a

> quantity that tells you when the desired accuracy is reached). Is

> there a function in M that can do double sums numerically and also

> checks for convergence?

>

> I appreciate any help, because this - I would guess - would make the

> calculation faster and the code more readable.

>

> Thanks for any suggestions

>

>

> Markus Huber

>

>

> PS: Don't be misled by the seeming simplicity of Appell's Function F4

> when you look it up. In all probability you will only find the

> standard series representation. There are also other representations

> with other regions of convergence that exist of 5 single

> hypergeometric series. Esp. those I need.

>

>

```

• Prev by Date: Re: minmum of a function
• Next by Date: MeijerG function
• Previous by thread: Re: Numerical calculation of a double sum (Appell's function F4)
• Next by thread: "hard" simplification