Matrix-exponential discretization

Computes exact linear sampled models under zero-order-hold assumptions.