The Q-D Method


This program approximates the roots of a polynomial. The approximate roots are helpful as guesses for the Bairstow method.


Introduction

The Q-D method program approximates the roots of a polynomial whose coefficients you enter. It is able to identify the roots without an initial estimate but requires many iterations to achieve accuracy comparable to the Bairstow method. Because the Bairstow method can fail to converge for unrealistic initial guesses, the Q-D method is recommended for identifying initial guesses for use in the Bairstow program.

The program can usually determine the roots with sufficient accuracy for the Bairstow program in 20 to 30 iterations. Finding roots with four-digit accuracy can take more than 100 iterations.

Before running the program, you must partition the calculator for at least 125 data registers, which is the default partitioning.

Because the calculato has built-in solutions for quadratic and cubic polynomials, the program requires you enter a polynomial of fourth order or higher. The largest polynomial the program can use is 25th order.

The polynomial you enter has the form

0 = anxn + … a3x3 + a2x2 + a1x + a0


The program generates values labeled e and q. To correctly interpret these values, you should be familiar with the principles of the Q-D method of finding roots.

Reference

Applied Numerical Analysis, Curtis F. Gerald, Addison-Wesley, 1970.

Starting the Program

To start the program that uses the Q-D method:
  1. Select { ZRO } from the MATHEMATICS menu.
    The calculator displays the FUNCTION ZEROS menu.
  2. Press { Q-D }.
    The program displays:
    D_Cd1NwV
  3. Enter the number of iterations you want performed and press { #it }.
  4. Press { EOD }.
    The program displays:
    D_dqBliN
  5. Select the appropriate options.
  • To enter a new polynomial, press { YES } and proceed to "Entering the Coefficients" on the next page.
  • To use a polynomial you have already entered (using either the Bairstow program or a previous run of this program), press { NO } and proceed to the section "Displaying the Approximate Roots".

Entering the Coefficients

When you select { YES } from the ENTER POLY? menu, the program displays
D_KFjTxu
  1. Enter the degree of the polynomial and press { n }.
  2. Press { EOD }.
    The program prompts you to enter the coefficient of each term, beginning with the first term, for example, for a 4th order polynomial, the following would be displayed.
    D_BzDOOG
  3. Enter the coefficient and press { ENT }.
  4. Repeat step 3 until all coefficients are entered.
    The program then displays the menu shown in the next section.

The Edit Menu

After you enter all the coefficients, the program displays an EDIT? menu that allows you to change any of the values.
D_ynZfG5
  • If you do not want to edit, press { NO } and go to the section "Displaying the Approximate Roots" below. The computation time for the approximate roots in some cases can be several minutes.
  • If you want to edit, press { YES } and use the procedure in the next section below.

Editing the Coefficients

When you select { YES } from the EDIT? menu, the program displays the same menu that you used to enter the coefficients, for example:
D_BzDOOG
  • To display the current value of the coefficient, press the [ CE ] key.
  • To accept the current value and proceed to the next coefficient, press { ENT }.
  • To edit the current value, enter the correct value and press { ENT }.
When you press { ENT } to accept or edit the value of the last coefficient, the program returns to the EDIT? menu.

Displaying the Approximate Roots

After you select { NO } from the either the ENTER POLY? or the EDIT? menu, the program displays the value of the first e. (If you are not using a printer and you plan to compute values for r and s, keep a written record of the e's and q's.)
D_7PChkK
  1. Press { NXT }.
    The program displays the corresponding value of q. Each q is an approximate root.
    D_o3qhOv
  2. Press { NXT }.
    • If there are remaining values for e and q, the program displays the next value of e.
    • If all values of e and q have been displayed, the program displays the menu shown in the next section "Selecting Options".

To correctly interpret the e's and q's, you must examine them after an initial number of iteration, and again after more iterations. From the trends evident in checking these values twice (or more), you can infer whether the roots are real or complex.

If all the e's are approaching zero, it is likely that all the roots are real, and the q's are the approximate roots. If an e is not approaching zero, it implies that a pair of complex roots is present.

Selecting Options

After the program has displayed all the values of e and q, it displays:
D_oxJGY5
{ Q-D } - Lets you perform additional iterations.
{ r/s } - Lets you compute the r and s values for use in the Bairstow program.
{ ESC } - Returns to the Q-D METHOD menu.

Performing Additional Iterations

To perform additional iterations using the Q-D method:
  1. Select the { Q-D } option.
    The program displays:
    D_1oM66K
  2. Enter the number of additional iterations you want performed and press { #it }.
  3. Press { EOD }.
    The program displays the approximate roots, as described in the section "Displaying the Approximate Roots" above.

Computing the r and s values

If you are using the Q-D program to compute r and s values for use in the Bairstow program, follow these steps.
  1. Choose a pair of roots from the list of e and q values.
  2. Select the { r/s } option.
    The program displays:
    D_CRBSLm
  3. Enter the number of one of the roots you have chose and press { 1st }. For example, if the root is the third in the list, enter 3.
  4. Enter the number of the other root and press { 2nd }.
  5. Press { EOD } to display the r value, for example:
    D_wy5kNr
  6. Press { NXT } to see the s value.
  7. Press { NXT } again to return to the COMPUTE r, s menu.


Back