From some point of view, after Riccati type equations, the simplest first order ordinary differential equations (ODEs) are those having as right hand side (RHS) a third degree polynomial in the dependent variable, also called Abel type ODEs:

PDEtools[declare](y(x), f(x), prime=x);  # turn ON the enhanced DE display

$y\left(x\right)\mathrm{will\; now\; be\; displayed\; as}y$

$f\left(x\right)\mathrm{will\; now\; be\; displayed\; as}f$

$\mathrm{derivatives\; with\; respect\; to}x\mathrm{of\; functions\; of\; one\; variable\; will\; now\; be\; displayed\; with\; \text{'}}$

diff(y(x),x) = f[3](x)*y(x)^3 + f[2](x)*y(x)^2 + f[1](x)*y(x) + f[0](x);

$\mathrm{y\text{'}}=f_3{y}^3+f_2{y}^2+f_{1}y+f_0$

where { f[1], f[2] f[3] } are arbitrary functions of x. Abel equations appear in the reduction of order process related to finding exact solutions to many second and higher order ODE families, and hence are frequently found in the modelling of real problems in varied areas. That has for a long time motivated people to study their integrable cases and related solving methods.

The implementation of methods for Abel ODEs dsolve follows the presentation done in Abel Equations: Equivalence and New Integrable Classes, by E.S. Cheb-Terrab and A.D. Roche, Computer Physics Communications 130 (2000).

A general "exact integration" strategy for these ODEs was first formulated by Liouville, in then 19th century, and is based on the concepts of classes, invariants and the solving of the equivalence problem. Generally speaking, two Abel ODEs belong to the same equivalence class if and only if one can be obtained from the other by means of a transformation of the form

TR := {x = F(t), y(x)= P(t)*u(t) + Q(t)};

$\mathrm{TR}≔\left\{x=F\left(t\right)\,y=P\left(t\right)u\left(t\right)+Q\left(t\right)\right\}$

where t and u(t) are respectively the new independent and dependent variables, and F, P, Q are arbitrary functions of t satisfying

diff(F(x),x)<>0, P <> 0;

$\mathrm{F\text{'}}\ne 0,P\ne 0$

Each class has infinitely many member and there are infinitely many classes. To each class there corresponds a different set of values of the invariants (built with the coefficients {f[1], f[2], f[3] } and their derivatives), and actually any one of them (we shall pick one and call it the invariant) is enough to characterize a class.

A simple integrable case happens when the invariant is constant; the solution to the ODE then follows straightforwardly in terms of quadratures, as explained in textbooks and in odeadvisor, Abel (the method for constant invariant Abel ODEs was implemented in Maple Release 5).

On the contrary, when the invariant is not constant, just a few integrable cases are known and the formulation of solving strategies is entirely based on the equivalence between two such Abel ODEs (one of which is integrable) under the transformation TR.

When a given Abel ODE belongs to one of a set of the solvable classes collected in the aforementioned CPC paper, dsolve's routines first determine this fact, without solving any differential equations, and use it to return a closed form solution without requiring further participation from the user. Apart from new solvable Abel classes, the ODE families that are covered include, as particular cases, all the Abel solvable cases presented in Kamke's and Murphy's books.

infolevel[dsolve] := 4;            # turn ON the display of userinfos

${\mathrm{infolevel}}_{\mathrm{dsolve}}≔4$

ode[185] := x^7*diff(y(x),x)+2*(x^2+1)*y(x)^3+5*x^3*y(x)^2 = 0;

${\mathrm{ode}}_{185}≔x^7\mathrm{y\text{'}}+2\left(x^2+1\right)y^3+5x^3{y}^2=0$

sol[185] := dsolve((6), implicit, useInt);

Methods for first order ODEs:
--- Trying classification methods ---
trying a quadrature
trying 1st order linear
trying Bernoulli
trying separable
trying inverse linear
trying homogeneous types:
trying Chini
differential order: 1; looking for linear symmetries
trying exact
trying Abel
The relative invariant s3 is: 10/27/x^12*(9*x^2+2)
The first absolute invariant s5^3/s3^5 is: -729/100*(90*x^4-15*x^2-14)^3/(9*x^2+2)^5
The second absolute invariant s3*s7/s5^2 is: 5/3*(9*x^2+2)*(972*x^6-324*x^4-15*x^2+98)/(90*x^4-15*x^2-14)^2
...checking Abel class AIL (45)
...checking Abel class AIL (310)
...checking Abel class AIR (36)
...checking Abel class AIL (301)
...checking Abel class AIL (1000)
...checking Abel class AIL (42)
...checking Abel class AIL (185)
inverse of the transformation solving the problem is: {t = x, u(t) = y(x)}
<- Abel successful

${\mathrm{sol}}_{185}≔\mathrm{c\_\_1}+\frac{x}{{\left({\left(\frac{1}{x}+\frac{x^2}{y}\right)}^2+1\right)}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}}+\frac{\left({\int }_{\phantom{\mathrm{`\; `}}}^{\frac{1}{x}+\frac{x^2}{y}}\frac{1}{{\left({\mathrm{\_a}}^2+1\right)}^{\raisebox{1ex}{$5$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}}\&DifferentialD;\mathrm{\_a}\right)}{2}=0$

odetest( sol[185], ode[185] );

$0$

tr := {x=F(t),y(x)=P(t)*u(t)+Q(t)};

$\mathrm{tr}≔\left\{x=F\left(t\right)\&comma;y=P\left(t\right)u\left(t\right)+Q\left(t\right)\right\}$

subs([t=x, u=y], PDEtools[dchange](tr,ode[185],[u(t),t]) );

$\frac{{F\left(x\right)}^7\left(\mathrm{P\text{'}}y+P\left(x\right)\mathrm{y\text{'}}+\mathrm{Q\text{'}}\right)}{\mathrm{F\text{'}}}+2\left({F\left(x\right)}^2+1\right){\left(P\left(x\right)y+Q\left(x\right)\right)}^3+5{F\left(x\right)}^3{\left(P\left(x\right)y+Q\left(x\right)\right)}^2=0$

collect(  isolate((10),diff(y(x),x)),  y, normal);

$\mathrm{y\text{'}}=-\frac{2\left({F\left(x\right)}^2+1\right){P\left(x\right)}^2\mathrm{F\text{'}}{y}^3}{{F\left(x\right)}^7}-\frac{P\left(x\right)\left(5{F\left(x\right)}^3+6{F\left(x\right)}^{2}Q\left(x\right)+6Q\left(x\right)\right)\mathrm{F\text{'}}{y}^2}{{F\left(x\right)}^7}-\frac{\left(\mathrm{P\text{'}}{F\left(x\right)}^7+10{F\left(x\right)}^{3}P\left(x\right)\mathrm{F\text{'}}Q\left(x\right)+6{F\left(x\right)}^{2}P\left(x\right)\mathrm{F\text{'}}{Q\left(x\right)}^2+6P\left(x\right)\mathrm{F\text{'}}{Q\left(x\right)}^2\right)y}{{F\left(x\right)}^{7}P\left(x\right)}-\frac{\mathrm{Q\text{'}}{F\left(x\right)}^7+5{F\left(x\right)}^3\mathrm{F\text{'}}{Q\left(x\right)}^2+2{F\left(x\right)}^2\mathrm{F\text{'}}{Q\left(x\right)}^3+2\mathrm{F\text{'}}{Q\left(x\right)}^3}{{F\left(x\right)}^{7}P\left(x\right)}$

dsolve((11), y(x), implicit, useInt);

Methods for first order ODEs:
--- Trying classification methods ---
trying a quadrature
trying 1st order linear
trying Bernoulli
trying separable
trying inverse linear
trying homogeneous types:
trying Chini
differential order: 1; looking for linear symmetries
trying exact
trying Abel
The first invariant is non-rational: -729/100*(90*F(x)^4-15*F(x)^2-14)^3/(9*F(x)^2+2)^5
-> Searching for a convenient change of variables...
<- Unable to rationalize the invariant
The relative invariant s3 is: 10/27*(9*F(x)^2+2)/F(x)^12*diff(F(x),x)^3*P(x)^3
The first absolute invariant s5^3/s3^5 is: -729/100*(90*F(x)^4-15*F(x)^2-14)^3/(9*F(x)^2+2)^5
The second absolute invariant s3*s7/s5^2 is: 5/3*(9*F(x)^2+2)*(972*F(x)^6-324*F(x)^4-15*F(x)^2+98)/(90*F(x)^4-15*F(x)^2-14)^2
...checking Abel class AIL (45)
...checking Abel class AIL (310)
...checking Abel class AIR (36)
...checking Abel class AIL (301)
...checking Abel class AIL (1000)
...checking Abel class AIL (42)
...checking Abel class AIL (185)
inverse of the transformation solving the problem is: {t = F(x), u(t) = P(x)*y(x)+Q(x)}
<- Abel successful

$\mathrm{c\_\_1}+\frac{F\left(x\right)}{{\left({\left(\frac{1}{F\left(x\right)}+\frac{{F\left(x\right)}^2}{P\left(x\right)y+Q\left(x\right)}\right)}^2+1\right)}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}}+\frac{\left({\int }_{\phantom{\mathrm{`\; `}}}^{\frac{1}{F\left(x\right)}+\frac{{F\left(x\right)}^2}{P\left(x\right)y+Q\left(x\right)}}\frac{1}{{\left({\mathrm{\_a}}^2+1\right)}^{\raisebox{1ex}{$5$}\!\left/ \!\raisebox{-1ex}{$4$}\right.}}\&DifferentialD;\mathrm{\_a}\right)}{2}=0$

ode[257] := x*(x*y(x)+x^4-1)*diff(y(x),x)-y(x)*(x*y(x)-x^4-1) = 0;

${\mathrm{ode}}_{257}≔x\left(xy+x^4-1\right)\mathrm{y\text{'}}-y\left(xy-x^4-1\right)=0$

dsolve(ode[257], implicit);

Methods for first order ODEs:
--- Trying classification methods ---
trying a quadrature
trying 1st order linear
trying Bernoulli
trying separable
trying inverse linear
trying homogeneous types:
trying Chini
differential order: 1; looking for linear symmetries
trying exact
trying Abel
The equivalent Abel ODE of 1st kind is: diff(u(x),x) = -2*x*(x^2+1)*(x^2-1)*u(x)^3-2/x^2*u(x)^2-1/x*u(x)
The relative invariant s3 is: -8/27*(27*x^8-9*x^4+2)/x^6
The first absolute invariant s5^3/s3^5 is: 729*(135*x^16-36*x^12+54*x^8-15*x^4+2)^3/(27*x^8-9*x^4+2)^5
The second absolute invariant s3*s7/s5^2 is: 1/3*(27*x^8-9*x^4+2)*(2835*x^24+243*x^20+2349*x^16-927*x^12+495*x^8-105*x^4+10)/(135*x^16-36*x^12+54*x^8-15*x^4+2)^2
...checking Abel class AIL (45)
...checking Abel class AIL (310)
...checking Abel class AIR (36)
...checking Abel class AIL (301)
...checking Abel class AIL (1000)
...checking Abel class AIL (42)
...checking Abel class AIL (185)
...checking Abel class AIA (by Halphen)
...checking Abel class AIL (205)
...checking Abel class AIA (147)
...checking Abel class AIL (581)
...checking Abel class AIL (200)
...checking Abel class AIL (257)
inverse of the transformation solving the problem is: {t = 1/x^2, u(t) = x*y(x)}
<- Abel successful

$\mathrm{c\_\_1}+\left(-xy+1\right){\&ExponentialE;}^{\frac{y\left(2x^3+y\right)}{2x^2}}=0$

Class[257] := diff(y(x),x) = -(x-1)*(x+1)*y(x)^3/(x^3)+y(x)^2;

${\mathrm{Class}}_{257}≔\mathrm{y\text{'}}=-\frac{\left(x-1\right)\left(x+1\right)y^3}{x^3}+y^2$

Abel_FirstKind_257 := convert(ode[257], FirstKind, 'keep'); # 'keep' the names (x, y(x))

$\mathrm{Abel\_FirstKind\_257}≔\mathrm{y\text{'}}=\left(-2x^5+2x\right)y^3-\frac{2y^2}{x^2}-\frac{y}{x},\left\{\mathrm{\_a}\left(x\right)=\frac{x^2+\left(-x^5+x\right)y}{y{x}^2}\right\}$

tr := {y(x) = u(t)/(sqrt(1/t)), x = sqrt(1/t)};

$\mathrm{tr}≔\left\{x=\sqrt{\frac{1}{t}}\&comma;y=\frac{u\left(t\right)}{\sqrt{\frac{1}{t}}}\right\}$

subs( [t=x, u=y], PDEtools[dchange](tr, Abel_FirstKind_257[1], normal));

$-\left(2x\mathrm{y\text{'}}+y\right)x=\frac{y\left(2x^2{y}^2-2y{x}^3-2y^2-x^2\right)}{x}$

collect( isolate( (18), diff(y(x),x) ), y(x), factor );

$\mathrm{y\text{'}}=-\frac{\left(x-1\right)\left(x+1\right)y^3}{x^3}+y^2$

ode[43] := diff(y(x),x)+(3*x^2*a+4*x*a^2+b)*y(x)^3+3*x*y(x)^2 = 0;

${\mathrm{ode}}_{43}≔\mathrm{y\text{'}}+\left(4x{a}^2+3x^{2}a+b\right)y^3+3y^{2}x=0$

Class[B] := diff(y(x),x) = 2*(x^2-C)*y(x)^3+2*(x+1)*y(x)^2;

${\mathrm{Class}}_B≔\mathrm{y\text{'}}=2\left(x^2-C\right)y^3+2\left(x+1\right)y^2$

tr := {x = -1/2*(2*a+3*t)/a, y(x) = -2/3*a^2*u(t)};

$\mathrm{tr}≔\left\{x=-\frac{2a+3t}{2a}\&comma;y=-\frac{2a^{2}u\left(t\right)}{3}\right\}$

eq_C := C = 1/4*(-3*b+4*a^3)/a^3;

$\mathrm{eq\_C}≔C=\frac{4a^3-3b}{4a^3}$

subs( [t=x, u=y, eq_C], PDEtools[dchange](tr, Class[B], [t,u(t)], normal));

isolate((24),diff(y(x),x));    # isolate y'

collect((25),y(x),factor);

infolevel[dsolve] := 0: # turn OFF userinfos

dsolve( ode[43] );

`dsolve/Abel_seed`(numbers);

$\left[1.1\&comma;1.2\&comma;1.3\&comma;1.4\&comma;1.5\&comma;1.51\&comma;1.52\&comma;1.53\&comma;1.6\&comma;1.8\&comma;1.9\&comma;33\&comma;36\&comma;42\&comma;45\&comma;147\&comma;185\&comma;200\&comma;201\&comma;205\&comma;257\&comma;301\&comma;310\&comma;400\&comma;515\&comma;581\&comma;815\&comma;1000\&comma;1001\right]$

nops((28));

$29$

`dsolve/Abel_seed`(33, y(x));  # instead of x and y may be any other letters

$\mathrm{y\text{'}}=\frac{\left(x-3\right)y^3}{8x}-\frac{5y^2}{4x}-\frac{3y}{8x}$