Solve equation with order higher than 32

I have a very high order equation, something like 0=x^50 + a*x^49 + b*x^48 +...+ 1 I tried to use fzero or program Newton Raphson solver to get the solution but matlab seems to have problem dealing with such high power equations.

Any ideas? Many thanks.

Answer 1

[Walter Roberson]

Except under limited conditions, any power that large is probably going to lead to numeric nonsense if you try to solve it that way.

For an equation such as that you should try roots().

Better yet, switch to the Symbolic Toolbox and solve() it and double() the results. The symbolic toolbox operates at higher internal precision, especially if you push up the digits setting.