I am trying to find all real solutions of the system equations
sol := solve(And(g'(x)=0,g''(x)<>0),x) assuming real;
I tried
restart; fprime := x-> x^6-(3/2)*x^5+2*x^4+(5/2)*x^3-7*x^2+2: f := unapply(simplify(int(fprime(x), x)), x): g := unapply(expand(f(x^2+2*x)), x): sol := solve(And(g'(x)=0,g''(x)<>0),x) assuming real; evalf(sol);
I don't get only real solutions. How can I get only real solutions?