뷔퐁의 바늘이란?
수학의 세계에는 직관에 반하는 놀라운 사실들과, 언뜻 무관해 보이는 현상이 아름답게 연결되는 정리들이 많이 존재합니다. 그중에서도 특히 유명하고 매력적인 문제 중 하나가 “뷔퐁의 바늘”(Buffon’s needle problem)입니다.
이 문제는 18세기 프랑스의 박물학자이자 수학자였던 조르주루이 르클레르 뷔퐁 백작(Georges-Louis Leclerc, Comte de Buffon)이 1733년에 제기하고, 1777년에 처음으로 풀이했습니다.
놀랍게도 이 문제는 “바닥에 무작위로 바늘을 떨어뜨린다"라는 극히 물리적이고 무작위적인 행위를 통해, 수학에서 가장 중요한 상수 중 하나인 원주율 $\pi$ 를 구할 수 있다는 것입니다. 이는 기하 확률론(Geometric probability)에서 가장 초기의 문제 중 하나로 알려져 있으며, 이후 몬테카를로 방법(Monte Carlo method)의 선구자라 할 수 있는 획기적인 발견이었습니다.
본 글에서는 뷔퐁의 바늘 의 문제 설정부터 수학적 증명, 그리고 현대 컴퓨터를 이용한 시뮬레이션을 통한 원주율 추정까지 상세하고 알기 쉽게 해설합니다.
문제의 기본 설정
뷔퐁의 바늘의 문제 설정은 매우 간단합니다.
- 평평한 바닥 위에 간격 $d$로 평행한 직선이 여러 개 그려져 있습니다.
- 길이 $l$인 바늘 한 개를 준비합니다.
- 이 바늘을 바닥을 향해 무작위로 떨어뜨립니다.
이때, “떨어진 바늘이 바닥에 그려진 평행선 중 하나와 교차할 확률은 얼마인가?” 라는 것이 뷔퐁이 제기한 문제입니다.
아래 다이어그램은 이 실험의 개념적 흐름을 보여줍니다.
graph TD
A["바늘을 무작위로 떨어뜨림"] -->|"중력과 공기 저항"| B["바닥에 착지"]
B --> C{"선과 교차했는가?"}
C -->|"예"| D["교차 횟수 카운트 (+1)"]
C -->|"아니오"| E["교차하지 않음"]
D --> F["확률 P 계산"]
E --> F
F -->|"P = 2l / (πd)"| G["원주율 π를 역산"]
여기서는 문제를 간단히 하기 위해 바늘의 길이 $l$이 평행선의 간격 $d$ 이하인($l \le d$) 짧은 바늘 의 경우를 고려합니다. 이 조건에서는 바늘이 동시에 2개 이상의 직선과 교차하는 일이 없습니다.
수학적 모델링과 확률 도출
이 문제를 수학적으로 풀기 위해서는 바늘의 상태를 수치화(매개변수화)할 필요가 있습니다. 바늘이 바닥에 떨어졌을 때, 그 위치와 방향은 완전히 무작위라고 가정합니다.
바늘의 위치를 결정하기 위해 다음 두 변수를 정의합니다.
- $x$: 바늘의 중심에서 가장 가까운 평행선까지의 수직 거리.
- $\theta$: 바늘과 평행선이 이루는 예각(또는 직각).
변수의 범위
먼저, 각 변수가 어떤 값을 취할 수 있는지 생각해 봅시다.
- 거리 $x$: 바늘의 중심은 인접한 두 평행선 사이의 어딘가에 떨어집니다. 가장 가까운 선까지의 거리를 고려하므로, $x$의 최솟값은 $0$(바늘의 중심이 선 위에 있을 때), 최댓값은 $\frac{d}{2}$(바늘의 중심이 두 선의 정확히 중간에 있을 때)입니다. 즉, $0 \le x \le \frac{d}{2}$입니다. 바늘을 무작위로 떨어뜨리므로 $x$는 이 범위에서 균일분포 를 따릅니다. 확률밀도함수는 $\frac{2}{d}$입니다.
- 각도 $\theta$: 바늘과 평행선이 이루는 각도는, 바늘이 직선과 평행할 때의 $0$에서 수직일 때의 $\frac{\pi}{2}$(90도)까지의 값을 취합니다. 대칭성으로 인해 이 이상의 각도를 고려할 필요가 없습니다. 따라서 $0 \le \theta \le \frac{\pi}{2}$입니다. 바늘의 방향도 무작위이므로 $\theta$도 이 범위에서 균일분포 를 따릅니다. 확률밀도함수는 $\frac{2}{\pi}$입니다.
변수 $x$와 $\theta$는 서로 독립이므로, 특정 쌍 $(x, \theta)$을 취하는 결합확률밀도함수 $f(x, \theta)$는 각각의 확률밀도함수의 곱으로 나타낼 수 있습니다.
$$ f(x, \theta) = \frac{2}{d} \times \frac{2}{\pi} = \frac{4}{d\pi} $$교차 조건
다음으로, 바늘이 직선과 교차하기 위한 조건을 생각합니다. 바늘이 직선과 교차하는 것은, 바늘의 중심에서 끝까지의 수직 방향 길이가 가장 가까운 선까지의 거리 $x$ 이상인 경우입니다.
바늘의 길이가 $l$이므로 중심에서 끝까지의 길이는 $\frac{l}{2}$입니다. 각도가 $\theta$일 때, 이 반쪽 바늘이 수직 방향으로 차지하는 거리(투영 길이)는 $\frac{l}{2} \sin \theta$입니다.
따라서 바늘이 직선과 교차하는 조건은 다음 부등식으로 나타낼 수 있습니다.
$$ x \le \frac{l}{2} \sin \theta $$확률의 계산
바늘이 선과 교차할 확률 $P$는, 교차 조건을 만족하는 영역에 대해 결합확률밀도함수 $f(x, \theta)$를 적분하여 구합니다.
$$ P = \iint_{\text{교차 영역}} f(x, \theta) \, dx \, d\theta $$구체적인 적분 범위는, $\theta$가 $0$에서 $\frac{\pi}{2}$까지 변하고, $x$가 $0$에서 교차 한계값 $\frac{l}{2} \sin \theta$까지 변하는 범위입니다.
$$ P = \int_{0}^{\frac{\pi}{2}} \int_{0}^{\frac{l}{2} \sin \theta} \frac{4}{d\pi} \, dx \, d\theta $$먼저, 안쪽 $x$에 대한 적분을 계산합니다.
$$ \int_{0}^{\frac{l}{2} \sin \theta} \frac{4}{d\pi} \, dx = \frac{4}{d\pi} \left[ x \right]_{0}^{\frac{l}{2} \sin \theta} = \frac{4}{d\pi} \left( \frac{l}{2} \sin \theta - 0 \right) = \frac{2l}{d\pi} \sin \theta $$다음으로, 바깥쪽 $\theta$에 대한 적분을 계산합니다.
$$ P = \int_{0}^{\frac{\pi}{2}} \frac{2l}{d\pi} \sin \theta \, d\theta = \frac{2l}{d\pi} \int_{0}^{\frac{\pi}{2}} \sin \theta \, d\theta $$$\sin \theta$의 적분은 $-\cos \theta$이므로,
$$ \int_{0}^{\frac{\pi}{2}} \sin \theta \, d\theta = \left[ -\cos \theta \right]_{0}^{\frac{\pi}{2}} = (-\cos \frac{\pi}{2}) - (-\cos 0) = -0 - (-1) = 1 $$따라서 구하려는 확률 $P$는 다음과 같습니다.
$$ P = \frac{2l}{d\pi} \times 1 = \frac{2l}{\pi d} $$이것이 뷔퐁의 바늘 의 기본 공식입니다. 바늘이 선과 교차할 확률은, 바늘의 길이 $l$의 2배를 원주율 $\pi$와 선의 간격 $d$의 곱으로 나눈 값입니다.
원주율 π의 추정 (몬테카를로 방법)
도출된 공식 $P = \frac{2l}{\pi d}$에는 훌륭하게 $\pi$가 포함되어 있습니다. 이를 $\pi$에 대해 풀면 다음과 같습니다.
$$ \pi = \frac{2l}{P d} $$이 식은 확률 $P$만 알면 원주율 $\pi$를 계산할 수 있음을 의미합니다. 물론 진정한 확률 $P$는 무한 번의 시행을 해야만 알 수 있지만, 실제 실험에서 바늘을 여러 번 떨어뜨림으로써 $P$의 근사값을 얻을 수 있습니다.
바늘을 떨어뜨린 총 횟수를 $N$, 그중 바늘이 선과 교차한 횟수를 $C$라 합시다. 시행 횟수 $N$이 충분히 크면, 큰 수의 법칙에 의해 경험적 확률 $\frac{C}{N}$은 이론적 확률 $P$에 가까워집니다.
$$ P \approx \frac{C}{N} $$이를 앞서의 식에 대입하면 원주율 $\pi$의 근사값을 구하는 공식을 얻습니다.
$$ \pi \approx \frac{2l \cdot N}{C \cdot d} $$계산이 가장 간단한 것은 바늘의 길이 $l$과 선의 간격 $d$가 같을 때($l = d$)입니다. 이 경우 식은 더욱 단순해집니다.
$$ \pi \approx \frac{2N}{C} $$즉, “바늘을 떨어뜨린 횟수"의 2배를 “교차한 횟수"로 나누기만 하면 원주율을 구할 수 있습니다!
Python 시뮬레이션
실제로 손으로 바늘을 수천 번이나 떨어뜨리는 것은 매우 고된 작업입니다(역사상 실제로 수천 번의 실험을 수행한 수학자도 존재하지만). 현대에는 컴퓨터를 사용하여 이 실험을 간단하게 시뮬레이션할 수 있습니다.
아래에 Python을 사용하여 뷔퐁의 바늘 실험을 시뮬레이션하고 원주율을 추정하는 간단한 코드 예시를 보여드립니다.
| |
이 코드를 실행하면, 난수를 사용하여 대량의 가상 바늘이 떨어지고, 매우 높은 정밀도로 $3.1415...$ 라는 원주율의 근사값을 얻을 수 있음을 확인할 수 있습니다. 이처럼 난수를 사용하여 확률적 문제의 근사해를 구하는 기법을 몬테카를로 방법 이라고 합니다.
마무리
뷔퐁의 바늘은 언뜻 보면 단순한 물리적 우연의 놀이처럼 보이지만, 그 뒤에는 확고한 수학적 이론이 존재합니다. 무작위 사건(확률)과 기하학적 형상(직선과 선분), 그리고 궁극의 무리수인 $\pi$가 하나의 간단한 수식 안에서 융합되는 모습은 수학의 아름다움을 체현한다고 할 수 있겠습니다.
또한, 이 문제는 현대 과학 기술에 필수적인 몬테카를로 방법의 원점으로서 역사적 중요성도 지니고 있습니다. 복잡한 계의 시뮬레이션이나 해석적으로 풀기 어려운 적분의 계산 등, 오늘날에도 뷔퐁의 아이디어는 다양한 형태로 우리의 세계를 지탱하고 있습니다.
종이와 펜, 그리고 이쑤시개 몇 개를 준비하여 집에서 이 위대한 수학 역사의 한 부분을 직접 체험해 보시는 것은 어떨까요?