std::expint, std::expintf, std::expintl
From cppreference.com
                    
                                        
                    < cpp | numeric | special math
                    
                                                            
                    | double      expint( double arg ); double      expint( float arg ); | (1) | (since C++17) | 
| double      expint( Integral arg ); | (2) | (since C++17) | 
4) A set of overloads or a function template accepting an argument of any integral type. Equivalent to (1) after casting the argument to double.
| Contents | 
[edit] Parameters
| arg | - | value of a floating-point or Integral type | 
[edit] Return value
If no errors occur, value of the exponential integral ofarg, that is -∫∞-arg
| e-t | 
| t | 
[edit] Error handling
Errors may be reported as specified in math_errhandling
- If the argument is NaN, NaN is returned and domain error is not reported
- If the argument is ±0, -∞ is returned
[edit] Notes
Implementations that do not support C++17, but support TR 29124, provide this function if __STDCPP_MATH_SPEC_FUNCS__ is defined by the implementation to a value at least 201003L and if the user defines __STDCPP_WANT_MATH_SPEC_FUNCS__ before including any standard library headers.
Implementations that do not support TR 29124 but support TR 19768, provide this function in the header tr1/cmath and namespace std::tr1
An implementation of this function is also available in boost.math
[edit] Example
Run this code
Output:
Ei(0) = -inf Ei(1) = 1.89512 Gompetz constant = 0.596347
[edit] External links
Weisstein, Eric W. "Exponential Integral." From MathWorld--A Wolfram Web Resource.