Segmented tau approximation for a parametric nonlinear neutral differential equation (Q2383819)

The nonlinear differential equation \[ y'(t)=ry(t)(a+y(t-\tau)+cy'(t-\tau)), t\geq 0, \] \[ y(t)=\psi(t), t\leq 0 \] is approximated by a sequence of equations of the same form in intervals \( I_{k}=[k\tau,(k+1)\tau], k\geq 0.\) For each \(k\) the solution is denoted \(y_{k}(t)\) and is found in polynomial form with the condition \(y_{k}(0) =y_{k-1}(1)\). The numerical results of this method are compared for three examples with results of the Runge-Kutta method.