The exponential matrix technique is used to solve the point kinetic equations of a reactor. This uses a 7th-order accurate Pade(3,3) approximation. This is stable and independent of the time-step size.
[[MELCOR]] – accident analysis code that uses this technique
el-genkPhysicsbasedDynamicModel