El siguiente comando

syms x real;
f = @(x) log(x^2)*exp(-1/(x^2));
fp(x) = diff(f(x),x);
fpp(x) = diff(fp(x),x);

y

solve(fpp(x)>0,x,'Real',true)

devolver el resultado

resolver ([0.0 <(8.0 * exp (-1.0 / x ^ 2)) / x ^ 4 - (2.0 * exp (-1.0 / x ^ 2)) / x ^ 2 - (6.0 * log (x ^ 2) * exp (-1.0 / x ^ 2)) / x ^ 4 + (4.0 * log (x ^ 2) * exp (-1.0 / x ^ 2)) / x ^ 6], [x == RD_NINF..RD_INF ])

que no es lo que espero

La primera pregunta:¿Es posible forzar a Matlabsolve devolver el conjunto de todas las soluciones?

(Esto está relacionado conesta pregunta.) Además, cuando trato de resolver la ecuación

solve(fpp(x)==0,x,'Real',true)

que vuelve

ans =
-1.5056100417680902125994180096313

No estoy satisfecho ya que no se devuelven todas las soluciones (son aproximadamente -1.5056, 1.5056, -0.5663 y 0.5663 obtenidas de WolframAlpha).

Yo sé esovpasolve con alguna conjetura inicial puede manejar esto. Pero no tengo ideacómo generalmente puedo encontrar valores iniciales adivinados para obtener todas las soluciones, que es mi segunda pregunta.

Se agradecen otras soluciones o sugerencias para resolver estos problemas.