Tutorial-2-Nonlinear 16/08/24, 12:11 PM
Tutorial-2
Problem 1
Use the Newton-Raphson method and a four-function calculator (+ − ×÷) to compute
√75
3
with four significant figure accuracy.
Problem 2
a. Show that if A is any positive number, then the sequence defined by
1 A
xn = xn−1 + , for n ≥ 1
2 2xn−1
converges to √A whenever x0 > 0.
b. What happens if x0 < 0? Explain.
Problem 3
Find an approximation to √3 correct to within 10−3 using the Bisection method.
Problem 4
The function
f(x) = x3 − 1.2x2 − 8.19x + 13.23
has a double root close to x = 2.
1. Find that root for ϵs = 10−7 using regular Newton-Raphson python code. How many
iterations are needed?
2. Modify the code to implement:
f(xi )
xi+1 = xi − m
f ′ (xi )
where m = 2 (corresponding to double roots). How many iterations are needed to
have ϵs = 10−7 .
Problem 5
http://localhost:8888/nbconvert/html/Root/Tutorial-2-Nonlinear.ipynb?download=false Page 1 of 2
Tutorial-2-Nonlinear 16/08/24, 12:11 PM
A cable is suspended as shown in the figure. Its length s and sag h are related to the
span L by
2 1
and h = (cosh
λL λL
s= sinh − 1) ,
λ 2 λ 2
where
λ = w0 /T0 (1)
w0 = weight of cable per unit length (2)
T0 = cable tension at O (3)
Compute s for L = 160 m and h = 15 m. Use python code for Ridder method.
Problem 6
Bernoulli's equation for fluid flow in an open channel with a small bump is
Q2 Q2
+ h0 = +h+H
2gb2 h20 2gb2 h2
where
Q = 1.2 m3 /s = volume of rate flow
g = 9.81 m/s2 = gravitational acceleration
b = 1.8 m = width of the channel
h0 = 0.6 m = upstream water level
H = 0.075 m = height of bump
h = water level above the bump
Determine h. Use brentq function from scipy library.
http://localhost:8888/nbconvert/html/Root/Tutorial-2-Nonlinear.ipynb?download=false Page 2 of 2