이 계산기의 기능
이 도구는 제1종 구면 베셀 함수 \(j_v(x)\), 제2종 구면 베셀 함수 \(y_v(x)\), 그리고 각각의 1차 도함수 \(j'_v(x)\)와 \(y'_v(x)\)를 계산합니다. 이 함수들은 구면 좌표계에서 파동 방정식과 헬름홀츠 방정식의 동경(radial) 방향 해이며, 산란 이론, 전자기파·음향 방사, 양자역학(자유 입자의 부분파) 등 물리학 전반에서 등장합니다. 새롭게 구현된 이 버전은 정수 차수 \(v \ge 0\)과 실수 인수 \(x > 0\)을 처리합니다.
사용 방법
차수 \(v\)(0, 1, 2와 같은 음이 아닌 정수)와 인수 \(x\)(양의 실수)를 입력하세요. 계산 버튼을 누르면 네 가지 값을 모두 얻을 수 있습니다. 다만 \(y_v(x)\)와 \(y'_v(x)\)는 \(x\)가 0에 가까워질수록 발산하므로, \(x = 0\)에서는 무한대로 표시됩니다. \(j_0(0)\)은 극한값으로 1이 됩니다.
공식 설명
이 함수들은 미분방정식 \(x^2 w'' + 2x w' + (x^2 - v(v+1))w = 0\)을 만족합니다. 닫힌 형태인 \(j_0 = \sin(x)/x\)와 \(y_0 = -\cos(x)/x\)에서 출발하여, 더 높은 차수는 삼항 점화식 $$f_{n+1} = \frac{2n+1}{x}\,f_n - f_{n-1}$$을 따릅니다. 상향 점화는 \(y_v\)에 대해서는 안정적이지만, \(j_v\)의 경우 \(n > x\)일 때 불안정합니다. 그래서 밀러(Miller)의 하향 점화법을 사용합니다. 즉 충분히 높은 차수에서 \(f\)를 각각 0과 1로 설정해 시작한 뒤 차수를 낮춰 가며 점화시키고, 0차 항이 \(\sin(x)/x\)와 일치하도록 모든 값을 다시 스케일링합니다. 도함수는 \(j'_v = j_{v-1} - \frac{v+1}{x}j_v\) 공식을 사용합니다.
계산 예시 (v = 0, x = 2)
$$j_0(2) = \frac{\sin(2)}{2} = 0.4546487134, \qquad y_0(2) = -\frac{\cos(2)}{2} = 0.2080734183$$입니다. \(j'_0 = -j_1\)이므로 \(j'_0(2) = -0.4353977750\)이 되고, \(y'_0 = -y_1\)이므로 \(y'_0(2) = 0.3506120043\)이 됩니다.
자주 묻는 질문
복소수 \(x\)도 지원하나요? 아니요. 원래 페이지는 복소수 인수를 받지만, 이 버전은 명확성과 속도를 위해 실수 \(x > 0\)으로 제한했습니다.
왜 \(y_v\)는 \(x = 0\)에서 무한대인가요? 제2종 함수는 원점에 극(pole)을 가지므로, \(x\)가 0에 가까워질수록 값이 한없이 커집니다.
정확도는 어느 정도인가요? 배정밀도(double precision)로 계산하여 약 15자리의 유효숫자를 제공하므로, 일반적인 공학·물리 작업에는 충분합니다.
