์ด ๊ณ์ฐ๊ธฐ์ ๊ธฐ๋ฅ
์ด ๋๊ตฌ๋ ์ฌ์ฉ์๊ฐ ์ง์ ํ x ๊ฐ ๋ฒ์์ ๋ํด ์๊ณก ์ฌ์ธ ์ ๋ถ Shi(x)์ ์๊ณก ์ฝ์ฌ์ธ ์ ๋ถ Chi(x)๋ฅผ ํ๋ก ์ ๋ฆฌํ๊ณ , ๋ ๊ณก์ ์ ํ๋์ ๊ทธ๋ํ์ ํจ๊ป ๊ทธ๋ ค์ค๋๋ค. ์ด ๋ ํจ์๋ ์ผ๊ฐ ํจ์์ ์ฌ์ธ ์ ๋ถ Si(x), ์ฝ์ฌ์ธ ์ ๋ถ Ci(x)์ ๋์ํ๋ ์๊ณก ๋ฒ์ ์ผ๋ก, ์ด์ ๋, ์ ํธ ํด์, ํน์ ํจ์์ ์ ๊ทผ ์ ๊ฐ ๋ฑ์์ ๋ฑ์ฅํฉ๋๋ค.
์ฌ์ฉ ๋ฐฉ๋ฒ
์ซ์ ์ธ ๊ฐ๋ฅผ ์ ๋ ฅํฉ๋๋ค. x์ ์ด๊น๊ฐ(์ฒซ ๋ฒ์งธ ํ), ์ฐ์๋ ํ ์ฌ์ด์ ์ฆ๊ฐ๋(์คํ ), ๊ทธ๋ฆฌ๊ณ ๋ฐ๋ณต ํ์(์์ฑํ ํ์ ๊ฐ์)์ ๋๋ค. ํ๋ \(i = 0\)๋ถํฐ count\(-1\)๊น์ง \(x_i = \text{startX} + i \cdot \text{stepX}\) ๊ท์น์ ๋ฐ๋ผ ๊ณ์ฐ๋ฉ๋๋ค. ์๋ฅผ ๋ค์ด ์ด๊น๊ฐ 0, ์คํ 0.5, ํ์ 3์ ์ ๋ ฅํ๋ฉด \(x = 0, 0.5, 1.0\)์์์ ๊ฐ์ด ํ์ผ๋ก ๋ง๋ค์ด์ง๋๋ค.
๊ณต์ ์ค๋ช
์ ์์ ๋ฐ๋ฅด๋ฉด $$\operatorname{Shi}(x) = \int_0^x \frac{\sinh(t)}{t}\,dt$$ ์ด๊ณ , $$\operatorname{Chi}(x) = \gamma + \ln|x| + \int_0^x \frac{\cosh t - 1}{t}\,dt$$ ์ ๋๋ค. ์ฌ๊ธฐ์ \(\gamma \approx 0.5772156649\) ๋ ์ค์ผ๋ฌ-๋ง์ค์ผ๋ก๋ ์์์ ๋๋ค. ์ด ๊ณ์ฐ๊ธฐ๋ ๋น ๋ฅด๊ฒ ์๋ ดํ๋ ๋ค์์ ๋ฉฑ๊ธ์๋ฅผ ์ฌ์ฉํด ๊ฐ์ ํ๊ฐํฉ๋๋ค. $$\operatorname{Shi}(x) = \sum_{k=0}^{\infty} \frac{x^{2k+1}}{(2k+1)\,(2k+1)!}$$ ์ด๊ณ $$\operatorname{Chi}(x) = \gamma + \ln x + \sum_{k=1}^{\infty} \frac{x^{2k}}{(2k)\,(2k)!}$$ ์ ๋๋ค. ๊ฐ ํญ์ ๋น์จ ๊ฐฑ์ ๋ฐฉ์์ผ๋ก ๋์ ๋์ด ํฉํ ๋ฆฌ์ผ ์ค๋ฒํ๋ก๋ฅผ ํผํ๋ฉฐ, ํญ์ด ๋ฌด์ํ ๋งํผ ์์์ง๋ฉด ๊ณ์ฐ์ ๋ฉ์ถฅ๋๋ค.
๊ณ์ฐ ์์
\(x = 1\)์ผ ๋: $$\operatorname{Shi}(1) = 1 + \frac{1}{18} + \frac{1}{600} + \frac{1}{35280} + \cdots \approx 1.0572509.$$ $$\operatorname{Chi}(1) = 0.5772157 + \ln 1 + \frac{1}{4} + \frac{1}{96} + \frac{1}{4320} + \cdots \approx 0.8378695.$$
์์ฃผ ๋ฌป๋ ์ง๋ฌธ
Chi(0)์ด ์ ์ ์๋์ง ์์์ผ๋ก ํ์๋๋์? Chi(x)์๋ \(\ln x\) ํญ์ด ๋ค์ด ์๋๋ฐ, \(x \rightarrow 0\)์ผ ๋ ์ด ๊ฐ์ด \(-\infty\)๋ก ๋ฐ์ฐํฉ๋๋ค. ๋ฐ๋ผ์ 0์์ Chi๋ ์ ํํ ๊ฐ์ ๊ฐ์ง์ง ์์ต๋๋ค.
์์ x๋ ์ด๋ป๊ฒ ๋๋์? Shi๋ ๊ธฐํจ์์ด๋ฏ๋ก \(\operatorname{Shi}(-x) = -\operatorname{Shi}(x)\)์ด๋ฉฐ ์ ์์ ์ผ๋ก ๊ณ์ฐ๋ฉ๋๋ค. ๋ฐ๋ฉด Chi(x)๋ \(x > 0\)์์๋ง ์ค์ซ๊ฐ์ ๊ฐ์ง๋ฉฐ(\(x < 0\)์ด๋ฉด ํ์๋ถ \(-i\pi\)๊ฐ ๋ํด์ง๋๋ค), ๋ฐ๋ผ์ ํ์์๋ \(x \le 0\)์ผ ๋ Chi๋ฅผ ์ ์๋์ง ์์์ผ๋ก ํ์ํฉ๋๋ค.
์ ํ๋๋ ์ด๋ ์ ๋์ธ๊ฐ์? \(|x|\)๊ฐ ์ ๋นํ ๋ฒ์(๋๋ต 10 ์ดํ)์์๋ ์ด ๊ธ์๊ฐ ๋ฐฐ์ ๋ฐ๋(double) ์ ์ฒด ์๋ฆฟ์๊น์ง ์ ํํ ๊ฐ์ ์ค๋๋ค. ๋ฐ๋ณต์ ๋ณดํต 20~40๊ฐ ํญ์์ ์๋ ดํฉ๋๋ค.