MS QP
MS QP
Important Instructions:
• Show at least 3 decimal places (floating point representations) for all results and intermediate
steps. You do not need to round-off after every floating-point operation.
• You may use a scientific calculator. Programmable calculators are not allowed.
• You need to show how you arrived at a number using the algorithm. Only showing the final
numbers will be interpreted as being copied from a neighbour and no marks will be awarded.
• This paper contains total 4 questions on two sides of the pages.
5 1 3
𝐴 = [2 7 6]
3 5 3
b) Using the matrices Q and R computed in (a), back-substitution algorithm and common matrix
operations (i.e., addition, subtraction, multiplication and transpose), compute the inverse of matrix A.
Do not use any other algorithm to compute inverse at any stage. 10
a) Compute the values of a and b such that the square of the Euclidean norm of the error (between the
true and the approximate function) is minimum. 10
b) Compute the values of a and b such that the square of the Euclidean semi-norm of the error
(between the true and the approximate function) in the net t = {0, 0.5, 1, 1.5, 2} is minimum. 10
c) If the maximum norm of the error occurs at t = 2.0 for the polynomial obtained in (a) and at t = 1.5
for the polynomial obtained in (b), compute the maximum true relative errors (in %) for these two
polynomials at the points of their respective maximum norm of the errors. From the result, comment
on which polynomial is a better fit to the function with respect to the maximum norm of error. 5
(b) Find an expression for the remainder term for your algorithm of computing bi’s if (x3 – α1x2 – α2x –
α3) is not an exact factor. 5
(c) Use the algorithm derived in (a) to estimate the other two roots of the polynomial shown below if
three of the roots are 1, –1 and 2. 5
x5 – 4x4 + 5x3 – 6x + 4 = 0
4x3 – 27xy2 + 25 = 0
4x2y – 3y3 – 1 = 0
b) The equation ex – x2 – 2x – 2 = 0 has a root near x = 2.5. In order to compute the root using the
Fixed Point method with the initial guess x0 = 2.5, which of the following iteration scheme(s) can be
used? If there are more than one usable scheme(s), which one will converge faster? Justify your
answer. 13
𝑒 𝑥𝑘 −2
(iii) 𝑥𝑘+1 = 𝑥𝑘 +2
(iv) 𝑥𝑘+1 = 𝑙𝑛(𝑥𝑘 2 + 2𝑥𝑘 + 2)
Following three algorithms are given which may be useful to solve the questions:
Algorithm 1
𝑏𝑛 𝑏𝑖 − ∑𝑛𝑗=𝑖+1 𝑎𝑖𝑗 𝑥𝑗
𝑥𝑛 = ; 𝑥𝑖 = ; 𝑖 = (𝑛 − 1), (𝑛 − 2), … 3, 2, 1
𝑎𝑛𝑛 𝑎𝑖𝑖
Algorithm 2
Algorithm 3
𝑘
(1) 𝒂(1) (𝑘+1) (𝑘+1) 𝑇 𝒛(𝑘+1)
𝒒 = ; 𝒛 =𝒂 − ∑(𝒂(𝑘+1) 𝒒(𝑖) )𝒒(𝑖) ; 𝒒(𝑘+1) =
‖𝒂(1) ‖2 𝑖=1
‖𝒛(𝑘+1) ‖2
𝑇
𝑟𝑖𝑗 = 𝒒(𝑖) 𝒂(𝑗)