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