MAT1252 Partial Differential Equations Final Exam Study Guide June 2025

Examination Overview: MAT1252 - Partial Differential Equations

The Final Examination for the Even Semester of the 2024/2025 academic year for the course MAT1252 - Partial Differential Equations (Persamaan Diferensial Parsial) was held on Wednesday, June 18, 2025. The examination took place from 10:30 to 12:30. The testing protocol was strictly defined as Closed Notes (Catatan Tertutup) and the use of calculators was prohibited (Tanpa Kalkulator). The document credits Ikhsan Rio Afriadi (@ikhsanrio16) and covers fundamental topics in partial differential equations, including wave equations, Fourier series, Laplace equations, and Sturm-Liouville problems.

Derivation of the One-Dimensional Wave Equation

The first requirement of the examination is to detail the steps for the derivation of the one-dimensional wave equation, expressed in the form utt=αuxxu_{tt} = \alpha u_{xx}. This derivation typically involves analyzing the transverse vibrations of a flexible string. The process begins by considering a small segment of a string with constant linear density ρ\rho under a tension force TT. By applying Newton's Second Law (\mathbf{F} = m·\mathbf{a}) to the vertical components of the tension forces at the endpoints of a segment [x,x+Δx][x, x + \Delta x], one identifies the net force acting on the segment. Assuming small angles of displacement, the vertical forces can be approximated using the slopes of the string, ux(x,t)u_x(x, t).

The mass of the segment is represented by ρΔx\rho \Delta x, and its vertical acceleration is uttu_{tt}. Equating the mass-acceleration product to the difference in the vertical components of tension results in the equation ρΔxuttT[ux(x+Δx,t)ux(x,t)]\rho \Delta x u_{tt} \approx T [u_x(x + \Delta x, t) - u_x(x, t)]. Dividing by Δx\Delta x and taking the limit as Δx0\Delta x \rightarrow 0 leads to the second-order partial differential equation utt=Tρuxxu_{tt} = \frac{T}{\rho} u_{xx}, where α=Tρ\alpha = \frac{T}{\rho} represents the square of the wave propagation speed, often denoted as c2c^2.

Determination of Fourier Coefficients for Periodic Functions

The second problem focuses on a periodic function f(x)f(x) with a period of 2π2\pi. The function is defined piecewise as follows: f(x)=1f(x) = -1 for πx<0-\pi \leq x < 0, and f(x)=x1f(x) = x - 1 for 0x<π0 \leq x < \pi. The goal is to determine the Fourier coefficient ana_n for the general Fourier series representation f(x)=a02+n=1(ancos(nx)+bnsin(nx))f(x) = \frac{a_0}{2} + \sum_{n=1}^{\infty} (a_n \cos(nx) + b_n \sin(nx)).

To find ana_n, one must evaluate the integral an=1πππf(x)cos(nx)dxa_n = \frac{1}{\pi} \int_{-\pi}^{\pi} f(x) \cos(nx) dx. This integral must be split into two parts according to the definition of f(x)f(x): an=1π[π0(1)cos(nx)dx+0π(x1)cos(nx)dx]a_n = \frac{1}{\pi} [ \int_{-\pi}^{0} (-1) \cos(nx) dx + \int_{0}^{\pi} (x - 1) \cos(nx) dx ].

The first part of the integral involves the basic integration of a cosine function, while the second part requires integration by parts for the term xcos(nx)x \cos(nx). The integration must be performed for n=1,2,3,n = 1, 2, 3, \dots to determine the general formula for the coefficients.

Fourier Cosine Series of a Sine Function

The third problem asks for the Fourier Cosine Series of the function f(x)=sin(x)f(x) = \sin(x) on the interval 0<x<π0 < x < \pi, assuming a period of 2π2\pi. Because the cosine series is requested, the function must be treated as an even extension of itself over the interval [π,π][-\pi, \pi]. Even though the function itself is a sine, the cosine series representation for an even extension on the interval [0,L][0, L] (where L=πL = \pi) is given by: f(x)=a02+n=1ancos(nπxL)f(x) = \frac{a_0}{2} + \sum_{n=1}^{\infty} a_n \cos(\frac{n \pi x}{L}).

For this specific case, L=πL = \pi, so the series becomes f(x)=a02+n=1ancos(nx)f(x) = \frac{a_0}{2} + \sum_{n=1}^{\infty} a_n \cos(nx). The coefficients are calculated using the formula an=2π0πsin(x)cos(nx)dxa_n = \frac{2}{\pi} \int_{0}^{\pi} \sin(x) \cos(nx) dx. This integral is typically solved using trigonometric product-to-sum identities: sin(x)cos(nx)=12[sin((1+n)x)+sin((1n)x)]\sin(x) \cos(nx) = \frac{1}{2} [ \sin((1+n)x) + \sin((1-n)x) ].

Solving the One-Dimensional Wave Equation with Specific Boundaries

The fourth problem presents a specific initial-boundary value problem (IBVP) for the wave equation: utt=4uxxu_{tt} = 4u_{xx} for the spatial domain 0<x<300 < x < 30 and time t>0t > 0.

The Boundary Conditions (SB) are homogeneous Dirichlet conditions: u(0,t)=0u(0, t) = 0 and u(30,t)=0u(30, t) = 0 for all t>0t > 0.

The Initial Conditions (SA) are given by:

  • The initial displacement u(x,0)=f(x)u(x, 0) = f(x).
  • The initial velocity ut(x,0)=0u_t(x, 0) = 0.

The function f(x)f(x) is defined as a triangle-like distribution: f(x)=x10f(x) = \frac{x}{10} for 0<x100 < x \leq 10. f(x)=30x20f(x) = \frac{30 - x}{20} for 10<x3010 < x \leq 30.

Since the initial velocity is zero, the solution takes the form of a separation of variables solution: u(x,t)=n=1Ansin(nπx30)cos(2nπt30)u(x, t) = \sum_{n=1}^{\infty} A_n \sin(\frac{n \pi x}{30}) \cos(\frac{2 n \pi t}{30}). The coefficients AnA_n must be determined by the Fourier Sine Series of the initial displacement function f(x)f(x) on the interval [0,30][0, 30], calculated as An=230030f(x)sin(nπx30)dxA_n = \frac{2}{30} \int_{0}^{30} f(x) \sin(\frac{n \pi x}{30}) dx.

Laplace's Equation on a Semi-Infinite Rectangular Plate

The fifth problem requires solving Laplace's equation u(x,y)u(x, y) on a rectangular plate defined by the boundaries 0<x<500 < x < 50 and 0<y<0 < y < \infty (a semi-infinite strip). The boundary conditions are:

  • u(0,y)=0u(0, y) = 0 and u(50,y)=0u(50, y) = 0 for 0<y<0 < y < \infty (lateral boundaries held at zero).
  • u(x,)=0u(x, \infty) = 0 (implied by the requirement for a physically bounded solution as yy approaches infinity).
  • The bottom boundary condition is given by u(x,0)=f(x)u(x, 0) = f(x), where: f(x)=xf(x) = x for 0<x250 < x \leq 25. f(x)=50xf(x) = 50 - x for 25x5025 \leq x \leq 50.

The general solution for Laplace's equation in this geometry is u(x,y)=n=1Cnsin(nπx50)enπy50u(x, y) = \sum_{n=1}^{\infty} C_n \sin(\frac{n \pi x}{50}) e^{-\frac{n \pi y}{50}}. The coefficients CnC_n are found by evaluating the Fourier Sine Series of f(x)f(x) on the interval [0,50][0, 50]: Cn=250050f(x)sin(nπx50)dxC_n = \frac{2}{50} \int_{0}^{50} f(x) \sin(\frac{n \pi x}{50}) dx.

Eigenvalues and Eigenfunctions of a Sturm-Liouville Problem

The final problem investigates a Sturm-Liouville problem defined by the differential equation y+λy=0y'' + \lambda y = 0 on the interval 0 \leq x \leq 1. The boundary conditions are mixed:

  • At x=0x = 0, the derivative is zero: y(0)=0y'(0) = 0.
  • At x=1x = 1, there is a boundary condition involving both the function and its derivative: y(1)+y(1)=0y'(1) + y(1) = 0.

To solve this, one must analyze three cases for the eigenvalue λ\lambda.

  1. If λ<0\lambda < 0 (let λ=k2\lambda = -k^2), the solution follows the form y(x)=c1cosh(kx)+c2sinh(kx)y(x) = c_1 \cosh(kx) + c_2 \sinh(kx). Setting y(0)=0y'(0) = 0 leads to c2=0c_2 = 0. The second condition leads to a transcendental equation that must be checked for non-trivial solutions.
  2. If λ=0\lambda = 0, the solution is linear: y(x)=c1x+c2y(x) = c_1 x + c_2. Setting y(0)=0y'(0) = 0 implies c1=0c_1 = 0, meaning y(x)=c2y(x) = c_2. The second condition y(1)+y(1)=0y'(1) + y(1) = 0 results in 0+c2=00 + c_2 = 0, yielding only the trivial solution.
  3. If λ>0\lambda > 0 (let λ=k2\lambda = k^2), the solution is y(x)=c1cos(kx)+c2sin(kx)y(x) = c_1 \cos(kx) + c_2 \sin(kx). The condition y(0)=0y'(0) = 0 forces c2=0c_2 = 0. Applying the second condition, y(1)+y(1)=0y'(1) + y(1) = 0, results in the equation ksin(k)+cos(k)=0-k \sin(k) + \cos(k) = 0, or tan(k)=1k\tan(k) = \frac{1}{k}. The values of kk that satisfy this transcendental equation define the eigenvalues λn=kn2\lambda_n = k_n^2, and the corresponding eigenfunctions are yn(x)=cos(knx)y_n(x) = \cos(k_n x).