# modelans - Natural Sciences Tripos Part IB Numerical...

This preview shows pages 1–3. Sign up to view the full content.

Natural Sciences Tripos Part IB Numerical Methods Model answers to some questions 6.1 Roots of an integrated quantity Consider the integral (29 (29 Fx f a x = ξξ d , where f ( x ) is an easily evaluated function. (a) Show how F ( x ) may be calculated from f ( x ) for some arbitrary value of x using the Trapezium Rule with the interval [ a,x ] subdivided into n– 1 subintervals. (b) Derive the error term in this approximation and show how Romberg Integration may be used to improve the accuracy of the solution given estimates of F ( x ) obtained from n– 1 and 2 n– 1 subintervals. Show that this estimate is equivalent to calculating the integral using Simpson’s Rule. (c) Suppose we wish to find the value of x such that F ( x ) = 1 . Using the Newton-Raphson method, predict the location of the root using a single iteration. You may assume F ( x ) is a monotonically increasing function of x from x = a to the neighbourhood of the root x = x* . Why is this assumption important? (d) Describe how the second and subsequent iterations may be calculated. Discuss how precisely the integral F ( x ) should be evaluated for each iteration and suggest an appropriate method for achieving this. Answer: (a) Let G n ( x ) be our approximation to F ( x ) using the Trapezium Rule with n- 1 subintervals each of a size ∆ξ = ( x-a ) / ( n- 1) . [Suggestion: draw a sketch of some function f( x ) showing how G n ( x ) is obtained using the Trapezium Rule. You may wish to state the integral from ξ to ξ + ∆ξ and then sum for the compound rule, or simply quote the compound rule as here] Then (29 ( 29 [] G x fa fx n =+ + + + + + - + ∆ξ ∆ξ ∆ξ ∆ξ 2 22 2 2 K (b) To obtain the error in G n ( x ) , consider first the error in one of the subintervals from ξ to ξ + ∆ξ , say, using a Taylor Series expansion of f we can see ( 29 (29 (29 (29 ( 29 ff f f f O f f O ′′ + + ′′′ + + + + + ∫∫ ξ ξ ξ ξ ξ ξ ξ ξ ξ dd ∆ξ ∆ξ ∆ξ ∆ξ ∆ξ ∆ξ ∆ξ ∆ξ ∆ξ ∆ξ ∆ξ 23 4 0 4 5 26 262 4 Comparing this with the Trapezium Rule approximation for this step

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
(29 (29 ( 29 [] (29 (29 (29 ff f f f f O
This is the end of the preview. Sign up to access the rest of the document.

## This note was uploaded on 05/01/2011 for the course CHBE 2120 taught by Professor Gallivan during the Spring '07 term at Georgia Institute of Technology.

### Page1 / 4

modelans - Natural Sciences Tripos Part IB Numerical...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document
Ask a homework question - tutors are online