4.16.8 ax2+by(x)+y(x)2=0

ODE
ax2+by(x)+y(x)2=0 ODE Classification

[[_homogeneous, `class G`]]

Book solution method
No Missing Variables ODE, Solve for y

Mathematica
cpu = 599.992 (sec), leaf count = 0 , timed out

$Aborted

Maple
cpu = 0.249 (sec), leaf count = 280

{_bx1_a2aby(x)(_a2aby(x)_a2y(x))1d_a+y(x)(ax2_fbx2_f)1_bx1(_fb1_a2a_fb2_a2a_fb)(_a2a_fb_a2_f)2d_ad_f+_C1=0,_bx1_a2aby(x)(_a2aby(x)_a+2y(x))1d_a+y(x)(ax2_fbx+2_f)1_bx1(_fb1_a2a_fb+2_a2a_fb)(_a2a_fb_a+2_f)2d_ad_f+_C1=0} Mathematica raw input

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

Mathematica raw output

$Aborted

Maple raw input

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

Maple raw output

Int(-1/((-_a^2*a-b*y(x))^(1/2)*_a-2*y(x))*(-_a^2*a-b*y(x))^(1/2),_a = _b .. x)+I
ntat(1/((-a*x^2-_f*b)^(1/2)*x-2*_f)-Int((-b*_f/(-_a^2*a-_f*b)^(1/2)-2*(-_a^2*a-_
f*b)^(1/2))/((-_a^2*a-_f*b)^(1/2)*_a-2*_f)^2,_a = _b .. x),_f = y(x))+_C1 = 0, I
nt(-1/((-_a^2*a-b*y(x))^(1/2)*_a+2*y(x))*(-_a^2*a-b*y(x))^(1/2),_a = _b .. x)+In
tat(-1/((-a*x^2-_f*b)^(1/2)*x+2*_f)-Int((b*_f/(-_a^2*a-_f*b)^(1/2)+2*(-_a^2*a-_f
*b)^(1/2))/((-_a^2*a-_f*b)^(1/2)*_a+2*_f)^2,_a = _b .. x),_f = y(x))+_C1 = 0