Add _eval_nseries method for expint()
Some tests for expint series
Merge pull request #2505 from skirpichev/2958-series-for-error_functions.py

Implement missing series methods for exponential integrals (first part of the issue 2958)