ガウス-ロバート数値積分
有限区間（a,b）の積分をガウス-ロバート求積法で計算します。
\(\normalsize Gauss-Lobatto\ quadrature\\ {\large\int_{\small -1}^{\small 1}}f(x)dx\simeq{\large\frac{2}{n(n-1)}}(f(-1)+f(1))+{\large\displaystyle \sum_{\small i=2}^{n-1}}w_{i}f(x_i)\\{\large\int_a^{b}}f(x)dx\simeq{\large\frac{b-a}{2}[{\large\frac{2(f(a)+f(b))}{n(n-1)}}+\displaystyle \sum_{\small i=1}^{n}}w_{i}f({\large\frac{b-a}{2}}x_i+{\large\frac{b+a}{2}})]\\\)
被積分関数f(x)は、解析的であることと周期関数でないことを前提としています。
\(\normalsize Gaussian\ quadrature\\ \hspace{20px} {\large\int_{\small a}^{\small b}}w(x)f(x)dx\simeq{\large\displaystyle \sum_{\small i=1}^{n}}w_{i}f(x_i)\\ Gauss-Lobatto\ quadrature\\ \hspace{30px} interval(a,b):\hspace{20px} [-1,\ 1]\\ \hspace{30px} w(x):\hspace{80px} 1\\ \hspace{30px} polynomials:\hspace{10px} P'_{n-1} (x)\\ \)
\(\normalsize Gaussian\ quadrature\\ \hspace{20px} {\large\int_{\small a}^{\small b}}w(x)f(x)dx\simeq{\large\displaystyle \sum_{\small i=1}^{n}}w_{i}f(x_i)\\ Gauss-Lobatto\ quadrature\\ \hspace{30px} interval(a,b):\hspace{20px} [-1,\ 1]\\ \hspace{30px} w(x):\hspace{80px} 1\\ \hspace{30px} polynomials:\hspace{10px} P'_{n-1} (x)\\ \)
この計算についてお客様の声はまだありません。