Folowing Changes:

Checks for unused variables
Removes unwanted return statements
Returns NonElementaryIntegral instead of None
Nitpicks in cds.py
Docstring for coupled_differential_system