Improved series expansion computation, and exposure to API.

The functions now switch between two different series expansion, depending
on the given t. This switching based on Navarro & Fuss (2009).

The series expansion is implemented for both symmetric and asymmetric
bounds. Only the symmetric case is available in the MATLAB interface.