Highly Accurate Method for Boundary Value Problems with Robin Boundary Conditions

The main aim of the current paper is to construct a numerical algorithm for the numerical solutions of second-order linear and nonlinear differential equations subject to Robin boundary conditions. A basis function in terms of the shifted Chebyshev polynomials of the first kind that satisfy the homogeneous Robin boundary conditions is constructed. It has established operational matrices for derivatives of the constructed polynomials. The obtained solutions are spectral and are consequences of the application of collocation method. This method converts the problem governed by their boundary conditions into systems of linear or nonlinear algebraic equations, which can be solved by any convenient numerical solver. The theoretical convergence and error estimates are discussed. Finally, we support the presented theoretical study by presenting seven examples to ensure the accuracy, efficiency, and applicability of the constructed algorithm. The obtained numerical results are compared with the exact solutions and results from other methods. The method produces highly accurate agreement between the approximate and exact solutions, which are displayed in tables and figures.