티스토리 뷰

mathe

1/√(1-x^2)의 적분 = pi/2

게으른 the lazy 2022. 10. 22. 12:46

 

 

아래 부정적분은,

 

10dx1x2

 

x=sinu로 치환하여 풀 수 있다.

 

10dx1x2=π/20cosuducosu=π2

 

왜 뜬금없이 원주율이 나올까?

 

 

 

P(Px,Py): y=1x2 위의 한 점

A(0,Py), P에서 y축에 내린 수선의 발

B(Px,0), P에서 x축에 내린 수선의 발

C(Px,1/1P2x): BP의 연장선이 y=1/1x2와 만나는 점

 

각 선분의 길이는 아래와 같다.

 

PA=PxPB=PyPC=1PyPy=1P2yPy=P2xPy

 

이제 다음과 같은 미소선분과 면적을 생각하자.

 

 

 

x=0에서 x=1까지 dAx를 적분하면 사분원의 면적이므로 π/4가 된다. dAy도 마찬가지이다.

 

x=1x=0dAx=y=1y=0dAy=π4

 

dAxdAy의 면적 비율은 다음과 같다.

 

dAydAx=PxdyPydx

 

dAydAf의 면적 비율은 다음과 같다.

 

dAfdAy=(P2x/Py)dxPxdy=PxdxPydy

 

그런데

 

PyPx=dxdy

 

이므로 결과적으로 dAf=dAy이다. 따라서 y=0에서 y=1까지 dAy를 적분한 것과 x=0에서 x=1까지 dAf를 적분한 것은 그 값이 같을 것이다. 정리하면,

 

x=1x=0dAx=y=1y=0dAy=x=1x=0dAf=π4

 

이다. 이것을 그림으로 표현하면 아래와 같다.

 

 

 

두 색칠된 영역의 면적은 π/4로 같으며, 원래 문제였던

 

10dx1x2

 

는 두 영역의 합과 같으므로 값은 π/2가 된다.

 

위 그림을 그린 매트랩 코드

 

x = linspace(0, 1-1e-3);
y1 = sqrt(1-x.^2);
y2 = 1./sqrt(1-x.^2);

figure, hold on,
plot(x, y1)
plot(x, y2)
xline(1, 'k--')
yline(1, 'k--')

patch([x, flip(x)], [zeros(size(x)), flip(y2)], 'b', FaceAlpha=0.1)
patch([x, flip(x)], [y1, flip(y2)], 'g', FaceAlpha=0.1)

axis equal
xlim([0, 1.2])
ylim([0, 2])

 

'mathe' 카테고리의 다른 글

입델 그런거 왜 함?  (2) 2022.11.25
[proof without words] determinant == area?  (0) 2022.11.25
벡터가 뭐지?  (2) 2022.11.23
공허참이란?  (2) 2022.11.22
[수학문제] |75-(5627)^(1/2)|와 |7-(341)^(1/3)| 중 어느 것이 더 클까?  (2) 2022.10.30
댓글