4.2.7 y(x)=a+bx+cy(x)2

ODE
y(x)=a+bx+cy(x)2 ODE Classification

[_Riccati]

Book solution method
Riccati ODE, Generalized ODE

Mathematica
cpu = 0.0220407 (sec), leaf count = 93

{{y(x)b(c1Ai(c(a+bx)(bc)2/3)+Bi(c(a+bx)(bc)2/3))(bc)2/3(c1Ai(c(a+bx)(bc)2/3)+Bi(c(a+bx)(bc)2/3))}}

Maple
cpu = 0.107 (sec), leaf count = 85

{y(x)=1b1c3(Ai(1)((bx+a)(b1c)23)_C1+Bi(1)((bx+a)(b1c)23))1c(Ai((bx+a)(b1c)23)_C1+Bi((bx+a)(b1c)23))1} Mathematica raw input

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

Mathematica raw output

{{y[x] -> (b*(AiryBiPrime[-((c*(a + b*x))/(-(b*c))^(2/3))] + AiryAiPrime[-((c*(a
 + b*x))/(-(b*c))^(2/3))]*C[1]))/((-(b*c))^(2/3)*(AiryBi[-((c*(a + b*x))/(-(b*c)
)^(2/3))] + AiryAi[-((c*(a + b*x))/(-(b*c))^(2/3))]*C[1]))}}

Maple raw input

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

Maple raw output

y(x) = (b/c^(1/2))^(1/3)*(AiryAi(1,-(b*x+a)/(b/c^(1/2))^(2/3))*_C1+AiryBi(1,-(b*
x+a)/(b/c^(1/2))^(2/3)))/c^(1/2)/(AiryAi(-(b*x+a)/(b/c^(1/2))^(2/3))*_C1+AiryBi(
-(b*x+a)/(b/c^(1/2))^(2/3)))