Reference: Ci

Ci(x) is the cosine integral.

Ci(x) = γ + ln(x) + ∫(cos(t) − 1)/t) dt from 0 to x where

γ = 0.57721 56649 01532... is the Euler–Mascheroni constant.