Scribd
Upload
68 views16 pages

NMpolynom 2

The document describes an algorithm for adapting Newton's method to find roots of polynomials. It presents the steps of the algorithm, which take in the polynomial coefficients, an initial guess, and other parameters and iteratively find better approximations to the roots. It then shows an implementation of the algorithm in Mathematica. Finally, it provides 6 examples of polynomials and applies Newton's method as implemented in Mathematica to find the real roots of each polynomial.

Uploaded by

Vasea Gurduza
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
68 views16 pages

NMpolynom 2

The document describes an algorithm for adapting Newton's method to find roots of polynomials. It presents the steps of the algorithm, which take in the polynomial coefficients, an initial guess, and other parameters and iteratively find better approximations to the roots. It then shows an implementation of the algorithm in Mathematica. Finally, it provides 6 examples of polynomials and applies Newton's method as implemented in Mathematica to find the real roots of each polynomial.

Uploaded by

Vasea Gurduza
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 16

Adaption of Newton’s method to

polynomial rootfinding

E. Cojuhari Mathematics in Computer Sciences


Algorithm

Polynew (a, n, x0 ,, itmax, root, b, ier)


1 Remark: a is a vector of coefficients, itmax the maximum
number of iterates to be computed, b the vector of
coefficients for the deflated polynomial, and ier an error
indicator.
2 itnum := 1
3 z := x0 , bn := c := an
4 For k = n − 1, . . . , 1, bk := ak + z · bk +1 c := bk + z · c
5 b0 := a0 + z · b1
6 If c = 0, ier := 2 and exit.
7 x1 := x0 − b0 /c
8 If |x1 − x0 | ≤ , then ier := 0, root := x1 , and exit.
9 If itnum = itmax, then ier := 1 and exit.
10 Otherwise, itnum := itnum + 1, x0 := x1 , and go to step 3.

E. Cojuhari Mathematics in Computer Sciences


Wolfram Mathematica algorithm realization

E. Cojuhari Mathematics in Computer Sciences


[Atkinson], p. 124, 125, N33(a, b), N38(a, b), N40 (a, b).

Use Newton’s method to calculate the real roots of the


following polynomials as accurately as possible. Estimate the
multiplicity of each root, and then if necessary, try an alternative
way of improving your calculated values.
1 x 4 − 3.2x 3 + 0.96x 2 + 4.608x − 3.456;
2 x 5 + 0.9x 4 − 1.62x 3 − 1.458x 2 + 0.6561x + 0.59049;
3 x 4 − 5.7x 3 − 0.47x 2 + 29.865x − 26.1602;
4 x 4 − 5.4x 3 + 10.56x 2 − 8.954x + 2.7951;
5 676039x 12 − 1939938x 10 + 2078505x 8 − 1021020x 6 +
225225x 4 − 18918x 2 + 231;
6 x 4 − 4.096152422706631x 3 + 3.284232335022705x 2 +
4.703847577293368x − 5.715767664977249.

E. Cojuhari Mathematics in Computer Sciences


x 4 − 3.2x 3 + 0.96x 2 + 4.608x − 3.456

E. Cojuhari Mathematics in Computer Sciences


x 4 − 3.2x 3 + 0.96x 2 + 4.608x − 3.456

E. Cojuhari Mathematics in Computer Sciences


x 5 + 0.9x 4 − 1.62x 3 − 1.458x 2 + 0.6561x + 0.59049

E. Cojuhari Mathematics in Computer Sciences


x 5 + 0.9x 4 − 1.62x 3 − 1.458x 2 + 0.6561x + 0.59049

E. Cojuhari Mathematics in Computer Sciences


x 4 − 5.7x 3 − 0.47x 2 + 29.865x − 26.1602

E. Cojuhari Mathematics in Computer Sciences


x 4 − 5.7x 3 − 0.47x 2 + 29.865x − 26.1602

E. Cojuhari Mathematics in Computer Sciences


x 4 − 5.4x 3 + 10.56x 2 − 8.954x + 2.7951

E. Cojuhari Mathematics in Computer Sciences


x 4 − 5.4x 3 + 10.56x 2 − 8.954x + 2.7951

E. Cojuhari Mathematics in Computer Sciences


676039x 12 − 1939938x 10 + 2078505x 8 − 1021020x 6 + 225225x 4 −
18918x 2 + 231

E. Cojuhari Mathematics in Computer Sciences


676039x 12 − 1939938x 10 + 2078505x 8 − 1021020x 6 + 225225x 4 −
18918x 2 + 231

E. Cojuhari Mathematics in Computer Sciences


x 4 − 4.096152422706631x 3 + 3.284232335022705x 2 +
4.703847577293368x − 5.715767664977249

E. Cojuhari Mathematics in Computer Sciences


x 4 − 4.096152422706631x 3 + 3.284232335022705x 2 +
4.703847577293368x − 5.715767664977249

E. Cojuhari Mathematics in Computer Sciences

You might also like