How to find the roots of a polynomial of n^(th)nth degree?
2 Answers
It gets more complicated fast as
Explanation:
Linear
ax+b = 0ax+b=0
Solution
Quadratic
ax^2+bx+c = 0ax2+bx+c=0
Solutions
Cubic
ax^3+bx^2+cx+d = 0ax3+bx2+cx+d=0
Substitute
at^3+b_1t+c_1 = 0at3+b1t+c1=0
If this has exactly one Real root, then I like to use Cardano's method, substituting
au^3+av^3+(3auv+b_1)(u+v)+c_1 = 0au3+av3+(3auv+b1)(u+v)+c1=0
then adding the constraintv = -b_1/(3au)v=−b13au to eliminate the(u+v)(u+v) term. Multiplying through byu^3u3 results in a quadratic inu^3u3 , whose two rootsr_1r1 andr_2r2 areu^3u3 andv^3v3 .
Hence:
x_1 = -b/(3a)+root(3)(r_1)+root(3)(r_2)x1=−b3a+3√r1+3√r2
is the Real root, and the Complex roots are given by:
x_2 = -b/(3a)+omega root(3)(r_1)+omega^2 root(3)(r_2)x2=−b3a+ω3√r1+ω23√r2
x_3 = -b/(3a)+omega^2 root(3)(r_1) + omega root(3)(r_2)x3=−b3a+ω23√r1+ω3√r2
where
There are other methods for solving cubics. If the cubic has
cos(3 theta) = 4cos^3(theta)-3cos(theta)cos(3θ)=4cos3(θ)−3cos(θ) .
Of course this is not quite as algebraically nice as expressions in terms of
Quartic
ax^4+bx^3+cx^2+dx+e = 0ax4+bx3+cx2+dx+e=0
Substitute
t^4+a_1t^2+b_1t+c_1 = 0t4+a1t2+b1t+c1=0
Next consider possible quadratic factors of this:
t^4+a_1t^2+b_1t+c_1 = (t^2+At+B)(t^2-At+C)t4+a1t2+b1t+c1=(t2+At+B)(t2−At+C)
=t^4+(B+C-A^2)t^2+A(C-B)t+BC=t4+(B+C−A2)t2+A(C−B)t+BC
Hence:
B+C = a_1+A^2B+C=a1+A2
C-B = b_1/AC−B=b1A
BC= c_1BC=c1
So:
(a_1+A^2)^2 = (B+C)^2 = (C-B)^2 + 4BC = b_1^2/(A^2) + 4c_1(a1+A2)2=(B+C)2=(C−B)2+4BC=b21A2+4c1
Hence:
A^2(a_1+A^2)^2-4c_1A^2-b_1^2 = 0A2(a1+A2)2−4c1A2−b21=0
If you multiply this out, you get a cubic in
Quintic
ax^5+bx^4+cx^3+dx^2+ex+f = 0ax5+bx4+cx3+dx2+ex+f=0
In general this has no algebraic solution in terms of
For example, the roots of
You can find approximations for the roots of
Explanation:
Let
Given a first approximation
a_(i+1) = a_i - f(a_i)/(f'(a_i))
to find successively better approximations.
This works with both Real and Complex zeros.
To find different roots, start with different initial approximations.
To find reasonable initial approximations, look at the behaviour of the function a little, perhaps evaluating