A new numerical method using generalized eigenvector decomposition for the solution of pulsed activation is described and tested. This method allows for the treatment of multiple occurences of isotopes in the decay scheme including the treatment of an arbitrary number of loops. Preliminary comparisons with existing activation codes show good agreement and confirm that the method produces accurate results.