Move all matrix generation functionality into one place #57
Description
Right now, the functions to generate the stability function and matrix representation for SDC are copy&pasted at different places (tests, plot_stability). They should all me moved to one place.