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