dsolve(x*diff(diff(y(x),x),x)+(2*a*x*ln(x)+1)*diff(y(x),x)+(a^2*x*ln(x)^2+a*ln(x)+a)*y(x)=0,y(x))