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)