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