A new investigation of the extended Krylov subspace method for matrix function evaluations