Exploring the topological physics of phonons is fundamentally important for understanding various practical applications. Here, we present a density-functional perturbation theory and finite displacement supercell based Python 3 software package called PH-NODE for efficiently computing phonon nodes present in real material through a first-principles approach. The present version of the code is interfaced with the WIEN2k, Elk, and ABINIT packages. In order to benchmark the code, six different types of materials are considered, which include (i) FeSi, a well-known double-Weyl point; (ii) LiCaAs, a half-Heusler single-type-I Weyl topological phonon (TP); and (iii) ScZn, coexisting nodal-line and nodal-ring TPs; (iv) TiS, six pairs of bulk Weyl nodes; (v) CdTe, type-II Weyl phonons; (vi) CsTe, coexisting TP and quadratic contact TP. In FeSi, the node points are found at Γ(0,0,0) and R(0.5,0.5,0.5) high symmetric points. Also, there are 21 energy values at which the node points are situated, corresponding to the full Brillouin zone. For LiCaAs, the previously reported literature claims that there is a node point along the W-X high symmetry direction between the highest longitudinal acoustic and the lowest transverse optical branch, while in our DFT calculations, a gap of 0.17 meV is found. Furthermore, ScZn hosts six nodal-ring TPs phonons at the boundary planes of the Brillouin zone in the vicinity of the M high-symmetric point. In addition to this, straight-line TPs are also found along the Γ-X and Γ-R high symmetric directions. Moreover, for TiS, six Weyl node points (WP1, WP2, WP3, WP4, WP5 and WP6) are found along H-K high-symmetric direction. In CdTe, it is found that Weyl points are located along the X-W high-symmetry direction. In the case of CsTe, a TP and a quadratic contact TP are found along the Γ-X direction and at the R high-symmetry point, respectively. The results obtained from the PH-NODE code are in good agreement with the experimentally and theoretically reported data for each material. Program summaryProgram title: PH-NODECPC Library link to program files:https://doi.org/10.17632/sjydzn49nw.1Licensing provisions: GNU General Public License 3.0Programming language: Python 3External routines/libraries: Math, Time, NumPy, SciPyNature of problem: Searching for the phonon-node points corresponding to the given number of phonon-branch using Nelder-Mead's simplex approach.Solution method: We present a density-functional perturbation theory and finite displacement supercell based Python 3 software package called PH-NODE for efficiently computing phonon nodes present in real material through a first-principles approach.