CHAPTER 2 EXTRACTION OF THE QUADRATICS FROM REAL ALGEBRAIC POLYNOMIAL
|
|
- Randell Hancock
- 5 years ago
- Views:
Transcription
1 24 CHAPTER 2 EXTRACTION OF THE QUADRATICS FROM REAL ALGEBRAIC POLYNOMIAL 2.1 INTRODUCTION Polynomial factorization is a mathematical problem, which is often encountered in applied sciences and many of the engineering disciplines. Development of numerical methods for root-extraction of algebraic polynomials is very much essential for higher order polynomial. In general, the root extraction is based on an initial approximation of one or two roots. The convergence on the result depends mainly on the nearness of the approximate root to the actual one. The problem of convergence on inaccurate results or divergence from actual roots due to improper initial approximation which is applied to numerical algorithms is addressed in this chapter. The numerical technique incorporating the methodologies for obtaining close initial approximate quadratic factor and error correcting technique are presented in the subsequent sections. The application of this technique to solve different class of polynomials is the central theme of this research work.
2 Nature of the Problem The literature study reveals that the iterative algorithms for obtaining the roots of an algebraic equation essentially begin with initial approximation and then provide a sequence of iterations converging to a root in the limit. This initial approximation should be sufficiently close to one of the roots to assure the convergence to accurate results. Consider a monic polynomial, P n (x) = a n x n + a n-1 x n a 1 x + a 0 (2.1) where n is degree and a i s are real coefficients with a n =1. The roots of the polynomial equation can be determined by extracting the quadratic factors from the polynomial. The roots of a quadratic polynomial can be determined using a known closed form formula. Let D 2 (x) be an approximate quadratic factor of the polynomial P n (x). If P n (x) is divided by this factor D 2 (x), then a quotient polynomial Q n-2 (x) of degree (n-2) and a remainder polynomial R 1 (x) of degree one are obtained, following the remainder theorem in appendix (A1.2.7). Thus, the polynomial P n (x) can be written as P n (x) = D 2 (x) Q n-2 (x) + R 1 (x) (2.2)
3 26 where D 2 (x) = x 2 + px + q (2.3) Q n-2 (x) = b n-2 x n-2 + b n-3 x n b 0 ; b n-2 =1 (2.4) R 1 (x) = e 1 x + e 2 (2.5) The coefficients e 1 and e 2 in the remainder polynomial R 1 (x) are the errors caused due to the approximation of the coefficients p and q in the divisor D 2 (x). If D 2 (x) is an exact factor of P n (x), then e 1 = e 2 =0. This chapter concentrates on determining a numerical method to arrive at the approximate quadratic factor D 2 (x) and develop an error correction technique which is applied in the successive iterations to refine this D 2 (x) such that e 1 and e 2 in R 1 (x) approach zero, conditioned by a specified error tolerance. 2.2 METHOD FOR THE EXTRACTION OF INITIAL QUADRATIC FACTOR The initial approximate quadratic factor D 2 (x) can be extracted from the coefficients of the polynomial P n (x), a n, a n-1, a n-2 a 1, a 0 by a (n+1) level triangular formation as shown in the Table 2.1. Let A n+1,n+1 denote the two dimensional triangular formation for extracting quadratic factor D 2 (x) from P n (x). a i,j A ( 0 ( 0 and the column index j denotes the coefficients of the polynomial in decending order.
4 27 An upper triangle is formed from 0 th level to n th level of rows by reducing column size by one at every level. The 0 th level of the triangular formation represents the coefficients of P n (x). This brings out the following relation; a 0,n = a n, a 0,(n-1) = a (n-1),, a 0,k = a k,, a 0,0 = a 0 where a 0,k is a k th element in the 0 th level in the triangular formation and a k is the coefficient of x k term in P n (x). The n th level represents the coefficient of P 0 (n). The 0 th level of triangular formation is assigned with the coefficients of the polynomial P n (x) in Equation (2.1) as a 0,j = a j ; j = 0, 1,, n (2.6) Each element in the 1 st level of the triangular formation is determined using the following relation a 1,j = a 0,(j+1) + a 0,j ; j = 0, 1,, n-1 (2.7) This reduces the degree of the polynomial to (n-1). The elements in the 2 nd level are determined using the following relation a 2,j = a 1,(j+1) + a 1,j ; j = 0, 1,, n-2 (2.8) Similarly, the elements in the third row are determined using the following relation a 3,j = a 2,(j+1) + a 2,j ; j = 0, 1,, n-3 (2.9)
5 28 And so on, until n th row of the triangular formation. In general, every element in each row in the triangular formation is generated using the following recurrence relation a (i+1),j = a i,j + a i,(j+1) ;the row index i = 0,1,, n ; the coefficient index j = 0, 1,, n-i (2.10) Table 2.1 Extraction of Initial Approximate Quadratic Factor Level Degree The Coefficients of P n (x) 0 n a 0,n a 0,(n-1) a 0,(n-2) a 0,(n-3) a 0,1 a 0,0 1 n-1 a 1,(n-1) a 1,(n-2) a 1,(n-3) a 1,1 a 1,0 2 n-2 a 2,(n-2) a 2,(n-3) a 2,1 a 2,0 : : : : : : n-2 2 a (n-2),2 a (n-2),1 a (n-2),0 n-1 1 a (n-1),1 a (n-1),0 n 0 a n,0 The initial approximate quadratic factor is extracted at the (n-2) th row. The elements a (n-2),2, a (n-2),1 and a (n-2),0 form the coefficients of the extracted quadratic factor as D 2 (x) = a (n-2),2 x 2 + a (n-2),1 x + a (n-2),0 (2.11) The above Equation (2.11) is made monic as D 2 (x) = x 2 + (a (n-2),1 / a (n-2).2 )x + (a (n-2),0 / a (n-2),2 ) (2.12)
6 29 can be written as In general, the extracted initial approximate quadratic factor D 2 (x) = x 2 +p 0 x+q 0 (2.13) where p 0 = (a (n-2),1 / a (n-2),2 ) and q 0 =(a (n-2),0 / a (n-2),2 ) Quadratic Synthetic Division (Paolo Ruffini 1809) The division of the polynomial Equation (2.1) by the initial quadratic factor Equation (2.13) is carried out by the standard quadratic synthetic division scheme as b k = a k - pb k-1 - qb k-2 ; k = 0, 1,, n (2.14) The quotient polynomial Q n-2 (x) = b n x n-2 +b n-1 x n-3 + +b 2 and the remainder polynomial R 1 (x) = b 1 x+b 0 are obtained from the above Equation (2.14). The remainder coefficients b 1 and b 0 are errors caused due to the approximation of p 0 and q 0 values in D 2 (x). 2.3 ERROR CORRECTION AND CONVERGENCE This iterative numerical technique begins with the approximate coefficients p 0 and q 0 of the extracted initial quadratic factor D 2 (x) in Equation (2.13), which are then successively corrected iteration by iteration such that the error coefficients e 1 =b 1 and e 2 =b 0 in the remainder polynomial obtained from Equation (2.14) should approach zero, within a specified accuracy limit. This can be achieved by applying the following error correcting technique. The bisection principle is also used to converge the quadratic to actual results.
7 Error Correction in the Initial Quadratic Factor Let (p t, q t ) be the true values of p and q in Equation (2.3) and Then p t t (2.15) The computation of the error corrections to the coefficients of initial approximant in Equation (2.3) is formulated as 1 2/n. (2.16) where n is the degree of P n (x) and e 1 and e 2 are the error coefficients in the remainder R 1 (x). 0, q 0 be the initial values of p and q, then the improved values for the first iteration are p 1 = p 0 1 = q 0 1 and q 1 are evaluated, the quadratic synthetic division is repeated, and the next improved p 2 and q 2 are determined from the relation p 2 = p 1 and q 2 = q 1 In general, the improved coefficients for the k th iteration are p k =p k-1 q k =q k-1 (2.17) k-1 and q=q k-1. accuracy.
8 Stopping Criterion Let p k and q k be the coefficients of D 2 (x) at k th iteration and p k+1 and q k+1 be the coefficients of D 2 (x) at (k+1) th iteration. In this error correcting technique, the change in the values from p k to p k+1 and q k to q k+1 in the two consecutive iterations are measured as the absolute errors. The process of iteration is terminated when desired level of accuracy is obtained. An objective criterion for terminating the iteration process is formulated as p k+1 p k a and q k+1 q k a a is an absolute error or Iterations iterations. a considered in this chapter is 10-6 for the printing purpose. Hence the obtained results are rounded up to the accuracy limit of However, the error tolerance can be extended to any higher limit based on the storage class in the implementation Convergence using Bisection Principle Any sign change from p k to p k+1 or q k to q k+1 in the successive iterations observed in Equation (2.17) indicates that the iterative process crosses the actual roots of the quadratic D 2 (x), following intermediate value property in A In order to bracket the roots within the coefficient intervals (p k to p k+1 ) and (q k to q k+1 ), the bisection technique is employed to enhance the convergence rate.
9 32 The midpoints are found to update the coefficients of D 2 (x) for the next iteration as p k+2 = (p k+1 + p k ) / 2 q k+2 = (q k+1 +q k ) / 2 (2.18) The repetition of error correcting technique and bisection principle is to be terminated when p and q have been obtained to the desired accuracy. Equation (2.17) as The convergence rate for p and q are determined from k = p k+1 - p k k = q k+1 - q k Let p k+1 = g(p k ), k = 0,1,2, q k+1 = g(q k ), k = 0,1,2, Let g(p k ) = p k + (e 1 / n) g(q k ) = q k + ( e 2 / n) Therefore g'(p k ) = g'(q k ) = 1. Hence, the scheme converges linearly. Since the convergence is linear to achieve a high degree of accuracy, a large number of iterations may be needed for some class of polynomial. However, the bisection technique is guaranteed to converge.
10 The Deflated Polynomial Let the extracted Quadratic factor be D 2 (x) = (x 2 + px + q) (2.19) The polynomial P n (x) can be divided by the Equation (2.18) using standard synthetic division, following the factor theorem property in A1.2.6 as Q n-2 (x) = P n (x) / D 2 (x) = b n x n-2 +b n-1 x n-3 + +b 3 x+b 2 ; b n =1 (2.20) where Q n-2 (x) is called as the deflated polynomial. In this case, the remainder polynomial R 1 (x) = 0 as the coefficients b 1 and b 0 in Equation (2.20) approach zero. The degree of the polynomial is reduced by 2. The next quadratic factor can be obtained in the similar process from the deflated polynomial. The application of this procedure is repeated until the degree of the polynomial becomes less than or equal to ALGORITHM FOR EXTRACTING ALL QUADRATICS In this section, the computational algorithms for the numerical technique incorporating the methodologies described in sections 2.2 and 2.3 to extract the quadratics of the polynomial P n (x) are presented. These algorithms are tested and analysed for a given P n (x) in the next section.
11 Algorithm for Initial Quadratic Factor The computational algorithm for the methodology formulated in section 2.2 to extract the initial quadratic is presented in this section. The computation involves only two rows in the formation of initial approximate quadratic factor D 2 (x) from the Table 2.1. The second row is determined from the first row and the first row is updated with the second row. The column size is reduced by 1 in successive iteration. The pseudo code of the algorithm for this scheme is deduced as: Algorithm 2.1 (Extracting the initial approximate quadratic D 2 (x) using Triangular Method): This algorithm extracts the initial approximate quadratic factor D 2 (x) from the given polynomial P n (x) = a n x n + a n-1 x n-1 + a n-2 x n-2 + +a 1 x+a 0, where n is the degree and a n, a n-1, a n-2,,a 1 and a 0 are real coefficients. Algorithm Initial_ApproxQuad(P n (x), n) /* The input parameters are the coefficients a n, a n-1,, a 1, a 0 and the degree n of P n (x). */ Step 1: // Triangular formation to extract initial quadratic D 2 (x) For i = 1 to n-2 do // Level of triangular formation For j = 0 to n-i do // Coefficient index Compute bj = a j + a j+1 ; Set a j = b j ; EndFor EndFor Step 2: // Coefficients of initial quadratic D 2 (x) Compute p = a 1 / a 2 and q = a 0 / a 2 ; Step 3: Write p, q; End Initial_ApproxQuad. Figure 2.1 Extracting the Initial Approximate Quadratic D 2 (x) using Triangular Method
12 35 The 0 th level of triangular formation is considered to be the coefficients of P n (x). The step 2 in the above algorithm shown in Figure 2.1, computes the relation for extracting initial quadratic in Equation (2.10). This step performs n(n+1)/2 additions only. The procedure is terminated at (n-2) th row to obtain the three elements to form the initial quadratic, thus saves three addition operations. Therefore, the total addition operations performed by this procedure is ((n(n+1)/2)-3); assume addition operation takes unit time. The step 3 calculates the initial coefficients p and q of D 2 (x) Algorithm for Extracting all Quadratics This algorithm invokes the algorithm 2.1 for extracting initial quadratic factor. The algorithm applies the error correction and convergence technique presented in section 2.3 to refine the initial quadratic to converge to actual D 2 (x) and also determines the deflated polynomial using standard quadratic synthetic division method. Algorithm 2.2 (Error Correction and Convergence of initial D 2 (x) in Triangular Method): This algorithm extracts an actual quadratic factor D 2 (x) = x 2 +px+q from a polynomial P n (x) = a n x n +a n-1 x n-1 + +a 1 x+a 0 of degree n and determines the deflated polynomial Q n-2 (x)=b n x n-2 +b n-1 x n b 3 x +b 2.
13 36 Algorithm Quad_Poly(P n (x),n) /* The input parameters are the coefficients a n, a n-1,, a 1, a 0 and the degree n of P n (x). */ Step 1: Set -6 ; // Accuracy limit Step 2: // Invoke algorithm 2.2 to extract the initial quadratic D 2 (x) Call Initial_ApproxQuad(P n (x),n); Step 3: // Error correction in the initial quadratic D 2 (x) For k = (n-1) to 0 do Compute b k = a k - pb k-1 - qb k-2 ; EndFor Compute 1 0/n; Compute p new new Step 4: // Convergence using Bisection principle If (sign(p) new) or sign(q) new)) Then Compute p new = ( p new + p ) / 2 ; Compute q new = ( q new + q ) / 2 ; EndIf Step 5: // Stopping criterion If ( ( p new or ( q new Then Set p = p new ; and q = q new ; goto 2; Else Write The Extracted Quadratic Factor is, x 2 + p new x +q new ; Write The coefficients of the deflated polynomial Q n-2 (x) are, b n, b n-1,, b 2; EndIf Step 6: // Update P n (x) after extracting D 2 (x) Set n = n-2; Set P n (x) = Q n-2 (x); End Quad_Poly. Figure 2.2 Error Correction and Convergence of Initial D 2 (x) in Triangular Method
14 37 In step 3 in the above algorithm shown in Figure 2.2, the error coefficients b 1 and b 0 in R 1 (x) are obtained by synthetic division of P n (x) by D 2 (x). These values are used for refining the approximate coefficients p new = b 1 /n and q new = b 1 /n in D 2 (x). For the next iteration, the updated new p new and q new of D 2 (x) are used for synthetic division. This process takes O(n) computations. Step 4 checks whether the iteration is crossing the actual results. If so, it applies the bisection technique to converge the iteration to the actual results. The absolute errors are measured with the coefficients p and q for two consecutive iterations as p new p and q new q in step 5. When the errors reach the tolerance limit of 10-6, the process of iteration is terminated and the obtained results are rounded up to this accuracy limit Generalized Algorithm for Solving Algebraic Polynomial The techniques presented in sections 2.2 and 2.3 can also be applied for ill-conditioned polynomial. This section presents a procedure to handle ill-conditioned polynomial and the procedure is incorporated with the methodologies described in sections 2.2 and 2.3 to develop a generalized algorithm for solving a given of polynomial Handling Ill-conditioned Polynomial If any coefficient a i (0 n(x) in Equation (2.1) is found negative, then the following procedure is to be applied before starting the step 2 in the algorithm2.2 Quad_Poly. Let S = n i 1 a (2.21) i where S is the sum of the coefficients of P n (x) except a 0 in Equation (2.1).
15 38 The condition a 0 the ill-conditioning of P n (x). This section handles this condition by inversing the polynomial. Now the Equation (2.1) becomes P n (1/x) = a n (1/x) n+1 + a n-1 (1/x) n a 0 (2.22) The Equation (2.21) can be written as P n (1/x) = a 0 x n + a 1 x n a n-1 x + a n (2.23) Scale the Equation (2.22) by a 0 to convert it into monic polynomial as P n (1/x) = x n + (a 1 /a 0 )x n (a n-1 / a 0 )x + (a n / a 0 ) (2.24) The equation (2.24) can be easily determined from equation (2.23) by the following procedure. Algorithm 2.3 (Handling Ill-conditioned polynomial P n (x) in Triangular Method): This algorithm inverses the ill-conditioned polynomial P n (x) and scales it to convert into monic polynomial, shown in Figure 2.3. Algorithm Inverse_Poly(P n (x),n) /* The input parameters are the coefficients a n, a n-1,, a 1, a 0 and the degree n of P n (x). */ Step 1: For i=0 to n do Set b[n-i] = a[i]; EndFor Step 2: For i= 0 to n do Compute a[i] = b[i] / b[n]; EndFor End Inverse_Poly Figure 2.3 Handling Ill-conditioned Polynomial P n (x) in Triangular Method
16 39 Thus, this section develops a generalized algorithm which extracts all quadratics of any of polynomial of degree n with real coefficients. The algorithm is deduced as Algorithm 2.4 (Generalised Algorithm for Extracting (n/2) quadratics D 2 (x) from P n (x) by Triangular method): This algorithm solves a polynomial P n (x) = a n x n + a n-1 x n-1 + a n-2 x n a 1 x + a 0 of degree n by extracting (n/2) quadratic factors x 2 +px+q. Algorithm Solve_Poly(P n (x),n) /* The input parameters are the coefficients a n, a n-1,, a 1, a 0 and the degree n of P n (x). */ Step 1: // Criterion and Handling the ill-condition of P n (x) If a i < 0 (0 Then Compute S = EndIf If a 0 > S Then n i 1 a ; Call Inverse_Poly(P n (x),n); EndIf Step 2: While (n > 2) Call Quad_Poly(P n (x),n); EndWhile End Solve_Poly. i Figure 2.4 Generalised Algorithm for Extracting (n/2) Quadratics D 2 (x) from P n (x) by Triangular Method
17 40 The step1 in the algorithm shown in Figure 2.4, handles the illconditioning of the polynomial P n (x). If the polynomial is found illconditioned, it is handled in by inversing it. The roots obtained from the ill-conditioned polynomial must be inverted back to get the actual roots. The algorithms 2.1, 2.2, 2.3 and 2.4 presented in this section are simple and straightforward to be implemented using any programming language. 2.5 NUMERICAL ILLUSTRATIONS The numerical procedure presented in this chapter for generating the initial quadratic factor and extracting all quadratics is tested on different class of polynomials in Equation (2.1). The numerical algorithms are implemented in C language and executed in a computing machine with the configuration Intel Core2 Duo CPU, 1.2GHz and 2GB RAM. In this section, three different illustrative test polynomials are taken such as polynomial having complex roots only, polynomial with cluster roots and an ill-conditioned polynomial with roots distributed in all the quadrants of complex plane. The illustrative examples show the performance efficiency of the methodologies present in this chapter Illustration 1 Consider the given polynomial to be: P 8 (x) = x 8 +9x 7 +39x x x x x x+60 (2.25)
18 41 The degree of the above polynomial is even and all the coefficients are positive integers. All the quadratics are extracted, thus the polynomial is solved using the algorithm Quad_Poly(P n (x),n) as explained in The results are tabulated in the Table 2.2. The execution time is measured for solving the polynomial Equation (2.25) in terms of milliseconds. Table 2.2 Extraction of Quadratics of the Polynomial in Equation (2.25) using Triangular Method Order of P n (x) Initial Quadratic D 2(x) Deflated Polynomial Q n-2 (x) Extracted Quadratic D 2 (x) Roots of D 2 (x) Exe. Time (ms) 8 x x x6 + 6x x x x x + 15 x 2 +3x ± j x x x 4 +4x 3 +6x 2 +4x+5 x 2 +2x+3-1 ± j x x x 2 +4x+5 x ± j x 2 +4x+5-2 ± j Illustration 2 Consider the polynomial P 4 (x) = x x x x+4.62 (2.26) The equation has even degree and positive real coefficients. The polynomial is solved using the algorithm 2.3 Quad_Poly(P n (x),n) and the results are tabulated in Table 2.3. The two quadratic factors and two pairs of cluster roots are obtained in milliseconds using the algorithm 2.3. The proposed procedure for solving P 4 (x) in Equation (2.26) is illustrated in Appendix 3.
19 42 Table 2.3 Extraction of Quadratics of the Polynomial in Equation (2.26) using Triangular Method Order of P n (x) Initial Quadratic D 2 (x) Deflated Polynomial Q n-2 (x) Extracted Quadratic D 2 (x) Roots of D 2 (x) Exe. Time (ms) 4 x x x x+4.2 x x+1.1-1, x x+4.2-2, Illustration 3 Consider the polynomial P(x) = x 4 + 4x 3-7x 2-22x + 24 (2.27) The order of the polynomial equation (2.27) is even and it has negative coefficients. It is found that the coefficient a 0 is greater than the sum of other coefficients i.e. a 4 +a 3 +a 2 +a 1. In this case, the algorithm 2.3 Inverse_Poly is applied to solve this above polynomial (2.27). The inverted polynomial is obtained from algorithm 2.3 as P(1/x) = x x x x (2.28) Since the polynomial (2.27) is inverted and the roots of the inverted extracted quadratics are obtained. The roots are again inverted to get the actual roots. This task is accomplished in 0.171milliseconds. And the results are tabulated in the Table 2.4.
20 43 Table 2.4 Extraction of Quadratics of the Polynomial in Equation (2.27) using Triangular Method Order of P n (x) Initial Quadratic D 2 (x) 4 x x Deflated Polynomial Q n-2 (x) x 2-1.5x Extracted Quadratic D 2 (x) x x x 2-1.5x+0.5 Roots of D 2 (x) Exe. Time (ms) -0.25, Inverted Roots: -4, , 0.5 Inverted Roots: 1, RESULTS AND DISCUSSIONS The observations made from the test illustrations 2.5.1, and are discussed here. The execution times measured in illustrations 2.5.1, and are comparable with Bairstow method. The initial approximate quadratic is obtained from the original given polynomial by the methodology discussed in section 2.2. Then quadratic synthetic division is carried out and error correction is applied to refine the original guess quadratic as discussed in section 2.3. Thus, the actual quadratic factor and the deflated polynomial are derived. Again the initial approximate quadratic is obtained from the deflated polynomial and the process is continued until all the quadratics is extracted. This guarantees that the method converges on accurate results. The bisection technique deployed during error correction process enhances the convergence rate.
21 44 In test illustration 2.5.3, it is found that the polynomial has negative coefficients. Then the coefficient value a 0 i.e. 24 is compared with the sum of other coefficients which is -24. In this case, a 0 is not lesser than the sum of remaining coefficients. This shows the illconditioning of the polynomial Equation (2.27). Hence, the polynomial is inverted and the proposed scheme is applied. In order to get the actual roots, the obtained roots are again inverted. 2.7 APPLICATION OF EXTRACTED INITIAL QUADRATIC TO BAIRSTOW METHOD Bairstow (1920) proposed a numerical technique which begins with an initial approximate quadratic to extract the actual quadratics of the polynomial P n (x). Lins (1922) suggested an initial approximant quadratic for Bairstow method. Birtwistle (1967) and Arthur (1972) examined the possibility of improving Bairstow's algorithm for finding the roots of polynomials with real coefficients. However, the Lins- Bairstow method is found divergent or converge to inaccurate results for certain types of polynomials. Hence, the initial approximate quadratic D 2 (x) extracted in section 2.1 can be suitably applied to Bairstow method. This application ensures Bairstow method converges on actual results. The relations in Equations (2.31) and (2.32) follow the Bairstow method to compute the error corrections in the approximate coefficients p and q in D 2 (x). Let the relation to determine c k using b k from Equation (2.4) be
22 45 c k = b k pc k-1 qc k-2 ;k=1,2,,n-1 and c 0 =1, c -1 = 0 (2.29) Let Det = (c n-2 ) 2 - c n-3 (c n-1 b n-1 ) (2.30) and (2.30) are computed as 2.29) nc n-3 - b n-1 c n-2 ) / Det (2.31) and n-1(c n-1 b n-1 ) b n c n-2 ) / Det (2.32) (2.32) are substituted in Equation (2.17) to improve the approximate values p and q of the initial quadratic factor D 2 (x). In step 3 of algorithm 2.2, the relations (2.29) to (2.32) are computed to find new p and q of D 2 (x) for the next iteration. The Bairstow method with the initial approximant D 2 (x) extracted in section 2.2 is tested for different class of polynomial. Three test polynomials are undertaken in this section such as a polynomial with many missing terms and having only complex roots, an illconditioned polynomial having roots distributed in all the quadrants of complex plane and a polynomial with very small cluster roots. In this section, along with the initial approximant D 2 (x), two other approximate quadratics are considered to analyse the performance of Bairstow method in every test. Let L 2 (x) be a Lins approximate which is formed from P n (x). This initial quadratic is written as L 2 (x) = x 2 + l 1 x+ l 2 where l 1 = (a 1 /a 2 ) and l 2 = (a 0 /a 2 ) ; a 2, a 1 and a 0 are the coefficients in P n (x). Let C 2 (x) = x 2 + c 1 x + c 2 be an approximate quadratic where the coefficients c 1 and c 2
23 46 are assigned with constant values i.e c 1 =0.5 and c 2 =0.5. The performances of Bairstow method with the initial approximant D 2 (x), Lins approximation L 2 (x) and a constant approximation C 2 (x) are compared in these tests Illustration 4 Consider the given polynomial to be: P 12 (x) = x 12 +x 8 +x (2.33) The degree of the above polynomial is even and all the coefficients are positive. The given polynomial is found ill-conditioned as many of the terms are zero. The roots of P 12 (x) are complex and distributed in both quadrants: ±0.565j, ±1.5281j, ±2.0756j, ±0.565j, ±1.5281j and ±2.0756j. The results of Bairstow method with initial approximant D 2 (x), Lins approximation and a constant approximation are shown in the Table 2.5. Bairstow method finds all the roots of the polynomial (2.33) in 74 iterations using the initial approximant D 2 (x) whereas it takes 157 iterations with the approximation C 2 (x) and it does not converge with Lins approximation Illustration 5 Consider the polynomial P(x) = x 5 -x 4-4x 3 +4x 2-5x-75 (2.34) The polynomial Equation (2.34) is ill-conditioned as a 0 > a n +a n-1 + +a 1. The roots are real and complex and distributed in
24 47 both quadrants of complex plane: 3, 1±2j and -2±1j. The initial approximant D 2 (x) is used to solve the polynomial. The Bairstow solves the polynomial (2.34) in 27 iterations with D 2 (x) whereas it does not converge for other approximations C 2 (x) and L 2 (x). The results are tabulated in Table 2.6. Table 2.5 Application of Triangular Approximation to Bairstow Method in Equation (2.33) and Comparison with Other Approximations Order P n (x) D 2 (x) C 2 (x) L 2 (x) D 2 (x) C 2 (x) L 2 (x) D 2 (x) C 2 (x) L 2 (x) D 2 (x) C 2 (x) L 2 (x) D 2 (x) C 2 (x) L 2 (x) D 2 (x) C 2 (x) L 2 (x) Initial Quadratic p = q = c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = p = q = c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = p = q = c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = p = q = c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = p = q = c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = p = q = c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = Extracted Factor D 2 (x) x x x x Does not converge x x x x Does not converge x x x x Does not converge x x x x Does not converge x x x x Does not converge x x x x Does not converge Roots of D 2 (x) ±j ±j ±j ±j ± j ± j ± j ± j ± j ± ± j ± j No. of Iterations
25 48 Table 2.6 Application of Triangular Approximation to Bairstow Method in Equation (2.34) and Comparison with other Approximations Order P n (x) Initial Quadratic Extracted Factor D 2 (x) Roots of D 2 (x) No. of Iterations 5 3 p = 0.6 q = 8.2 c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = p = -1.5 q = 1 c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = x x Does not converge Does not converge x x Does not converge Does not converge ±j ±j x Illustration 6 Consider the polynomial P(x) = x x 7 +x 6 +2x 5-2x 4-3x 3 +3x 2 +2x+1 (2.35) The order of the polynomial Equation (2.35) is even and it has negative coefficients. It is found that the polynomial is ill-conditioned as
26 49 the coefficient a 7 is greater than the sum of other coefficients i.e. (a 8 +a 6 +a 5 +a 4 +a 3 +a 2 +a 1 +a 0 ). The polynomial (2.35) has both real and complex roots. These roots are far away distributed in the complex plane: , , ± j , ± j and ± j Matlab does not show the real part of complex roots as considered to be very small compared to other roots. Bairstow method solves this polynomial with the initial approximant D 2 (x) using lesser iterations compared to other approximations. The results are shown in the Table Results and Discussion The initial approximate quadratic D 2 (x) is obtained by the scheme presented in 2.2 from the given polynomial. Then Bairstow procedure is carried out and error correction is applied to refine the original guess quadratic. Thus, the actual quadratic factor and the deflated polynomial are derived. Again the initial approximate quadratic is obtained from the deflated polynomial and the process is continued until all the quadratics are extracted. The Bairstow method is repeated for other initial guesses such as Lins approximation L 2 (x) and constant approximation C 2 (x). Even though Bairstow method has high convergence rate, in some cases it is observed that it may either diverge or converge on inaccurate results. But, the illustrative examples show that the triangular method proposed in the section2.2 guarantees the performance of Bairstow method to converge on accurate results for any class of polynomial whether it is well conditioned or ill-conditioned.
27 50 Table 2.7 Application of Triangular Approximation to Bairstow Method in Equation (2.35) and Comparison with other Approximations Order P n (x) Initial Quadratic p = q = c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = p = q = c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = p = q = c 1 = 0.5 c 2 = 0.5 (a 1 /a 2 ) = (a 0 /a 2 ) = Extracted Factor D 2 (x) x x x x x x x x x x x x x x x x x x x x x x x x Roots of D 2 (x) ±j ±j ±j ±j ±j ±j ±j ±j ±j No. of Iterations In test illustration 4, initial approximant D 2 (x) accelerates the Bairstow method to converge 53% faster than the approximation C 2 (x), whereas Bairstow method does not converge for Lins approximation L 2 (x). Bairstow method converges only for the initial approximant D 2 (x)
28 51 and diverges for other approximations C 2 (x) and L 2 (x) in illustration 5. In this case, the initial approximant D 2 (x) proves 100% better than other approximations. In illustration 6, the quadratic D 2 (x) is 18.46% faster than the approximation C 2 (x) and 94% better than the Lins approximation L 2 (x). Thus, the methodology presented in 2.2 can be treated as a global method for generating initial approximations for solving algebraic polynomial of any class. 2.8 CONCLUSION In this chapter, a numerical technique incorporating the methodologies for extracting initial approximant and correcting residue error to obtain the actual quadratics of an algebraic polynomial is presented in the sections 2.2 and 2.3. The appropriate algorithms are developed in the section 2.4. The error correcting process is enhanced, since the initial quadratic approximant is extracted from the given polynomial. The proposed method is reliable in finding both real and complex roots of a given polynomial up to a chosen accuracy. The method is also suitable for ill-conditioned polynomial. For the ill-conditioned polynomial such as polynomial having very large coefficients, negative coefficients and missing terms, it is suggested in the section that the polynomial can be inverted and scaled, and then the present method can be used to extract all the quadratics. The illustrative examples shown in the section 2.5 prove that the method is simple to be implemented using any programming language. As shown in section 2.7 the proposed procedure for initial approximation can be
29 52 treated as the first step in Bairstow method. The same procedure can be applied to other root-finding numerical procedures. However, the procedure is found slower for certain ill-conditioned polynomials. The reason for the slow convergence is attributed to the fact that the error correction process does not depend upon the relative location of the roots of the polynomial. In chapter 3, a method is proposed in which the extraction of initial approximation and error correction technique follow the characteristic of the roots of the given polynomial. This scheme accelerates the convergence.
Scientific Computing. Roots of Equations
ECE257 Numerical Methods and Scientific Computing Roots of Equations Today s s class: Roots of Equations Polynomials Polynomials A polynomial is of the form: ( x) = a 0 + a 1 x + a 2 x 2 +L+ a n x n f
More information3.3 Dividing Polynomials. Copyright Cengage Learning. All rights reserved.
3.3 Dividing Polynomials Copyright Cengage Learning. All rights reserved. Objectives Long Division of Polynomials Synthetic Division The Remainder and Factor Theorems 2 Dividing Polynomials In this section
More informationZEROS OF POLYNOMIAL FUNCTIONS ALL I HAVE TO KNOW ABOUT POLYNOMIAL FUNCTIONS
ZEROS OF POLYNOMIAL FUNCTIONS ALL I HAVE TO KNOW ABOUT POLYNOMIAL FUNCTIONS TOOLS IN FINDING ZEROS OF POLYNOMIAL FUNCTIONS Synthetic Division and Remainder Theorem (Compressed Synthetic Division) Fundamental
More informationAppendix: Synthetic Division
Appendix: Synthetic Division AP Learning Objectives In this section, we will learn how to: 1. Divide polynomials using synthetic division. Synthetic division is a short form of long division with polynomials.
More informationRoots of Equations Roots of Polynomials
Roots of Equations Roots of Polynomials Content Roots of Polynomials Birge Vieta Method Lin Bairstow Method Roots of Polynomials The methods used to get all roots (real and complex) of a polynomial are:
More informationChapter 6. Nonlinear Equations. 6.1 The Problem of Nonlinear Root-finding. 6.2 Rate of Convergence
Chapter 6 Nonlinear Equations 6. The Problem of Nonlinear Root-finding In this module we consider the problem of using numerical techniques to find the roots of nonlinear equations, f () =. Initially we
More informationMore Polynomial Equations Section 6.4
MATH 11009: More Polynomial Equations Section 6.4 Dividend: The number or expression you are dividing into. Divisor: The number or expression you are dividing by. Synthetic division: Synthetic division
More informationSection 4.3. Polynomial Division; The Remainder Theorem and the Factor Theorem
Section 4.3 Polynomial Division; The Remainder Theorem and the Factor Theorem Polynomial Long Division Let s compute 823 5 : Example of Long Division of Numbers Example of Long Division of Numbers Let
More informationSolution of Algebric & Transcendental Equations
Page15 Solution of Algebric & Transcendental Equations Contents: o Introduction o Evaluation of Polynomials by Horner s Method o Methods of solving non linear equations o Bracketing Methods o Bisection
More informationWarm-Up. Simplify the following terms:
Warm-Up Simplify the following terms: 81 40 20 i 3 i 16 i 82 TEST Our Ch. 9 Test will be on 5/29/14 Complex Number Operations Learning Targets Adding Complex Numbers Multiplying Complex Numbers Rules for
More informationGoals for This Lecture:
Goals for This Lecture: Learn the Newton-Raphson method for finding real roots of real functions Learn the Bisection method for finding real roots of a real function Look at efficient implementations of
More informationProcess Model Formulation and Solution, 3E4
Process Model Formulation and Solution, 3E4 Section B: Linear Algebraic Equations Instructor: Kevin Dunn dunnkg@mcmasterca Department of Chemical Engineering Course notes: Dr Benoît Chachuat 06 October
More informationChapter 3: Root Finding. September 26, 2005
Chapter 3: Root Finding September 26, 2005 Outline 1 Root Finding 2 3.1 The Bisection Method 3 3.2 Newton s Method: Derivation and Examples 4 3.3 How To Stop Newton s Method 5 3.4 Application: Division
More informationLinear System of Equations
Linear System of Equations Linear systems are perhaps the most widely applied numerical procedures when real-world situation are to be simulated. Example: computing the forces in a TRUSS. F F 5. 77F F.
More informationThe purpose of computing is insight, not numbers. Richard Wesley Hamming
Systems of Linear Equations The purpose of computing is insight, not numbers. Richard Wesley Hamming Fall 2010 1 Topics to Be Discussed This is a long unit and will include the following important topics:
More informationM.SC. PHYSICS - II YEAR
MANONMANIAM SUNDARANAR UNIVERSITY DIRECTORATE OF DISTANCE & CONTINUING EDUCATION TIRUNELVELI 627012, TAMIL NADU M.SC. PHYSICS - II YEAR DKP26 - NUMERICAL METHODS (From the academic year 2016-17) Most Student
More informationSection 3.1: Characteristics of Polynomial Functions
Chapter 3: Polynomial Functions Section 3.1: Characteristics of Polynomial Functions pg 107 Polynomial Function: a function of the form f(x) = a n x n + a n 1 x n 1 +a n 2 x n 2 +...+a 2 x 2 +a 1 x+a 0
More informationROOT FINDING REVIEW MICHELLE FENG
ROOT FINDING REVIEW MICHELLE FENG 1.1. Bisection Method. 1. Root Finding Methods (1) Very naive approach based on the Intermediate Value Theorem (2) You need to be looking in an interval with only one
More informationDividing Polynomials: Remainder and Factor Theorems
Dividing Polynomials: Remainder and Factor Theorems When we divide one polynomial by another, we obtain a quotient and a remainder. If the remainder is zero, then the divisor is a factor of the dividend.
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More information6x 3 12x 2 7x 2 +16x 7x 2 +14x 2x 4
2.3 Real Zeros of Polynomial Functions Name: Pre-calculus. Date: Block: 1. Long Division of Polynomials. We have factored polynomials of degree 2 and some specific types of polynomials of degree 3 using
More informationPART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435
PART I Lecture Notes on Numerical Solution of Root Finding Problems MATH 435 Professor Biswa Nath Datta Department of Mathematical Sciences Northern Illinois University DeKalb, IL. 60115 USA E mail: dattab@math.niu.edu
More informationSECTION 2.3: LONG AND SYNTHETIC POLYNOMIAL DIVISION
2.25 SECTION 2.3: LONG AND SYNTHETIC POLYNOMIAL DIVISION PART A: LONG DIVISION Ancient Example with Integers 2 4 9 8 1 In general: dividend, f divisor, d We can say: 9 4 = 2 + 1 4 By multiplying both sides
More informationx 2 + 6x 18 x + 2 Name: Class: Date: 1. Find the coordinates of the local extreme of the function y = x 2 4 x.
1. Find the coordinates of the local extreme of the function y = x 2 4 x. 2. How many local maxima and minima does the polynomial y = 8 x 2 + 7 x + 7 have? 3. How many local maxima and minima does the
More informationNumerical Solution of f(x) = 0
Numerical Solution of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@pdx.edu ME 350: Finding roots of f(x) = 0 Overview Topics covered in these slides
More informationMHF4U Unit 2 Polynomial Equation and Inequalities
MHF4U Unit 2 Polynomial Equation and Inequalities Section Pages Questions Prereq Skills 82-83 # 1ac, 2ace, 3adf, 4, 5, 6ace, 7ac, 8ace, 9ac 2.1 91 93 #1, 2, 3bdf, 4ac, 5, 6, 7ab, 8c, 9ad, 10, 12, 15a,
More informationNumerical Methods. Root Finding
Numerical Methods Solving Non Linear 1-Dimensional Equations Root Finding Given a real valued function f of one variable (say ), the idea is to find an such that: f() 0 1 Root Finding Eamples Find real
More informationLesson 7.1 Polynomial Degree and Finite Differences
Lesson 7.1 Polynomial Degree and Finite Differences 1. Identify the degree of each polynomial. a. 3x 4 2x 3 3x 2 x 7 b. x 1 c. 0.2x 1.x 2 3.2x 3 d. 20 16x 2 20x e. x x 2 x 3 x 4 x f. x 2 6x 2x 6 3x 4 8
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 5. Nonlinear Equations
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T Heath Chapter 5 Nonlinear Equations Copyright c 2001 Reproduction permitted only for noncommercial, educational
More informationCOURSE Iterative methods for solving linear systems
COURSE 0 4.3. Iterative methods for solving linear systems Because of round-off errors, direct methods become less efficient than iterative methods for large systems (>00 000 variables). An iterative scheme
More informationAdvanced Math Quiz Review Name: Dec Use Synthetic Division to divide the first polynomial by the second polynomial.
Advanced Math Quiz 3.1-3.2 Review Name: Dec. 2014 Use Synthetic Division to divide the first polynomial by the second polynomial. 1. 5x 3 + 6x 2 8 x + 1, x 5 1. Quotient: 2. x 5 10x 3 + 5 x 1, x + 4 2.
More informationNotes for Numerical Analysis Math 5465 by S. Adjerid Virginia Polytechnic Institute and State University. (A Rough Draft)
Notes for Numerical Analysis Math 5465 by S. Adjerid Virginia Polytechnic Institute and State University (A Rough Draft) 1 2 Contents 1 Error Analysis 5 2 Nonlinear Algebraic Equations 7 2.1 Convergence
More informationComputational Methods. Solving Equations
Computational Methods Solving Equations Manfred Huber 2010 1 Solving Equations Solving scalar equations is an elemental task that arises in a wide range of applications Corresponds to finding parameters
More informationComplex Numbers: Definition: A complex number is a number of the form: z = a + bi where a, b are real numbers and i is a symbol with the property: i
Complex Numbers: Definition: A complex number is a number of the form: z = a + bi where a, b are real numbers and i is a symbol with the property: i 2 = 1 Sometimes we like to think of i = 1 We can treat
More informationIntro to Scientific Computing: How long does it take to find a needle in a haystack?
Intro to Scientific Computing: How long does it take to find a needle in a haystack? Dr. David M. Goulet Intro Binary Sorting Suppose that you have a detector that can tell you if a needle is in a haystack,
More informationFinding the Roots of f(x) = 0. Gerald W. Recktenwald Department of Mechanical Engineering Portland State University
Finding the Roots of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab:
More informationFinding the Roots of f(x) = 0
Finding the Roots of f(x) = 0 Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@me.pdx.edu These slides are a supplement to the book Numerical Methods with Matlab:
More informationLesson 7.1 Polynomial Degree and Finite Differences
Lesson 7.1 Polynomial Degree and Finite Differences 1. Identify the degree of each polynomial. a. 1 b. 0.2 1. 2 3.2 3 c. 20 16 2 20 2. Determine which of the epressions are polynomials. For each polynomial,
More informationL1 2.1 Long Division of Polynomials and The Remainder Theorem Lesson MHF4U Jensen
L1 2.1 Long Division of Polynomials and The Remainder Theorem Lesson MHF4U Jensen In this section you will apply the method of long division to divide a polynomial by a binomial. You will also learn to
More informationRoots & Zeros of Polynomials. How the roots, solutions, zeros, x-intercepts and factors of a polynomial function are related.
Roots & Zeros of Polynomials How the roots, solutions, zeros, x-intercepts and factors of a polynomial function are related. A number a is a zero or root of a function y = f (x) if and only if f (a) =
More information5.3 The Power Method Approximation of the Eigenvalue of Largest Module
192 5 Approximation of Eigenvalues and Eigenvectors 5.3 The Power Method The power method is very good at approximating the extremal eigenvalues of the matrix, that is, the eigenvalues having largest and
More informationLong and Synthetic Division of Polynomials
Long and Synthetic Division of Polynomials Long and synthetic division are two ways to divide one polynomial (the dividend) by another polynomial (the divisor). These methods are useful when both polynomials
More informationBisection Method. and compute f (p 1 ). repeat with p 2 = a 2+b 2
Bisection Method Given continuous function f (x) on the interval [a, b] with f (a) f (b) < 0, there must be a root in (a, b). To find a root: set [a 1, b 1 ] = [a, b]. set p 1 = a 1+b 1 2 and compute f
More informationLinear System of Equations
Linear System of Equations Linear systems are perhaps the most widely applied numerical procedures when real-world situation are to be simulated. Example: computing the forces in a TRUSS. F F 5. 77F F.
More informationCOURSE Numerical methods for solving linear systems. Practical solving of many problems eventually leads to solving linear systems.
COURSE 9 4 Numerical methods for solving linear systems Practical solving of many problems eventually leads to solving linear systems Classification of the methods: - direct methods - with low number of
More informationExperience in Factoring Large Integers Using Quadratic Sieve
Experience in Factoring Large Integers Using Quadratic Sieve D. J. Guan Department of Computer Science, National Sun Yat-Sen University, Kaohsiung, Taiwan 80424 guan@cse.nsysu.edu.tw April 19, 2005 Abstract
More informationInverses. Stephen Boyd. EE103 Stanford University. October 28, 2017
Inverses Stephen Boyd EE103 Stanford University October 28, 2017 Outline Left and right inverses Inverse Solving linear equations Examples Pseudo-inverse Left and right inverses 2 Left inverses a number
More informationUnit 1: Polynomial Functions SuggestedTime:14 hours
Unit 1: Polynomial Functions SuggestedTime:14 hours (Chapter 3 of the text) Prerequisite Skills Do the following: #1,3,4,5, 6a)c)d)f), 7a)b)c),8a)b), 9 Polynomial Functions A polynomial function is an
More informationL1 2.1 Long Division of Polynomials and The Remainder Theorem Lesson MHF4U Jensen
L1 2.1 Long Division of Polynomials and The Remainder Theorem Lesson MHF4U Jensen In this section you will apply the method of long division to divide a polynomial by a binomial. You will also learn to
More information7.3 The Jacobi and Gauss-Siedel Iterative Techniques. Problem: To solve Ax = b for A R n n. Methodology: Iteratively approximate solution x. No GEPP.
7.3 The Jacobi and Gauss-Siedel Iterative Techniques Problem: To solve Ax = b for A R n n. Methodology: Iteratively approximate solution x. No GEPP. 7.3 The Jacobi and Gauss-Siedel Iterative Techniques
More information3.1 Introduction. Solve non-linear real equation f(x) = 0 for real root or zero x. E.g. x x 1.5 =0, tan x x =0.
3.1 Introduction Solve non-linear real equation f(x) = 0 for real root or zero x. E.g. x 3 +1.5x 1.5 =0, tan x x =0. Practical existence test for roots: by intermediate value theorem, f C[a, b] & f(a)f(b)
More informationMath 3 Variable Manipulation Part 3 Polynomials A
Math 3 Variable Manipulation Part 3 Polynomials A 1 MATH 1 & 2 REVIEW: VOCABULARY Constant: A term that does not have a variable is called a constant. Example: the number 5 is a constant because it does
More informationCS 450 Numerical Analysis. Chapter 5: Nonlinear Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 3 Lecture 3 3.1 General remarks March 4, 2018 This
More informationProcedure for Graphing Polynomial Functions
Procedure for Graphing Polynomial Functions P(x) = a nx n + a n-1x n-1 + + a 1x + a 0 To graph P(x): As an example, we will examine the following polynomial function: P(x) = 2x 3 3x 2 23x + 12 1. Determine
More informationMatrix notation. A nm : n m : size of the matrix. m : no of columns, n: no of rows. Row matrix n=1 [b 1, b 2, b 3,. b m ] Column matrix m=1
Matrix notation A nm : n m : size of the matrix m : no of columns, n: no of rows Row matrix n=1 [b 1, b 2, b 3,. b m ] Column matrix m=1 n = m square matrix Symmetric matrix Upper triangular matrix: matrix
More informationDay 6: 6.4 Solving Polynomial Equations Warm Up: Factor. 1. x 2-2x x 2-9x x 2 + 6x + 5
Day 6: 6.4 Solving Polynomial Equations Warm Up: Factor. 1. x 2-2x - 15 2. x 2-9x + 14 3. x 2 + 6x + 5 Solving Equations by Factoring Recall the factoring pattern: Difference of Squares:...... Note: There
More informationAppendix 8 Numerical Methods for Solving Nonlinear Equations 1
Appendix 8 Numerical Methods for Solving Nonlinear Equations 1 An equation is said to be nonlinear when it involves terms of degree higher than 1 in the unknown quantity. These terms may be polynomial
More informationChapter 1. Root Finding Methods. 1.1 Bisection method
Chapter 1 Root Finding Methods We begin by considering numerical solutions to the problem f(x) = 0 (1.1) Although the problem above is simple to state it is not always easy to solve analytically. This
More informationRobotics 1 Inverse kinematics
Robotics 1 Inverse kinematics Prof. Alessandro De Luca Robotics 1 1 Inverse kinematics what are we looking for? direct kinematics is always unique; how about inverse kinematics for this 6R robot? Robotics
More informationNumerical Analysis using Maple and Matlab
Numerical Analysis using Maple and Matlab Dr. Seongjai Kim Professor of Mathematics Department of Mathematics and Statistics Mississippi State University Mississippi State, MS 39762 skim@math.msstate.edu
More informationChapter 2. Solving Systems of Equations. 2.1 Gaussian elimination
Chapter 2 Solving Systems of Equations A large number of real life applications which are resolved through mathematical modeling will end up taking the form of the following very simple looking matrix
More informationOutline. Scientific Computing: An Introductory Survey. Nonlinear Equations. Nonlinear Equations. Examples: Nonlinear Equations
Methods for Systems of Methods for Systems of Outline Scientific Computing: An Introductory Survey Chapter 5 1 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... Introduce the topic of numerical methods Consider the Error analysis and sources of errors Introduction A numerical method which
More informationEIGENVALUE PROBLEMS. EIGENVALUE PROBLEMS p. 1/4
EIGENVALUE PROBLEMS EIGENVALUE PROBLEMS p. 1/4 EIGENVALUE PROBLEMS p. 2/4 Eigenvalues and eigenvectors Let A C n n. Suppose Ax = λx, x 0, then x is a (right) eigenvector of A, corresponding to the eigenvalue
More information15 Nonlinear Equations and Zero-Finders
15 Nonlinear Equations and Zero-Finders This lecture describes several methods for the solution of nonlinear equations. In particular, we will discuss the computation of zeros of nonlinear functions f(x).
More informationChapter 2 Formulas and Definitions:
Chapter 2 Formulas and Definitions: (from 2.1) Definition of Polynomial Function: Let n be a nonnegative integer and let a n,a n 1,...,a 2,a 1,a 0 be real numbers with a n 0. The function given by f (x)
More information6.5 Dividing Polynomials
Name Class Date 6.5 Dividing Polynomials Essential Question: What are some ways to divide polynomials, and how do you know when the divisor is a factor of the dividend? Explore Evaluating a Polynomial
More informationWe say that a polynomial is in the standard form if it is written in the order of decreasing exponents of x. Operations on polynomials:
R.4 Polynomials in one variable A monomial: an algebraic expression of the form ax n, where a is a real number, x is a variable and n is a nonnegative integer. : x,, 7 A binomial is the sum (or difference)
More informationCollege Algebra Notes
Metropolitan Community College Contents Introduction 2 Unit 1 3 Rational Expressions........................................... 3 Quadratic Equations........................................... 9 Polynomial,
More informationCourse Notes: Week 1
Course Notes: Week 1 Math 270C: Applied Numerical Linear Algebra 1 Lecture 1: Introduction (3/28/11) We will focus on iterative methods for solving linear systems of equations (and some discussion of eigenvalues
More informationChapter 2 Solutions of Equations of One Variable
Chapter 2 Solutions of Equations of One Variable 2.1 Bisection Method In this chapter we consider one of the most basic problems of numerical approximation, the root-finding problem. This process involves
More informationOutline. Math Numerical Analysis. Intermediate Value Theorem. Lecture Notes Zeros and Roots. Joseph M. Mahaffy,
Outline Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Lecture Notes Zeros and Roots Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationChapter 4 No. 4.0 Answer True or False to the following. Give reasons for your answers.
MATH 434/534 Theoretical Assignment 3 Solution Chapter 4 No 40 Answer True or False to the following Give reasons for your answers If a backward stable algorithm is applied to a computational problem,
More informationS56 (5.1) Polynomials.notebook August 25, 2016
Q1. Simplify Daily Practice 28.6.2016 Q2. Evaluate Today we will be learning about Polynomials. Q3. Write in completed square form x 2 + 4x + 7 Q4. State the equation of the line joining (0, 3) and (4,
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Taylor s Theorem Can often approximate a function by a polynomial The error in the approximation
More informationCLASS NOTES Models, Algorithms and Data: Introduction to computing 2018
CLASS NOTES Models, Algorithms and Data: Introduction to computing 2018 Petros Koumoutsakos, Jens Honore Walther (Last update: April 16, 2018) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More informationDividing Polynomials
3-3 3-3 Dividing Polynomials Warm Up Lesson Presentation Lesson Quiz Algebra 2 Warm Up Divide using long division. 1. 161 7 2. 12.18 2.1 23 5.8 Divide. 3. 4. 6x + 15y 3 7a 2 ab a 2x + 5y 7a b Objective
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationPOLYNOMIALS. Maths 4 th ESO José Jaime Noguera
POLYNOMIALS Maths 4 th ESO José Jaime Noguera 1 Algebraic expressions Book, page 26 YOUR TURN: exercises 1, 2, 3. Exercise: Find the numerical value of the algebraic expression xy 2 8x + y, knowing that
More informationPolynomial and Rational Functions. Copyright Cengage Learning. All rights reserved.
2 Polynomial and Rational Functions Copyright Cengage Learning. All rights reserved. 2.3 Real Zeros of Polynomial Functions Copyright Cengage Learning. All rights reserved. What You Should Learn Use long
More informationLetting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc.
1 Polynomial Matrices 1.1 Polynomials Letting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc., n ws ( ) as a
More informationThe minimal polynomial
The minimal polynomial Michael H Mertens October 22, 2015 Introduction In these short notes we explain some of the important features of the minimal polynomial of a square matrix A and recall some basic
More informationA Fast Algorithm for Computing High-dimensional Risk Parity Portfolios
A Fast Algorithm for Computing High-dimensional Risk Parity Portfolios Théophile Griveau-Billion Quantitative Research Lyxor Asset Management, Paris theophile.griveau-billion@lyxor.com Jean-Charles Richard
More informationLesson 2.1: Quadratic Functions
Quadratic Functions: Lesson 2.1: Quadratic Functions Standard form (vertex form) of a quadratic function: Vertex: (h, k) Algebraically: *Use completing the square to convert a quadratic equation into standard
More informationChapter 3: Polynomial and Rational Functions
Chapter 3: Polynomial and Rational Functions 3.1 Polynomial Functions and Their Graphs A polynomial function of degree n is a function of the form P (x) = a n x n + a n 1 x n 1 + + a 1 x + a 0 The numbers
More information5.1 Monomials. Algebra 2
. Monomials Algebra Goal : A..: Add, subtract, multiply, and simplify polynomials and rational expressions (e.g., multiply (x ) ( x + ); simplify 9x x. x Goal : Write numbers in scientific notation. Scientific
More informationETNA Kent State University
C 8 Electronic Transactions on Numerical Analysis. Volume 17, pp. 76-2, 2004. Copyright 2004,. ISSN 1068-613. etnamcs.kent.edu STRONG RANK REVEALING CHOLESKY FACTORIZATION M. GU AND L. MIRANIAN Abstract.
More informationComputational Methods for Engineers Programming in Engineering Problem Solving
Computational Methods for Engineers Programming in Engineering Problem Solving Abu Hasan Abdullah January 6, 2009 Abu Hasan Abdullah 2009 An Engineering Problem Problem Statement: The length of a belt
More informationx 3 2x = (x 2) (x 2 2x + 1) + (x 2) x 2 2x + 1 = (x 4) (x + 2) + 9 (x + 2) = ( 1 9 x ) (9) + 0
1. (a) i. State and prove Wilson's Theorem. ii. Show that, if p is a prime number congruent to 1 modulo 4, then there exists a solution to the congruence x 2 1 mod p. (b) i. Let p(x), q(x) be polynomials
More informationCLASS NOTES Computational Methods for Engineering Applications I Spring 2015
CLASS NOTES Computational Methods for Engineering Applications I Spring 2015 Petros Koumoutsakos Gerardo Tauriello (Last update: July 2, 2015) IMPORTANT DISCLAIMERS 1. REFERENCES: Much of the material
More informationScientific Computing. Roots of Equations
ECE257 Numerical Methods and Scientific Computing Roots of Equations Today s s class: Roots of Equations Bracketing Methods Roots of Equations Given a function f(x), the roots are those values of x that
More informationGENG2140, S2, 2012 Week 7: Curve fitting
GENG2140, S2, 2012 Week 7: Curve fitting Curve fitting is the process of constructing a curve, or mathematical function, f(x) that has the best fit to a series of data points Involves fitting lines and
More informationNumerical Analysis Fall. Roots: Open Methods
Numerical Analysis 2015 Fall Roots: Open Methods Open Methods Open methods differ from bracketing methods, in that they require only a single starting value or two starting values that do not necessarily
More informationSynthetic Substitution
Write your questions and thoughts here! 7.2 Synthetic and Long Polynomial Division 1 Use direct substitution to evaluate t a = a! a! + a + 5 when a = - 2 Synthetic Substitution Synthetic substitution is
More informationChapter 8. Exploring Polynomial Functions. Jennifer Huss
Chapter 8 Exploring Polynomial Functions Jennifer Huss 8-1 Polynomial Functions The degree of a polynomial is determined by the greatest exponent when there is only one variable (x) in the polynomial Polynomial
More informationLearning Objectives. Zeroes. The Real Zeros of a Polynomial Function
The Real Zeros of a Polynomial Function 1 Learning Objectives 1. Use the Remainder and Factor Theorems 2. Use the Rational Zeros Theorem to list the potential rational zeros of a polynomial function 3.
More information1) Synthetic Division: The Process. (Ruffini's rule) 2) Remainder Theorem 3) Factor Theorem
J.F. Antona 1 Maths Dep. I.E.S. Jovellanos 1) Synthetic Division: The Process (Ruffini's rule) 2) Remainder Theorem 3) Factor Theorem 1) Synthetic division. Ruffini s rule Synthetic division (Ruffini s
More informationON THE GLOBAL KRYLOV SUBSPACE METHODS FOR SOLVING GENERAL COUPLED MATRIX EQUATIONS
ON THE GLOBAL KRYLOV SUBSPACE METHODS FOR SOLVING GENERAL COUPLED MATRIX EQUATIONS Fatemeh Panjeh Ali Beik and Davod Khojasteh Salkuyeh, Department of Mathematics, Vali-e-Asr University of Rafsanjan, Rafsanjan,
More information