The purpose of this paper is to study solvability of the higher order nonlinear neutral delay differential equation d n d t n [ x ( t ) + c ( t ) x ( t − τ ) ] + ( − 1 ) n + 1 f ( t , x ( σ 1 ( t ) ) , x ( σ 2 ( t ) ) , … , x ( σ k ( t ) ) ) = g ( t ) , t ≥ t 0 , $$\begin{aligned}& \frac{d^{n}}{dt^{n}}\bigl[x(t)+c(t)x(t-\tau)\bigr]+(-1)^{n+1}f\bigl(t,x \bigl(\sigma _{1}(t)\bigr),x\bigl(\sigma_{2}(t)\bigr), \ldots,x\bigl(\sigma_{k}(t)\bigr)\bigr) \\ & \quad =g(t),\quad t\geq t_{0}, \end{aligned}$$ where n and k are positive integers, τ > 0 , t 0 ∈ R , f ∈ C ( [ t 0 , + ∞ ) × R k , R ) $f\in C ([t_{0},+\infty)\times {\mathbb{R}}^{k},{\mathbb{R}} ) $ , c , g , σ i ∈ C ( [ t 0 , + ∞ ) , R ) $c,g,\sigma_{i}\in C([t_{0},+\infty),{\mathbb{R}})$ and lim t → + ∞ σ i ( t ) = + ∞ $\lim_{t\rightarrow +\infty}\sigma_{i}(t)=+ \infty$ for i ∈ { 1 , 2 , … , k } . Under suitable conditions, several existence results of uncountably many nonoscillatory solutions and convergence of Mann iterative approximations for the above equation are shown. Three nontrivial examples are given to demonstrate the advantage of our results over the existing ones in the literature.