MCP๋กœ ์—ฐ๊ฒฐ โ†’

๊ณ„์‚ฐ ์ž…๋ ฅ

๊ณต์‹

๊ด‘๊ณ 

๊ฒฐ๊ณผ

Simpson's rule result (n = 64)
3.1415926535892149
์ •์ ๋ถ„ ๊ทผ์‚ฌ๊ฐ’
์‚ฌ๋‹ค๋ฆฌ๊ผด ๊ณต์‹
3.1415519634856541
์ค‘์  ๊ณต์‹
3.1416129986418473
์‹ฌํ”„์Šจ ๊ณต์‹
3.1415926535892149
n ์‚ฌ๋‹ค๋ฆฌ๊ผด ์ค‘์  ์‹ฌํ”„์Šจ
2 3.1000000000000001 3.1623529411764704 3.1333333333333333
4 3.1311764705882354 3.1468005183939427 3.1415686274509804
8 3.1389884944910889 3.1428947295916889 3.1415925024587068
16 3.1409416120413889 3.1419181743085600 3.1415926512248218
32 3.1414298931749749 3.1416740337963374 3.1415926535528360
64 3.1415519634856541 3.1416129986418473 3.1415926535892149

์ด ๊ณ„์‚ฐ๊ธฐ๋กœ ํ•  ์ˆ˜ ์žˆ๋Š” ๊ฒƒ

์ด ๋„๊ตฌ๋Š” ์œ ํ•œ ๊ตฌ๊ฐ„ [a, b]์—์„œ ํ•จ์ˆ˜ f(x)์˜ ์ •์ ๋ถ„์„ ์„ธ ๊ฐ€์ง€ ๊ณ ์ „์ ์ธ ํ•ฉ์„ฑ ์ ๋ถ„ ๊ณต์‹์œผ๋กœ ๊ทผ์‚ฌํ•ฉ๋‹ˆ๋‹ค. ๋ฐ”๋กœ ์‚ฌ๋‹ค๋ฆฌ๊ผด ๊ณต์‹(Trapezoidal Rule), ์ค‘์  ๊ณต์‹(Midpoint Rule), ์‹ฌํ”„์Šจ ๊ณต์‹(Simpson's Rule)์ž…๋‹ˆ๋‹ค. ๋‹จ์ˆœํžˆ ๊ฐ’ ํ•˜๋‚˜๋งŒ ๋‚ด๋†“๋Š” ๋Œ€์‹ , ๋ถ„ํ•  ์ˆ˜๋ฅผ 2, 4, 8, 16, โ€ฆ์ฒ˜๋Ÿผ 2๋ฐฐ์”ฉ ๋Š˜๋ ค๊ฐ€๋ฉฐ ์ง€์ •ํ•œ ์ตœ๋Œ€๊ฐ’ N๊นŒ์ง€ ๊ฐ ๊ณต์‹์„ ๊ณ„์‚ฐํ•˜๊ณ , ๊ทธ ๊ฒฐ๊ณผ๋ฅผ ์ˆ˜๋ ดํ‘œ๋กœ ๋ณด์—ฌ์ค๋‹ˆ๋‹ค. ์ถ”์ •๊ฐ’์ด ์ ์  ์•ˆ์ •๋˜๋Š” ๋ชจ์Šต์„ ์ง์ ‘ ๋ˆˆ์œผ๋กœ ํ™•์ธํ•˜๋ฉฐ ์ •ํ™•๋„๋ฅผ ๊ฐ€๋Š ํ•  ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค.

a์™€ b ์‚ฌ์ด์˜ ์ˆ˜์ง ๋ ๋กœ ๊ทผ์‚ฌํ•œ ๊ณก์„  ์•„๋ž˜ ๋ฉด์ 
์ˆ˜์น˜ ์ ๋ถ„์€ ๊ฐ„๋‹จํ•œ ๋ ๋ฅผ ํ•ฉ์‚ฐํ•˜์—ฌ a๋ถ€ํ„ฐ b๊นŒ์ง€ f(x) ์•„๋ž˜ ๋ฉด์ ์„ ์ถ”์ •ํ•ฉ๋‹ˆ๋‹ค.

์‚ฌ์šฉ ๋ฐฉ๋ฒ•

ํ”ผ์ ๋ถ„ ํ•จ์ˆ˜๋ฅผ ๋ณ€์ˆ˜ x์— ๋Œ€ํ•œ ์ˆ˜์‹์œผ๋กœ ์ž…๋ ฅํ•˜์„ธ์š”(์˜ˆ: 4/(1+x^2) ๋˜๋Š” sin(x)*exp(-x)). ์‚ฌ์šฉ ๊ฐ€๋Šฅํ•œ ์—ฐ์‚ฐ์ž๋Š” ๊ด„ํ˜ธ์™€ ํ•จ๊ป˜ + - * / ^ ์ด๋ฉฐ, sin, cos, tan, exp, log/ln, sqrt, abs ๊ฐ™์€ ์ž์ฃผ ์“ฐ๋Š” ํ•จ์ˆ˜์™€ ์ƒ์ˆ˜ pi, e๋„ ์ง€์›ํ•ฉ๋‹ˆ๋‹ค. ํ•˜ํ•œ a์™€ ์ƒํ•œ b๋ฅผ ์„ค์ •ํ•˜๊ณ , ์ตœ๋Œ€ ๋ถ„ํ•  ์ˆ˜ N(2์˜ ๊ฑฐ๋“ญ์ œ๊ณฑ)์„ ๊ณ ๋ฅธ ๋’ค, ํ™”๋ฉด์— ํ‘œ์‹œํ•  ์ž๋ฆฟ์ˆ˜๋ฅผ ์„ ํƒํ•˜์„ธ์š”. ๋Œ€ํ‘œ ๊ฒฐ๊ณผ๋กœ๋Š” \(n = N\)์ผ ๋•Œ์˜ ์‹ฌํ”„์Šจ ์ถ”์ •๊ฐ’์ด ํ‘œ์‹œ๋˜๋Š”๋ฐ, ๋ณดํ†ต ๊ฐ€์žฅ ๋น ๋ฅด๊ฒŒ ์ˆ˜๋ ดํ•˜๊ธฐ ๋•Œ๋ฌธ์ž…๋‹ˆ๋‹ค.

๊ณต์‹ ์ž์„ธํžˆ ์•Œ์•„๋ณด๊ธฐ

๊ตฌ๊ฐ„์„ n๊ฐœ๋กœ ๋‚˜๋ˆ„๋ฉด ๊ฐ„๊ฒฉ์€ \(h = (b - a)/n\) ์ด๊ณ  ๋ถ„์ ์€ \(x_i = a + i \cdot h\) ์ž…๋‹ˆ๋‹ค. ์‚ฌ๋‹ค๋ฆฌ๊ผด ๊ณต์‹์€ ํ•จ์ˆ˜๊ฐ’์„ ๋”ํ•˜๋˜ ์–‘ ๋์ ์— 1/2์˜ ๊ฐ€์ค‘์น˜๋ฅผ ๋‘ก๋‹ˆ๋‹ค. ์ค‘์  ๊ณต์‹์€ ๊ฐ ์†Œ๊ตฌ๊ฐ„์˜ ์ค‘์‹ฌ์—์„œ ํ•จ์ˆ˜๊ฐ’์„ ์ƒ˜ํ”Œ๋งํ•ฉ๋‹ˆ๋‹ค. ์‹ฌํ”„์Šจ ๊ณต์‹์€ ๋‘ ๋ฐฉ์‹์„ 1, 4, 2, 4, โ€ฆ, 4, 1์˜ ๊ฐ€์ค‘์น˜๋กœ ์„ž์œผ๋ฉฐ, 3์ฐจ ๋‹คํ•ญ์‹๊นŒ์ง€ ์ •ํ™•ํ•˜๊ฒŒ ์ ๋ถ„ํ•ฉ๋‹ˆ๋‹ค.

$$\int_{a}^{b} f(x)\,dx \approx \frac{h}{3}\left[ f(x_0) + 4\sum_{i\,\text{odd}} f(x_i) + 2\sum_{i\,\text{even}} f(x_i) + f(x_N) \right]$$

๊ทธ ๊ฒฐ๊ณผ ์˜ค์ฐจ๊ฐ€ \(h^4\) ์ฐจ์ˆ˜๋กœ ์ค„์–ด๋“œ๋Š”๋ฐ, ์ด๋Š” ๋‹ค๋ฅธ ๋‘ ๊ณต์‹์˜ \(h^2\)๋ณด๋‹ค ํ›จ์”ฌ ๋น ๋ฆ…๋‹ˆ๋‹ค.

๊ณก์„  ์•„๋ž˜ ์‚ฌ๋‹ค๋ฆฌ๊ผดยท์ค‘์ ยท์‹ฌํ”„์Šจ ๋  ๋ชจ์–‘์„ ๋น„๊ตํ•˜๋Š” ์„ธ ๊ฐœ์˜ ํŒจ๋„
์„ธ ๊ทœ์น™์€ ๊ฐ ๋ ๊ฐ€ ๊ณก์„ ์„ ๊ทผ์‚ฌํ•˜๋Š” ๋ฐฉ์‹์ด ๋‹ค๋ฆ…๋‹ˆ๋‹ค: ์ง์„  ํ˜„, ์ค‘์ ์˜ ํ‰ํ‰ํ•œ ์ง์‚ฌ๊ฐํ˜•, ๋งž์ถ˜ ํฌ๋ฌผ์„ .

๊ณ„์‚ฐ ์˜ˆ์‹œ

๊ตฌ๊ฐ„ [0, 1]์—์„œ \(f(x) = 4/(1+x^2)\)์˜ ์ •ํ™•ํ•œ ์ ๋ถ„๊ฐ’์€ \(\pi = 3.14159265\ldots\) ์ž…๋‹ˆ๋‹ค. \(n = 4\), \(h = 0.25\)์ผ ๋•Œ ์‚ฌ๋‹ค๋ฆฌ๊ผด ์ถ”์ •๊ฐ’์€ ์•ฝ \(3.131176\), ์ค‘์  ์ถ”์ •๊ฐ’์€ ์•ฝ \(3.146801\), ์‹ฌํ”„์Šจ ์ถ”์ •๊ฐ’์€ ์•ฝ \(3.141569\)๋กœ, ์ด๋ฏธ ์†Œ์ˆ˜์  ๋‹ค์„ฏ ์ž๋ฆฌ๊นŒ์ง€ ์ •ํ™•ํ•ฉ๋‹ˆ๋‹ค. n์„ ๋Š˜๋ฆด์ˆ˜๋ก ์„ธ ๊ฐ’ ๋ชจ๋‘ \(\pi\)์— ๊ฐ€๊นŒ์›Œ์ง‘๋‹ˆ๋‹ค.

์ž์ฃผ ๋ฌป๋Š” ์งˆ๋ฌธ

์™œ N์ด 2์˜ ๊ฑฐ๋“ญ์ œ๊ณฑ์ด์–ด์•ผ ํ•˜๋‚˜์š”? ๋ถ„ํ•  ์ˆ˜๋ฅผ 2๋ฐฐ์”ฉ ๋Š˜๋ฆฌ๋ฉด ํ‘œ์˜ ์ธ์ ‘ํ•œ ํ–‰๋ผ๋ฆฌ ๋น„๊ตํ•˜๊ธฐ ์‰ฝ๊ณ , ์‹ฌํ”„์Šจ ๊ณต์‹์— ํ•„์š”ํ•œ ์ง์ˆ˜ ๋ถ„ํ• ๋„ ์ž๋™์œผ๋กœ ๋ณด์žฅ๋ฉ๋‹ˆ๋‹ค.

์ ์šฉํ•  ์ˆ˜ ์—†๋Š” ํ•จ์ˆ˜๋Š” ๋ฌด์—‡์ธ๊ฐ€์š”? ์ˆ˜๋ ดํ‘œ๋Š” ํ”ผ์ ๋ถ„ ํ•จ์ˆ˜๊ฐ€ ๋งค๋„๋Ÿฝ๊ณ (ํ•ด์„์ ) ์ฃผ๊ธฐ์ ์ด์ง€ ์•Š๋‹ค๊ณ  ๊ฐ€์ •ํ•ฉ๋‹ˆ๋‹ค. 0์„ ์ง€๋‚˜๋Š” \(1/x\)์ฒ˜๋Ÿผ ๊ตฌ๊ฐ„ [a, b] ๋‚ด๋ถ€์— ํŠน์ด์ ์ด ์žˆ๋Š” ํ•จ์ˆ˜๋Š” ๋ฌดํ•œ๋Œ€๋‚˜ ์—‰๋šฑํ•œ ๊ฐ’์„ ๋‚ด๋†“์Šต๋‹ˆ๋‹ค.

a์™€ b๊ฐ€ ๊ฐ™๊ฑฐ๋‚˜ a๊ฐ€ b๋ณด๋‹ค ํฌ๋ฉด ์–ด๋–ป๊ฒŒ ๋˜๋‚˜์š”? \(a = b\)์ด๋ฉด ์ ๋ถ„๊ฐ’์€ 0์ž…๋‹ˆ๋‹ค. \(a > b\)์ด๋ฉด ๊ฒฐ๊ณผ๋Š” [b, a] ๊ตฌ๊ฐ„ ์ ๋ถ„๊ฐ’์˜ ๋ถ€ํ˜ธ๋ฅผ ๋’ค์ง‘์€(์Œ์ˆ˜) ๊ฐ’์ด ๋ฉ๋‹ˆ๋‹ค.

์ตœ์ข… ์—…๋ฐ์ดํŠธ: