Linear time-invariant ODEs with nil potent matrices yield polynomial solutions. It might be faster to calculate those instead of invoking a matrix exponential. - View it on GitHub