print(solve(x=default(realprecision,99),2,x^3-x-1))
\q