4.2.10 y(x)=a0+a1y(x)+a2y(x)2

ODE
y(x)=a0+a1y(x)+a2y(x)2 ODE Classification

[_quadrature]

Book solution method
Separable ODE, Independent variable missing

Mathematica
cpu = 0.0369065 (sec), leaf count = 54

{{y(x)4a0a2a12tan(12(c1+x)4a0a2a12)a12a2}}

Maple
cpu = 0.012 (sec), leaf count = 41

{x214a0a2a12arctan(2a2y(x)+a14a0a2a12)+_C1=0} Mathematica raw input

DSolve[y'[x] == a0 + a1*y[x] + a2*y[x]^2,y[x],x]

Mathematica raw output

{{y[x] -> (-a1 + Sqrt[-a1^2 + 4*a0*a2]*Tan[(Sqrt[-a1^2 + 4*a0*a2]*(x + C[1]))/2]
)/(2*a2)}}

Maple raw input

dsolve(diff(y(x),x) = a0+a1*y(x)+a2*y(x)^2, y(x),'implicit')

Maple raw output

x-2/(4*a0*a2-a1^2)^(1/2)*arctan((2*a2*y(x)+a1)/(4*a0*a2-a1^2)^(1/2))+_C1 = 0