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