Unwrapping the phase response functions for non-linear systems