OFFSET
0,1
COMMENTS
The other real (negative) root is -A060007.
One of the pair of complex conjugate roots is obtained by negating sqrt(2*u) and sqrt(u) in the formula for r below, giving 0.248126062... - 1.033982060...*i.
FORMULA
r = (-sqrt(2)*u + sqrt(sqrt(2*u) - 2*u^2))/(2*sqrt(u)), with u = (Ap^(1/3) + ep*Am^(1/3))/3, where Ap = (3/16)*(9 + sqrt(3*283)), Am = (3/16)*(9 - sqrt(3*283)) and ep = (-1 + sqrt(3)*i)/2, with i = sqrt(-1). For the trigonometric version set u = (2/3)*sqrt(3)*sinh((1/3)*arcsinh((3/16)* sqrt(3))).
EXAMPLE
r = 0.724491959000515611588372282187036565786494481350011017270...
MATHEMATICA
First[RealDigits[x/.N[{x->Root[-1+#1+#1^4 &, 2, 0]}, 75]]] (* Stefano Spezia, Aug 27 2022 *)
PROG
(PARI) solve(x=0, 1, x^4 + x - 1) \\ Michel Marcus, Aug 28 2022
CROSSREFS
KEYWORD
AUTHOR
Wolfdieter Lang, Aug 27 2022
STATUS
approved