Kürkçü, Ömür Kıvanç2023-06-162023-06-1620202008-13592251-7456https://doi.org/10.1007/s40096-019-00314-8https://hdl.handle.net/20.500.14365/1019In this paper, we introduce a numerical method to obtain an accurate approximate solution of the integro-differential delay equations with state-dependent bounds. The method is based basically on the generalized Mott polynomial with the parameter-beta Chebyshev-Lobatto collocation points and matrix structures. These matrices are gathered under a unique matrix equation and then solved algebraically, which produce the desired solution. We discuss the behavior of the solutions, controlling their parameterized form via beta and so we monitor the effectiveness of the method. We improve the obtained solutions by employing the Mott-residual error estimation. In addition to comparing the results in tables, we also illustrate the solutions in figures, which are made up of the phase plane, logarithmic and standard scales. All results indicate that the present method is simple-structured, reliable and straightforward to write a computer program module on any mathematical software.eninfo:eu-repo/semantics/openAccessCollocation pointsError estimationMatrix methodMott polynomialDelayDifferential EquationCollocation MethodDicksonTaylorModelArticle10.1007/s40096-019-00314-82-s2.0-85112664858