Hello
load(to_poly_solver)$
eq:(2*x-3)^(1/2)=(x-2)^(1/2)$
sol:to_poly_solve(eq,x)$
sol:subst(%union="[",sol)$
map(rhs,flatten(sol));
gives the output
[1]
The solution 1 leads to non-real values on the left and the right sides of the equation ( e.g. (2*1-3)^(1/2) can't be computed in the real numbers ).
I'd like to eliminate such solutions. So the solution set of the equation above should be [ ] (empty set).
Can you help me, please?
By the way, what is the solution set of x*(2*x-3)^(1/2)=x*(x-2)^(1/2) if non-real values are not allowed? What would you prefer?