4.6.16 \(x^2 y'(x)=a+b x^n+x^2 y(x)^2\)

ODE
\[ x^2 y'(x)=a+b x^n+x^2 y(x)^2 \] ODE Classification

[_rational, _Riccati]

Book solution method
Riccati ODE, Generalized ODE

Mathematica
cpu = 0.135075 (sec), leaf count = 1260

\[\left \{\left \{y(x)\to \frac {-n^{\frac {2 \sqrt {(1-4 a) n^2}}{n^2}+1} \left (x^n\right )^{\frac {i \sqrt {4 a-1}}{n}+1} J_{\frac {\sqrt {(1-4 a) n^2}}{n^2}-1}\left (\frac {2 \sqrt {b} \sqrt {x^n}}{n}\right ) \Gamma \left (\frac {n+\sqrt {1-4 a}}{n}\right ) b^{\frac {i \sqrt {4 a-1}}{n}+\frac {1}{2}}+n^{\frac {2 \sqrt {(1-4 a) n^2}}{n^2}+1} \left (x^n\right )^{\frac {i \sqrt {4 a-1}}{n}+1} J_{\frac {\sqrt {(1-4 a) n^2}}{n^2}+1}\left (\frac {2 \sqrt {b} \sqrt {x^n}}{n}\right ) \Gamma \left (\frac {n+\sqrt {1-4 a}}{n}\right ) b^{\frac {i \sqrt {4 a-1}}{n}+\frac {1}{2}}-i \sqrt {4 a-1} n^{\frac {2 \sqrt {(1-4 a) n^2}}{n^2}+1} \left (x^n\right )^{\frac {i \sqrt {4 a-1}}{n}+\frac {1}{2}} J_{\frac {\sqrt {(1-4 a) n^2}}{n^2}}\left (\frac {2 \sqrt {b} \sqrt {x^n}}{n}\right ) \Gamma \left (\frac {n+\sqrt {1-4 a}}{n}\right ) b^{\frac {i \sqrt {4 a-1}}{n}}-n^{\frac {2 \sqrt {(1-4 a) n^2}}{n^2}+1} \left (x^n\right )^{\frac {i \sqrt {4 a-1}}{n}+\frac {1}{2}} J_{\frac {\sqrt {(1-4 a) n^2}}{n^2}}\left (\frac {2 \sqrt {b} \sqrt {x^n}}{n}\right ) \Gamma \left (\frac {n+\sqrt {1-4 a}}{n}\right ) b^{\frac {i \sqrt {4 a-1}}{n}}+n^{\frac {2 \sqrt {(1-4 a) n^2}}{n^2}} \sqrt {(1-4 a) n^2} \left (x^n\right )^{\frac {i \sqrt {4 a-1}}{n}+\frac {1}{2}} J_{\frac {\sqrt {(1-4 a) n^2}}{n^2}}\left (\frac {2 \sqrt {b} \sqrt {x^n}}{n}\right ) \Gamma \left (\frac {n+\sqrt {1-4 a}}{n}\right ) b^{\frac {i \sqrt {4 a-1}}{n}}-n^{\frac {2 i \sqrt {4 a-1}}{n}} \left (-i \sqrt {4 a-1} n+n+\sqrt {(1-4 a) n^2}\right ) \left (x^n\right )^{\frac {\sqrt {(1-4 a) n^2}}{n^2}+\frac {1}{2}} J_{-\frac {\sqrt {(1-4 a) n^2}}{n^2}}\left (\frac {2 \sqrt {b} \sqrt {x^n}}{n}\right ) c_1 \Gamma \left (1-\frac {\sqrt {1-4 a}}{n}\right ) b^{\frac {\sqrt {(1-4 a) n^2}}{n^2}}-n^{\frac {2 i \sqrt {4 a-1}}{n}+1} \left (x^n\right )^{\frac {\sqrt {(1-4 a) n^2}}{n^2}+1} J_{-\frac {\sqrt {(1-4 a) n^2}}{n^2}-1}\left (\frac {2 \sqrt {b} \sqrt {x^n}}{n}\right ) c_1 \Gamma \left (1-\frac {\sqrt {1-4 a}}{n}\right ) b^{\frac {\sqrt {(1-4 a) n^2}}{n^2}+\frac {1}{2}}+n^{\frac {2 i \sqrt {4 a-1}}{n}+1} \left (x^n\right )^{\frac {\sqrt {(1-4 a) n^2}}{n^2}+1} J_{1-\frac {\sqrt {(1-4 a) n^2}}{n^2}}\left (\frac {2 \sqrt {b} \sqrt {x^n}}{n}\right ) c_1 \Gamma \left (1-\frac {\sqrt {1-4 a}}{n}\right ) b^{\frac {\sqrt {(1-4 a) n^2}}{n^2}+\frac {1}{2}}}{2 n x \sqrt {x^n} \left (b^{\frac {i \sqrt {4 a-1}}{n}} n^{\frac {2 \sqrt {(1-4 a) n^2}}{n^2}} J_{\frac {\sqrt {(1-4 a) n^2}}{n^2}}\left (\frac {2 \sqrt {b} \sqrt {x^n}}{n}\right ) \Gamma \left (\frac {n+\sqrt {1-4 a}}{n}\right ) \left (x^n\right )^{\frac {i \sqrt {4 a-1}}{n}}+b^{\frac {\sqrt {(1-4 a) n^2}}{n^2}} n^{\frac {2 i \sqrt {4 a-1}}{n}} J_{-\frac {\sqrt {(1-4 a) n^2}}{n^2}}\left (\frac {2 \sqrt {b} \sqrt {x^n}}{n}\right ) c_1 \Gamma \left (1-\frac {\sqrt {1-4 a}}{n}\right ) \left (x^n\right )^{\frac {\sqrt {(1-4 a) n^2}}{n^2}}\right )}\right \}\right \}\]

Maple
cpu = 0.099 (sec), leaf count = 196

\[ \left \{ y \left ( x \right ) ={\frac {1}{2\,x} \left ( 2\,\sqrt {b} \left ( {{\sl Y}_{{\frac {\sqrt {1-4\,a}+n}{n}}}\left (2\,{\frac {\sqrt {b}{x}^{n/2}}{n}}\right )}{\it \_C1}+{{\sl J}_{{\frac {\sqrt {1-4\,a}+n}{n}}}\left (2\,{\frac {\sqrt {b}{x}^{n/2}}{n}}\right )} \right ) {x}^{n/2}- \left ( \sqrt {1-4\,a}+1 \right ) \left ( {{\sl Y}_{{\frac {1}{n}\sqrt {1-4\,a}}}\left (2\,{\frac {\sqrt {b}{x}^{n/2}}{n}}\right )}{\it \_C1}+{{\sl J}_{{\frac {1}{n}\sqrt {1-4\,a}}}\left (2\,{\frac {\sqrt {b}{x}^{n/2}}{n}}\right )} \right ) \right ) \left ( {{\sl Y}_{{\frac {1}{n}\sqrt {1-4\,a}}}\left (2\,{\frac {\sqrt {b}{x}^{n/2}}{n}}\right )}{\it \_C1}+{{\sl J}_{{\frac {1}{n}\sqrt {1-4\,a}}}\left (2\,{\frac {\sqrt {b}{x}^{n/2}}{n}}\right )} \right ) ^{-1}} \right \} \] Mathematica raw input

DSolve[x^2*y'[x] == a + b*x^n + x^2*y[x]^2,y[x],x]

Mathematica raw output

{{y[x] -> (-(b^(Sqrt[(1 - 4*a)*n^2]/n^2)*n^(((2*I)*Sqrt[-1 + 4*a])/n)*(n - I*Sqr
t[-1 + 4*a]*n + Sqrt[(1 - 4*a)*n^2])*(x^n)^(1/2 + Sqrt[(1 - 4*a)*n^2]/n^2)*Besse
lJ[-(Sqrt[(1 - 4*a)*n^2]/n^2), (2*Sqrt[b]*Sqrt[x^n])/n]*C[1]*Gamma[1 - Sqrt[1 - 
4*a]/n]) - b^(1/2 + Sqrt[(1 - 4*a)*n^2]/n^2)*n^(1 + ((2*I)*Sqrt[-1 + 4*a])/n)*(x
^n)^(1 + Sqrt[(1 - 4*a)*n^2]/n^2)*BesselJ[-1 - Sqrt[(1 - 4*a)*n^2]/n^2, (2*Sqrt[
b]*Sqrt[x^n])/n]*C[1]*Gamma[1 - Sqrt[1 - 4*a]/n] + b^(1/2 + Sqrt[(1 - 4*a)*n^2]/
n^2)*n^(1 + ((2*I)*Sqrt[-1 + 4*a])/n)*(x^n)^(1 + Sqrt[(1 - 4*a)*n^2]/n^2)*Bessel
J[1 - Sqrt[(1 - 4*a)*n^2]/n^2, (2*Sqrt[b]*Sqrt[x^n])/n]*C[1]*Gamma[1 - Sqrt[1 - 
4*a]/n] - b^((I*Sqrt[-1 + 4*a])/n)*n^(1 + (2*Sqrt[(1 - 4*a)*n^2])/n^2)*(x^n)^(1/
2 + (I*Sqrt[-1 + 4*a])/n)*BesselJ[Sqrt[(1 - 4*a)*n^2]/n^2, (2*Sqrt[b]*Sqrt[x^n])
/n]*Gamma[(Sqrt[1 - 4*a] + n)/n] - I*Sqrt[-1 + 4*a]*b^((I*Sqrt[-1 + 4*a])/n)*n^(
1 + (2*Sqrt[(1 - 4*a)*n^2])/n^2)*(x^n)^(1/2 + (I*Sqrt[-1 + 4*a])/n)*BesselJ[Sqrt
[(1 - 4*a)*n^2]/n^2, (2*Sqrt[b]*Sqrt[x^n])/n]*Gamma[(Sqrt[1 - 4*a] + n)/n] + b^(
(I*Sqrt[-1 + 4*a])/n)*n^((2*Sqrt[(1 - 4*a)*n^2])/n^2)*Sqrt[(1 - 4*a)*n^2]*(x^n)^
(1/2 + (I*Sqrt[-1 + 4*a])/n)*BesselJ[Sqrt[(1 - 4*a)*n^2]/n^2, (2*Sqrt[b]*Sqrt[x^
n])/n]*Gamma[(Sqrt[1 - 4*a] + n)/n] - b^(1/2 + (I*Sqrt[-1 + 4*a])/n)*n^(1 + (2*S
qrt[(1 - 4*a)*n^2])/n^2)*(x^n)^(1 + (I*Sqrt[-1 + 4*a])/n)*BesselJ[-1 + Sqrt[(1 -
 4*a)*n^2]/n^2, (2*Sqrt[b]*Sqrt[x^n])/n]*Gamma[(Sqrt[1 - 4*a] + n)/n] + b^(1/2 +
 (I*Sqrt[-1 + 4*a])/n)*n^(1 + (2*Sqrt[(1 - 4*a)*n^2])/n^2)*(x^n)^(1 + (I*Sqrt[-1
 + 4*a])/n)*BesselJ[1 + Sqrt[(1 - 4*a)*n^2]/n^2, (2*Sqrt[b]*Sqrt[x^n])/n]*Gamma[
(Sqrt[1 - 4*a] + n)/n])/(2*n*x*Sqrt[x^n]*(b^(Sqrt[(1 - 4*a)*n^2]/n^2)*n^(((2*I)*
Sqrt[-1 + 4*a])/n)*(x^n)^(Sqrt[(1 - 4*a)*n^2]/n^2)*BesselJ[-(Sqrt[(1 - 4*a)*n^2]
/n^2), (2*Sqrt[b]*Sqrt[x^n])/n]*C[1]*Gamma[1 - Sqrt[1 - 4*a]/n] + b^((I*Sqrt[-1 
+ 4*a])/n)*n^((2*Sqrt[(1 - 4*a)*n^2])/n^2)*(x^n)^((I*Sqrt[-1 + 4*a])/n)*BesselJ[
Sqrt[(1 - 4*a)*n^2]/n^2, (2*Sqrt[b]*Sqrt[x^n])/n]*Gamma[(Sqrt[1 - 4*a] + n)/n]))
}}

Maple raw input

dsolve(x^2*diff(y(x),x) = a+b*x^n+x^2*y(x)^2, y(x),'implicit')

Maple raw output

y(x) = 1/2*(2*b^(1/2)*(BesselY(((1-4*a)^(1/2)+n)/n,2*b^(1/2)/n*x^(1/2*n))*_C1+Be
sselJ(((1-4*a)^(1/2)+n)/n,2*b^(1/2)/n*x^(1/2*n)))*x^(1/2*n)-((1-4*a)^(1/2)+1)*(B
esselY((1-4*a)^(1/2)/n,2*b^(1/2)/n*x^(1/2*n))*_C1+BesselJ((1-4*a)^(1/2)/n,2*b^(1
/2)/n*x^(1/2*n))))/x/(BesselY((1-4*a)^(1/2)/n,2*b^(1/2)/n*x^(1/2*n))*_C1+BesselJ
((1-4*a)^(1/2)/n,2*b^(1/2)/n*x^(1/2*n)))