GTU Mathematics 2 (3110015) Semester Semester 1 & 2 Summer 2025 Solved Paper
Complete step-by-step worked solutions for GTU Mathematics - 2 (3110015) Summer 2025 examination paper covering vector calculus, differential equations, series solutions, Fourier integrals, and Laplace transforms.
Introduction
This document contains fully worked, student-friendly solutions for the Gujarat Technological University (GTU) Mathematics - 2 (3110015) examination of Summer 2025. It provides detailed derivations, systematic steps for differential equations, clear integration limits for vector theorems, and rigorous methods for Laplace transforms and series solutions.
Paper Information
- University
- Gujarat Technological University
- Department
- First Year Engineering (All Branches)
- Semester
- Semester 1 & 2
- Subject
- Mathematics - 2
- Subject Code
- 3110015
- Exam
- Summer 2025
- Year
- 2025
Questions and Solutions
Q.1 (a)
Question
Show that F = (y² - z² + 3yz - 2x)i + (3xz + 2xy)j + (3xy - 2xz + 2z)k is both solenoidal and irrotational. [03 Marks]
Solution
Definition
- A vector field F is solenoidal if its divergence is zero:
div F = ∇ · F = 0. - A vector field F is irrotational if its curl is zero:
curl F = ∇ × F = 0.
---
Step 1: Check Solenoidal (∇ · F)
Given:
- F₁ = y² - z² + 3yz - 2x
- F₂ = 3xz + 2xy
- F₃ = 3xy - 2xz + 2z
Compute the partial derivatives:
- ∂F₁/∂x = ∂/∂x (y² - z² + 3yz - 2x) = -2
- ∂F₂/∂y = ∂/∂y (3xz + 2xy) = 2x
- ∂F₃/∂z = ∂/∂z (3xy - 2xz + 2z) = -2x + 2
Add them together: ∇ · F = ∂F₁/∂x + ∂F₂/∂y + ∂F₃/∂z ∇ · F = -2 + 2x + (-2x + 2) = 0
Since ∇ · F = 0, the vector field F is solenoidal.
---
Step 2: Check Irrotational (∇ × F)
Compute the curl using the determinant:
∇ × F = | i j k |
| ∂/∂x ∂/∂y ∂/∂z |
| F₁ F₂ F₃ |
Calculate each component:
- i-component:
∂F₃/∂y - ∂F₂/∂z = ∂/∂y (3xy - 2xz + 2z) - ∂/∂z (3xz + 2xy)
= 3x - 3x = 0
- j-component:
∂F₁/∂z - ∂F₃/∂x = ∂/∂z (y² - z² + 3yz - 2x) - ∂/∂x (3xy - 2xz + 2z)
= (-2z + 3y) - (3y - 2z) = -2z + 3y - 3y + 2z = 0
- k-component:
∂F₂/∂x - ∂F₁/∂y = ∂/∂x (3xz + 2xy) - ∂/∂y (y² - z² + 3yz - 2x)
= (3z + 2y) - (2y + 3z) = 3z + 2y - 2y - 3z = 0
Therefore: ∇ × F = 0i + 0j + 0k = 0
Since ∇ × F = 0, the vector field F is irrotational.
Conclusion: The vector field F is both solenoidal and irrotational.
Q.1 (b)
Question
Find Fourier integral representation of function f(x) = {1, |x| < 1; 0, |x| > 1}. Hence, Evaluate ∫ [sin(ω) cos(ωx) / ω] dω. [04 Marks]
Solution
Step 1: Identify Nature of Function
The given function is:
- f(x) = 1 for -1 < x < 1
- f(x) = 0 for |x| > 1
Notice that f(-x) = f(x), which means f(x) is an even function.
---
Step 2: Fourier Cosine Integral Formulation
For an even function, the Fourier sine component B(ω) = 0, and the Fourier integral representation simplifies to the Fourier Cosine Integral: f(x) = (2/π) ∫₀^∞ A(ω) cos(ωx) dω where: A(ω) = ∫₀^∞ f(t) cos(ωt) dt
Substitute f(t) = 1 for 0 ≤ t < 1 and f(t) = 0 for t > 1: A(ω) = ∫₀¹ 1 · cos(ωt) dt = [sin(ωt) / ω]₀¹ = sin(ω) / ω
Therefore, the Fourier integral representation is: f(x) = (2/π) ∫₀^∞ [sin(ω) / ω] cos(ωx) dω
---
Step 3: Evaluate the Integral
Rearranging the equation: ∫₀^∞ [sin(ω) cos(ωx) / ω] dω = (π/2) f(x)
Now consider the values based on x:
- For |x| < 1:
f(x) = 1 => ∫₀^∞ [sin(ω) cos(ωx) / ω] dω = π/2
- For |x| > 1:
f(x) = 0 => ∫₀^∞ [sin(ω) cos(ωx) / ω] dω = 0
- At the points of discontinuity x = ±1:
By Dirichlet's condition, the integral equals (π/2) · [f(1⁺) + f(1⁻)] / 2 = (π/2) · (0 + 1)/2 = π/4.
Final Answer: ∫₀^∞ [sin(ω) cos(ωx) / ω] dω =
- π/2, for |x| < 1
- π/4, for x = ±1
- 0, for |x| > 1
Q.1 (c)
Question
Verify Green’s theorem in the plane for ∮_C [(3x² - 8y²) dx + (4y - 6xy) dy], where C is the boundary of the region bounded by x = 0, y = 0, x + y = 1. [07 Marks]
Solution
Statement of Green's Theorem
Green's theorem relates a line integral around a simple closed curve C to a double integral over the plane region R bounded by C: ∮C (M dx + N dy) = ∬R (∂N/∂x - ∂M/∂y) dx dy
Here:
- M = 3x² - 8y²
- N = 4y - 6xy
---
Part 1: Double Integral (Right-Hand Side)
Compute the partial derivatives:
- ∂M/∂y = ∂/∂y (3x² - 8y²) = -16y
- ∂N/∂x = ∂/∂x (4y - 6xy) = -6y
So: ∂N/∂x - ∂M/∂y = -6y - (-16y) = 10y
The region R is the triangle bounded by x = 0, y = 0, and y = 1 - x (from x = 0 to x = 1): ∬_R 10y dx dy = ∫₀¹ [ ∫₀^(1-x) 10y dy ] dx
Evaluate the inner integral: ∫₀^(1-x) 10y dy = [5y²]₀^(1-x) = 5(1 - x)² = 5(1 - 2x + x²)
Evaluate the outer integral: ∫₀¹ 5(1 - 2x + x²) dx = 5 [x - x² + x³/3]₀¹ = 5 [1 - 1 + 1/3] = 5/3
RHS = 5/3
---
Part 2: Line Integral (Left-Hand Side)
The boundary curve C consists of three straight line segments oriented counter-clockwise:
- C₁: Along y = 0 from (0, 0) to (1, 0)
- y = 0, dy = 0
- x goes from 0 to 1
- ∫_C₁ = ∫₀¹ (3x² - 0) dx + 0 = [x³]₀¹ = 1
- C₂: Along x + y = 1 from (1, 0) to (0, 1)
- y = 1 - x => dy = -dx
- x goes from 1 to 0
- M = 3x² - 8(1 - x)² = 3x² - 8(1 - 2x + x²) = -5x² + 16x - 8
- N = 4(1 - x) - 6x(1 - x) = 4 - 4x - 6x + 6x² = 6x² - 10x + 4
- Integrand: M dx + N dy = (-5x² + 16x - 8) dx + (6x² - 10x + 4)(-dx)
= (-5x² + 16x - 8 - 6x² + 10x - 4) dx
= (-11x² + 26x - 12) dx - ∫_C₂ = ∫₁⁰ (-11x² + 26x - 12) dx
= [-11x³/3 + 13x² - 12x]₁⁰
= 0 - [-11/3 + 13 - 12]
= - [-11/3 + 1] = - [-8/3] = 8/3
- C₃: Along x = 0 from (0, 1) to (0, 0)
- x = 0, dx = 0
- y goes from 1 to 0
- ∫_C₃ = ∫₁⁰ (4y - 0) dy = [2y²]₁⁰ = 0 - 2 = -2
---
Sum of Line Integrals:
∮C (M dx + N dy) = ∫C₁ + ∫C₂ + ∫C₃ = 1 + 8/3 - 2 = -1 + 8/3 = 5/3
LHS = 5/3
Conclusion
Since LHS = RHS = 5/3, Green's theorem is verified.
Q.2 (a)
Question
Solve y e^x dx + (2y + e^x) dy = 0. [03 Marks]
Solution
Step 1: Test for Exactness
The equation is of the form M dx + N dy = 0, where:
- M = y e^x
- N = 2y + e^x
Find partial derivatives:
- ∂M/∂y = e^x
- ∂N/∂x = e^x
Since ∂M/∂y = ∂N/∂x = e^x, the differential equation is exact.
---
Step 2: Integrate to Find the General Solution
The solution is given by: ∫ M dx (treating y as constant) + ∫ (terms of N independent of x) dy = C
- Term 1: ∫ y e^x dx = y e^x
- Term 2: Terms of N not containing x is
2y. So: ∫ 2y dy = y²
Therefore, the general solution is: y e^x + y² = C where C is an arbitrary constant.
Q.2 (b)
Question
Find the directional Derivative of ∅ = xy² + yz² at the point (2, -1, 1) in the direction of the vector i + 2j + 2k. [04 Marks]
Solution
Step 1: Find the Gradient ∇∅
The gradient is defined as: ∇∅ = (∂∅/∂x) i + (∂∅/∂y) j + (∂∅/∂z) k
Given ∅ = xy² + yz²:
- ∂∅/∂x = y²
- ∂∅/∂y = 2xy + z²
- ∂∅/∂z = 2yz
At the point P(2, -1, 1):
- ∂∅/∂x = (-1)² = 1
- ∂∅/∂y = 2(2)(-1) + (1)² = -4 + 1 = -3
- ∂∅/∂z = 2(-1)(1) = -2
Thus: ∇∅ at (2, -1, 1) = 1i - 3j - 2k
---
Step 2: Find the Unit Vector along the Given Direction
Given vector: a = i + 2j + 2k Magnitude of a: |a| = √(1² + 2² + 2²) = √(1 + 4 + 4) = √9 = 3
Unit vector û: û = a / |a| = (i + 2j + 2k) / 3
---
Step 3: Compute the Directional Derivative
The directional derivative is the dot product of ∇∅ and û: Directional Derivative = ∇∅ · û = (1i - 3j - 2k) · (1/3)(i + 2j + 2k) = (1/3) [ (1)(1) + (-3)(2) + (-2)(2) ] = (1/3) [ 1 - 6 - 4 ] = -9 / 3 = -3
Final Answer: The directional derivative is -3.
Q.2 (c)
Question
Find the power series solution of d²y/dx² + xy = 0. [07 Marks]
Solution
Step 1: Nature of the Ordinary Point
The equation is y'' + xy = 0. Here, x = 0 is an ordinary point because the coefficient of y'' is 1 (non-zero everywhere).
Let the power series solution about x = 0 be: y = ∑(m=0)^∞ am x^m = a₀ + a₁x + a₂x² + a₃x³ + ...
Differentiating term by term: y' = ∑(m=1)^∞ m am x^(m-1) y'' = ∑(m=2)^∞ m(m-1) am x^(m-2)
---
Step 2: Substitute into Differential Equation
y'' + xy = 0 ∑(m=2)^∞ m(m-1) am x^(m-2) + x ∑(m=0)^∞ am x^m = 0 ∑(m=2)^∞ m(m-1) am x^(m-2) + ∑(m=0)^∞ am x^(m+1) = 0
Write out initial terms to equalize powers of x:
- For the first sum, at m = 2: 2(1) a₂ = 2a₂
- Shift index of first sum by letting m-2 = k => m = k+2:
∑(k=0)^∞ (k+2)(k+1) a(k+2) x^k
= 2a₂ + ∑(k=1)^∞ (k+2)(k+1) a(k+2) x^k - For the second sum, shift index by letting m+1 = k => m = k-1:
∑(k=1)^∞ a(k-1) x^k
Combining: 2a₂ + ∑(k=1)^∞ [ (k+2)(k+1) a(k+2) + a_(k-1) ] x^k = 0
---
Step 3: Equating Coefficients to Zero
- Constant term (x⁰): 2a₂ = 0 => a₂ = 0
- Recurrence relation for k ≥ 1:
a(k+2) = -a(k-1) / [(k+2)(k+1)]
---
Step 4: Calculate Coefficients in Terms of a₀ and a₁
- For k = 1:
a₃ = -a₀ / (3 · 2) = -a₀ / 6 - For k = 2:
a₄ = -a₁ / (4 · 3) = -a₁ / 12 - For k = 3:
a₅ = -a₂ / (5 · 4) = 0 (since a₂ = 0) - For k = 4:
a₆ = -a₃ / (6 · 5) = -(-a₀ / 6) / 30 = a₀ / 180 - For k = 5:
a₇ = -a₄ / (7 · 6) = -(-a₁ / 12) / 42 = a₁ / 504
---
Step 5: Write the General Solution
Grouping terms with a₀ and a₁: y = a₀ [ 1 - x³/6 + x⁶/180 - ... ] + a₁ [ x - x⁴/12 + x⁷/504 - ... ] where a₀ and a₁ are arbitrary constants.
Q.2 (c) OR
Question
Find the series solution of (1 - x²)y'' - 2xy' + 2y = 0. [07 Marks]
Solution
Step 1: Standard Legendre Form
Notice that this is Legendre's differential equation: (1 - x²)y'' - 2xy' + n(n+1)y = 0 with n = 1, since n(n+1) = 1(2) = 2. Here, x = 0 is an ordinary point.
Let the series solution be: y = ∑(m=0)^∞ am x^m y' = ∑(m=1)^∞ m am x^(m-1) y'' = ∑(m=2)^∞ m(m-1) am x^(m-2)
---
Step 2: Substitute into the Equation
(1 - x²) ∑(m=2)^∞ m(m-1) am x^(m-2) - 2x ∑(m=1)^∞ m am x^(m-1) + 2 ∑(m=0)^∞ am x^m = 0
Expand the terms: ∑(m=2)^∞ m(m-1) am x^(m-2) - ∑(m=2)^∞ m(m-1) am x^m - ∑(m=1)^∞ 2m am x^m + ∑(m=0)^∞ 2 am x^m = 0
Shift index in the first summation by setting m - 2 = k => m = k + 2: ∑(k=0)^∞ (k+2)(k+1) a(k+2) x^k - ∑(k=2)^∞ k(k-1) ak x^k - ∑(k=1)^∞ 2k ak x^k + ∑(k=0)^∞ 2 ak x^k = 0
---
Step 3: Equate Coefficients of Powers of x to Zero
- For x⁰ (k = 0):
(2)(1) a₂ + 2 a₀ = 0 => 2a₂ + 2a₀ = 0 => a₂ = -a₀
- For x¹ (k = 1):
(3)(2) a₃ - 2(1) a₁ + 2 a₁ = 0 => 6a₃ = 0 => a₃ = 0
- For x^k (k ≥ 2):
(k+2)(k+1) a(k+2) - [k(k-1) + 2k - 2] ak = 0
k(k-1) + 2k - 2 = k² - k + 2k - 2 = k² + k - 2 = (k+2)(k-1)
So: (k+2)(k+1) a(k+2) = (k+2)(k-1) ak a(k+2) = [(k - 1) / (k + 1)] ak
---
Step 4: Calculate Subsequent Coefficients
- For k = 2:
a₄ = (1/3) a₂ = (1/3)(-a₀) = -a₀/3 - For k = 3:
a₅ = (2/4) a₃ = 0 (since a₃ = 0) - For k = 4:
a₆ = (3/5) a₄ = (3/5)(-a₀/3) = -a₀/5 - For all odd k ≥ 3:
a_k = 0 because a₃ = 0.
---
Step 5: Final Series Solution
y = a₀ [ 1 - x² - x⁴/3 - x⁶/5 - ... ] + a₁ x where a₀ and a₁ are arbitrary constants. Notice that y₁ = x is a polynomial solution (the Legendre polynomial P₁(x) = x).
Q.3 (a)
Question
Solve (4D² - 4D + 1)y = e^(x/2). [03 Marks]
Solution
Step 1: Complementary Function (y_c)
Auxiliary Equation (A.E.): 4m² - 4m + 1 = 0 (2m - 1)² = 0 => m = 1/2, 1/2 (real and repeated roots)
Therefore, the Complementary Function is: y_c = (c₁ + c₂ x) e^(x/2)
---
Step 2: Particular Integral (y_p)
y_p = [1 / (4D² - 4D + 1)] e^(x/2) = [1 / (2D - 1)²] e^(x/2)
Since replacing D by 1/2 makes the denominator zero (case of failure with multiplicity 2): Using the standard formula: [1 / (D - a)^n] e^(ax) = (x^n / n!) e^(ax)
Rewrite: [1 / 4(D - 1/2)²] e^(x/2) = (1/4) · (x² / 2!) e^(x/2) = (x² / 8) e^(x/2)
---
Step 3: General Solution
y = yc + yp y = (c₁ + c₂ x) e^(x/2) + (x² / 8) e^(x/2)
Q.3 (b)
Question
Solve (D² + 6D + 8)y = cos 2x. [04 Marks]
Solution
Step 1: Complementary Function (y_c)
Auxiliary Equation: m² + 6m + 8 = 0 (m + 2)(m + 4) = 0 => m = -2, -4
y_c = c₁ e^(-2x) + c₂ e^(-4x)
---
Step 2: Particular Integral (y_p)
y_p = [1 / (D² + 6D + 8)] cos 2x
Substitute D² = -2² = -4: y_p = [1 / (-4 + 6D + 8)] cos 2x = [1 / (6D + 4)] cos 2x = [1 / (2(3D + 2))] cos 2x
Multiply numerator and denominator by conjugate (3D - 2): y_p = (3D - 2) / [2(9D² - 4)] cos 2x
Again substitute D² = -4: Denominator = 2 [ 9(-4) - 4 ] = 2 [ -36 - 4 ] = 2(-40) = -80
Numerator: (3D - 2) cos 2x = 3 D(cos 2x) - 2 cos 2x = 3(-2 sin 2x) - 2 cos 2x = -6 sin 2x - 2 cos 2x
So: y_p = (-6 sin 2x - 2 cos 2x) / (-80) = (3 sin 2x + cos 2x) / 40
---
Step 3: General Solution
y = yc + yp y = c₁ e^(-2x) + c₂ e^(-4x) + (1/40) (3 sin 2x + cos 2x)
Q.3 (c)
Question
Solve (D² - 2D + 3)y = x³ + sin x by using the method of undetermined coefficients. [07 Marks]
Solution
Step 1: Complementary Function (y_c)
Auxiliary equation: m² - 2m + 3 = 0 m = [2 ± √(4 - 12)] / 2 = [2 ± √(-8)] / 2 = 1 ± i√2
y_c = e^x [ c₁ cos(√2 x) + c₂ sin(√2 x) ]
---
Step 2: Form of the Particular Solution (y_p)
The right side is R(x) = x³ + sin x. Since no term in R(x) appears in yc, the trial solution is: yp = yp1 + yp2 where:
- y_p1 = A x³ + B x² + C x + D
- y_p2 = E cos x + F sin x
---
Step 3: Determine y_p1 for R₁(x) = x³
yp1 = A x³ + B x² + C x + D yp1' = 3A x² + 2B x + C y_p1'' = 6A x + 2B
Substitute into (D² - 2D + 3)y = x³: (6A x + 2B) - 2(3A x² + 2B x + C) + 3(A x³ + B x² + C x + D) = x³ 3A x³ + (3B - 6A) x² + (3C - 4B + 6A) x + (3D - 2C + 2B) = x³
Equating coefficients:
- x³: 3A = 1 => A = 1/3
- x²: 3B - 6A = 0 => 3B = 6(1/3) = 2 => B = 2/3
- x¹: 3C - 4B + 6A = 0 => 3C = 4(2/3) - 6(1/3) = 8/3 - 2 = 2/3 => C = 2/9
- x⁰: 3D - 2C + 2B = 0 => 3D = 2(2/9) - 2(2/3) = 4/9 - 4/3 = -8/9 => D = -8/27
So: y_p1 = (1/3)x³ + (2/3)x² + (2/9)x - 8/27
---
Step 4: Determine y_p2 for R₂(x) = sin x
yp2 = E cos x + F sin x yp2' = -E sin x + F cos x y_p2'' = -E cos x - F sin x
Substitute into (D² - 2D + 3)y = sin x: (-E cos x - F sin x) - 2(-E sin x + F cos x) + 3(E cos x + F sin x) = sin x (-E - 2F + 3E) cos x + (-F + 2E + 3F) sin x = sin x (2E - 2F) cos x + (2E + 2F) sin x = sin x
Equating coefficients:
- cos x: 2E - 2F = 0 => E = F
- sin x: 2E + 2F = 1 => 4E = 1 => E = 1/4, F = 1/4
So: y_p2 = (1/4) cos x + (1/4) sin x
---
Step 5: General Solution
y = yc + yp1 + y_p2 y = e^x [ c₁ cos(√2 x) + c₂ sin(√2 x) ] + (1/3)x³ + (2/3)x² + (2/9)x - 8/27 + (1/4)(cos x + sin x)
Q.3 (a) OR
Question
Solve (D² + a²)y = cos ax. [03 Marks]
Solution
Step 1: Complementary Function (y_c)
Auxiliary Equation: m² + a² = 0 => m = ±ia y_c = c₁ cos(ax) + c₂ sin(ax)
---
Step 2: Particular Integral (y_p)
y_p = [1 / (D² + a²)] cos ax
Since replacing D² with -a² gives zero in the denominator, this is a case of failure. We use the standard formula: [1 / (D² + a²)] cos ax = (x / 2a) sin ax
Derivation: [1 / (D² + a²)] cos ax = Real Part of [1 / (D² + a²)] e^(iax) = Real Part of e^(iax) [1 / ((D + ia)² + a²)] 1 = Real Part of e^(iax) [1 / (D² + 2iaD)] 1 = Real Part of e^(iax) (1 / 2ia) [1 / D] 1 = Real Part of e^(iax) [x / (2ia)] = Real Part of (-i x / 2a) (cos ax + i sin ax) = (x / 2a) sin ax
---
Step 3: General Solution
y = yc + yp y = c₁ cos(ax) + c₂ sin(ax) + (x / 2a) sin ax
Q.3 (b) OR
Question
Solve d²y/dx² + 2 dy/dx + 2y = e^x sin x + 7. [04 Marks]
Solution
Step 1: Complementary Function (y_c)
Auxiliary Equation: m² + 2m + 2 = 0 m = [-2 ± √(4 - 8)] / 2 = -1 ± i
y_c = e^(-x) [ c₁ cos x + c₂ sin x ]
---
Step 2: Particular Integral y_p1 (for e^x sin x)
y_p1 = [1 / (D² + 2D + 2)] (e^x sin x) Using the exponential shift property: [1 / f(D)] (e^(ax) V) = e^(ax) [1 / f(D + a)] V
Replace D with D + 1: y_p1 = e^x [1 / ((D + 1)² + 2(D + 1) + 2)] sin x = e^x [1 / (D² + 2D + 1 + 2D + 2 + 2)] sin x = e^x [1 / (D² + 4D + 5)] sin x
Substitute D² = -1² = -1: y_p1 = e^x [1 / (-1 + 4D + 5)] sin x = e^x [1 / (4D + 4)] sin x = (e^x / 4) [1 / (D + 1)] sin x
Multiply numerator and denominator by (D - 1): y_p1 = (e^x / 4) [(D - 1) / (D² - 1)] sin x = (e^x / 4) [(D - 1) / (-1 - 1)] sin x = (e^x / 4) [(D - 1) / (-2)] sin x = - (e^x / 8) [D(sin x) - sin x] = - (e^x / 8) [cos x - sin x] = (e^x / 8) [sin x - cos x]
---
Step 3: Particular Integral y_p2 (for 7)
y_p2 = [1 / (D² + 2D + 2)] (7 e^(0x)) = 7 / (0 + 0 + 2) = 7/2
---
Step 4: General Solution
y = yc + yp1 + y_p2 y = e^(-x) [ c₁ cos x + c₂ sin x ] + (e^x / 8) (sin x - cos x) + 7/2
Q.3 (c) OR
Question
Solve x² d²y/dx² + x dy/dx + y = log x · sin(log x) by Cauchy’s linear equation. [07 Marks]
Solution
Step 1: Transform Cauchy-Euler Equation into Constant Coefficients
Let x = e^z, so that z = log x. Define θ = d/dz. Then:
- x (dy/dx) = θ y
- x² (d²y/dx²) = θ(θ - 1) y
Substitute into the differential equation: [ θ(θ - 1) + θ + 1 ] y = z sin z [ θ² - θ + θ + 1 ] y = z sin z (θ² + 1) y = z sin z
---
Step 2: Complementary Function (y_c)
Auxiliary equation: m² + 1 = 0 => m = ±i y_c = c₁ cos z + c₂ sin z
---
Step 3: Particular Integral (y_p)
y_p = [1 / (θ² + 1)] (z sin z)
We know that z sin z = Imaginary Part of (z e^(iz)): y_p = Im { [1 / (θ² + 1)] (z e^(iz)) }
Using exponential shift property: = Im { e^(iz) [1 / ((θ + i)² + 1)] z } = Im { e^(iz) [1 / (θ² + 2iθ - 1 + 1)] z } = Im { e^(iz) [1 / (θ(θ + 2i))] z } = Im { e^(iz) [1 / (2iθ (1 + θ / 2i))] z } = Im { [e^(iz) / (2i)] (1/θ) [ 1 - θ/(2i) ] z } = Im { [-i e^(iz) / 2] (1/θ) [ z - 1/(2i) ] } = Im { [-i e^(iz) / 2] (1/θ) [ z + i/2 ] }
Integrating with respect to z: (1/θ) [ z + i/2 ] = z²/2 + iz/2
Now multiply: [-i e^(iz) / 2] [ z²/2 + iz/2 ] = (1/4) e^(iz) [ -i z² + z ] = (1/4) (cos z + i sin z) (z - i z²) = (1/4) [ z cos z + z² sin z + i (z sin z - z² cos z) ]
Taking the Imaginary Part: y_p = (1/4) [ z sin z - z² cos z ]
---
Step 4: Express General Solution in Terms of x
Substitute z = log x into yc and yp: y = c₁ cos(log x) + c₂ sin(log x) + (1/4) [ log x · sin(log x) - (log x)² · cos(log x) ]
Q.4 (a)
Question
Find the Laplace Transform of f(t) = sinh at. [03 Marks]
Solution
Definition
The Laplace transform of a function f(t) is defined as: L{f(t)} = ∫₀^∞ e^(-st) f(t) dt
---
Derivation
By definition, sinh at = (e^(at) - e^(-at)) / 2.
Applying linearity of the Laplace transform: L{sinh at} = L{ (e^(at) - e^(-at)) / 2 } = (1/2) [ L{e^(at)} - L{e^(-at)} ]
We know that: L{e^(at)} = 1 / (s - a), for s > a L{e^(-at)} = 1 / (s + a), for s > -a
Therefore: L{sinh at} = (1/2) [ 1/(s - a) - 1/(s + a) ] = (1/2) [ ( (s + a) - (s - a) ) / ((s - a)(s + a)) ] = (1/2) [ (2a) / (s² - a²) ] = a / (s² - a²), for s > |a|
Final Answer: L{sinh at} = a / (s² - a²).
Q.4 (b)
Question
Find the Laplace Transform of t e^(4t) cos 2t. [04 Marks]
Solution
Step 1: Laplace Transform of cos 2t
L{cos 2t} = s / (s² + 2²) = s / (s² + 4)
---
Step 2: Multiplication by t Property
By property: L{t f(t)} = - d/ds [ L{f(t)} ]
L{t cos 2t} = - d/ds [ s / (s² + 4) ] = - [ ( (s² + 4)(1) - s(2s) ) / (s² + 4)² ] = - [ (s² + 4 - 2s²) / (s² + 4)² ] = - [ (4 - s²) / (s² + 4)² ] = (s² - 4) / (s² + 4)²
---
Step 3: First Shifting Theorem
By the first shifting theorem: L{e^(at) g(t)} = G(s - a) Here, a = 4. Replace s with s - 4:
L{e^(4t) · (t cos 2t)} = [ (s - 4)² - 4 ] / [ (s - 4)² + 4 ]² = [ s² - 8s + 16 - 4 ] / [ s² - 8s + 16 + 4 ]² = (s² - 8s + 12) / (s² - 8s + 20)²
Final Answer: (s² - 8s + 12) / (s² - 8s + 20)²
Q.4 (c)
Question
Solve the differential equation using initial value problem: y'' + 6y' = 1, y(0) = 2, y'(0) = 0. [07 Marks]
Solution
Step 1: Apply Laplace Transform to Both Sides
L{y'' + 6y'} = L{1} L{y''} + 6 L{y'} = 1/s
Using initial value formulas:
- L{y'} = s Y(s) - y(0)
- L{y''} = s² Y(s) - s y(0) - y'(0)
Substitute y(0) = 2 and y'(0) = 0: [ s² Y(s) - 2s - 0 ] + 6 [ s Y(s) - 2 ] = 1/s (s² + 6s) Y(s) - 2s - 12 = 1/s s(s + 6) Y(s) = 2s + 12 + 1/s = (2s² + 12s + 1) / s
Y(s) = (2s² + 12s + 1) / [ s²(s + 6) ]
---
Step 2: Partial Fraction Decomposition
Let: (2s² + 12s + 1) / [ s²(s + 6) ] = A/s + B/s² + C/(s + 6)
Multiply through by s²(s + 6): 2s² + 12s + 1 = A s (s + 6) + B (s + 6) + C s²
- Put s = 0:
1 = B(6) => B = 1/6
- Put s = -6:
2(-6)² + 12(-6) + 1 = C(-6)²
2(36) - 72 + 1 = 36C
72 - 72 + 1 = 36C => C = 1/36
- Equate coefficients of s²:
2 = A + C
A = 2 - 1/36 = 71/36
Thus: Y(s) = (71/36)(1/s) + (1/6)(1/s²) + (1/36)(1/(s + 6))
---
Step 3: Taking Inverse Laplace Transform
y(t) = L⁻¹{Y(s)} = (71/36) L⁻¹{1/s} + (1/6) L⁻¹{1/s²} + (1/36) L⁻¹{1/(s + 6)} = (71/36)(1) + (1/6) t + (1/36) e^(-6t)
Final Answer: y(t) = 71/36 + (1/6)t + (1/36)e^(-6t)
Q.4 (a) OR
Question
Find the inverse Laplace transform of log((s² + b²) / (s² + a²)). [03 Marks]
Solution
Step 1: Use Property of Derivative of Laplace Transform
Let F(s) = log((s² + b²) / (s² + a²)) = log(s² + b²) - log(s² + a²).
We know that: L⁻¹{ F'(s) } = -t f(t) => f(t) = - (1/t) L⁻¹{ F'(s) }
---
Step 2: Differentiate F(s) with Respect to s
F'(s) = d/ds [ log(s² + b²) - log(s² + a²) ] = 2s / (s² + b²) - 2s / (s² + a²)
---
Step 3: Take Inverse Laplace Transform
L⁻¹{ F'(s) } = 2 L⁻¹{ s / (s² + b²) } - 2 L⁻¹{ s / (s² + a²) } = 2 cos bt - 2 cos at
---
Step 4: Compute f(t)
Since -t f(t) = 2(cos bt - cos at): f(t) = - [ 2(cos bt - cos at) ] / t = 2 (cos at - cos bt) / t
Final Answer: L⁻¹{ log((s² + b²) / (s² + a²)) } = 2(cos at - cos bt) / t
Q.4 (b) OR
Question
Find the Laplace Transform of [e^(-2t) sin 2t cosh t] / t. [04 Marks]
Solution
Step 1: Simplify the Numerator Function
Let g(t) = sin 2t cosh t. Using definition: cosh t = (e^t + e^(-t)) / 2 g(t) = (1/2) [ e^t sin 2t + e^(-t) sin 2t ]
We know that L{sin 2t} = 2 / (s² + 4). By shifting theorem: L{g(t)} = (1/2) [ 2 / ((s - 1)² + 4) + 2 / ((s + 1)² + 4) ] = 1 / ((s - 1)² + 4) + 1 / ((s + 1)² + 4)
---
Step 2: Apply Division by t Property
L{ g(t) / t } = ∫s^∞ L{g(t)} du = ∫s^∞ [ 1 / ((u - 1)² + 4) + 1 / ((u + 1)² + 4) ] du
Using standard integral ∫ 1 / (x² + a²) dx = (1/a) arctan(x/a): = [ (1/2) arctan((u - 1)/2) + (1/2) arctan((u + 1)/2) ]_s^∞ = (1/2) [ (π/2 + π/2) - (arctan((s - 1)/2) + arctan((s + 1)/2)) ] = (1/2) [ π - arctan((s - 1)/2) - arctan((s + 1)/2) ]
Since π/2 - arctan(x) = arccot(x): = (1/2) [ arccot((s - 1)/2) + arccot((s + 1)/2) ]
---
Step 3: Apply First Shifting for e^(-2t)
By shifting property, replace s with s + 2: L{ e^(-2t) [sin 2t cosh t / t] } = (1/2) [ arccot((s + 1)/2) + arccot((s + 3)/2) ]
Or equivalently in terms of arctan: = (1/2) [ π - arctan((s + 1)/2) - arctan((s + 3)/2) ]
Q.4 (c) OR
Question
Find the Inverse Laplace Transform of (5s + 3) / [ (s - 1)(s² + 2s + 5) ]. [07 Marks]
Solution
Step 1: Partial Fraction Decomposition
Let: (5s + 3) / [ (s - 1)(s² + 2s + 5) ] = A / (s - 1) + (Bs + C) / (s² + 2s + 5)
Multiplying by (s - 1)(s² + 2s + 5): 5s + 3 = A(s² + 2s + 5) + (Bs + C)(s - 1)
- Set s = 1:
5(1) + 3 = A(1 + 2 + 5) + 0
8 = 8A => A = 1
- Equate coefficients of s²:
0 = A + B => B = -A => B = -1
- Equate constant terms:
3 = 5A - C
3 = 5(1) - C => C = 5 - 3 => C = 2
So the expression becomes: 1 / (s - 1) + (-s + 2) / (s² + 2s + 5) = 1 / (s - 1) - (s - 2) / (s² + 2s + 5)
---
Step 2: Complete the Square in Denominator
s² + 2s + 5 = (s + 1)² + 4 = (s + 1)² + 2²
Rewrite numerator in terms of (s + 1): s - 2 = (s + 1) - 3
Therefore: (s - 2) / [ (s + 1)² + 2² ] = (s + 1) / [ (s + 1)² + 2² ] - 3 / [ (s + 1)² + 2² ] = (s + 1) / [ (s + 1)² + 2² ] - (3/2) · 2 / [ (s + 1)² + 2² ]
---
Step 3: Invert Term by Term
- L⁻¹{ 1 / (s - 1) } = e^t
- L⁻¹{ (s + 1) / [ (s + 1)² + 2² ] } = e^(-t) cos 2t
- L⁻¹{ 2 / [ (s + 1)² + 2² ] } = e^(-t) sin 2t
Combining everything: f(t) = L⁻¹{ 1/(s - 1) } - [ L⁻¹{ (s + 1)/((s + 1)² + 2²) } - (3/2) L⁻¹{ 2/((s + 1)² + 2²) } ] = e^t - e^(-t) cos 2t + (3/2) e^(-t) sin 2t
Final Answer: y(t) = e^t + e^(-t) [ (3/2) sin 2t - cos 2t ]
Q.5 (a)
Question
Solve dy/dx + y cot x = 2 cos x. [03 Marks]
Solution
Step 1: Identify Equation Type
This is a standard first-order linear differential equation: dy/dx + P(x) y = Q(x) where:
- P(x) = cot x
- Q(x) = 2 cos x
---
Step 2: Integrating Factor (I.F.)
I.F. = e^(∫ P(x) dx) = e^(∫ cot x dx) = e^(log(sin x)) = sin x
---
Step 3: General Solution
The general solution is given by: y · (I.F.) = ∫ Q(x) · (I.F.) dx + C
y · sin x = ∫ (2 cos x) (sin x) dx + C y · sin x = ∫ sin 2x dx + C y · sin x = - (cos 2x) / 2 + C
Alternatively, since 2 sin x cos x = d(sin² x)/dx: y · sin x = sin² x + C y = sin x + C csc x
Q.5 (b)
Question
Solve 3x⁴p² - xp - y = 0. [04 Marks]
Solution
Step 1: Express in Solvable Form for y
The given differential equation is of first order and higher degree (where p = dy/dx): 3x⁴p² - xp - y = 0 => y = 3x⁴p² - xp
---
Step 2: Substitution Method
Let X = 1/x and Y = y. (Alternatively, let u = 1/x, which transforms Clairaut-like equations). Let us substitute: u = x² or try transforming coordinates. Notice that 3x⁴p² = 3(x²p)². Let v = x² and w = y. Then: dy/dx = (dw/dv)(dv/dx) = P · 2x => p = 2x P (where P = dw/dv). Substitute into the equation: 3x⁴(2x P)² - x(2x P) - y = 0 => not Clairaut's form directly.
Let's rewrite the equation by putting y = 3(x²p)² - (xp): Let u = 1/x, then: p = dy/dx = (dy/du)(du/dx) = (-1/x²) (dy/du) = -u² Pu, where Pu = dy/du. Then x p = (1/u)(-u² Pu) = -u Pu. And x⁴ p² = (x² p)² = [ (1/u²) (-u² Pu) ]² = Pu².
Substitute into the original equation: 3(Pu)² - (-u Pu) - y = 0 => y = u Pu + 3 Pu²
---
Step 3: Clairaut’s Form
This is in Clairaut's form: y = u Pu + f(Pu), where f(Pu) = 3 Pu².
The general solution of Clairaut's equation is obtained directly by replacing P_u with an arbitrary constant c: y = c u + 3 c²
---
Step 4: Substitute Back u = 1/x
y = c (1/x) + 3c² xy = c + 3c² x or 3c²x - xy + c = 0 where c is an arbitrary constant.
Q.5 (c)
Question
Using the method variation of parameter solve the differential equation (D² + 1)y = x sin x. [07 Marks]
Solution
Step 1: Complementary Function (y_c)
(D² + 1)y = 0 => m² + 1 = 0 => m = ±i y_c = c₁ cos x + c₂ sin x
Here:
- y₁ = cos x
- y₂ = sin x
- R(x) = x sin x
---
Step 2: Wronskian (W)
W(y₁, y₂) = | y₁ y₂ | | y₁' y₂' | = | cos x sin x | | -sin x cos x | = cos² x - (-sin² x) = cos² x + sin² x = 1
---
Step 3: Method of Variation of Parameters
The particular solution is: y_p = u(x) y₁ + v(x) y₂ where: u(x) = - ∫ [y₂ R(x) / W] dx = - ∫ sin x (x sin x) / 1 dx = - ∫ x sin² x dx v(x) = ∫ [y₁ R(x) / W] dx = ∫ cos x (x sin x) / 1 dx = ∫ x sin x cos x dx = (1/2) ∫ x sin 2x dx
---
Step 4: Evaluate u(x)
Since sin² x = (1 - cos 2x) / 2: u(x) = - (1/2) ∫ x (1 - cos 2x) dx = - (1/2) [ ∫ x dx - ∫ x cos 2x dx ] = - (1/2) [ x²/2 - ( x (sin 2x / 2) - ∫ (sin 2x / 2) dx ) ] = - (1/2) [ x²/2 - (x sin 2x)/2 - (cos 2x)/4 ] = - x²/4 + (x sin 2x)/4 + (cos 2x)/8
---
Step 5: Evaluate v(x)
v(x) = (1/2) ∫ x sin 2x dx Using integration by parts: = (1/2) [ x (-cos 2x / 2) - ∫ 1 · (-cos 2x / 2) dx ] = (1/2) [ - (x cos 2x)/2 + (sin 2x)/4 ] = - (x cos 2x)/4 + (sin 2x)/8
---
Step 6: Form Particular Integral y_p
y_p = u(x) cos x + v(x) sin x = [ -x²/4 + (x sin 2x)/4 + (cos 2x)/8 ] cos x + [ - (x cos 2x)/4 + (sin 2x)/8 ] sin x = - (x²/4) cos x + (x/4) [ sin 2x cos x - cos 2x sin x ] + (1/8) [ cos 2x cos x + sin 2x sin x ]
Using trigonometric identities:
- sin 2x cos x - cos 2x sin x = sin(2x - x) = sin x
- cos 2x cos x + sin 2x sin x = cos(2x - x) = cos x
So: y_p = - (x²/4) cos x + (x/4) sin x + (1/8) cos x
Notice that (1/8) cos x can be absorbed into the complementary function c₁ cos x.
---
Step 7: General Solution
y = c₁ cos x + c₂ sin x - (x²/4) cos x + (x/4) sin x
Q.5 (a) OR
Question
Solve (xy - 2y²) dx - (x² - 3xy) dy = 0. [03 Marks]
Solution
Step 1: Check Homogeneity
Rewrite in the form M dx + N dy = 0:
- M = xy - 2y² (degree 2)
- N = -(x² - 3xy) = -x² + 3xy (degree 2)
Since both M and N are homogeneous functions of degree 2, this is a homogeneous differential equation.
---
Step 2: Substitution Method
Put y = v x => dy = v dx + x dv
Substitute into the equation: [ x(vx) - 2(vx)² ] dx - [ x² - 3x(vx) ] (v dx + x dv) = 0 [ vx² - 2v²x² ] dx - [ x² - 3vx² ] (v dx + x dv) = 0
Divide by x² (for x ≠ 0): (v - 2v²) dx - (1 - 3v)(v dx + x dv) = 0 [ (v - 2v²) - (v - 3v²) ] dx - x (1 - 3v) dv = 0 [ v - 2v² - v + 3v² ] dx - x (1 - 3v) dv = 0 v² dx - x (1 - 3v) dv = 0
---
Step 3: Separate Variables and Integrate
v² dx = x (1 - 3v) dv dx / x = [ (1 - 3v) / v² ] dv dx / x = [ 1/v² - 3/v ] dv
Integrate both sides: ∫ dx / x = ∫ v^(-2) dv - 3 ∫ (1/v) dv log|x| = - 1/v - 3 log|v| + C log|x| + 3 log|v| + 1/v = C log|x · v³| + 1/v = C
---
Step 4: Substitute Back v = y/x
log| x · (y/x)³ | + x/y = C log| y³ / x² | + x/y = C
Q.5 (b) OR
Question
Solve dy/dx + (1/x) y = y² / x². [04 Marks]
Solution
Step 1: Recognize Bernoulli’s Differential Equation
The equation is of the form: dy/dx + P(x) y = Q(x) yⁿ with n = 2, P(x) = 1/x, and Q(x) = 1/x².
---
Step 2: Divide by y²
(1/y²) dy/dx + (1/x)(1/y) = 1/x²
---
Step 3: Substitution
Let z = 1/y = y⁻¹ Then dz/dx = - y⁻² (dy/dx) = - (1/y²) (dy/dx) So: (1/y²) dy/dx = - dz/dx
Substitute into the equation:
- dz/dx + (1/x) z = 1/x²
dz/dx - (1/x) z = - 1/x²
This is a standard linear differential equation in z.
---
Step 4: Integrating Factor (I.F.)
I.F. = e^(∫ - (1/x) dx) = e^(-log x) = 1/x
---
Step 5: Solve for z
z · (I.F.) = ∫ Q(x) · (I.F.) dx + C z · (1/x) = ∫ (-1/x²)(1/x) dx + C z / x = - ∫ x⁻³ dx + C z / x = - [ x⁻² / (-2) ] + C = 1 / (2x²) + C
Multiply through by x: z = 1 / (2x) + Cx
---
Step 6: Substitute Back z = 1/y
1/y = 1 / (2x) + Cx = (1 + 2Cx²) / (2x)
y = 2x / (1 + 2Cx²) where C is an arbitrary constant.
Q.5 (c) OR
Question
Solve (D⁴ - 16)y = e^(2x) + x⁴. [07 Marks]
Solution
Step 1: Complementary Function (y_c)
Auxiliary Equation: m⁴ - 16 = 0 (m² - 4)(m² + 4) = 0 (m - 2)(m + 2)(m² + 4) = 0 Roots are: m = 2, -2, ±2i
y_c = c₁ e^(2x) + c₂ e^(-2x) + c₃ cos 2x + c₄ sin 2x
---
Step 2: Particular Integral y_p1 (for e^(2x))
y_p1 = [1 / (D⁴ - 16)] e^(2x)
Since m = 2 is a single root of the auxiliary equation, substituting D = 2 gives 0 (case of failure). Factor out the failing factor: 1 / (D⁴ - 16) = 1 / [ (D - 2)(D + 2)(D² + 4) ]
Substitute D = 2 into the non-zero terms: y_p1 = [1 / (D - 2)] [ 1 / ((2 + 2)(2² + 4)) ] e^(2x) = [1 / (D - 2)] [ 1 / (4 · 8) ] e^(2x) = (1/32) [1 / (D - 2)] e^(2x)
Using the standard formula [1 / (D - a)] e^(ax) = x e^(ax): y_p1 = (x / 32) e^(2x)
---
Step 3: Particular Integral y_p2 (for x⁴)
y_p2 = [1 / (D⁴ - 16)] x⁴ = [1 / -16(1 - D⁴/16)] x⁴ = - (1/16) [ 1 - D⁴/16 ]⁻¹ x⁴
Expanding binomially: (1 - u)⁻¹ = 1 + u + u² + ...: = - (1/16) [ 1 + D⁴/16 + ... ] x⁴
Apply the derivatives to x⁴:
- D(x⁴) = 4x³
- D²(x⁴) = 12x²
- D³(x⁴) = 24x
- D⁴(x⁴) = 24
So: y_p2 = - (1/16) [ x⁴ + 24/16 ] = - (1/16) [ x⁴ + 3/2 ] = - x⁴/16 - 3/32
---
Step 4: General Solution
y = yc + yp1 + y_p2 y = c₁ e^(2x) + c₂ e^(-2x) + c₃ cos 2x + c₄ sin 2x + (x / 32) e^(2x) - x⁴/16 - 3/32
Frequently Asked Questions
What is the condition for a vector field to be both solenoidal and irrotational?
A vector field F is solenoidal if its divergence is zero (∇ · F = 0), meaning there are no sources or sinks. It is irrotational if its curl is zero (∇ × F = 0), meaning it has no rotational tendency and can be expressed as the gradient of a scalar potential.
How do you handle a case of failure when finding a particular integral with D² + a² and cos(ax)?
When substituting D² = -a² into 1/(D² + a²) cos(ax), the denominator becomes zero. To solve this, convert cos(ax) into the real part of e^(iax), apply the exponential shift property, and integrate to get the standard shortcut formula: [1 / (D² + a²)] cos(ax) = (x / 2a) sin(ax).
When does Green's Theorem apply?
Green's Theorem applies to a positively oriented, piecewise-smooth, simple closed curve C in the xy-plane that bounds a region R, provided the functions M(x, y) and N(x, y) have continuous partial derivatives in an open region containing R.