What's the best way to solve the following with $x\sim \text{Normal}(0,1)$?
$$E_x[\log(1-2 a x^2 + a^2 x^4)]=0$$
Idea 1: FindRoot
a0 = 2.42;
term[a_] :=
TransformedDistribution[Log[1 - 2 a x^2 + a^2 x^4],
x \[Distributed] NormalDistribution[]];
FindRoot[NExpectation[t, t \[Distributed] term[a]], {a, a0}]
(* NIntegrate failed to converge to prescribed accuracy *)
Idea 2: increase MaxRecursion
, AccuracyGoal
a0 = 2.42;
method = {"NIntegrate", {MinRecursion -> 12, MaxRecursion -> 30,
WorkingPrecision -> 30, AccuracyGoal -> 12}};
term[a_] :=
TransformedDistribution[Log[1 - 2 a x^2 + a^2 x^4],
x \[Distributed] NormalDistribution[]];
FindRoot[
NExpectation[t, t \[Distributed] term[a], Method -> method], {a, a0}]
(* {a -> 2.42125 - 2.16855*10^-9 I} *)
Idea 3: look at density to find tricks
a0 = 2.42;
func = With[{a = a0},
PDF[TransformedDistribution[Log[1 - 2 a x^2 + a^2 x^4],
x \[Distributed] NormalDistribution[]]]];
Plot[func[x], {x, -3, 3}]
Now it seems an integral splitting technique would work, but I can't get PDF for symbolic $a$, any ideas?
Background: this gives largest $a$ such that $w=w-a w x^2$ converges almost surely for any starting $w$