Any order imaginary time propagation method for solving the Schrödinger equation