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