Construction of Exact Parametric or Closed Form Solutions of Some Unsolvable Classes of Nonlinear ODEs ( Abel ’ s Nonlinear ODEs of the First Kind and Relative Degenerate Equations )

We provide a new mathematical technique leading to the construction of the exact parametric or closed form solutions of the classes of Abel’s nonlinear differential equations ODEs of the first kind. These solutions are given implicitly in terms of Bessel functions of the first and the second kind Neumann functions , as well as of the free member of the considered ODE; the parameter ν being introduced furnishes the order of the above Bessel functions and defines also the desired solutions of the considered ODE as one-parameter family of surfaces. The nonlinear initial or boundary value problems are also investigated. Finally, introducing a relative mathematical methodology, we construct the exact parametric or closed form solutions for several degenerate Abel’s equation of the first kind.


Introduction
One of the most known nonlinear differential equations ODEs is the Abel equation of the first kind.There are admissible functional transformations that can reduce the general first kind nonlinear Abel ODE u ξ f 3 ξ u 3 f 2 ξ u 2 f 1 ξ u f 0 ξ to an Abel equation of the first kind of the normal form y x y 3 F x .None of these types of equations admits general solutions in terms of known tabulated functions, since only very special cases depending on the form of the free member F x or on f i ξ , i 0, 1, 2, 3 can be solved in parametric form 1-3 .We underline here that adding a cubic term in a Riccati ODE , an Abel nonlinear ODE of the first kind arises.
Our motivation of studying this problem was to provide a mathematical methodology through which the construction of exact parametric or closed form solutions of an Abel International Journal of Mathematics and Mathematical Sciences first kind nonlinear ODE form is possible.According to this method, separating properly the terms of the considered ODE, one leads to an equivalent system of two equations.The first one is a Riccati ODE, and the second one is a cubic equation that must be simultaneously hold true.In other words, any solution which satisfies these two equations the Riccati ODE and the cubic one satisfies also the Abel equation of the first kind of the normal form.Using basic results given by Kamke 1 , the above Riccati equation is transformed to a second-order linear ODE of variable and determinable coefficients.We succeed in defining the coefficients of the above second-order linear ODE, so that the prescribed cubic equation hold true, and in the sequel in solving exactly this secondorder linear ODE.This introducing procedure guides to the construction of the abovementioned exact parametric or closed form solutions of the initial Abel equation.Here, we underline that all the above introduced functional transformations are admissible and that several basic results concerning the solutions of such types of equations are included in Zarmpoutis 4 .A nonlinear initial or boundary value problem is investigated under the assumption that the constructed solution of the Abel equation is unique inside a main interval of the independent variable.Finally, relative to the above mathematical methodology, we succeed in constructing the exact or closed form solutions for several degenerate Abel's nonlinear ODEs of the first kind special cases which are introduced in 5-7 .

Analysis-Construction of the General Solutions
An Abel equation of the first kind has the general form Here, the notation y x dy/dx, y xx d 2 y/dx 2 , . . . is used for the total derivatives.The admissible functional transformations 3, Chapter 0.1.5-2: brings the general form 2.1 in the canonical normal form:

2.3b
International Journal of Mathematics and Mathematical Sciences 3 The first step of our new mathematical method is to separate properly the terms of the Abel equation 2.3a in order to find a more convenient system of equations.In more details, 2.3a can be rewritten as such that the solution of 2.3a identically satisfies the equation Introducing a subsidiary, smooth function Q x , one splits 2.5 into the following two equivalent equations:

2.6b
The first of them 2.6a is a Riccati ODE , while the second one 2.6b is a cubic equation.In fact, the elimination of Q x between 2.6a and 2.6b leads to 2.4a .This means that any solution of the Riccati ODE 2.6a , that satisfies the restriction 2.6b , is also a solution of the Abel equation 2.4a and 2.4b .It is well known that any Riccati equation can be transformed to a second-order linear ODE of variable coefficients 1 .
The second step of our method is to construct a solvable second-order linear ODE of variable coefficients, which is equivalent to the Riccati 2.6a that satisfies the restriction 2.6b .A kind of solvable Riccati equation of the normal form 3, pp 105; type 38 is of the type which admits a particular solution of the form where f * and g * are arbitrary functions that must be determined.

International Journal of Mathematics and Mathematical Sciences
Supposing that the Riccati equation 2.6a can be written in the form 2.7 , in order to solve 2.6a , it is necessary to define the above-mentioned f * and g * .Thus, one reads so that 2.6a admits a particular solution of the type 2.8 .From now on, differentiating 2.9b , we extract while combining the last equation with 2.9a and 2.9b , one gets By now, the system 2.9a and 2.9b becomes Introducing a second subsidiary differentiable function Ω x such that 13 system 2.12 can be rewritten as

2.14b
According to 1, p 441; type 15 , the Bessel second-order linear ODE where G x is an arbitrary smooth function and μ, ν are arbitrary parameters, admits the two partial linearly independent solutions are the Bessel functions of the first and the second kind Neumann function of order v, respectively.Thus, the general solution of 2.15 is given by Comparing now 2.14a with 2.15 , one finds with

2.17c
Here, y x is the general solution of the Riccati equation 2.6a under the restriction 2.6b .Integration of 2.17a furnishes while inserting the result 2.18 into the second equation 2.17b , one finds By the assumption μ 0, the Bessel equation 2.15 together with 2.14a , 2.13 , 2.6b , and 2.18 becomes Here, g * p 1 , g * p 2 are two partial linearly independent solutions of the Bessel ODE 2.14a .

International Journal of Mathematics and Mathematical Sciences
By now, using equations 2.9b , 2.8 , and 2.18 , the restriction 2.6b as well as the properties of the Bessel functions 3, 8, 9 , we get the estimation as in equation 2.19 for μ 0.

2.21b
Thus, the two partial linearly independent solutions of the Riccati equation 2.6a y p 1 , y p 2 under the restriction 2.6b become where

2.23
Furthermore, according to 3 , the general solution of the above Riccati equation 2.6a under the restriction 2.6b is given by the formula C integration constant.

2.24b
Solving the first of them for U and taking the derivatives in both 2.21a and 2.21b , we obtain the result

2.25
Since any solution of the Riccati equation 2.6a under the restriction 2.6b is also a solution of the Abel equation 2.3a , we conclude that the general solution 2.24a of the Riccati 2.6a under the restriction 2.6b is also the general solution of the Abel equation 2.3a .Taking into account the relations where F x is the free member of the considered Abel equation, while C is an integration constant and v is a free parameter denoting the order of Bessel functions being introduced.
From the above analysis, one concludes that the implicit functional relation 2.28 expressing the general solution of the Abel equation of the first kind of the normal form 2.3a may not be unique inside a main interval of the values x 0 , x 1 of the independent variable x, because of various types that function G x 2.29 admits; in other words, it can be divided into several branches of solutions valid inside consecutive subintervals.In this case, matching of the corresponding solutions including various integration constants C as well as various values of the parameter v must be followed in all the consecutive subintervals that solution 2.28 changes.

Nonlinear Initial or Boundary Value Problem
We suppose that solution 2.28 remains unchanged inside the first main interval x 0 , x 1 .Following the initial data x 1 x 0 1 0, y 1 y 0 1 , we provide by 2.29 where the first upper right index denotes the value of the corresponding quantity at the origin x 0 0, y y 0 .Thus, for the origin 0, y 0 1 , we dispose the following nonlinear algebraictranscendental system of equations Equation 2.28

3.3
The set of equations 3.2a and 3.2b is sufficient for the estimation of the constant of integration C 0 as well as for the parameter v 0 , so that the one-parameter family of surfaces solution of the original Abel equation takes the concrete form with given v v 0 and C C 0 .

Applications and Discussion
In this section, we will examine some special cases of degenerate Abel's nonlinear ODEs of the first kind that cover all the examined applications 5-7 .A modified mathematical methodology is introduced that simplify the mathematical methodology and technique being introduced in previous sections.We will investigate and discuss two degenerate type of equations, that is, y x f x y 3 g x y 2 h x or y x f x y 3 g x y 2 .

4.1
Consider the first of these equations introducing h x 0 .From the second of these, setting which, further, using the functional transformation 3, page 50 International Journal of Mathematics and Mathematical Sciences 9 results in the simpler form This new equation is an Abel equation of the second kind that, further, by the substitution 3, page 40 becomes an Abel equation of the second kind of the normal form, that is, This last equation by the admissible substitution reduces to the following Abel degenerate nonlinear ODE of the first kind: 4.9 Using the final admissible functional transformation 1 , tξ t / 0, 4.10 performing we transform 4.9 to the new nonlinear ODE of the second kind of the Emden-Fowler type, namely, where F ξ is a known smooth function of the variable ξ.

International Journal of Mathematics and Mathematical Sciences
By now, one uses the mathematical methodology and the results are included in 10 .In fact, introducing the functional transformation and supposing that G z ≡ F ξ P 3 /ξ − 2zP , we derive the following nonlinear ODE when G z is a function which must be determined: Thus, the ODE 4.14 becomes the following first-order Abel's nonlinear type of equations −z 2 P G P z P ; P P z .

4.17
We apply now the Julia construction 1, page 27 , and we derive the following functional relation between its coefficients: λ, μ free parameters constants of integration .

4.23
Furthermore, the dependent variable ξ can be estimated by the combination of 4.20 together with the second of 4.16 , namely, by means of which one gets Since F ξ is a known function of ξ, the last equation permits us to evaluate ξ with respect to the independent variable t and thus functions w ξ and y ξ through relations 4.10 and 4.8 , respectively.Following the inverse course, one finally defines y x by the way of 4.1 .Summarizing, we underline that, in this section, by using the results of 10 , we proposed a new mathematical methodology and we succeed in giving new mathematical techniques in constructing analytic or closed form solutions of a class of degenerate Abel's ODEs in mathematical physics and nonlinear mechanics.The methodology being introduced simplifies that of Section 2, which refers to the solution of an Abel equation of the first kind.For strengthening, we refer to some categories of degenerate Abel's ODEs of the first kind which are included in 5-7 , namely, i ax b 2 y x ax b y 3 cy 3 0; ii y x cy 2 0; iii x 2 y x ay 3 − ax 2 y 2 0; iv y x 3ax 2 4x b y 3 3xy 2 0, and so forth.

Conclusion
A new mathematical technique leading to the construction of the exact parametric or closed form solutions of the classes of Abel's nonlinear differential equations ODEs of the first kind is presented.The solutions are given implicitly using Bessel's functions of the first and second kind Neumann functions , the free member of the considered ODE as well as the subindex v, that defines the order of the Bessel functions.The nonlinear initial or boundary value problems are also investigated.Finally, through a new simpler mathematical methodology, we succeed in giving the exact parametric solutions of some degenerate cases of Abel's equation of the first kind included and investigated in 1, 5, 7 .
μ free parameters constants of integrations .International Journal of Mathematics and Mathematical Sciences that, because F ξ is a known function of ξ, furnishes ξ as function of t.From now on, one follows the inverse course.