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

๊ณ„์‚ฐ ์ž…๋ ฅ

This rebuilt version computes integer order v โ‰ฅ 0 for real x > 0. (The original tool also supports complex x and arbitrary real order.)

๊ณต์‹

๊ณต์‹: ๊ตฌ๋ฉด ๋ฒ ์…€ ํ•จ์ˆ˜ ๊ณ„์‚ฐ๊ธฐ

๊ด‘๊ณ 

๊ฒฐ๊ณผ

๊ตฌ๋ฉด ๋ฒ ์…€ jv(x)
0.4546487134
์ œ1์ข…
ํ•จ์ˆ˜ ๊ฐ’
jv(x) โ€” first kind 0.4546487134
yv(x) โ€” second kind 0.2080734183
j'v(x) โ€” derivative -0.435397775
y'v(x) โ€” derivative 0.3506120043

์ด ๊ณ„์‚ฐ๊ธฐ์˜ ๊ธฐ๋Šฅ

์ด ๋„๊ตฌ๋Š” ์ œ1์ข… ๊ตฌ๋ฉด ๋ฒ ์…€ ํ•จ์ˆ˜ \(j_v(x)\), ์ œ2์ข… ๊ตฌ๋ฉด ๋ฒ ์…€ ํ•จ์ˆ˜ \(y_v(x)\), ๊ทธ๋ฆฌ๊ณ  ๊ฐ๊ฐ์˜ 1์ฐจ ๋„ํ•จ์ˆ˜ \(j'_v(x)\)์™€ \(y'_v(x)\)๋ฅผ ๊ณ„์‚ฐํ•ฉ๋‹ˆ๋‹ค. ์ด ํ•จ์ˆ˜๋“ค์€ ๊ตฌ๋ฉด ์ขŒํ‘œ๊ณ„์—์„œ ํŒŒ๋™ ๋ฐฉ์ •์‹๊ณผ ํ—ฌ๋ฆ„ํ™€์ธ  ๋ฐฉ์ •์‹์˜ ๋™๊ฒฝ(radial) ๋ฐฉํ–ฅ ํ•ด์ด๋ฉฐ, ์‚ฐ๋ž€ ์ด๋ก , ์ „์ž๊ธฐํŒŒยท์Œํ–ฅ ๋ฐฉ์‚ฌ, ์–‘์ž์—ญํ•™(์ž์œ  ์ž…์ž์˜ ๋ถ€๋ถ„ํŒŒ) ๋“ฑ ๋ฌผ๋ฆฌํ•™ ์ „๋ฐ˜์—์„œ ๋“ฑ์žฅํ•ฉ๋‹ˆ๋‹ค. ์ƒˆ๋กญ๊ฒŒ ๊ตฌํ˜„๋œ ์ด ๋ฒ„์ „์€ ์ •์ˆ˜ ์ฐจ์ˆ˜ \(v \ge 0\)๊ณผ ์‹ค์ˆ˜ ์ธ์ˆ˜ \(x > 0\)์„ ์ฒ˜๋ฆฌํ•ฉ๋‹ˆ๋‹ค.

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

์ฐจ์ˆ˜ \(v\)(0, 1, 2์™€ ๊ฐ™์€ ์Œ์ด ์•„๋‹Œ ์ •์ˆ˜)์™€ ์ธ์ˆ˜ \(x\)(์–‘์˜ ์‹ค์ˆ˜)๋ฅผ ์ž…๋ ฅํ•˜์„ธ์š”. ๊ณ„์‚ฐ ๋ฒ„ํŠผ์„ ๋ˆ„๋ฅด๋ฉด ๋„ค ๊ฐ€์ง€ ๊ฐ’์„ ๋ชจ๋‘ ์–ป์„ ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค. ๋‹ค๋งŒ \(y_v(x)\)์™€ \(y'_v(x)\)๋Š” \(x\)๊ฐ€ 0์— ๊ฐ€๊นŒ์›Œ์งˆ์ˆ˜๋ก ๋ฐœ์‚ฐํ•˜๋ฏ€๋กœ, \(x = 0\)์—์„œ๋Š” ๋ฌดํ•œ๋Œ€๋กœ ํ‘œ์‹œ๋ฉ๋‹ˆ๋‹ค. \(j_0(0)\)์€ ๊ทนํ•œ๊ฐ’์œผ๋กœ 1์ด ๋ฉ๋‹ˆ๋‹ค.

๊ณต์‹ ์„ค๋ช…

์ด ํ•จ์ˆ˜๋“ค์€ ๋ฏธ๋ถ„๋ฐฉ์ •์‹ \(x^2 w'' + 2x w' + (x^2 - v(v+1))w = 0\)์„ ๋งŒ์กฑํ•ฉ๋‹ˆ๋‹ค. ๋‹ซํžŒ ํ˜•ํƒœ์ธ \(j_0 = \sin(x)/x\)์™€ \(y_0 = -\cos(x)/x\)์—์„œ ์ถœ๋ฐœํ•˜์—ฌ, ๋” ๋†’์€ ์ฐจ์ˆ˜๋Š” ์‚ผํ•ญ ์ ํ™”์‹ $$f_{n+1} = \frac{2n+1}{x}\,f_n - f_{n-1}$$์„ ๋”ฐ๋ฆ…๋‹ˆ๋‹ค. ์ƒํ–ฅ ์ ํ™”๋Š” \(y_v\)์— ๋Œ€ํ•ด์„œ๋Š” ์•ˆ์ •์ ์ด์ง€๋งŒ, \(j_v\)์˜ ๊ฒฝ์šฐ \(n > x\)์ผ ๋•Œ ๋ถˆ์•ˆ์ •ํ•ฉ๋‹ˆ๋‹ค. ๊ทธ๋ž˜์„œ ๋ฐ€๋Ÿฌ(Miller)์˜ ํ•˜ํ–ฅ ์ ํ™”๋ฒ•์„ ์‚ฌ์šฉํ•ฉ๋‹ˆ๋‹ค. ์ฆ‰ ์ถฉ๋ถ„ํžˆ ๋†’์€ ์ฐจ์ˆ˜์—์„œ \(f\)๋ฅผ ๊ฐ๊ฐ 0๊ณผ 1๋กœ ์„ค์ •ํ•ด ์‹œ์ž‘ํ•œ ๋’ค ์ฐจ์ˆ˜๋ฅผ ๋‚ฎ์ถฐ ๊ฐ€๋ฉฐ ์ ํ™”์‹œํ‚ค๊ณ , 0์ฐจ ํ•ญ์ด \(\sin(x)/x\)์™€ ์ผ์น˜ํ•˜๋„๋ก ๋ชจ๋“  ๊ฐ’์„ ๋‹ค์‹œ ์Šค์ผ€์ผ๋งํ•ฉ๋‹ˆ๋‹ค. ๋„ํ•จ์ˆ˜๋Š” \(j'_v = j_{v-1} - \frac{v+1}{x}j_v\) ๊ณต์‹์„ ์‚ฌ์šฉํ•ฉ๋‹ˆ๋‹ค.

๊ฐ์‡ ํ•˜๋Š” ์ง„๋™์œผ๋กœ ๋‚˜ํƒ€๋‚ธ ์ œ1์ข… ๊ตฌ๋ฉด ๋ฒ ์…€ ํ•จ์ˆ˜ ๊ทธ๋ž˜ํ”„
์ œ1์ข… ๊ตฌ๋ฉด ๋ฒ ์…€ ํ•จ์ˆ˜ \(j_v(x)\)๋Š” \(x\)๊ฐ€ ์ปค์งˆ์ˆ˜๋ก ์ง„๋™ํ•˜๋ฉฐ ๊ฐ์‡ ํ•œ๋‹ค.

๊ณ„์‚ฐ ์˜ˆ์‹œ (v = 0, x = 2)

$$j_0(2) = \frac{\sin(2)}{2} = 0.4546487134, \qquad y_0(2) = -\frac{\cos(2)}{2} = 0.2080734183$$์ž…๋‹ˆ๋‹ค. \(j'_0 = -j_1\)์ด๋ฏ€๋กœ \(j'_0(2) = -0.4353977750\)์ด ๋˜๊ณ , \(y'_0 = -y_1\)์ด๋ฏ€๋กœ \(y'_0(2) = 0.3506120043\)์ด ๋ฉ๋‹ˆ๋‹ค.

x=2 ๋ถ€๊ทผ์— ์ ์ด ํ‘œ์‹œ๋œ j0(x) = sin x / x ๊ทธ๋ž˜ํ”„
0์ฐจ ๊ณก์„  \(j_0(x) = \sin(x)/x\), ์˜ˆ์ œ์—์„œ \(x = 2\)์—์„œ ๊ณ„์‚ฐํ•จ.

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

๋ณต์†Œ์ˆ˜ \(x\)๋„ ์ง€์›ํ•˜๋‚˜์š”? ์•„๋‹ˆ์š”. ์›๋ž˜ ํŽ˜์ด์ง€๋Š” ๋ณต์†Œ์ˆ˜ ์ธ์ˆ˜๋ฅผ ๋ฐ›์ง€๋งŒ, ์ด ๋ฒ„์ „์€ ๋ช…ํ™•์„ฑ๊ณผ ์†๋„๋ฅผ ์œ„ํ•ด ์‹ค์ˆ˜ \(x > 0\)์œผ๋กœ ์ œํ•œํ–ˆ์Šต๋‹ˆ๋‹ค.

์™œ \(y_v\)๋Š” \(x = 0\)์—์„œ ๋ฌดํ•œ๋Œ€์ธ๊ฐ€์š”? ์ œ2์ข… ํ•จ์ˆ˜๋Š” ์›์ ์— ๊ทน(pole)์„ ๊ฐ€์ง€๋ฏ€๋กœ, \(x\)๊ฐ€ 0์— ๊ฐ€๊นŒ์›Œ์งˆ์ˆ˜๋ก ๊ฐ’์ด ํ•œ์—†์ด ์ปค์ง‘๋‹ˆ๋‹ค.

์ •ํ™•๋„๋Š” ์–ด๋А ์ •๋„์ธ๊ฐ€์š”? ๋ฐฐ์ •๋ฐ€๋„(double precision)๋กœ ๊ณ„์‚ฐํ•˜์—ฌ ์•ฝ 15์ž๋ฆฌ์˜ ์œ ํšจ์ˆซ์ž๋ฅผ ์ œ๊ณตํ•˜๋ฏ€๋กœ, ์ผ๋ฐ˜์ ์ธ ๊ณตํ•™ยท๋ฌผ๋ฆฌ ์ž‘์—…์—๋Š” ์ถฉ๋ถ„ํ•ฉ๋‹ˆ๋‹ค.

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