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