지난시간에 배운 내용을 잠깐 리뷰해봅시다.
오일러는 12!=√π212!=√π2 라는 결과에 영감을 받아 팩토리얼이 적분과 관련이 있을 것이라고 생각하게 되고, 아래 적분을 떠올립니다.
∫10xe(1−x)ndx∫10xe(1−x)ndx
이 적분식을 지난시간에 아래와 같이 변형했습니다.
∫10xe(1−x)ndx=n!(e+1)(e+2)⋯(e+n)(e+n+1)∫10xe(1−x)ndx=n!(e+1)(e+2)⋯(e+n)(e+n+1)
팩토리얼과 적분을 연결한 식이 유도되었습니다. 오일러는 위 수식을 변형해서 감마함수 적분형을 유도합니다. 유도해보겠습니다.
먼저 ee 를 fgfg 로 치환합시다.
∫10xfg(1−x)ndx=n!(fg+1)(fg+2)⋯(fg+n)(fg+n+1)∫10xfg(1−x)ndx=n!(fg+1)(fg+2)⋯(fg+n)(fg+n+1)
아래와 같이 우변 분모의 각 항을 통분해줍니다.
∫10xfg(1−x)ndx=n!(f+gg)(f+2gg)⋯(f+ngg)(f+(n+1)gg)∫10xfg(1−x)ndx=n!(f+gg)(f+2gg)⋯(f+ngg)(f+(n+1)gg)
우변을 아래와 같이 변형합시다.
∫10xfg(1−x)ndx=n!(f+g)(f+2g)⋯(f+ng)(f+(n+1)g)gn+1∫10xfg(1−x)ndx=n!(f+g)(f+2g)⋯(f+ng)(f+(n+1)g)gn+1
아래와 같이 우변을 둘로 분리해줍니다.
∫10xfg(1−x)ndx=gn+1f+(n+1)g⋅n!(f+g)(f+2g)⋯(f+ng)∫10xfg(1−x)ndx=gn+1f+(n+1)g⋅n!(f+g)(f+2g)⋯(f+ng)
아래와 같이 우변의 첫항을 좌변으로 옮겨줍니다.
f+(n+1)ggn+1⋅∫10xfg(1−x)ndx=n!(f+g)(f+2g)⋯(f+ng)f+(n+1)ggn+1⋅∫10xfg(1−x)ndx=n!(f+g)(f+2g)⋯(f+ng)
좌우 변을 바꿔줍니다.
n!(f+g)(f+2g)⋯(f+ng)=f+(n+1)ggn+1⋅∫10xfg(1−x)ndxn!(f+g)(f+2g)⋯(f+ng)=f+(n+1)ggn+1⋅∫10xfg(1−x)ndx
치환을 할겁니다. 아래와 같이 치환합니다.
x=sgf+gx=sgf+g
dx=gf+gs−ff+gdsdx=gf+gs−ff+gds
위 식에 대입합시다.
n!(f+g)(f+2g)⋯(f+ng)=f+(n+1)ggn+1⋅∫10sff+g(1−sgf+g)ngf+gs−ff+gdsn!(f+g)(f+2g)⋯(f+ng)=f+(n+1)ggn+1⋅∫10sff+g(1−sgf+g)ngf+gs−ff+gds
아래와 같이 소거합니다.
n!(f+g)(f+2g)⋯(f+ng)=f+(n+1)ggn+1⋅∫10gf+g(1−sgf+g)ndsn!(f+g)(f+2g)⋯(f+ng)=f+(n+1)ggn+1⋅∫10gf+g(1−sgf+g)nds
아래와 같이 변형합니다.
n!(f+g)(f+2g)⋯(f+ng)=f+(n+1)ggn+1⋅∫10(gf+g)n+1(gf+g)n(1−sgf+g)ndsn!(f+g)(f+2g)⋯(f+ng)=f+(n+1)ggn+1⋅∫10(gf+g)n+1(gf+g)n(1−sgf+g)nds
아래와 같이 변형합니다. 적분식 안의 첫 항의 분자는 적분식 바깥의 식과 곱해주고, 분모는 적분식 안의 두번째 항과 합쳐준 것입니다.
n!(f+g)(f+2g)⋯(f+ng)=f+(n+1)g(f+g)n+1⋅∫10(1−sgf+gg/f+g)ndsn!(f+g)(f+2g)⋯(f+ng)=f+(n+1)g(f+g)n+1⋅∫10(1−sgf+gg/f+g)nds
이제 f를 1로 보내고, g를 0으로 보낼 것입니다.
limf→1,g→0n!(f+g)(f+2g)⋯(f+ng)=limf→1,g→0f+(n+1)g(f+g)n+1⋅∫10(1−sgf+gg/f+g)ndslimf→1,g→0n!(f+g)(f+2g)⋯(f+ng)=limf→1,g→0f+(n+1)g(f+g)n+1⋅∫10(1−sgf+gg/f+g)nds
적분안의 항만 따로 계산해봅시다. 아래와 같습니다.
limf→1,g→01−sgf+gg/f+glimf→1,g→01−sgf+gg/f+g
계산을 편하게 하기 위해 gf+ggf+g 를 z로 치환합시다.
limz→01−szzlimz→01−szz
로피탈 정리를 이용하여 극한값을 구하겠습니다. 분자와 분모를 z로 미분합니다. (ax)′=lna⋅ax 를 이용합니다.
limz→01−szz=limz→0−lns⋅sz1=−lns
우리가 계산하던 식을 다시 가져옵시다.
limf→1,g→0n!(f+g)(f+2g)⋯(f+ng)=limf→1,g→0f+(n+1)g(f+g)n+1⋅∫10(1−sgf+gg/f+g)nds
극한값은 아래와 같습니다.
n!=∫10(−lns)ndt
치환을 하겠습니다. −lns를 t로 치환합시다.
−lns=ts=e−tds=−e−tdt
적분범위를 구하기 위해 그래프를 그리면 아래와 같습니다.

s가 0에서 1로 변할 때, t는 무한대에서 0으로 변합니다.
n!=∫0∞tn(−e−t)dt
아래와 같이 변형합시다. 적분구간을 뒤집으면 부호가 바뀝니다.
n!=∫∞0tne−tdt
감마함수는 팩토리얼함수에서 정의역을 확장한 것이므로 n대신 n-1을 넣어서 팩토리얼 함수로 만들어줍니다.
(n−1)!=∫∞0tn−1e−tdt
n을 실수 영역으로 확장하면 감마함수가 됩니다.
Γ(z)=∫∞0tz−1e−tdt
'@ 필수과목 > 손으로 푸는 통계' 카테고리의 다른 글
[손으로 푸는 통계 ver1.0] 66. 표본분산의 분포 유도 (31) 감마함수 수렴성 증명 #2 (0) | 2021.09.18 |
---|---|
[손으로 푸는 통계 ver1.0] 65. 표본분산의 분포 유도 (30) 감마함수 수렴성 증명 #1 (0) | 2021.08.01 |
[손으로 푸는 통계 ver1.0] 64. 표본분산의 분포 유도 (29) 감마 1/2 계산하기 (0) | 2021.08.01 |
[손으로 푸는 통계 ver1.0] 63. 표본분산의 분포 유도 (28) 감마함수 적분형의 재귀적 성질 (0) | 2021.08.01 |
[손으로 푸는 통계 ver1.0] 61. 표본분산의 분포 유도 (26) 팩토리얼과 적분의 연결 (0) | 2021.06.19 |
[손으로 푸는 통계 ver1.0] 60. 표본분산의 분포 유도 (25) 이분의일 팩토리얼이 이분의 루트 파이임을 증명 (0) | 2021.06.19 |
[손으로 푸는 통계 ver1.0] 59. 표본분산의 분포 유도 (24) 왈리스 공식 유도3 (Wallis product) (0) | 2021.06.19 |
[손으로 푸는 통계 ver1.0] 58. 표본분산의 분포 유도 (23) 왈리스 공식 유도2 (Wallis product) (0) | 2021.06.12 |
댓글
bigpicture님의
글이 좋았다면 응원을 보내주세요!
이 글이 도움이 됐다면, 응원 댓글을 써보세요. 블로거에게 지급되는 응원금은 새로운 창작의 큰 힘이 됩니다.
응원 댓글은 만 14세 이상 카카오계정 이용자라면 누구나 편하게 작성, 결제할 수 있습니다.
글 본문, 댓글 목록 등을 통해 응원한 팬과 응원 댓글, 응원금을 강조해 보여줍니다.
응원금은 앱에서는 인앱결제, 웹에서는 카카오페이 및 신용카드로 결제할 수 있습니다.