[[_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
Maple ✓
cpu = 0.646 (sec), leaf count = 849
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