SymPy currently only supports an extension of 0 (you can simulate other endpoints by performing a shift). It does not support expansion at infinity, which is used in algorithmic analysis.
But it would be a good base package for it, and if you implement it, we will gladly accept the patch (nb: I am the main developer of SymPy).
Remember that the whole problem is complex, especially if you have two variables or even symbolic constants. It is also difficult if you want to support oscillation functions. EDIT: If you're interested in oscillating features, this SymPy mailing list discussion contains interesting articles.
EDIT 2: And I would recommend not trying to build it yourself from scratch, without using a computer algebra system. You will eventually have to write your own computer algebra system, which is a lot of work and even more work if you want to do it right and not have it slow. There are already many existing systems, including many that can act as libraries for code to be built on top of them (for example, SymPy).
source share