일본 와산(和算) 학자 공식으로 원주율(π) 구하기

계산 입력

계산 공식 다케베 가타히로 1722년 (π²/9 급수) 마츠나가 요시스케 1739년 (π/3 급수)

항의 개수

표시 자릿수 6 10 14 18 22 26 30 34 38 42 46 50

공식

Show calculation steps (1)

1. Takebe Katahiro (1722)

   

   Series converging to pi^2/9; the k-th term has numerator (k!)^2 over a product of consecutive integers from 3. pi = 3 times the square root of the sum.

결과

π 근삿값

3.141592653589794

summed over 100 terms

급수 합산값	1.047197551196598
항의 개수	100
요청한 표시 자릿수	26

이 계산은 배정밀도 부동소수점을 사용하므로, 표시 자릿수를 어떻게 설정하든 정확도는 유효 자릿수 약 15~16자리에서 한계에 다다릅니다. 다케베 가타히로와 마츠나가 요시스케가 이뤄낸 역사적인 41자리·52자리 결과를 재현하려면 임의 정밀도 연산이 필요합니다.

이 계산기로 할 수 있는 것

이 도구는 일본의 전통 수학인 '와산(和算)'을 연구하던 학자들이 고안한 두 가지 역사적 급수 중 하나를 더해 원주율 π를 근사 계산합니다. \(\pi^2/9\)로 수렴하는 다케베 가타히로(竹部賢弘)의 급수(1722년)와 \(\pi/3\)으로 수렴하는 마츠나가 요시스케(松永良弼)의 급수(1739년) 가운데 원하는 것을 고를 수 있습니다. 역사적 배경은 일본의 것이지만, 그 바탕이 되는 급수 자체는 순수하고 보편적인 수학이므로 어디에서나 동일하게 π로 수렴합니다.

사용 방법

먼저 드롭다운에서 공식을 선택하고, 합산할 항의 개수 \(N\)을 입력합니다(항이 많을수록 결과가 더 정확해집니다). 그다음 표시할 자릿수를 지정하세요. 계산기는 π의 근삿값과 함께 급수 합의 원시값(마츠나가는 \(\pi/3\), 다케베는 \(\pi^2/9\))을 함께 보여 주므로 중간 계산값을 직접 확인할 수 있습니다.

공식 풀이

마츠나가 급수에서는 각 항이 직전 항에 \((2k-1)^2\)을 곱하고 \(4k\times(4k+2)\)로 나눈 값이 됩니다. 누적 합 \(S\)로부터 \(\pi = 3S\)가 됩니다. 다케베 급수에서는 각 항에 \(k^2\)을 곱하고 \((2k+1)(2k+2)\)로 나눕니다. 이때 누적 합으로부터 \(\pi = 3 \times \sqrt{S}\)가 됩니다. 이렇게 항을 점진적으로 갱신하는 점화식을 사용하면 거대한 계승(팩토리얼)을 직접 계산할 필요가 없어 오버플로를 방지할 수 있습니다.

$$\frac{\pi}{3} = 1 + \frac{1^2}{4\cdot6} + \frac{1^2\cdot3^2}{4\cdot6\cdot8\cdot10} + \cdots$$$$\frac{\pi}{3} = \sum_{k=0}^{N-1} \frac{((2k-1)!!)^2}{4\cdot6\cdots(4k+2)}$$$$\frac{\pi^2}{9} = \sum_{k=0}^{N-1} \frac{(k!)^2}{3\cdot4\cdots(2k+2)}$$

와산 급수에서 더하는 각 항이 점점 작아지므로 부분합이 파이로 수렴합니다.

계산 예시

마츠나가 급수를 \(N = 4\)개 항으로 계산하면: $$1 + \frac{1}{24} + \frac{9}{1920} + \frac{225}{322560} = 1.0470517113$$이 되고, 따라서 π는 약 \(3 \times 1.0470517113 = 3.1411551340\)입니다. 항을 100개까지 더하면 배정밀도(double)가 표현할 수 있는 한계 정확도, 즉 약 \(3.14159265358979\)에 도달합니다.

부분합(1항, 2항, 3항…)이 파이 값으로 점점 올라갑니다.

자주 묻는 질문

자릿수를 늘려도 정확도가 더 좋아지지 않는 이유는? 이 계산은 배정밀도 부동소수점 연산으로 수행되며, 유효 자릿수가 약 15~16자리로 제한됩니다. 다케베와 마츠나가가 이뤄낸 역사적인 41자리·52자리 계산을 재현하려면 임의 정밀도(BigDecimal) 연산이 필요합니다.

항을 1개만 입력하면 어떻게 되나요? 급수가 맨 앞의 1만 반환하므로 두 공식 모두 π를 약 3으로 계산합니다.

어느 급수가 더 빨리 수렴하나요? 두 급수 모두 꾸준히 수렴하며, 실제로는 수백 개 항만 더해도 배정밀도의 한계에 도달하기에 충분합니다.