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