Template:cpp/numeric/special functions/dsc expint