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