“…Therefore, their analytical and numerical solutions have been an interesting area of research for many years. Recently, several numerical methods such as Crank-Nicolson Hermite cubic orthogonal spline collocation method [1], polynomial-based mean weighted residuals methods [2], finite difference schemes [3][4][5][6][7], Laplace transform method [8], finite element method [9], wavelets method [10], Ritz-Galerkin method [11], spectral meshless radial point interpolation method [12], -method [13], implicit Euler method [14] and operational approach of the Tau method [15] have been proposed to solve systems including nonlocal boundary conditions. It is well known that, Bernstein polynomials (BPs) play an important role as basis functions for various numerical techniques for solving different mathematical systems.…”