4.21.3 a2y(x)2(b2(cxay(x))2)+c2(b2(cxay(x))2)+2ab2cy(x)=0

ODE
a2y(x)2(b2(cxay(x))2)+c2(b2(cxay(x))2)+2ab2cy(x)=0 ODE Classification

[[_homogeneous, `class C`], _dAlembert]

Book solution method
Clairaut’s equation and related types, d’Alembert’s equation (also call Lagrange’s)

Mathematica
cpu = 1.54616 (sec), leaf count = 71

{{y(x)cc1b2c2(xc1)2a},{y(x)b2c2(xc1)2+cc1a}}

Maple
cpu = 1.334 (sec), leaf count = 195

{a2(y(x))22y(x)acx+c2x22b2a2=0,xy(x)cxaa(2_a2a24b2)c(_a2a22b2a2_a2(_a2a22b2))d_a_C1=0,xy(x)cxaa(2_a2a24b2)c(_a2a22b2+a2_a2(_a2a22b2))d_a_C1=0} Mathematica raw input

DSolve[c^2*(b^2 - (c*x - a*y[x])^2) + 2*a*b^2*c*y'[x] + a^2*(b^2 - (c*x - a*y[x])^2)*y'[x]^2 == 0,y[x],x]

Mathematica raw output

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

Maple raw input

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

Maple raw output

(a^2*y(x)^2-2*y(x)*a*c*x+c^2*x^2-2*b^2)/a^2 = 0, x-Intat(-1/2*a*(_a^2*a^2-2*b^2+
(-a^2*_a^2*(_a^2*a^2-2*b^2))^(1/2))/(_a^2*a^2-2*b^2)/c,_a = y(x)-1/a*c*x)-_C1 = 
0, x-Intat(-1/2*a*(_a^2*a^2-2*b^2-(-a^2*_a^2*(_a^2*a^2-2*b^2))^(1/2))/(_a^2*a^2-
2*b^2)/c,_a = y(x)-1/a*c*x)-_C1 = 0