[_quadrature]
Book solution method
Missing Variables ODE, Dependent variable missing, Solve for
Mathematica ✓
cpu = 2.68298 (sec), leaf count = 828
Maple ✓
cpu = 0.225 (sec), leaf count = 166
DSolve[4*y'[x] + 4*y'[x]^3 == x,y[x],x]
Mathematica raw output
{{y[x] -> ((1 + I*Sqrt[3])*(9*Sqrt[3]*x + Sqrt[64 + 27*x^2])*(-9*x + Sqrt[192 +
81*x^2])^(1/3))/(96*3^(1/6)) + ((I/256)*3^(1/6)*(I + Sqrt[3])*Sqrt[64 + 27*x^2]*
(-512 - 576*x^2 - 243*x^4 + 32*x*Sqrt[192 + 81*x^2] + 27*x^3*Sqrt[192 + 81*x^2])
)/((-9*x + Sqrt[192 + 81*x^2])^(1/3)*(-64 - 27*x^2 + 3*x*Sqrt[192 + 81*x^2])) +
((1 - I*Sqrt[3])*(-9*x + Sqrt[192 + 81*x^2])^(2/3)*(-3584 - 5184*x^2 - 729*x^4 +
480*x*Sqrt[192 + 81*x^2] + 81*x^3*Sqrt[192 + 81*x^2]))/(1536*3^(1/3)*(-32 - 27*
x^2 + 3*x*Sqrt[192 + 81*x^2])) + C[1]}, {y[x] -> ((1 - I*Sqrt[3])*(9*Sqrt[3]*x +
Sqrt[64 + 27*x^2])*(-9*x + Sqrt[192 + 81*x^2])^(1/3))/(96*3^(1/6)) - ((I/256)*3
^(1/6)*(-I + Sqrt[3])*Sqrt[64 + 27*x^2]*(-512 - 576*x^2 - 243*x^4 + 32*x*Sqrt[19
2 + 81*x^2] + 27*x^3*Sqrt[192 + 81*x^2]))/((-9*x + Sqrt[192 + 81*x^2])^(1/3)*(-6
4 - 27*x^2 + 3*x*Sqrt[192 + 81*x^2])) + ((1 + I*Sqrt[3])*(-9*x + Sqrt[192 + 81*x
^2])^(2/3)*(-3584 - 5184*x^2 - 729*x^4 + 480*x*Sqrt[192 + 81*x^2] + 81*x^3*Sqrt[
192 + 81*x^2]))/(1536*3^(1/3)*(-32 - 27*x^2 + 3*x*Sqrt[192 + 81*x^2])) + C[1]},
{y[x] -> (-((9*Sqrt[3]*x + Sqrt[64 + 27*x^2])*(-9*x + Sqrt[192 + 81*x^2])^(1/3))
/(8*3^(1/6)) + ((-9*x + Sqrt[192 + 81*x^2])^(2/3)*(3584 + 5184*x^2 + 729*x^4 - 4
80*x*Sqrt[192 + 81*x^2] - 81*x^3*Sqrt[192 + 81*x^2]))/(128*3^(1/3)*(-32 - 27*x^2
+ 3*x*Sqrt[192 + 81*x^2])) + (3*3^(1/6)*Sqrt[64 + 27*x^2]*(-512 - 576*x^2 - 243
*x^4 + 32*x*Sqrt[192 + 81*x^2] + 27*x^3*Sqrt[192 + 81*x^2]))/(64*(-9*x + Sqrt[19
2 + 81*x^2])^(1/3)*(-64 - 27*x^2 + 3*x*Sqrt[192 + 81*x^2])))/6 + C[1]}}
Maple raw input
dsolve(4*diff(y(x),x)^3+4*diff(y(x),x) = x, y(x),'implicit')
Maple raw output
y(x) = Int(1/12*I/(27*x+3*(81*x^2+192)^(1/2))^(1/3)*((-3^(1/2)+I)*(27*x+3*(81*x^
2+192)^(1/2))^(2/3)-12*I-12*3^(1/2)),x)+_C1, y(x) = Int(1/12*I/(27*x+3*(81*x^2+1
92)^(1/2))^(1/3)*((3^(1/2)+I)*(27*x+3*(81*x^2+192)^(1/2))^(2/3)-12*I+12*3^(1/2))
,x)+_C1, y(x) = Int(1/6*((27*x+3*(81*x^2+192)^(1/2))^(2/3)-12)/(27*x+3*(81*x^2+1
92)^(1/2))^(1/3),x)+_C1