ME-261
Numerical Analysis
Topic – Solution of System of Linear and Nonlinear Equations
(Iterative methods), and Linear Regression Analysis
1. The following system of linear equations is generated by applying the Kirchhoff's law to the
circuit shown in Fig. 1,
−25𝐼1 − 4𝐼3 + 29𝐼4 = 100
55𝐼1 − 25𝐼4 = −200
−37𝐼3 − 4𝐼4 = −250
Figure 1: Circuit diagram
where the I’s are the current in the loop in Amperes. Solve the above system of linear equations
using the Gauss-Seidel method. Write numerical values with an accuracy of up to three decimal
places and perform at least five iterations. [for initial guess, take I1(0) = I3(0) = I4(0) = 1]
2. Determine the roots of the following system of non-linear equations using Newton-Raphson
methods:
𝑥2 + 𝑦2 = 4
𝑒𝑥 + 𝑦 = 1
Use initial guesses as [x(0), y(0)] = [1, −1.7]. Perform only one set of iteration.
3. An investigator has reported the data tabulated below for an experiment to determine the growth
rate of bacteria, k (per day), as a function of oxygen concentration, c (mg/L). It is known that
such data can be modeled by the following equation:
𝑘𝑚𝑎𝑥 𝑐 2
𝑘=
𝑐𝑠 + 𝑐 2
where cs and kmax are the parameters. Use linear regression analysis to determine these parameters,
and predict the growth rate at c = 2 mg/L.
C (mg/L) 0.5 0.8 1.5 2.5 4.0 6.0
k (per day) 1.1 2.4 5.3 7.6 8.9 10.9
ME-261
Numerical Analysis (Extra Problems)
Topic – Solution of System of Linear and Nonlinear Equations
(Iterative methods) and Linear Regression Analysis
Problem 1:
Three masses are suspended vertically by a series of identical springs where mass-1 is at the top
and mass-3 is at the bottom. From the free-body diagram of each mass, we get the following system
of linear equations,
𝑘1 𝑥1 + 𝑘2 (𝑥1 − 𝑥2 ) = 𝑚1 𝑔
𝑘2 (𝑥2 − 𝑥1 ) + 𝑘3 (𝑥2 − 𝑥3 ) = 𝑚2 𝑔
𝑘3 (𝑥3 − 𝑥2 ) = 𝑚3 𝑔
2
If g = 9.81 m/s , m1 = 2 kg, m2 = 3 kg, m3 = 2.5 kg, and the k’s = 10 N/m, where mi = Mass of
block i, k = spring constant and g = gravitational acceleration, then solve the above system of
linear equations using Gauss-Seidel method with Successive Over Relaxation (SOR). Perform
up to 10 iterations and take λ = 1.2. Write numerical values with an accuracy of up to three decimal
places. Assume any suitable initial guess.
Problem 2:
Determine the roots of the following non-linear system of equations using Newton-Raphson
method:
𝑥𝑦𝑧 − 𝑥 2 + 𝑦 2 = 1.34
𝑥𝑦 − 𝑧 2 = 0.09
𝑒 𝑥 − 𝑒 𝑦 + 𝑧 = 0.41
Assume any suitable initial guesses and perform two sets of iteration. Write numerical values with
an accuracy of up to three decimal places.