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