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