Fixed matrix exponential for filter design.

- The previous implementation did not work on
  non-diagonalizable matrices; the new one does.
- Also expanded the test of `cont2discrete`,
  and made it not require Scipy.