The numerical solution of the point kinetics equation using the matrix exponential method

► A numerical solution of the stiff differential equations using a matrix exponential method. ► Each sampling time interval is divided into two subintervals to be considered. ► Various input examples applied to the point kinetics equations are studied using. ► The method shows superior useful and effective activity.