Consider using different approach to calculating time-dependent quantities

[smith-garrett]

Currently using numerical solution of matrix exponential. Could also consider:

• Uniformization
• Erlangization (mandjes2022numerical)

Might be faster and/or better for automatic differentiation? More numerically stable?