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