Linear transformation method for multinuclide decay calculation

A linear transformation method for generic multinuclide decay calculations is presented together with its properties and implications. The method takes advantage of the linear form of the decay solution N(t)=F(t)N0, where N(t) is a column vector that represents the numbers of atoms of the radioactive nuclides in the decay chain, N0 is the initial value vector of N(t), and F(t) is a lower triangular matrix whose time-dependent elements are independent of the initial values of the system.