Metodos Numericos
Enviado por diegoquizhpi • 11 de Diciembre de 2014 • 1.440 Palabras (6 Páginas) • 137 Visitas
CHAPTER 22
22.1 Analytical solution:
The first iteration involves computing 1 and 2 segment trapezoidal rules and combining them as
and computing the approximate error as
The computation can be continues as in the following tableau until a < 0.5%.
1 2 3
n a 1.6908% 0.0098%
1 27.62500000 25.87500000 25.83456463
2 26.31250000 25.83709184
4 25.95594388
The true error of the final result can be computed as
22.2 Analytical solution:
1 2 3 4
t 5.8349% 0.1020% 0.0004%
n a 26.8579% 0.3579% 0.0015862%
1 90.38491615 43.57337260 41.21305531 41.17125852
2 55.27625849 41.36057514 41.17191160
4 44.83949598 41.18370307
8 42.09765130
22.3
1 2 3
n a 7.9715% 0.0997%
1 1.34376994 1.97282684 1.94183605
2 1.81556261 1.94377297
4 1.91172038
22.4 Change of variable:
Therefore, the transformed function is
Two-point formula:
Three-point formula:
Four-point formula:
22.5 Change of variable:
Therefore, the transformed function is
Two-point formula:
Three-point formula:
Four-point formula:
22.6 Change of variable:
Therefore, the transformed function is
Five-point formula:
22.7 Here is the Romberg tableau for this problem.
1 2 3
n a 5.5616% 0.0188%
1 224.36568786 288.56033084 289.43080513
2 272.51167009 289.37640049
4 285.16021789
Therefore, the estimate is 289.430805.
22.8 Change of variable:
Therefore, the transformed function is
Two-point formula:
The remaining formulas can be implemented with the results summarized in this table.
n Integral t
2 1.5 39.95%
3 3.1875 27.60%
4 2.189781 12.34%
5 2.671698 6.95%
6 2.411356 3.47%
Thus, the results are converging, but at a very slow rate. Insight into this behavior can be gained by looking at the function and its derivatives.
We can plot the second derivative as
The second and higher derivatives are large. Thus, the integral evaluation is inaccurate because the error is related to the magnitudes of the derivatives.
22.9 (a)
We can use 8 applications of the extended midpoint rule.
This result is close to the analytical solution
(b)
For the first part, we can use 4 applications of Simpson’s 1/3 rule
For the second part,
...