손으로 제곱근 계산하기

$\sqrt{2}=1.4142\cdots$임은 외우고 있을 것이다. 하지만 $\sqrt{30}$의 근삿값을 소수점 아래 다섯 자리까지 찾으라고 한다면 조금은 당황스러울 것이다. 이제 손으로 제곱근 값을 계산하는 방법을 알아보자.

바빌로니아 방법(Babylonian method)

$\sqrt{S}$의 근삿값 $x_n$을 찾고 다음으로 $x_{n+1}$을 아래와 같은 점화식으로 찾아 나가는 방식이다. wiki

$$x_{n+1}=\frac{1}{2}(x_n+\frac{S}{x_n})$$

$\sqrt{30}$은 $x_0 =5$로 놓고 풀면 아래와 같이 아주 빠르게 근삿값에 가까워짐을 알 수 있다.

$$x_1=\frac{1}{2}(5+\frac{30}{5})=5.5$$

$$x_2=\frac{1}{2}(5.5+\frac{30}{5.5})=5.477272727$$

$$x_3=\frac{1}{2}(5.477225575+\frac{30}{5.477225575})=5.477225575$$

$$x_4=\frac{1}{2}(5.477225575+\frac{30}{5.477225575})=5.477225575$$

이 방법은 아래와 같이 증명할 수 있다.

$\varepsilon >0$은 이미 구한 근삿값 $x_n$보다 아주 작은 수라고 가정하자. 

$$\sqrt{S}=x_n+\epsilon =x_{n+1}$$

$$S=(x_n+\epsilon {)}^2={x_n}^2+2\epsilon x_n+{\epsilon }^2={x_n}^2+\epsilon (2x_n+\epsilon )$$

정리하면

$$\epsilon =\frac{S-{x_n}^2}{2x_n+\epsilon }$$

$\varepsilon >0$은 아주 작은 수이므로

$$\epsilon =\frac{S-{x_n}^2}{2x_n}$$

라고 할 수 있다. 이제 새로운 근삿값 $x_{n+1}$을 구하면

$$x_{n+1}=x_n+\frac{S-{x_n}^2}{2x_n}=\frac{{x_n}^2+S}{2x_n}=\frac{1}{2}(x_n+\frac{S}{x_n})$$

이 방법은 방정식 $x^2 -S=0$의 해를 찾을 때, $f(x)=x^2 -S$의 $x$절편을 찾는 뉴튼 방법과 같다.

$(x_0 ,f(x_0))$에서 접선 방정식은 아래와 같다.  $$y-f(x_0)=f^{\prime }(x_0)(x-x_0)$$

한편, $f^{\prime}(x_0)=2x_0$이므로

$$0-f(x_0)=2x_0(x_1-x_0)$$

$$x_1=x_0-\frac{f(x_0)}{2x_0}=x_0-\frac{x_0^2-S}{2x_0}=\frac{1}{2}(x_0+\frac{S}{x_0})$$

일반화하면 같은 결과를 얻는다.

바빌로니아 수학이 참 놀랍다. 하지만 엄밀한 증명으로 나아가지 못하고 근삿값 계산에 머물렀기 때문에 더 발전하지 못한 것으로 생각한다. 역시 수학의 꽃은 증명이다.

연분수로 구하기

아래와 같이 연분수 계산으로 제곱근의 근삿값을 구할 수도 있다.

$$\sqrt{2}=1+\sqrt{2}-1=1+\frac{1}{1+\sqrt{2}}=1.5(\sqrt{2}\approx 1)$$

$$\sqrt{2}=1+\frac{1}{2+\frac{1}{1+\sqrt{2}}}=1+\frac{1}{2+\frac{1}{1+1}}=1+\frac{2}{5}=1.4$$

$$\sqrt{2}=1+\frac{1}{2+\frac{1}{1+\sqrt{2}}}=1+\frac{1}{2+\frac{1}{2+\frac{1}{2+\frac{1}{2+\cdots }}}}$$

와 같이 나타낼 수 있으므로 근삿값을 바라는 만큼 정확하게 구할 수 있다.

 연분수(Continued fraction) 자세히 알아보기

소수 표현으로 구하기

$\sqrt2 = a_0.a_1 a_2$라고 하자.

$$\sqrt{2}=a_0+a_1\times \frac{1}{10}+a_2\times \frac{1}{{10}^2}$$

$$2=(a_0+a_1\times \frac{1}{10}+a_2\times \frac{1}{{10}^2}{)}^2$$

$$2={a_0}^2+{a_1}^2\times \frac{1}{{10}^2}+{a_2}^2\times \frac{1}{{10}^4}+2a_0{a}_1\frac{1}{10}+2a_1{a}_2\frac{1}{{10}^3}+2a_3{a}_0\frac{1}{{10}^2}$$

$a_i =0,1,2,\cdots,9$이므로 $a_0 =1$임을 쉽게 알수 있다.

$$2-1={a_1}^2\times \frac{1}{{10}^2}+{a_2}^2\times \frac{1}{{10}^4}+2a_0{a}_1\times \frac{1}{10}+2a_1{a}_2\times \frac{1}{{10}^3}+2a_3{a}_0\times \frac{1}{{10}^2}$$

$${10}^2(2-1)={a_1}^2+{a_2}^2\times \frac{1}{{10}^2}+2a_0{a}_1\times 10+2a_1{a}_2\times \frac{1}{10}+2a_3{a}_0$$

$${10}^2(2-1)={a_1}^2+2a_0{a}_1\times 10+{a_2}^2\times \frac{1}{{10}^2}+2a_1{a}_2\times \frac{1}{10}+2a_3{a}_0$$

$${10}^2(2-1)-a_1(2a_0\times 10+a_1)={a_2}^2\times \frac{1}{{10}^2}+2a_1{a}_2\times \frac{1}{10}+2a_3{a}_0$$

${a_1}(2a_0 \times 10+a_1)$가 $100$을 넘지 않는 최대 정수이려면 $a_1 =4$임을 찾을 수 있다.

이것을 일반화하면 아래와 같은 제곱근 계산법을 얻을 수 있다.

펠 방정식으로 구하기

 

펠 방정식(Pell's equation)은 아래와 같은 디오판토스 방정식이다.

$$x^2-N{y}^2=1$$

$N=2$일 때, 이 방정식의 해를 좌표평면에 나타내면 쌍곡선이 된다. 이때, 방정식의 정수해를 찾는 것은 쌍곡선 위의 점 가운데 $x,y$좌표가 모두 정수인 것을 찾는 것과 같다. $(-1,0) ,\;\;(1,0)$은 당연한(trivial) 해이다. 쉽게 찾을 수 있는 6개의 해는 아래와 같다.

인도의 브라마굽타는 아래와 같은 방법으로 알려진 두 개의 해 $(x_1,y_1),\;\;(x_2,y_2)$로 다른 해를 찾았다고 한다.

${x_1}^2 -N{y_1}^2 =1,\;\;{x_2}^2 -N{y_2}^2 =1$라고 하면

$(x_1^2 - Ny_1^2)(x_2^2 - Ny_2^2) = (x_1x_2 + Ny_1y_2)^2 - N(x_1y_2 + x_2y_1)^2 = (x_1x_2 - Ny_1y_2)^2 - N(x_1y_2 - x_2y_1)^2=1$이므로 $(x_1x_2 + Ny_1y_2,x_1y_2 + x_2y_1)$와 $(x_1x_2 - Ny_1y_2,x_1y_2 - x_2y_1)$도 해가 된다.

방정식 $x^2 -2y^2 =1$의 해인 $(3,2)$로 새로운 해 $(3\cdot 3+2\cdot 2\cdot 2, \;\;3\cdot 2+3\cdot 2)=(17,\;12)$을 얻을 수 있다. 얻어낸 두 해로 같은 계산을 되풀이하여 $(99,70)$을 얻을 수 있다. 되풀이하면 무한히 많은 해를 얻을 수 있다. 한편 방정식을 정리하면

$$N=\frac{x^2-1}{y^2}=(\frac{x}{y}{)}^2-\frac{1}{y^2}$$

이다. $y\rightarrow \infty$라고 하면 $\displaystyle{\frac{x}{y} \rightarrow \sqrt{N}}$이다.(쌍곡선의 점근선 기울기를 생각하면 된다.) 그러므로 펠 방정식의 정수해를 구함으로써 제곱근의 근삿값을 구할 수 있다.

http://en.wikipedia.org/wiki/Pell%27s_equation

Pell's equation - Wikipedia