4.3.5 y(x)=y(x)a+by(x)

ODE
y(x)=y(x)a+by(x) ODE Classification

[_quadrature]

Book solution method
Separable ODE, Independent variable missing

Mathematica
cpu = 0.0310508 (sec), leaf count = 27

{{y(x)asech2(12a(c1+x))b}}

Maple
cpu = 0.015 (sec), leaf count = 24

{x+21aArtanh(a+by(x)a)+_C1=0} Mathematica raw input

DSolve[y'[x] == y[x]*Sqrt[a + b*y[x]],y[x],x]

Mathematica raw output

{{y[x] -> -((a*Sech[(Sqrt[a]*(x + C[1]))/2]^2)/b)}}

Maple raw input

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

Maple raw output

x+2/a^(1/2)*arctanh((a+b*y(x))^(1/2)/a^(1/2))+_C1 = 0