4.36.50 y(x)=a+4b2y(x)+3by(x)2+3y(x)y(x)

ODE
y(x)=a+4b2y(x)+3by(x)2+3y(x)y(x) ODE Classification

[[_2nd_order, _missing_x], [_2nd_order, _reducible, _mu_x_y1], [_2nd_order, _reducible, _mu_xy]]

Book solution method
TO DO

Mathematica
cpu = 10.979 (sec), leaf count = 2742

{{y(x)i4b33ab3/223(b3+4b63ab3)2b334b33a2b3/2a4b33a2b3/2b4b63ab3b3+3(e2bxa)4b33a2b3/2I4b63ab32b3(32ae2bxac1b)Γ(4b63ab32b3+1)c14b33a2b3/2+i4b33ab3/2234b63ab32b3+1234b33a2b3/2a4b33a2b3/2b4b63ab3b3+324b33a(e2bxa)4b33a2b3/2I4b63ab32b3(32ae2bxac1b)Γ(4b63ab32b3+1)c14b33a2b3/2i4b33ab3/2234b63ab32b3+1234b33a2b3/2a4b33a2b3/2b4b63ab3b34b63ab3(e2bxa)4b33a2b3/2I4b63ab32b3(32ae2bxac1b)Γ(4b63ab32b3+1)c14b33a2b3/2+234b33a2b3/2+1234b63ab32b3a4b63ab32b3b4b33ab3/2(2b34b33ab3/2+4b63ab3)(e2bxa)4b63ab32b3I4b63ab32b3(32ae2bxac1b)c2Γ(14b63ab32b3)c14b63ab32b3+i4b33ab3/2234b63ab32b3312(4b33ab3/2+1)a12(4b33ab3/2+1)b4b63ab3b3+2(e2bxa)12(4b33ab3/2+1)I4b63ab32b31(32ae2bxac1b)Γ(4b63ab32b3+1)c112(4b33ab3/2+1)+i4b33ab3/2234b63ab32b3312(4b33ab3/2+1)a12(4b33ab3/2+1)b4b63ab3b3+2(e2bxa)12(4b33ab3/2+1)I4b63ab32b3+1(32ae2bxac1b)Γ(4b63ab32b3+1)c112(4b33ab3/2+1)+234b33a2b3/23b3+4b63ab32b3ab3+4b63ab32b3b4b33ab3/2+2(e2bxa)b3+4b63ab32b3I4b63ab32b31(32ae2bxac1b)c2Γ(14b63ab32b3)c1b3+4b63ab32b3+234b33a2b3/23b3+4b63ab32b3ab3+4b63ab32b3b4b33ab3/2+2(e2bxa)b3+4b63ab32b3I14b63ab32b3(32ae2bxac1b)c2Γ(14b63ab32b3)c1b3+4b63ab32b332b2(i4b33ab3/2234b63ab32b334b33a2b3/2a4b33a2b3/2b4b63ab3b3(e2bxa)4b33a2b3/2I4b63ab32b3(32ae2bxac1b)Γ(4b63ab32b3+1)c14b33a2b3/2+234b33a2b3/234b63ab32b3a4b63ab32b3b4b33ab3/2(e2bxa)4b63ab32b3I4b63ab32b3(32ae2bxac1b)c2Γ(14b63ab32b3)c14b63ab32b3)}}

Maple
cpu = 0.646 (sec), leaf count = 849

{y(x)2b2(4_ab33b2_a2+(RootOf(6K1/24b3+3a4b43abb(1/23_Zb2)_C1b2+9K1/24b3+3a4b43abb(1/23_Zb2)_C1b_a34b43abK1/24b3+3a4b43abb(1/23_Zb2)_C1+3K1/24b3+24b43abb3a4b43abb(1/23_Zb2)3_C1_Z+6I1/24b3+3a4b43abb(1/23_Zb2)b2+9I1/24b3+3a4b43abb(1/23_Zb2)b_a3I1/24b3+24b43abb3a4b43abb(1/23_Zb2)3_Z34b43abI1/24b3+3a4b43abb(1/23_Zb2)))2ab)1d_ax_C2=0,y(x)2b2(4_ab33b2_a2+(RootOf(6K1/24b3+3a4b43abb(1/23_Zb2)_C1b2+9K1/24b3+3a4b43abb(1/23_Zb2)_C1b_a34b43abK1/24b3+3a4b43abb(1/23_Zb2)_C1+3K1/24b3+24b43abb3a4b43abb(1/23_Zb2)3_C1_Z+6I1/24b3+3a4b43abb(1/23_Zb2)b2+9I1/24b3+3a4b43abb(1/23_Zb2)b_a3I1/24b3+24b43abb3a4b43abb(1/23_Zb2)3_Z34b43abI1/24b3+3a4b43abb(1/23_Zb2)))2ab)1d_ax_C2=0} Mathematica raw input

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

Mathematica raw output

{{y[x] -> (2^(1/2 + (3*Sqrt[-3*a + 4*b^3])/(2*b^(3/2)))*3^(Sqrt[-3*a*b^3 + 4*b^6
]/(2*b^3))*a^(Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3))*b^(Sqrt[-3*a + 4*b^3]/b^(3/2))*(-2
*b^3 - b^(3/2)*Sqrt[-3*a + 4*b^3] + Sqrt[-3*a*b^3 + 4*b^6])*(1/(a*E^(2*b*x)))^(S
qrt[-3*a*b^3 + 4*b^6]/(2*b^3))*BesselI[-Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3), (Sqrt[3/
2]*Sqrt[a]*Sqrt[1/(a*E^(2*b*x))]*Sqrt[C[1]])/b]*C[1]^(Sqrt[-3*a*b^3 + 4*b^6]/(2*
b^3))*C[2]*Gamma[1 - Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3)] + 2^((3*Sqrt[-3*a + 4*b^3])
/(2*b^(3/2)))*3^((b^3 + Sqrt[-3*a*b^3 + 4*b^6])/(2*b^3))*a^((b^3 + Sqrt[-3*a*b^3
 + 4*b^6])/(2*b^3))*b^(2 + Sqrt[-3*a + 4*b^3]/b^(3/2))*(1/(a*E^(2*b*x)))^((b^3 +
 Sqrt[-3*a*b^3 + 4*b^6])/(2*b^3))*BesselI[-1 - Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3), (
Sqrt[3/2]*Sqrt[a]*Sqrt[1/(a*E^(2*b*x))]*Sqrt[C[1]])/b]*C[1]^((b^3 + Sqrt[-3*a*b^
3 + 4*b^6])/(2*b^3))*C[2]*Gamma[1 - Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3)] + 2^((3*Sqrt
[-3*a + 4*b^3])/(2*b^(3/2)))*3^((b^3 + Sqrt[-3*a*b^3 + 4*b^6])/(2*b^3))*a^((b^3 
+ Sqrt[-3*a*b^3 + 4*b^6])/(2*b^3))*b^(2 + Sqrt[-3*a + 4*b^3]/b^(3/2))*(1/(a*E^(2
*b*x)))^((b^3 + Sqrt[-3*a*b^3 + 4*b^6])/(2*b^3))*BesselI[1 - Sqrt[-3*a*b^3 + 4*b
^6]/(2*b^3), (Sqrt[3/2]*Sqrt[a]*Sqrt[1/(a*E^(2*b*x))]*Sqrt[C[1]])/b]*C[1]^((b^3 
+ Sqrt[-3*a*b^3 + 4*b^6])/(2*b^3))*C[2]*Gamma[1 - Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3)
] - I^(Sqrt[-3*a + 4*b^3]/b^(3/2))*2^((3*(b^3 + Sqrt[-3*a*b^3 + 4*b^6]))/(2*b^3)
)*3^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*a^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*b^(3 + S
qrt[-3*a*b^3 + 4*b^6]/b^3)*(1/(a*E^(2*b*x)))^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*Be
sselI[Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3), (Sqrt[3/2]*Sqrt[a]*Sqrt[1/(a*E^(2*b*x))]*S
qrt[C[1]])/b]*C[1]^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*Gamma[1 + Sqrt[-3*a*b^3 + 4*
b^6]/(2*b^3)] + I^(Sqrt[-3*a + 4*b^3]/b^(3/2))*2^(1/2 + (3*Sqrt[-3*a*b^3 + 4*b^6
])/(2*b^3))*3^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*a^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2))
)*b^(3/2 + Sqrt[-3*a*b^3 + 4*b^6]/b^3)*Sqrt[-3*a + 4*b^3]*(1/(a*E^(2*b*x)))^(Sqr
t[-3*a + 4*b^3]/(2*b^(3/2)))*BesselI[Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3), (Sqrt[3/2]*
Sqrt[a]*Sqrt[1/(a*E^(2*b*x))]*Sqrt[C[1]])/b]*C[1]^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)
))*Gamma[1 + Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3)] - I^(Sqrt[-3*a + 4*b^3]/b^(3/2))*2^
(1/2 + (3*Sqrt[-3*a*b^3 + 4*b^6])/(2*b^3))*3^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*a^
(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*b^(Sqrt[-3*a*b^3 + 4*b^6]/b^3)*Sqrt[-3*a*b^3 + 
4*b^6]*(1/(a*E^(2*b*x)))^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*BesselI[Sqrt[-3*a*b^3 
+ 4*b^6]/(2*b^3), (Sqrt[3/2]*Sqrt[a]*Sqrt[1/(a*E^(2*b*x))]*Sqrt[C[1]])/b]*C[1]^(
Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*Gamma[1 + Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3)] + I^(S
qrt[-3*a + 4*b^3]/b^(3/2))*2^((3*Sqrt[-3*a*b^3 + 4*b^6])/(2*b^3))*3^((1 + Sqrt[-
3*a + 4*b^3]/b^(3/2))/2)*a^((1 + Sqrt[-3*a + 4*b^3]/b^(3/2))/2)*b^(2 + Sqrt[-3*a
*b^3 + 4*b^6]/b^3)*(1/(a*E^(2*b*x)))^((1 + Sqrt[-3*a + 4*b^3]/b^(3/2))/2)*Bessel
I[-1 + Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3), (Sqrt[3/2]*Sqrt[a]*Sqrt[1/(a*E^(2*b*x))]*
Sqrt[C[1]])/b]*C[1]^((1 + Sqrt[-3*a + 4*b^3]/b^(3/2))/2)*Gamma[1 + Sqrt[-3*a*b^3
 + 4*b^6]/(2*b^3)] + I^(Sqrt[-3*a + 4*b^3]/b^(3/2))*2^((3*Sqrt[-3*a*b^3 + 4*b^6]
)/(2*b^3))*3^((1 + Sqrt[-3*a + 4*b^3]/b^(3/2))/2)*a^((1 + Sqrt[-3*a + 4*b^3]/b^(
3/2))/2)*b^(2 + Sqrt[-3*a*b^3 + 4*b^6]/b^3)*(1/(a*E^(2*b*x)))^((1 + Sqrt[-3*a + 
4*b^3]/b^(3/2))/2)*BesselI[1 + Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3), (Sqrt[3/2]*Sqrt[a
]*Sqrt[1/(a*E^(2*b*x))]*Sqrt[C[1]])/b]*C[1]^((1 + Sqrt[-3*a + 4*b^3]/b^(3/2))/2)
*Gamma[1 + Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3)])/(3*Sqrt[2]*b^2*(2^((3*Sqrt[-3*a + 4*
b^3])/(2*b^(3/2)))*3^(Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3))*a^(Sqrt[-3*a*b^3 + 4*b^6]/
(2*b^3))*b^(Sqrt[-3*a + 4*b^3]/b^(3/2))*(1/(a*E^(2*b*x)))^(Sqrt[-3*a*b^3 + 4*b^6
]/(2*b^3))*BesselI[-Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3), (Sqrt[3/2]*Sqrt[a]*Sqrt[1/(a
*E^(2*b*x))]*Sqrt[C[1]])/b]*C[1]^(Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3))*C[2]*Gamma[1 -
 Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3)] + I^(Sqrt[-3*a + 4*b^3]/b^(3/2))*2^((3*Sqrt[-3*
a*b^3 + 4*b^6])/(2*b^3))*3^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*a^(Sqrt[-3*a + 4*b^3
]/(2*b^(3/2)))*b^(Sqrt[-3*a*b^3 + 4*b^6]/b^3)*(1/(a*E^(2*b*x)))^(Sqrt[-3*a + 4*b
^3]/(2*b^(3/2)))*BesselI[Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3), (Sqrt[3/2]*Sqrt[a]*Sqrt
[1/(a*E^(2*b*x))]*Sqrt[C[1]])/b]*C[1]^(Sqrt[-3*a + 4*b^3]/(2*b^(3/2)))*Gamma[1 +
 Sqrt[-3*a*b^3 + 4*b^6]/(2*b^3)]))}}

Maple raw input

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

Maple raw output

Intat(-2/(-4*_a*b^3-3*b^2*_a^2+RootOf(6*BesselK(1/2*(-4*b^3+3*a)/(4*b^4-3*a*b)^(
1/2)/b,-1/2/b^2*3^(1/2)*_Z)*_C1*b^2+9*BesselK(1/2*(-4*b^3+3*a)/(4*b^4-3*a*b)^(1/
2)/b,-1/2/b^2*3^(1/2)*_Z)*_C1*b*_a-3*(4*b^4-3*a*b)^(1/2)*BesselK(1/2*(-4*b^3+3*a
)/(4*b^4-3*a*b)^(1/2)/b,-1/2/b^2*3^(1/2)*_Z)*_C1+3*BesselK(1/2*(4*b^3+2*(4*b^4-3
*a*b)^(1/2)*b-3*a)/(4*b^4-3*a*b)^(1/2)/b,-1/2/b^2*3^(1/2)*_Z)*3^(1/2)*_C1*_Z+6*B
esselI(-1/2*(-4*b^3+3*a)/(4*b^4-3*a*b)^(1/2)/b,-1/2/b^2*3^(1/2)*_Z)*b^2+9*Bessel
I(-1/2*(-4*b^3+3*a)/(4*b^4-3*a*b)^(1/2)/b,-1/2/b^2*3^(1/2)*_Z)*b*_a-3*BesselI(1/
2*(4*b^3+2*(4*b^4-3*a*b)^(1/2)*b-3*a)/(4*b^4-3*a*b)^(1/2)/b,-1/2/b^2*3^(1/2)*_Z)
*3^(1/2)*_Z-3*(4*b^4-3*a*b)^(1/2)*BesselI(-1/2*(-4*b^3+3*a)/(4*b^4-3*a*b)^(1/2)/
b,-1/2/b^2*3^(1/2)*_Z))^2-a*b)*b^2,_a = y(x))-x-_C2 = 0, Intat(-2/(-4*_a*b^3-3*b
^2*_a^2+RootOf(6*BesselK(1/2*(-4*b^3+3*a)/(4*b^4-3*a*b)^(1/2)/b,1/2/b^2*3^(1/2)*
_Z)*_C1*b^2+9*BesselK(1/2*(-4*b^3+3*a)/(4*b^4-3*a*b)^(1/2)/b,1/2/b^2*3^(1/2)*_Z)
*_C1*b*_a-3*(4*b^4-3*a*b)^(1/2)*BesselK(1/2*(-4*b^3+3*a)/(4*b^4-3*a*b)^(1/2)/b,1
/2/b^2*3^(1/2)*_Z)*_C1+3*BesselK(1/2*(4*b^3+2*(4*b^4-3*a*b)^(1/2)*b-3*a)/(4*b^4-
3*a*b)^(1/2)/b,1/2/b^2*3^(1/2)*_Z)*3^(1/2)*_C1*_Z+6*BesselI(-1/2*(-4*b^3+3*a)/(4
*b^4-3*a*b)^(1/2)/b,1/2/b^2*3^(1/2)*_Z)*b^2+9*BesselI(-1/2*(-4*b^3+3*a)/(4*b^4-3
*a*b)^(1/2)/b,1/2/b^2*3^(1/2)*_Z)*b*_a-3*BesselI(1/2*(4*b^3+2*(4*b^4-3*a*b)^(1/2
)*b-3*a)/(4*b^4-3*a*b)^(1/2)/b,1/2/b^2*3^(1/2)*_Z)*3^(1/2)*_Z-3*(4*b^4-3*a*b)^(1
/2)*BesselI(-1/2*(-4*b^3+3*a)/(4*b^4-3*a*b)^(1/2)/b,1/2/b^2*3^(1/2)*_Z))^2-a*b)*
b^2,_a = y(x))-x-_C2 = 0