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