이 계산기의 기능
이 도구는 n점 일반화 가우스-라게르 구적법(generalized Gauss-Laguerre quadrature)의 노드(가로좌표) \(x_i\)와 가중치 \(w_i\)를 계산합니다. 국가나 지역에 관계없이 동일하게 적용되는 순수 수학 기반의 수치적분 도구입니다. 이 구적 규칙은 가중함수 \(x^{\alpha} e^{-x}\)를 포함하는 반무한 구간 [0, ∞) 위의 적분을 근사합니다.
$$\int_{0}^{\infty} x^{\alpha} e^{-x} f(x)\, dx \approx \sum_{i=1}^{n} w_i\, f(x_i).$$노드는 일반화 라게르 다항식 \(L_n^{(\alpha)}(x)\)의 양의 근이며, f가 차수 \(2n-1\) 이하의 다항식일 때 이 규칙은 정확한 값을 줍니다.
사용 방법
차수 n(점의 개수, 2부터 100까지)을 선택하고, 지수 매개변수 α(−1보다 큰 임의의 실수; 고전적인 가우스-라게르 규칙은 \(\alpha = 0\)을 사용)를 입력한 뒤, 표시할 유효 자릿수를 지정하세요. 결과에는 각 노드와 그에 대응하는 가중치가 \(x_i\)의 오름차순으로 정리되어 나타나며, 자체 검산 항목도 함께 제공됩니다.
공식과 계산 방법
각 가중치는 닫힌 형태의 공식 $$w_i = \frac{\Gamma(n+\alpha+1)\cdot x_i}{n!\cdot\left[(n+1)L_{n+1}^{(\alpha)}(x_i)\right]^2}$$을 만족합니다. 내부적으로는 이와 동등하면서 수치적으로 안정적인 골루브-벨치(Golub-Welsch) 방법을 사용합니다. 즉, 대각 성분이 \(a_k = 2k+\alpha+1\)이고 비대각 성분이 \(b_k = \sqrt{k(k+\alpha)}\)인 대칭 삼중대각 야코비 행렬을 구성합니다. 이 행렬의 고윳값이 노드가 되고, 각 가중치는 \(\mu_0\cdot(\text{고유벡터의 첫 성분})^2\)과 같으며, 여기서 \(\mu_0 = \Gamma(\alpha+1)\)은 0차 모멘트입니다. 이 방식은 큰 팩토리얼에서 발생하는 오버플로를 방지해 줍니다.
계산 예시
\(n = 2\), \(\alpha = 0\)인 경우: \(L_2^{(0)}(x) = (x^2-4x+2)/2\)이므로 근은 \(x = 2 \pm \sqrt{2}\)이고, 따라서 \(x_1 = 0.5857864\), \(x_2 = 3.4142136\)이 됩니다. 가중치는 \(w_1 = (2+\sqrt{2})/4 = 0.8535534\), \(w_2 = (2-\sqrt{2})/4 = 0.1464466\)입니다. 두 가중치의 합은 \(1 = \Gamma(1)\)로, 결과가 올바름을 확인할 수 있습니다.
자주 묻는 질문
α는 어떤 역할을 하나요? α는 가중치 \(x^{\alpha}\)를 결정합니다. \(\alpha = 0\)이면 표준 가우스-라게르 규칙이 되고, \(\alpha > 0\)이면 원점에서 멀어진 영역에 더 큰 가중치를 둡니다. 반드시 −1보다 커야 합니다.
정확도는 어느 정도인가요? n점 규칙은 차수 \(2n-1\) 이하의 다항식을 정확히 적분하며, 매끄러운 함수에 대해서는 빠르게 수렴합니다.
결과를 어떻게 검산하나요? 모든 가중치의 합은 항상 \(\Gamma(\alpha+1)\)과 같으며, 이는 0차 모멘트 행에 표시됩니다.
