Functional differential equations of neutral type are a class of differential equations in which the derivative of the unknown functions depends on the history of the function and its derivative as well. Due to this nature the explicit solutions of these equations are not easy to compute and sometime even not possible. Therefore, one must use some numerical technique to find an approximate solution to these equations. In this paper, we used a spectral collocation method which is based on Bernstein polynomials to find the approximate solution. The disadvantage of using Bernstein polynomials is that they are not orthogonal and therefore one cannot use the properties of orthogonal polynomials for the efficient evaluation of differential equations. In order to avoid this issue and to fully use the properties of orthogonal polynomials, a change of basis transformation from Bernstein to Legendre polynomials is used. An error analysis in infinity norm is provided, followed by several numerical examples to justify the efficiency and accuracy of the proposed scheme.
Citation: Ishtiaq Ali. Bernstein collocation method for neutral type functional differential equation[J]. Mathematical Biosciences and Engineering, 2021, 18(3): 2764-2774. doi: 10.3934/mbe.2021140
Functional differential equations of neutral type are a class of differential equations in which the derivative of the unknown functions depends on the history of the function and its derivative as well. Due to this nature the explicit solutions of these equations are not easy to compute and sometime even not possible. Therefore, one must use some numerical technique to find an approximate solution to these equations. In this paper, we used a spectral collocation method which is based on Bernstein polynomials to find the approximate solution. The disadvantage of using Bernstein polynomials is that they are not orthogonal and therefore one cannot use the properties of orthogonal polynomials for the efficient evaluation of differential equations. In order to avoid this issue and to fully use the properties of orthogonal polynomials, a change of basis transformation from Bernstein to Legendre polynomials is used. An error analysis in infinity norm is provided, followed by several numerical examples to justify the efficiency and accuracy of the proposed scheme.
[1] | R. Farouki, Legendre-Bernstein basis transformations, J. Comput. Appl. Math., 119 (2000), 145–160. doi: 10.1016/S0377-0427(00)00376-9 |
[2] | R. Farouki, T. Goodman, T. Sauer, Construction of orthogonal bases for polynomials in Bernstein form on triangular and simplex domains, Comput. Aided Geom. Des., 20 (2003), 209–230. |
[3] | K. Höllig, J. Hörner, Approximation and Modeling with B-Splines, Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 2013. |
[4] | G. Farin, Curves and Surface for Computer Aided Geometric Design, Academic Press: Boston, MA, USA, (1993), 32–58. |
[5] | R. Farouki, V. Rajan, Algorithms for polynomials in Bernstein form, Comput. Aided Geom. Des., 5 (1988), 1–26. |
[6] | K. Höllig, J. Hörner, Approximation and Modelling with B-Splines, SIAM, Philadelphia, PA, USA, 132 (2013), 32–58. |
[7] | Y. Liu, Numerical solution of implicit neutral functional differential equations, SIAM J. Nume. Anal., 36 (1999), 516–528. doi: 10.1137/S003614299731867X |
[8] | I. Ali, H. Brunner, T. Tang, A spectral method for pantograph-type delay differential equations and its convergence analysis, J. Comput. Math., 27 (2009), 254–265. |
[9] | I. Ali, H. Brunner, T. Tang, Spectral methods for pantograph-type differential and integral equations with multiple delays, Front. Math. China, 4 (2009), 49–61. doi: 10.1007/s11464-009-0010-z |
[10] | C. Canuto, M. Hussaini, A. Quarteroni, T. A. Zang, Spectral Methods: Fundamentals in Single Domains, Springer, Berlin, 2006. |
[11] | J. Shen, T. Tang, Spectral and High-Order Methods with Applications, Science Press, Beijing, 2006. |
[12] | H. Brunner, Collocation Methods for Volterra Integral and Related Functional Equations, Cambridge University Press, Cambridge, 2004. |
[13] | H. Brunner, Q. Hu, Optimal superconvergence results for delay integro-differential equations of pantograph type, SIAM J. Numer. Anal., 45 (2007), 986–1004. doi: 10.1137/060660357 |
[14] | H. Brunner, Q. Hu, Q. Lin, Geometric meshes in Collocation Methods for Volterra Integral with proportional delay, IMA J. Numer. Anal., 21 (2001), 783–798. doi: 10.1093/imanum/21.4.783 |
[15] | A. Bellen, M. Zennaro, Numerical Methods for Delay Differentials Equations, Oxford University Press, Oxford, 2003. |
[16] | A. Bellen, Preservation of superconvergence in the numerical integration of delay differential equations with proportional delay, IMA J. Numer. Anal., 22 (2002), 529–536. doi: 10.1093/imanum/22.4.529 |
[17] | A. Bataineh, O. Işik, N. Aloushoush, et al., Bernstein operational matrix with error analysis for solving high order delay differential equations, Int. J. Appl. Comput. Math., 3 (2017), 1749–1762. doi: 10.1007/s40819-016-0212-5 |
[18] | P. Sahu, R. Saha, A new numerical approach for the solution of nonlinear Fredholm integral equations system of second kind by using Bernstein collocation method, Math. Methods Appl. Sci., 38 (2015), 274–280. doi: 10.1002/mma.3067 |
[19] | P. Sahu, R. Saha, Legendre spectral collocation method for the solution of the model describing biological species living together, J. Comput. Appl. Math., 296 (2016), 47–55. doi: 10.1016/j.cam.2015.09.011 |
[20] | M. Bhatti, P. Bracken, Solutions of differential equations in a Bernstein polynomial basis, J. Comput. Appl. Math., 205 (2007), 272–280. doi: 10.1016/j.cam.2006.05.002 |
[21] | G. Mastroianni, D. Occorsio, Optimal systems of nodes for Lagrange interpolation on bounded intervals: A survey, J. Comput. Appl. Math., 134 (2001), 325–341. doi: 10.1016/S0377-0427(00)00557-4 |
[22] | A. Iserles, On nonlinear delay differential equations, Trans. Amer. Math. Soc., 344 (1994), 441–447. doi: 10.1090/S0002-9947-1994-1225574-4 |
[23] | D. Trif, Direct operational tau method for pantograph-type equations, Appl. Math. Comput., 219 (2012), 2194–2203. |
[24] | E. Ishiwata, Y. Muroya, Rational approximation method for delay differential equations with proportional delay, Appl. Math. Comput., 187 (2007), 741–747. |
[25] | I. Ali, S. Khan, Analysis of stochastic delayed SIRS model with exponential birth and saturated incidence rate, Chaos, Solitons Fractals, 138 (2020), 110008. doi: 10.1016/j.chaos.2020.110008 |
[26] | S. Khan, I. Ali, Applications of Legendre spectral collocation method for solving system of time delay differential equations, Adv. Mech. Eng., 12 (2020), 1–13. |
[27] | S. Khan, M. Ali, I. Ali, A spectral collocation method for stochastic Volterra integro-differential equations and its error analysis, Adv. Differ. Equations, 1 (2019), 161. |
[28] | S. Khan, I. Ali, Application of Legendre spectral-collocation method to delay differential and stochastic delay differential equation, AIP Adv., 8 (2018), 035301. doi: 10.1063/1.5016680 |
[29] | O. Isik, Z. Güney, M. Sezer, Bernstein series solutions of pantograph equations using polynomial interpolation, J. Differ. Equations Appl., 18 (2012), 357–374. doi: 10.1080/10236198.2010.496456 |
[30] | A. Romero, P. Galvín, J. Cámara-Molina, A. Tadeu, On the formulation of a BEM in the Bezíer-Bernstein space for the solution of Helmholtz equation, Appl. Math. Modell., 74 (2019), 301–319. doi: 10.1016/j.apm.2019.05.001 |
[31] | A. Romero, P. Galvín, A. Tadeu, An accurate treatment of non-homogeneous boundary conditions for development of the BEM, Eng. Anal. Boundary Elem., 116 (2020), 93–101. doi: 10.1016/j.enganabound.2020.04.008 |