Find the complete solution of x3dx3d3y+3x2dx2d2y+xdxdy+y=xlogx.
Technique
This is a Cauchy–Euler equation. Substitute x=et (so t=logx) to convert it to a constant-coefficient ODE in t, using xkDk=θ(θ−1)⋯(θ−k+1) with θ=dtd.
Solution
With x=et, θ=dtd, the standard operator identities are
xD=θ,x2D2=θ(θ−1),x3D3=θ(θ−1)(θ−2).
The left side becomes
θ(θ−1)(θ−2)+3θ(θ−1)+θ+1.
Expand:
θ(θ−1)(θ−2)=θ3−3θ2+2θ,3θ(θ−1)=3θ2−3θ,sum=θ3−3θ2+2θ+3θ2−3θ+θ+1=θ3+1.
So the equation transforms to (RHS: xlogx=ett)
(θ3+1)y=tet.
Complementary function. Auxiliary equation m3+1=0⇒(m+1)(m2−m+1)=0, roots
m=−1,m=21±i3.
Thus
yc=C1e−t+et/2[C2cos23t+C3sin23t].
Back to x (et=x, t=logx):
yc=xC1+x[C2cos(23logx)+C3sin(23logx)].
Particular integral.yp=θ3+11tet. Shift: θ3+11ett=et(θ+1)3+11t. Now (θ+1)3+1=θ3+3θ2+3θ+2; applied to t keep terms up to first order (θ2t=θ3t=0):
2+3θ+⋯1t=211+23θ+⋯1t=21(1−23θ+⋯)t=21(t−23).
So
yp=et⋅21(t−23)=2x(logx−23)=2xlogx−43x.