4.41.45 y(x)(a0+a1sin2(y(x)))+a2y(x)(a1sin2(y(x))+a3)+a1y(x)2+a1y(x)2sin(y(x))cos(y(x))=0

ODE
y(x)(a0+a1sin2(y(x)))+a2y(x)(a1sin2(y(x))+a3)+a1y(x)2+a1y(x)2sin(y(x))cos(y(x))=0 ODE Classification

[[_2nd_order, _missing_x]]

Book solution method
TO DO

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

$Aborted

Maple
cpu = 3.671 (sec), leaf count = 818

{y(x)((a0+a1)(tan(_h))2+a0)ea1arctan((a0+a1)tan(_h)1a0(a0+a1))1a0(a0+a1)12((((tan(_h))2(sin(_h))2_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h+(sin(_h))2_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h)a1+((tan(_h))2_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h+_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h)a3)a2a0+a12a2(tan(_h))2(sin(_h))2_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h+a1a2a3(tan(_h))2_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h_C1/2)((tan(_h))2+1)((a0+a1)(tan(_h))2+a0)d_hx_C2=0,y(x)((a0+a1)(tan(_h))2+a0)ea1arctan((a0+a1)tan(_h)1a0(a0+a1))1a0(a0+a1)12((((tan(_h))2(sin(_h))2_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h+(sin(_h))2_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h)a1+((tan(_h))2_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h+_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h)a3)a2a0+a12a2(tan(_h))2(sin(_h))2_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h+a1a2a3(tan(_h))2_h((tan(_h))2+1)(a0+a1(sin(_h))2)(ea1a0(a0+a1)arctan((a0+a1)tan(_h)a0(a0+a1)))2d_h_C1/2)((tan(_h))2+1)((a0+a1)(tan(_h))2+a0)d_hx_C2=0} Mathematica raw input

DSolve[a2*(a3 + a1*Sin[y[x]]^2)*y[x] + a1*y'[x]^2 + a1*Cos[y[x]]*Sin[y[x]]*y'[x]^2 + (a0 + a1*Sin[y[x]]^2)*y''[x] == 0,y[x],x]

Mathematica raw output

$Aborted

Maple raw input

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

Maple raw output

Intat(exp(a1/(a0*(a0+a1))^(1/2)*arctan((a0+a1)*tan(_h)/(a0*(a0+a1))^(1/2)))/(-2*
(((Int(1/(tan(_h)^2+1)*exp(a1/(a0*(a0+a1))^(1/2)*arctan((a0+a1)*tan(_h)/(a0*(a0+
a1))^(1/2)))^2*tan(_h)^2/(a0+a1*sin(_h)^2)*sin(_h)^2*_h,_h)+Int(1/(tan(_h)^2+1)*
exp(a1/(a0*(a0+a1))^(1/2)*arctan((a0+a1)*tan(_h)/(a0*(a0+a1))^(1/2)))^2/(a0+a1*s
in(_h)^2)*sin(_h)^2*_h,_h))*a1+(Int(1/(tan(_h)^2+1)*exp(a1/(a0*(a0+a1))^(1/2)*ar
ctan((a0+a1)*tan(_h)/(a0*(a0+a1))^(1/2)))^2*tan(_h)^2/(a0+a1*sin(_h)^2)*_h,_h)+I
nt(1/(tan(_h)^2+1)*exp(a1/(a0*(a0+a1))^(1/2)*arctan((a0+a1)*tan(_h)/(a0*(a0+a1))
^(1/2)))^2/(a0+a1*sin(_h)^2)*_h,_h))*a3)*a2*a0+a1^2*a2*Int(1/(tan(_h)^2+1)*exp(a
1/(a0*(a0+a1))^(1/2)*arctan((a0+a1)*tan(_h)/(a0*(a0+a1))^(1/2)))^2*tan(_h)^2/(a0
+a1*sin(_h)^2)*sin(_h)^2*_h,_h)+a1*a2*a3*Int(1/(tan(_h)^2+1)*exp(a1/(a0*(a0+a1))
^(1/2)*arctan((a0+a1)*tan(_h)/(a0*(a0+a1))^(1/2)))^2*tan(_h)^2/(a0+a1*sin(_h)^2)
*_h,_h)-1/2*_C1)*(tan(_h)^2+1)*((a0+a1)*tan(_h)^2+a0))^(1/2)*((a0+a1)*tan(_h)^2+
a0),_h = y(x))-x-_C2 = 0, Intat(-exp(a1/(a0*(a0+a1))^(1/2)*arctan((a0+a1)*tan(_h
)/(a0*(a0+a1))^(1/2)))/(-2*(((Int(1/(tan(_h)^2+1)*exp(a1/(a0*(a0+a1))^(1/2)*arct
an((a0+a1)*tan(_h)/(a0*(a0+a1))^(1/2)))^2*tan(_h)^2/(a0+a1*sin(_h)^2)*sin(_h)^2*
_h,_h)+Int(1/(tan(_h)^2+1)*exp(a1/(a0*(a0+a1))^(1/2)*arctan((a0+a1)*tan(_h)/(a0*
(a0+a1))^(1/2)))^2/(a0+a1*sin(_h)^2)*sin(_h)^2*_h,_h))*a1+(Int(1/(tan(_h)^2+1)*e
xp(a1/(a0*(a0+a1))^(1/2)*arctan((a0+a1)*tan(_h)/(a0*(a0+a1))^(1/2)))^2*tan(_h)^2
/(a0+a1*sin(_h)^2)*_h,_h)+Int(1/(tan(_h)^2+1)*exp(a1/(a0*(a0+a1))^(1/2)*arctan((
a0+a1)*tan(_h)/(a0*(a0+a1))^(1/2)))^2/(a0+a1*sin(_h)^2)*_h,_h))*a3)*a2*a0+a1^2*a
2*Int(1/(tan(_h)^2+1)*exp(a1/(a0*(a0+a1))^(1/2)*arctan((a0+a1)*tan(_h)/(a0*(a0+a
1))^(1/2)))^2*tan(_h)^2/(a0+a1*sin(_h)^2)*sin(_h)^2*_h,_h)+a1*a2*a3*Int(1/(tan(_
h)^2+1)*exp(a1/(a0*(a0+a1))^(1/2)*arctan((a0+a1)*tan(_h)/(a0*(a0+a1))^(1/2)))^2*
tan(_h)^2/(a0+a1*sin(_h)^2)*_h,_h)-1/2*_C1)*(tan(_h)^2+1)*((a0+a1)*tan(_h)^2+a0)
)^(1/2)*((a0+a1)*tan(_h)^2+a0),_h = y(x))-x-_C2 = 0