Exercse-2014F-2-Answer.docVIP

Exercise-2 for Numerical Methods (P50(1))Determine rigorously if each function has a unique fixed point on the given interval. g(x)= 1-x2/4 on [0,1]; Solution: g(x)=1/x on [0.5,5.2]. Solution:(P85(1))Let f(x)= x2-x+2. Find the Newton-Raphson formula pk=g(pk-1). Solution: Start with p0=-1.5 and find p1,p2, and p3. Solution: Matlab Exercises: (P51(1))Use the program of fixed point iteration to approximate the fixed points (if any) of each function. Produce a graph of each function and the line y=x that clearly shows any fixed points. g(x)=x5-3x3-2x2+2; g(x)=cos(sin(x)); g(x)=x2-sin(x+0.15); g(x)=xx-cos(x) Solution: (a) See the programs in the M-files. The result is the following: P =0, 2, 2 k =3 p =2 err =0 Actually, at least there are two other roots for this equation as shown in the following plot: (b) See the programs in the M-files. The result is the following: P =0 1.0000 0.6664 0.8150 0.7467 0.7781 0.7636 0.7703 0.7672 0.7686 0.7680 0.7683 0.7681 0.7682 0.7682 0.7682 0.7682 k = 17 p = 0.7682 err =6.1856e-006(c) See the programs in the M-files. The result is the following: Choose p0=1.5, tol=0.00001, and maxNum=100. The iterations cannot converge: 1. We got an error message “maximum number of iterations exceeded”. 2. The last points are P =…. 0.3234 -0.3513 0.3234 -0.3513 0.3234 -0.3513 0.3234 -0.3513 0.3234 -0.3513 The data show that the iterations oscillate between 0.3234 and -0.3513. 3. err =0.6747Choosing p0=0.5, we got the similar result: k =100; p =-0.3513; err =0.6747; P =0.5000 -0.3552 0.3299 … -0.3513 0.3234 -0.3513 0.3234 -0.3513(d) See the programs in the M-files. The result is the following: Choose p0=1.2, tol=0.00001, and maxNum=20. It is convergent. k =17; p =1.0000; err =7.9231e-006 P =1.2000 1.1650 1.1248 1.0850 … 1.0033 1.0015 1.0007 1.0003 1.0002 1.0001 1.0000 1.0000 1.0000 (P87(22))Halley’s method is a way to speed up convergence of Newton-Raphson method. Do programming for the two methods and