Engineering Mathematics 2: cram pack

Formula sheets, the key method steps for every topic and 16 practice problems, with full solutions at the back.

New problemsBack to the course

In the print window, choose “Save as PDF” to keep a copy on your device.

1Formulas and methods

Differential equations

Exponential growth/decay
dy/dt = ky → y = y0ekt
Half-life
k = ln 2 / t½
Newton's law of cooling
T(t) = Ta + (T0 − Ta)e−kt
Integrating factor
y′ + p(x)y = q(x): μ = e∫p dx, (μy)′ = μq
Mixing tank
dm/dt = r·c − (r/V)·m

Matrices & determinants

2×2 determinant
det = ad − bc
2×2 inverse
A⁻¹ = (1/det)·[d −b; −c a]
Product
(AB)ij = Σk AikBkj
Eigenvalues
det(A − λI) = 0
Cramer’s rule
xi = det(Ai)/det(A)

Vectors: dot & cross product

Magnitude
|a| = √(ax² + ay² + az²)
Unit vector
û = a/|a|
Components from an angle
Fx = F·cos θ, Fy = F·sin θ
Dot product
a·b = axbx + ayby + azbz = |a||b|cos θ
Scalar projection of a onto b
a·b/|b|
Work
W = F·d
Cross product
a × b = (aybz − azby, azbx − axbz, axby − aybx)
Areas
parallelogram |a × b|, triangle ½|a × b|
Moment about O
MO = r × F, r = A − O
Volume of a parallelepiped
V = |a·(b × c)|

Method

  1. Write every vector in components; for a force given by a magnitude and an angle, use F·cos θ and F·sin θ.
  2. For lengths and directions use the magnitude √(x² + y² + z²) and the unit vector a/|a|.
  3. Use the dot product for angles, projections, perpendicularity and work: a·b = |a||b|cos θ.
  4. Use the cross product for perpendicular vectors, areas and moments: write the three components carefully, keeping the order a × b.

Integrals

Power rule
∫xⁿ dx = xⁿ⁺¹/(n + 1) + C
Definite integral
∫ab f dx = F(b) − F(a)
Common integrals
∫eˣ dx = eˣ, ∫1/x dx = ln|x|, ∫sin x dx = −cos x
Area between curves
∫ab (f − g) dx

Method

  1. Rewrite the integrand as a sum of simple terms where possible.
  2. Apply the power rule: ∫xⁿ dx = xⁿ⁺¹/(n + 1), and standard integrals for ex, 1/x, sin and cos.
  3. Use substitution when you see a function and its derivative together.
  4. For a definite integral, evaluate F(b) − F(a).

Multivariable calculus

Partial derivatives
fx = ∂f/∂x (y constant), fxy = fyx
Gradient
∇f = (fx, fy), largest slope |∇f|
Directional derivative
Duf = ∇f·u, u = v/|v|
Tangent plane
z = f(a, b) + fx(x − a) + fy(y − b)
Chain rule
dz/dt = fx·dx/dt + fy·dy/dt
Second-derivative test
D = fxxfyy − fxy²: D > 0, fxx > 0 minimum; D > 0, fxx < 0 maximum; D < 0 saddle
Double integral
∬R f dA = ∫ab∫g(x)h(x) f dy dx
Polar coordinates
x = r cos θ, y = r sin θ, dA = r dr dθ

Method

  1. For a partial derivative, treat the other variable as a constant and use the one-variable rules.
  2. For gradient questions compute ∇f = (fx, fy) at the point; use |∇f| for the steepest slope and ∇f·u with a unit vector u for a directional derivative.
  3. For extrema solve fx = 0 and fy = 0 together, then classify every point with D = fxx·fyy − fxy².
  4. For a double integral sketch the region, choose the order (or polar coordinates for discs) so the inner limits are the bounding curves, and integrate inside first.

Sequences & series

Arithmetic sequence
aₙ = a₁ + (n − 1)d, Sₙ = n(a₁ + aₙ)/2
Geometric sequence
aₙ = a₁rⁿ⁻¹, Sₙ = a₁(1 − rⁿ)/(1 − r)
Infinite geometric series
S = a₁/(1 − r), |r| < 1
p-series
Σ1/kᵖ converges only for p > 1
Ratio test
L = lim |ak+1/ak|: L < 1 converges, L > 1 diverges
Radius of convergence
R = lim |ck/ck+1|
Taylor polynomial
Pₙ(x) = Σ f⁽ᵏ⁾(a)(x − a)ᵏ/k!
Maclaurin series
eˣ = Σxᵏ/k!, sin x = x − x³/3! + …, cos x = 1 − x²/2! + …, ln(1 + x) = x − x²/2 + …
Lagrange error bound
|Rₙ| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!

Method

  1. Identify the type: arithmetic (constant difference), geometric (constant ratio), a standard series, a power series or a Taylor polynomial.
  2. For sums use Sₙ = n(a₁ + aₙ)/2, Sₙ = a₁(1 − rⁿ)/(1 − r), or S = a₁/(1 − r) when |r| < 1, with a₁ the first term of the sum.
  3. For convergence: n-th term test first, then compare with Σ1/kᵖ or a geometric series, and use the ratio or root test for factorials and k-th powers.
  4. For Taylor polynomials list f(a), f′(a), f″(a), …, divide the k-th derivative by k!, and bound the error with M·|x − a|ⁿ⁺¹/(n + 1)!.

Convolution & LTI systems

Discrete convolution
y[n] = Σ x[k]·h[n − k]
Continuous convolution
y(t) = ∫ x(τ)·h(t − τ) dτ
Length and sum
Ly = Lx + Lh − 1, Σy = Σx·Σh
Start, end and area
t1 + t3 … t2 + t4, ∫y = ∫x·∫h
Step into e^(−at)u(t)
y(t) = (A/a)(1 − e−at)
First-order system
y′ + a·y = b·x → h(t) = b·e−atu(t)
First-order recursion
y[n] = a·y[n − 1] + b·x[n] → h[n] = b·aⁿ
Sinusoid through h(t)
|H(jω)|, ∠H(jω); H(jω) = b/(a + jω)
Fourier series amplitudes
square ±A: 4A/(kπ), odd k; triangle ±A: 8A/(π²k²), odd k; pulse train: A·D and (2A/(kπ))|sin kπD|

Method

  1. Sketch both signals and note where each one is nonzero.
  2. Flip h and slide it: for each output time, find the range where both signals overlap.
  3. Multiply and add (discrete) or multiply and integrate (continuous) over that overlap only.
  4. Check the result: the output starts at the sum of the start times, ends at the sum of the end times, and its area or sum is the product of the input areas or sums.

2Practice problems

1.Differential equationsMedium
t (min)T (°C)46Ta​187 °C120 °C

A forged steel part at 440 °C is left to cool in a workshop at 20 °C. After 46 min its temperature is 187 °C. It follows Newton's law of cooling, dT/dt = −k·(T − Tₐ). How long after the start, t, does it take to reach 120 °C?

T0 = 440 °CTa = 20 °CT(46) = 187 °CT = 120 °C
t = min
2.Differential equationsHard
t (min)m (kg)t = ?16 kg

A tank holds 200 L of water in which 10 kg of salt is dissolved. Brine with 0.1 kg of salt per litre flows in at 10 L/min, and the well-mixed solution flows out at the same rate. After how many minutes t does the tank contain 16 kg of salt?

V = 200 Lm0 = 10 kgc = 0.1 kg/Lr = 10 L/min
t = min
3.Matrices & determinantsMedium
A =−10−41−31440

Find the determinant det A of the 3×3 matrix A shown.

A = [−1 0 −4; 1 −3 1; 4 4 0]
det A =
4.Matrices & determinantsHard
A =11−35

Find the larger eigenvalue λ of the matrix A shown.

A = [1 1; −3 5]
λ =
5.Vectors: dot & cross productMedium
xyzab

Find the dot product a·b of a = (6, −3, −3) and b = (−6, 4, 0).

a = (6, −3, −3)b = (−6, 4, 0)
a·b =
6.Vectors: dot & cross productHard
xyzabc

Find the volume V of the parallelepiped whose edges from one corner are a = (2, −3, −1), b = (2, 2, −4) and c = (2, −4, 4).

a = (2, −3, −1)b = (2, 2, −4)c = (2, −4, 4)
V =
7.IntegralsMedium
xy11.5y = x ln x

Evaluate I = ∫₁1.5 x·ln x dx.

∫₁1.5 x ln x dx
I =
8.IntegralsHard
xy1y = 3xy = 3x²

The region between the line y = 3x and the parabola y = 3x² is rotated about the x-axis. Find the volume V of the solid.

y = 3xy = 3x²
V =
9.Multivariable calculusMedium
xy−1−1P

Use the tangent plane of f(x, y) = sin(x − 2y) at (−1, −1) to approximate f(−1.2, −0.8). Give the linear approximation L.

f(x, y) = sin(x − 2y)(a, b) = (−1, −1)(x, y) = (−1.2, −0.8)
L =
10.Multivariable calculusHard
xy−20202(2, 2)(−2, 2)

f(x, y) = −x³ + 12x − y² + 4y has a critical point at (−2, 2). Compute the discriminant D = fxxfyy − fxy² there.

f(x, y) = −x³ + 12x − y² + 4yP = (−2, 2)
D =
11.Sequences & seriesMedium
AΣ (k + 4)/k²k = 1 … ∞BΣ 3/k³k = 1 … ∞CΣ (−1)ᵏ/kk = 1 … ∞DΣ 1/(k² + 1)k = 1 … ∞

Exactly one of the four series diverges. Click it.

A: Σ (k + 4)/k²B: Σ 3/k³C: Σ (−1)ᵏ/kD: Σ 1/(k² + 1)
  1. 1A: Σ (k + 4)/k²
  2. 2B: Σ 3/k³
  3. 3C: Σ (−1)ᵏ/k
  4. 4D: Σ 1/(k² + 1)
12.Sequences & seriesHard
xyfP₄1

Use the Maclaurin polynomial of degree 4 of f(x) = cos(2x) to approximate f(1). Give P₄(1).

f(x) = cos(2x)n = 4x = 1
P₄(1) =
13.Laplace transformsMedium
ReIm−2−4−6−8−10−12−14−16−18

F(s) = (2s + 32)/(s² + 5s + 6). Its inverse transform has the form f(t) = A·e−2t + B·e−3t. Find the coefficient A by partial fractions.

F(s) = (2s + 32)/(s² + 5s + 6)
A =
14.Laplace transformsHard
t (s)y(t)t = 0.6y(0.6) = ?

Solve y″ + 3y′ + 2y = 0 with y(0) = 1 and y′(0) = 13 using the Laplace transform, and find y(0.6).

y″ + 3y′ + 2y = 0y(0) = 1y′(0) = 13t = 0.6 s
y(0.6) =
15.Convolution & LTI systemsMedium
x(t)t0.51h(t)t1

The pulse x(t) = 1 for 0 ≤ t < 0.5 s (zero elsewhere) is applied to a system with h(t) = e−1tu(t). Find the output y(2) at t = 2 s, after the pulse has ended.

x(t) = 1, 0 ≤ t < 0.5 sh(t) = e−1tu(t)t = 2 s
y(2) =
16.Convolution & LTI systemsHard
u[n]3+y[n]0.8·z⁻¹+

A system obeys y[n] = 0.8·y[n − 1] + 3·x[n] and starts at rest. The input is the unit step x[n] = u[n]. Find y[3].

y[n] = 0.8y[n − 1] + 3x[n]x[n] = u[n]
y[3] =

3Solutions

1.Differential equations

Answer: t = 71.58 min

First-order differential equations: dy/dt = k·y has the solution y = y₀·ekt (growth for k > 0, decay for k < 0); a half-life T gives k = ln 2/T. Newton's law of cooling, dT/dt = −k·(T − Tₐ), gives T(t) = Tₐ + (T₀ − Tₐ)·e−kt. A linear equation y′ + p(x)·y = q(x) is solved with the integrating factor μ = e∫p dx: (μy)′ = μq, so μy = ∫μq dx + C.

Cooling constant

k = ln(420/167)/46 = 0.02005 min⁻¹

Time to reach the target

t = ln(420/100)/0.02005 = 71.58 min

Answer

t = 71.58 min
2.Differential equations

Answer: t = 18.33 min

First-order differential equations: dy/dt = k·y has the solution y = y₀·ekt (growth for k > 0, decay for k < 0); a half-life T gives k = ln 2/T. Newton's law of cooling, dT/dt = −k·(T − Tₐ), gives T(t) = Tₐ + (T₀ − Tₐ)·e−kt. A linear equation y′ + p(x)·y = q(x) is solved with the integrating factor μ = e∫p dx: (μy)′ = μq, so μy = ∫μq dx + C.

Result

m(t) = 20 + (−10)·e−t/20

Solve for t

e−t/20 = (20 − 16)/(20 − 10) = 0.4
t = 20·ln(2.5) = 18.33 min

Answer

t = 18.33 min
3.Matrices & determinants

Answer: det A = −60

Matrices: the determinant of a 2×2 matrix is ad − bc. For 3×3, expand along a row: each entry times the determinant of the 2×2 matrix left after removing its row and column, with alternating signs + − +. Eigenvalues λ solve det(A − λI) = 0; for 2×2 that is λ² − (trace)·λ + det = 0.

2×2 minors of the first row

M₁₁ = −4, M₁₂ = −4, M₁₃ = 16

Cofactor expansion

det A = (−1)·(−4) − 0·(−4) + (−4)·16 = −60

Answer

det A = −60
4.Matrices & determinants

Answer: λ = 4

Matrices: the determinant of a 2×2 matrix is ad − bc. For 3×3, expand along a row: each entry times the determinant of the 2×2 matrix left after removing its row and column, with alternating signs + − +. Eigenvalues λ solve det(A − λI) = 0; for 2×2 that is λ² − (trace)·λ + det = 0.

Characteristic equation

λ² − (6)λ + (8) = 0

Solve

λ = (6 ± √(4))/2 → λ = 2, 4
λ = 4

Answer

λ = 4
5.Vectors: dot & cross product

Answer: a·b = −48

Dot product: a·b = axbx + ayby + azbz = |a|·|b|·cos θ, a scalar. It gives the angle between two vectors, cos θ = a·b/(|a|·|b|), and it is zero when they are perpendicular. The scalar projection of a onto b is a·b/|b|, and a constant force F moving through a displacement d does work W = F·d.

Multiply matching components and add

a·b = 6·(−6) + (−3)·4 + (−3)·0 = (−36) + (−12) + 0 = −48

Answer

a·b = −48

Negative: the angle between the vectors is more than 90°.

6.Vectors: dot & cross product

Answer: V = 44

Cross product: a × b = (aybz − azby, azbx − axbz, axby − aybx) is perpendicular to both vectors, and |a × b| = |a|·|b|·sin θ is the area of the parallelogram they span (the triangle has half of it). Order matters: b × a = −(a × b). The moment of a force F applied at point A about point O is MO = r × F, with r = A − O.

Cross product b × c

(b × c)x = bycz − bzcy = 2·4 − (−4)·(−4) = −8
(b × c)y = bzcx − bxcz = (−4)·2 − 2·4 = −16
(b × c)z = bxcy − bycx = 2·(−4) − 2·2 = −12

Dot product with a

a·(b × c) = 2·(−8) + (−3)·(−16) + (−1)·(−12) = 44
V = |44| = 44

Answer

V = 44
7.Integrals

Answer: I = 0.1436

Integration by parts: ∫u·dv = u·v − ∫v·du, from the product rule. Choose u as the factor that gets simpler when you differentiate it (a power of x, or ln x), and dv as the factor you can integrate (eˣ, sin x, cos x, or x when u = ln x). For a definite integral, evaluate u·v at both limits: ∫ₐᵇ u dv = [u·v]ₐᵇ − ∫ₐᵇ v du. A factor x² needs two rounds.

Choose u and dv

u = ln x, du = dx/x; dv = x dx, v = x²/2

Antiderivative

∫x ln x dx = (x²/2)ln x − ∫(x/2) dx = (x²/2)ln x − x²/4

Limits

(1.125·ln 1.5 − 0.5625) − (0 − 0.25) = 0.1436

Answer

I = 0.1436
8.Integrals

Answer: V = 3.77

Volumes of revolution: rotating the region under y = f(x) from x = a to x = b about the x-axis makes a solid whose cross-section at x is a disc of radius f(x). Adding the discs gives V = π∫ₐᵇ f(x)² dx. When the region lies between an outer curve R(x) and an inner curve r(x), each slice is a washer: V = π∫ₐᵇ (R(x)² − r(x)²) dx. Square the radii before you subtract them.

Intersection points

3x = 3x² → x = 0, x = 1

Washer method

V = π∫₀1 ((3x)² − (3x²)²) dx

Integrate

V = π(9·1³/3 − 9·1⁵/5) = π·1.2 = 3.77

Answer

V = 3.77
9.Multivariable calculus

Answer: L = 0.5173

Gradient: ∇f = (fx, fy) points in the direction of steepest ascent, and |∇f| is the largest rate of change. The directional derivative along a direction v is Duf = ∇f·u with the unit vector u = v/|v|. The tangent plane at (a, b) is z = f(a, b) + fx(a, b)(x − a) + fy(a, b)(y − b); it gives the linear approximation of f near that point. Chain rule: dz/dt = fx·dx/dt + fy·dy/dt.

Value and slopes at the point

f(−1, −1) = 0.8415
fx = 1 cos(x − 2y) → fx(−1, −1) = 0.5403
fy = −2 cos(x − 2y) → fy(−1, −1) = −1.081

Tangent plane

z = 0.8415 + 0.5403(x + 1) − 1.081(y + 1)
z = 0.5403x − 1.081y + 0.3012

Linear approximation

L = 0.8415 + 0.5403·(−0.2) + (−1.081)·0.2 = 0.5173

Answer

L = 0.5173

The exact value is f = 0.389418; the plane is close because the step is small.

10.Multivariable calculus

Answer: D = −24

Critical points and double integrals: a critical point has fx = fy = 0. Second-derivative test with D = fxxfyy − fxy²: D > 0 and fxx > 0 is a local minimum, D > 0 and fxx < 0 a local maximum, D < 0 a saddle point, and D = 0 decides nothing. For a double integral over a general region, let the inner limits be the curves that bound the region; for discs and sectors use polar coordinates, x = r cos θ, y = r sin θ and dA = r dr dθ (do not forget the factor r).

Second derivatives

fxx = 12, fyy = −2, fxy = 0

Discriminant

D = 12·(−2) − 0² = −24
→ saddle point

Answer

D = −24
11.Sequences & series

Answer: A: Σ (k + 4)/k²

Convergence: if the terms do not tend to 0 the series diverges (n-th term test). Σ1/kᵖ converges only for p > 1, and Σrᵏ only for |r| < 1. Compare a series with one of these when its terms behave like them for large k. Ratio test: L = lim |ak+1/ak|; root test: L = lim |ak|1/k. L < 1 converges, L > 1 diverges, L = 1 decides nothing. A power series Σck(x − x₀)ᵏ converges for |x − x₀| < R with R = lim |ck/ck+1|; an alternating series with terms decreasing to 0 converges.

Test each series

A: Σ (k + 4)/k² → Behaves like 1/k for large k: diverges by comparison with the harmonic series.
B: Σ 3/k³ → p-series with p = 3 > 1: converges.
C: Σ (−1)ᵏ/k → Alternating, with 1/k decreasing to 0: converges (conditionally).
D: Σ 1/(k² + 1) → Smaller than 1/k², a convergent p-series: converges by comparison.

The answer

A: Σ (k + 4)/k²

Answer

A: Σ (k + 4)/k²
12.Sequences & series

Answer: P₄(1) = −0.3333

Taylor polynomials: near x = a a smooth function is approximated by Pₙ(x) = f(a) + f′(a)(x − a) + f″(a)(x − a)²/2! + … + f⁽ⁿ⁾(a)(x − a)ⁿ/n!; for a = 0 this is the Maclaurin polynomial. Standard series: eˣ = Σxᵏ/k!, sin x = x − x³/3! + x⁵/5! − …, cos x = 1 − x²/2! + x⁴/4! − …, ln(1 + x) = x − x²/2 + x³/3 − …, 1/(1 − x) = Σxᵏ. Lagrange error bound: |f(x) − Pₙ(x)| ≤ M·|x − a|ⁿ⁺¹/(n + 1)!, where M is the largest value of |f⁽ⁿ⁺¹⁾| between a and x.

Maclaurin polynomial

P₄(x) = 1 − 2x² + 0.6667x⁴

Substitute

P₄(1) = 1 − 2 + 0.666667 = −0.3333

Answer

P₄(1) = −0.3333

The exact value is f(1) = −0.416147; the error −0.0828135 shrinks when you take more terms or a point closer to 0.

13.Laplace transforms

Answer: A = 28

Laplace transform: F(s) = ∫₀∞ f(t)·e−st dt turns derivatives into multiplications by s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Standard pairs: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), and a shift s → s + σ multiplies the time function by e−σt. To invert a rational F(s), split it into partial fractions: the residue of a simple pole p is A = [(s − p)·F(s)] at s = p. Final value theorem: f(∞) = lims→0 s·F(s), valid only when all poles of s·F(s) lie in the left half-plane. Initial value theorem: f(0⁺) = lims→∞ s·F(s).

Factor the denominator

s² + 5s + 6 = (s + 2)(s + 3)
F(s) = A/(s + 2) + B/(s + 3)

Cover-up rule

Multiply by (s + 2) and set s = −2: the B-term vanishes.

A = (2s + 32)/(s + 3) |s = −2 = 28/1 = 28

Answer

A = 28
(B = −26, so f(t) = 28·e−2t − 26·e−3t)
14.Laplace transforms

Answer: y(0.6) = 4.015

Laplace transform: F(s) = ∫₀∞ f(t)·e−st dt turns derivatives into multiplications by s: ℒ{f′} = sF(s) − f(0), ℒ{f″} = s²F(s) − s·f(0) − f′(0). Standard pairs: 1 ↔ 1/s, e−at ↔ 1/(s + a), t·e−at ↔ 1/(s + a)², sin ωt ↔ ω/(s² + ω²), cos ωt ↔ s/(s² + ω²), and a shift s → s + σ multiplies the time function by e−σt. To invert a rational F(s), split it into partial fractions: the residue of a simple pole p is A = [(s − p)·F(s)] at s = p. Final value theorem: f(∞) = lims→0 s·F(s), valid only when all poles of s·F(s) lie in the left half-plane. Initial value theorem: f(0⁺) = lims→∞ s·F(s).

Transform the equation

s²Y − 1s − 13 + 3(sY − 1) + 2Y = 0
Y(s) = (s + 16)/((s + 1)(s + 2))

Partial fractions

A = 15/1 = 15
B = 14/(−1) = −14

Result

y(t) = 15·e−t − 14·e−2t
y(0.6) = 4.015

Answer

y(0.6) = 4.015

Check: y(0) = A + B = 1.

15.Convolution & LTI systems

Answer: y(2) = 0.08779

Continuous convolution: y(t) = ∫ x(τ)·h(t − τ) dτ. Flip h, slide it to position t and integrate the product of the overlapping parts. For causal signals the integral runs from 0 to t. Useful results: the area of y is the area of x times the area of h; if x lasts from t₁ to t₂ and h from t₃ to t₄, y lasts from t₁ + t₃ to t₂ + t₄. For h(t) = e−atu(t) and a step of height A, y(t) = (A/a)(1 − e−at).

Output at the end of the pulse

y(0.5) = (1/1)(1 − e−0.5) = 0.3935

Decay after the pulse

y(2) = 0.3935·e−1·1.5 = 0.08779

Answer

y(2) = 0.08779
16.Convolution & LTI systems

Answer: y[3] = 8.856

Impulse response and Fourier series: the impulse response h(t) describes an LTI system completely: for y′ + a·y = b·x it is h(t) = b·e−atu(t), and for y[n] = a·y[n − 1] + b·x[n] it is h[n] = b·aⁿu[n]. A sinusoid cos ωt passes through with gain |H(jω)| and phase ∠H(jω), where H is the Fourier transform of h. A periodic signal is a sum of harmonics at k/T: a square wave of ±A has only odd harmonics with amplitudes 4A/(kπ), a triangle wave of peak A has 8A/(π²k²), and a pulse train of height A and duty cycle D has the mean A·D and the harmonics (2A/(kπ))·|sin(kπD)|.

Iterate

y[0] = 3
y[1] = 5.4
y[2] = 7.32
y[3] = 8.856

Result

y[3] = 3·(1 − (0.8)4)/(1 − 0.8) = 8.856

Answer

y[3] = 8.856

ExamProblems · examproblems.com/courses/engineering-mathematics-2