dsolve(diff(y(x),x) = (-y(x)^2+4*a*x)^2/y(x),y(x))