xが小さい場合のクンマー関数U(a,b,x)の数値計算

• Item type: Journal(1)
• 公開日: 2021-04-15
• タイトル: xが小さい場合のクンマー関数U(a,b,x)の数値計算
• タイトル / 言語: en
• タイトル: Computation of Kummer Function U(a,b,x) for Small Argument x
• 言語: jpn
• キーワード / 主題Scheme: Other
• キーワード / 主題: [一般論文] クンマー関数，合流型超幾何関数，ガンマ関数の逆数のテイラー展開
• 資源タイプ / 資源タイプ識別子: http://purl.org/coar/resource_type/c_6501
• 資源タイプ: journal article
• ID登録: 10.20729/00210557
• ID登録 / ID登録タイプ: JaLC
• 著者所属: 中部大学
• 著者所属: 中部大学
• 著者所属: 中部大学
• 著者所属(英): en
• 著者所属(英): Chubu University
• 著者所属(英): en
• 著者所属(英): Chubu University
• 著者所属(英): en
• 著者所属(英): Chubu University
• 著者名: 吉田, 年雄; 奥居, 哲; 足達, 義則
• 著者名(英): Toshio, Yoshida; Satoshi, Okui; Yoshinori, Adachi
• 論文抄録 / 内容記述タイプ: Other
• 論文抄録 / 内容記述: $a>0$，$b\ge 0$で，$x$（$>0$）が小さい場合のクンマー関数$U(a,b,x)$の新しい数値計算法を述べている．$U(a,b,x)$は$U(a,b,x)=(\pi/\sin \pi b)\{M(a,b,x)/(\Gamma(a+1-b)\Gamma(b))-x^{1-b}M(a+1-b,2-b,x)/(\Gamma(a)\Gamma(2-b))\}$で定義される．上式の$\{\ \}$の前半部と後半部の減算では，$b$が整数に近くないときには，桁落ちは起こらないので，そのまま計算すればよいが，$b$が整数に近いときには，桁落ちが起こり，工夫が必要である．たとえば，$0\le b\le 0.5$の場合には，上式の$\{\ \}$の部分は$b/(\Gamma(1+b)\Gamma(a+1-b))$+$\sum_{k=1}^{\infty}\{A_{k}(a,b)-B_{k}(a,b)\}x^{k}$の形で表わすことができる．ここで，$A_{k}(a,b)$は2つの項の差で表され，特に$A_{1}(a,b)=a/(\Gamma(1+b)\Gamma(a+1-b))-1/(\Gamma(2-b)\Gamma(a))$である．この式は，$b$が0に近いときには桁落ちが生ずるが，ガンマ関数の逆数$1/\Gamma(1+b)$，$1/\Gamma(2-b)$，$1/\Gamma(a+1-b)$の$b=0$周りのテイラー展開を用いれば，桁落ちを回避することができる．そのため，本論文では，プログラム実行時に，$1/\Gamma(1+b)$や$1/\Gamma(2-b)$，さらに，変数$a$を含む$1/\Gamma(a+1-b)$の場合でも，必要な次数までの微分係数を能率的に求められる方法を提案している．$A_{k}(a,b)\ (k\ge 2)$は$A_{k-1}(a,b)$から漸化式を用いて計算できる．このようにして，桁落ちを回避することにより$U(a,b,x)$を求めている．
• 論文抄録(英) / 内容記述タイプ: Other
• 論文抄録(英) / 内容記述: Numerical method for computation of Kummer functions $U(a,b,x)$ with small argument $x$ ($>0$) $a>0$, $b\ge 0$）is proposed. The function $U(a,b,x)$ is defined as $U(a,b,x)=(\pi/\sin \pi b)\{M(a,b,x)/(\Gamma(a+1-b)\Gamma(b))-x^{1-b}M(a+1-b,2-b,x)/(\Gamma(a)\Gamma(2-b))\}$. When $b$ is not near to an integer, the cancellation of significant digits does not occur by subtraction of these two parts and so we can compute the above expression as it is. However, when $b$ is near to an integer, such cancellation occurs. For example, in case of $0\le b\le 0.5$, the part $\{\ \}$ of the above equation can be expressed as $b/(\Gamma(b+1)\Gamma(a+1-b))$+$\sum_{k=1}^{\infty}\{A_{k}(a,b)-B_{k}(a,b)\}x^{k}$, where $A_{k}(a,b)$ is a subtraction of two terms, especially $A_{1}(a,b)=a/(\Gamma(1+b)\Gamma(a+1-b))-1/(\Gamma(2-b)\Gamma(a))$. When $b$ is near to 0, the cancellation of significant digits occurs in this equation. By using Taylor series about $b=0$ for reciprocal gamma functions, $1/\Gamma(1+b)$, $1/\Gamma(2-b)$ and $1/\Gamma(a+1-b)$, such cancellation can be removed. In this paper, we describe an efficient method for the runtime computation of derivatives as much as necessary in case of $1/\Gamma(1+b)$, $1/\Gamma(2-b)$ and $1/\Gamma(a+1-b)$ which contains a variable $a$. By the recurrence, $A_{k}(a,b)\ (k\ge 2)$ is computed from $A_{k-1}(a,b)$. Removing the cancellation of significant digits in this way, we can compute $U(a,b,x)$.
• 書誌レコードID / 収録物識別子タイプ: NCID
• 書誌レコードID / 収録物識別子: AN00116647
• 書誌情報: 情報処理学会論文誌; 巻 62, 号 4, p. 1056-1064, 発行日 2021-04-15
• ISSN / 収録物識別子タイプ: ISSN
• ISSN / 収録物識別子: 1882-7764