Exam 1(a) Wed 20 Feb 2002

1. Consider the representation of floating-point numbers with base 10 and 2 digits in the fraction part. The values for the exponents are between -10 and +10.

   1. What is the machine precision in this number system?

   2. Represent the numbers 17 and 333 as floating point numbers and illustrate the calculation of 17+333, using rounding. What is the calculated sum?

                                                           +--------+
                                                           |    /15 |
                                                           +--------+
   

2. Below make the plot of g(x) = 10/(x-3). Starting at x_0 = 4, illustrate on the plot below how to produce four more points defined by x(k+1) = g(x(k)), k=0,1,...

   Compute the convergence rate of this fixed-point iteration

   1. for x=-2:

   2. for x=5:

   What conclusions can you make from the rates you computed above?

                                                           +--------+
                                                           |    /20 |
                                                           +--------+
   

3. Apply 3 steps of the golden section search method to find the minimum of the function f(x) = x^2 - 3x in the interval [0,5].

   Write the values for a, b, x_1, x_2, f(x_1), and f(x_2) in the table (4 decimal places):

   +------+---------+---------+---------+---------+----------+----------+
   | step |    a    |    b    |   x_1   |   x_2   |  f(x_1)  |  f(x_2)  |
   +------+---------+---------+---------+---------+----------+----------+
   |   0  |   0.000 |   5.000 |         |         |          |          |
   |   1  |         |         |         |         |          |          |
   |   2  |         |         |         |         |          |          |
   +------+---------+---------+---------+---------+----------+----------+
   

                                                           +--------+
                                                           |    /15 |
                                                           +--------+
   

4. Consider the matrix

              [  4.447E-01   7.919E+08  ]
          A = [                         ]
              [  6.154E-01   9.218E-01  ]
   
    with its inverse
   
              [ -1.891E-09   1.625E+00  ]
     A^(-1) = [                         ].
              [  1.263E-09  -9.124E-10  ]
   
   

   1. Compute the condition number of A using the norm ||.||_1.

   2. Suppose we wish to solve the system A x = b, using the matrix A from above. Assuming a relative error of 10^(-16) on the coefficients of the matrix, what is the bound on the relative error of the solution?

                                                           +--------+
                                                           |    /20 |
                                                           +--------+
   

5. Consider the matrix

             [   2.000     8.000     0.000   ]
         A = [   9.000     4.000     4.000   ].
             [   6.000     7.000     8.000   ]
   

   1. Compute the LU decomposition of A with partial pivoting.

      Calculate with four decimal places, using rounding: write the answer of every step rounded to four decimal places, and use the rounded number in the calculations of the next step.

   2. What is the determinant of A?

                                                           +--------+
                                                           |    /30 |
                                                           +--------+