Abstract Harmonic numbers and generalized harmonic numbers have been studied since the distant past and involved in a wide range of diverse fields such as analysis of algorithms in computer science, various branches of number theory, elementary particle physics and theoretical physics. Here we aim at presenting further interesting identities about certain finite or infinite series involving harmonic numbers and generalized harmonic numbers by applying an algorithmic method to a known summation formula for the hypergeometric function 5 F 4 (1).