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