A new class of polynomial functions equipped with a parameter

In this study, a new class of polynomial functions although equipped with a parameter is introduced. This class can be employed for computational solution of linear or non-linear functional equations, including ordinary differential equations or integral equations. The extra parameter permits us to obtain more accurate results. In the present paper, a number of numerical examples show the ability of this class of polynomial functions.

In many problems, utilizing polynomial functions can assist us in finding reasonable solutions. For example, in physics, engineering, and economics, we can find many computational methods applying polynomial functions. A list of polynomial functions and their applications can be found in [2]. In recent years, a great number of new families of polynomials are introduced, including q-analogue of Hermite polynomials [6], d-orthogonal polynomials by Cheikha and Romdhane [5], unified family of generalized Apostol-Bernoulli, and Euler and Genocchi polynomials by El-Desouky and Gomaa [7].
In numerical solution of an Ordinary Differential Equation (ODE), we usually put a linear combination of polynomial functions with unknown coefficients in the main equation and then find the unknowns. This can be done using collocation or, Galerkin or Tau methods. As we know, in this scheme, we can select a set of orthogonal or orthonormal polynomial functions for the solution.
In the current research, we aim to introduce a new class of polynomial functions equipped with a parameter that can be used for computational solution of linear or non-linear functional equations, including ordinary differential equations, partial differential equations, or integral equations. The extra parameter contributes to obtaining more accurate results. Numerical examples demonstrate the ability of this class of polynomial functions to obtain better solutions.

Polynomial functions
In this section, a new class of polynomial functions is introduced. This definition is based on the Chebyshev polynomial of the second kind, U n ðxÞ. These functions have already been employed. For example, Boubaker polynomials are defined as B n ðxÞ ¼ U n ð x 2 Þ þ 3U nÀ2 ð x 2 Þ, for n ! 2. Furthermore, the modified Boubaker polynomials [13] are defined asB n ðxÞ ¼ 6xU nÀ1 ðxÞ À 2U n ðxÞ, for n ! 1.
Definition 2.1 Suppose that a is an arbitrary constant and U n ðxÞ is the Chebyshev polynomial of the second kind. Let A 0 ðxÞ ¼ 1 for any n ! 1. The new class of polynomial functions is defined by the following: A n ðxÞ ¼ axU nÀ1 ðxÞ þ U n ðxÞ: We can find very easily that A 1 ðxÞ ¼ð2 þ aÞx; . . . A nþ1 ðxÞ ¼2xA n ðxÞ À A nÀ1 ðxÞ: We can see that for n ! 2, A n ðxÞ ¼ ð1 þ a 2 ÞU n ðxÞþ a 2 U nÀ2 ðxÞ: Proposition 2.1 Let T n ðxÞ be the Chebyshev polynomial of the first kind. Then, A n ðxÞ ¼ ða þ 1Þ x n T 0 n ðxÞ þ T n ðxÞ for n ! 1.

Proposition 2.2
Let b:c designate the floor function. Then, for n ! 1, The polynomial A n ðxÞ for n ! 0 satisfies the following equation: For n ! 0, the following equations are obtained: Proposition 2.5 By considering the inner product we have ðA n ðxÞ; A m ðxÞÞ ¼ Corollary 2.1 The system fA 4n ðxÞg 1 n¼0 is an orthogonal system of polynomial functions, with respect to the inner product defined in Proposition 2.5.
The mth-order non-linear differential equation In this section, the new class of polynomial functions is considered to find the numerical solutions of some nonlinear differential equations. The solutions obtained show the efficiency of this class of functions. In this class, we have an unknown parameter, and by selecting it, we can reduce the error of PES. In all of the following examples, the least square method is applied so as to find the unknown parameter. All iterative algorithms, which we needed for our initial guess, were initialized at zero.
The mth-order non-linear differential equation with the mixed boundary conditions X mÀ1 k¼0 À a j;k y ðkÞ ðaÞ þ b j;k y ðkÞ Á ¼ c j ; j ¼ 0; 1; . . .; m À 1; ð2Þ where y is an unknown function, g and P k;r are continuous functions, and a j;k , b j;k , and c j are real or complex constants, is considered in [13]. By PFS collocation method, we can solve this problem. In this method, we approximate the solution of (1) in the truncated series form: where c k , k ¼ 0; 1; . . .; N are the unknown coefficients. In this method, N ! m is selected in any positive integer number. The matrix form of the collocation method is discussed in [13]. In this method, a non-linear system with N þ 1 equations is constructed to find N þ 1 unknowns. However, here, we have N þ 2 unknowns, i.e., c k , k ¼ 0; 1; . . .; N and the unknown parameter of a, of the new class of polynomial functions of A k . By adding one collocation point to the set of points, we can construct a nonlinear system with N þ 2 equations. Via the least square method, we can minimize the L 2 norm of the residual in the augmented non-linear system to obtain a. Afterwards, we can continue the matrix collocation method to find other unknowns like [13].
Example 3.1 We consider the Riccati differential equation: y 0 ðxÞ ¼ yðxÞ À 2y 2 ðxÞ; 0 x 1; with the initial condition yð0Þ ¼ 1, [4,9,13], which is a special case of (1). The exact solution is yðxÞ ¼ 1=ð2 À e Àx Þ. Table 1 indicates the absolute error e N ðxÞ ¼ jyðxÞ À y N ðxÞj for different methods, including the Taylor method [4], the decomposition method [9], the PES by Bessel polynomials [13], and the present method. Table 2 demonstrates a comparison between the new class of polynomial functions in present method and Boubaker polynomials, pointing to the efficiency of the new class of polynomial functions.
Example 3.2 We consider the non-linear Abel differential equation of the second kind: yðxÞy 0 ðxÞ þ xyðxÞ þ y 2 ðxÞ þ x 2 y 3 ðxÞ ¼ xe Àx þ x 2 e À3x ; 0 x 1; with the initial condition yð0Þ ¼ 1, [8,13], which is special case of (1). The exact solution is yðxÞ ¼ e Àx . Table 3 shows the absolute error e N ðxÞ for different methods. The Taylor method [8] in first column, the PES by Bessel polynomials in second column, the PES by Boubaker polynomials in third column, and the present method in last column are shown. Table 4 shows a comparison with Boubaker polynomials, and so shows the efficiency of the new class of polynomial functions.

Conclusion
In this study, a new class of polynomial functions equipped with a parameter was introduced. We found that the extra parameter permits us to obtain more accurate results. By solving Riccati and Abel non-linear differential equations, the ability and efficiency of this class of polynomial functions were supported. The application of this class of polynomials in solving complicated ordinary differential equations and partial differential equations and integral equations is expected.