We describe an algorithm that uses arithmetic operations for computing the determinant of the matrix , where is the adjacency matrix of an order tree. Combining this algorithm with interpolation, we derive a simple algorithm requiring arithmetic operations to find the characteristic polynomial of the adjacency matrix of any tree. We apply our algorithm and recompute a 22-degree characteristic polynomial, which had been incorrectly reported in the quantum chemistry literature.