Solving point burnup equations by Magnus method

The burnup equation of nuclides is one of the most equations in nuclear reactor physics, which is generally coupled with transport calculations. The burnup equation describes the variation of the nuclides with time. Because of its very stiffness and the need for large time step, this equation is solved by special methods, for example transmutation trajectory analysis (TTA) or the matrix exponential methods where the matrix exponential is approximated by CRAM.