[[_3rd_order, _missing_y]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.0235126 (sec), leaf count = 42
Maple ✓
cpu = 0.029 (sec), leaf count = 34
DSolve[4*x^2*y'[x] - 4*x^3*y''[x] + 4*x^4*y'''[x] == 1,y[x],x]
Mathematica raw output
{{y[x] -> -1/(36*x) + (x^2*(2*C[1] - C[2]))/4 + C[3] + (x^2*C[2]*Log[x])/2}}
Maple raw input
dsolve(4*x^4*diff(diff(diff(y(x),x),x),x)-4*x^3*diff(diff(y(x),x),x)+4*x^2*diff(y(x),x) = 1, y(x),'implicit')
Maple raw output
y(x) = 1/36*(18*x^3*_C1*ln(x)-1+(-9*_C1+18*_C2)*x^3+36*_C3*x)/x