solve((b*y*u)^2=((b^2-a^2)*u+a*x)^2*(1-u^2),u)