Über die Näherung einer unendlichen Reihe

Angesichts

$$S_p=2\sum_{n=0}^p(-1)^n\frac {R^{(k+n)\alpha}y^{k+n}}{\Gamma(k)\,n!\,\theta^{k+n}(k+n)\,((k+n)\alpha+2)}$$ lassen$x=\frac{y }{\theta }R^{\alpha }$zu machen $$S_p=\frac 2{\Gamma(k) }\sum_{n=0}^p(-1)^n\frac {x^{(k+n)}}{\,n!\,(k+n)\,((k+n)\alpha+2)}$$ welche machen $$S_\infty=1-\frac{\Gamma (k,x)}{\Gamma (k)}+\frac{\Gamma \left(k+\frac{2}{\alpha },x\right)-\Gamma \left(k+\frac{2}{\alpha }\right)}{\Gamma (k)}\,x^{-2/\alpha }$$ Unter Verwendung Ihrer Nummern geben Sie dies an$S_\infty=0.89$(mit$55$nachlaufend$0$'S).

Jetzt,$S_p$schreiben in Bezug auf$\, _2F_2(.)$hypergeometrische Funktionen und wahrscheinlich ist dies die Schwierigkeit aus numerischer Sicht. Die Terme sind von Anfang an extrem schwankend (sie sind negativ für ungerade Werte von$p$und positiv für gerade Werte von$p$). Für die ersten

$$\left( \begin{array}{cc} p & S_p \\ 0 & +3.05072\times 10^{25} \\ 1 & -5.01200\times 10^{27} \\ 2 & +4.13536\times 10^{29} \\ 3 & -2.28392\times 10^{31} \\ 4 & +9.49546\times 10^{32} \\ 5 & -3.16900\times 10^{34} \\ 6 & +8.84121\times 10^{35} \\ 7 & -2.12041\times 10^{37} \\ 8 & +4.46180\times 10^{38} \\ 9 & -8.36644\times 10^{39} \\ 10 & +1.41526\times 10^{41} \end{array} \right)$$

Fahren Sie mit geraden Werten von fort$p$, werden wir feststellen, dass sie beginnen abzunehmen

$$\left( \begin{array}{cc} 10 & 1.41526\times 10^{41} \\ 20 & 4.48410\times 10^{51} \\ 30 & 9.94112\times 10^{59} \\ 40 & 8.48485\times 10^{66} \\ 50 & 6.33071\times 10^{72} \\ 60 & 6.71603\times 10^{77} \\ 70 & 1.39937\times 10^{82} \\ 80 & 7.21169\times 10^{85} \\ 90 & 1.09282\times 10^{89} \\ 100 & 5.57116\times 10^{91} \\ 110 & 1.06427\times 10^{94} \\ 120 & 8.32196\times 10^{95} \\ 130 & 2.86727\times 10^{97} \\ 140 & 4.63325\times 10^{98} \\ 150 & 3.70455\times 10^{99} \\ 160 & 1.53530\times 10^{100} \\ 170 & 3.43503\times 10^{100} \\ 180 & 4.30083\times 10^{100} \\ 190 & 3.11128\times 10^{100} \\ 200 & 1.33820\times 10^{100} \\ 300 & 1.19093\times 10^{86} \\ 400 & 2.52084\times 10^{57} \\ 500 & 6.78163\times 10^{17} \\ 600 & 0.89000 \end{array} \right)$$

Um das zu erklären, müssen wir für Ihre Zahlen rechnen$p$so dass

$$Q_{p+1}=\frac{\left(\frac{2000}{11}\right)^{p+21}}{20! (p+21) (p+22) (p+1)!} \leq \epsilon$$

Der$Q_{p+1}$Laufzeit geht durch ein Maximum bei$p=179.5$(siehe Tabelle !) und für diesen Wert$Q_{p+1}=8.61902\times 10^{100}$.

Für$\epsilon=10^{-20}$,wir brauchen$p=581$(sogar für$\epsilon=10^{-2}$,$p=544$).

Eine grobe Schätzung von$p$wird von gegeben

$$p=-\frac{\log (\epsilon )}{W\left(-\frac{11 }{2000 e}\log (\epsilon )\right)}$$ Wo$W(.)$ist die Lambert-Funktion.

Für den allgemeinsten Fall und eine bessere Näherung, indem wir Logarithmen nehmen und die Stirling-Näherung für den allerersten Term verwenden, beenden wir mit der Gleichung

$$\color{blue}{p \log (ex)- \left( p+\frac32\right) \log (p)=K} \qquad\text{where} \quad \color{blue}{K=\log \left(\frac{\sqrt{{\pi }} \,\epsilon \, x^{-(k+1)}\, \Gamma (k)}{\alpha \,\sqrt 2}\right)}$$ Seit$p$groß ist, könnten wir die Gleichung durch approximieren $$\color{blue}{\left(p+\frac{3}{2}\right) \log (e x)-\left(p+\frac{3}{2}\right) \log \left(p+\frac{3}{2}\right)=K}$$ für die die Lösung als Lambert-Funktion angegeben ist $$\color{blue}{p=-\frac{K}{W\left(-\frac{K}{e x}\right)}-\frac 32}$$ Angewendet auf den Arbeitsfall ergibt sich$p=599$für$\epsilon=10^{-20}$.

Ich denke, dass alles erklärt ist.

Bearbeiten

In Anbetracht des allgemeinen Begriffs

$$a_n=\frac {x^{(k+n)}}{\,n!\,(k+n)\,((k+n)\alpha+2)}$$ und unter Verwendung der logarithmischen Differenzierung haben wir$\frac{\partial a_n}{\partial n}=0$Wenn $$-2 \alpha (k+n)-(k+n) \psi (n+1) (\alpha (k+n)+2)+(k+n) \log (x) (\alpha (k+n)+2)-2=0$$ Vorausgesetzt, dass$n$groß ist, ist die Erweiterung der oben genannten Menge $$\alpha \log \left(\frac{x}{n}\right)n^2 + \left(2 (\alpha k+1) \log \left(\frac{x}{n}\right)-\frac{5 \alpha }{2}\right)n+\cdots$$ So$a_n$ist maximiert für$n \sim x$.

Unter Verwendung Ihrer Zahlen, einer rigorosen Maximierung von$a_n$zeigt, dass es bei passiert$n=179.5$(erinnere dich daran$x=\frac{2000}{11}=181.8$) wofür$a_n=5.24\times 10^{117}$.