[[_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
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