An investigation is made of a new approach to solving the multimode kinetics equations. These equations are cast into a matrix generalization, , of the point kinetics equations and successfully solved by applying an analytic inversion of the matrix [[I] - ε[A]] to both the Padé (2,0) and modified Padé (1,1) approximations to exp(Δt[A]).