Skip to main content

Chapter 4

Differential Equations

TU BSc past exam questions

Past questions and answers

17 questions set from this chapter, 1 of them more than once. Most repeated first.

  • Asked 2 times
  • 2082 · Differential Equation (MAT 202) (New Course)
  • 2081 · Differential Equation (MAT 202) · 8 marks

Consider the system x1′=x1−2x2, x1(0)=1x_1' = x_1 - 2x_2,\ x_1(0) = 1 x2′=3x1−4x2, x2(0)=2x_2' = 3x_1 - 4x_2,\ x_2(0) = 2 (a) Transform the above system into a single equation of second order. (b) Find x1x_1 and x2x_2 that also satisfy the given initial conditions. [4+4]

Answer

(a) Reduction to a single second-order equation

Given [ x_1'=x_1-2x_2,\qquad x_2'=3x_1-4x_2. ]

From the first equation, [ 2x_2=x_1-x_1' ] so [ x_2=\frac{x_1-x_1'}{2}. ]

Differentiate: [ x_2'=\frac{x_1'-x_1''}{2}. ]

Substitute this and (x_2=\frac{x_1-x_1'}{2}) into the second equation: [ \frac{x_1'-x_1''}{2}=3x_1-4\left(\frac{x_1-x_1'}{2}\right). ]

Simplify the right-hand side: [ 3x_1-2(x_1-x_1')=3x_1-2x_1+2x_1'=x_1+2x_1'. ]

Thus [ \frac{x_1'-x_1''}{2}=x_1+2x_1'. ]

Multiplying by (2): [ x_1'-x_1''=2x_1+4x_1', ] so [ x_1''+3x_1'+2x_1=0. ]

Hence the required second-order equation is [ \boxed{x_1''+3x_1'+2x_1=0}. ]

Also, from (x_2=\frac{x_1-x_1'}{2}), once (x_1) is known, (x_2) follows.


(b) Solution satisfying the initial conditions

The initial conditions are [ x_1(0)=1,\qquad x_2(0)=2. ]

From (x_1'=x_1-2x_2), [ x_1'(0)=x_1(0)-2x_2(0)=1-4=-3. ]

So we solve [ x_1''+3x_1'+2x_1=0,\qquad x_1(0)=1,\quad x_1'(0)=-3. ]

The characteristic equation is [ r^2+3r+2=0 ] [ (r+1)(r+2)=0. ]

Thus [ r=-1,,-2. ]

Therefore, [ x_1(t)=C_1e^{-t}+C_2e^{-2t}. ]

Using (x_1(0)=1): [ C_1+C_2=1. \quad (1) ]

Differentiate: [ x_1'(t)=-C_1e^{-t}-2C_2e^{-2t}. ]

Using (x_1'(0)=-3): [ -C_1-2C_2=-3 ] or [ C_1+2C_2=3. \quad (2) ]

Subtracting (1) from (2): [ C_2=2. ]

Then from (1): [ C_1=-1. ]

Hence [ \boxed{x_1(t)=-e^{-t}+2e^{-2t}}. ]

Now use [ x_2=\frac{x_1-x_1'}{2}. ]

We have [ x_1'=e^{-t}-4e^{-2t}. ]

Therefore [ x_2(t)=\frac{(-e^{-t}+2e^{-2t})-(e^{-t}-4e^{-2t})}{2} ] [ =\frac{-2e^{-t}+6e^{-2t}}{2} ] [ =-e^{-t}+3e^{-2t}. ]

So [ \boxed{x_2(t)=-e^{-t}+3e^{-2t}}. ]


Final answer: [ \boxed{x_1(t)=-e^{-t}+2e^{-2t}} ] [ \boxed{x_2(t)=-e^{-t}+3e^{-2t}} ]

These satisfy (x_1(0)=1) and (x_2(0)=2).

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2082 · Differential Equation (MAT 202) (New Course)

(b) If the interest rate on 30 years mortgage is fixed at 10% and if a monthly payment of $ 1000 is the maximum that the buyer can afford, what is the maximum mortgage loan that can be made under these conditions?

Answer

The maximum loan is the present value of an ordinary annuity of 1000 per month for 30 years.

Given: monthly payment R=1000R = 1000; monthly rate i=0.1012=0.0083333i = \dfrac{0.10}{12} = 0.0083333; number of payments n=30×12=360n = 30 \times 12 = 360.

PV=R[1−(1+i)−ni]=1000[1−(1.0083333)−3600.0083333]PV = R\left[\frac{1-(1+i)^{-n}}{i}\right] = 1000\left[\frac{1-(1.0083333)^{-360}}{0.0083333}\right]

Since (1.0083333)360≈19.837(1.0083333)^{360} \approx 19.837, we have (1.0083333)−360≈0.05041(1.0083333)^{-360} \approx 0.05041. Then

PV=1000×1−0.050410.0083333=1000×113.95≈113,951PV = 1000\times\frac{1-0.05041}{0.0083333} = 1000 \times 113.95 \approx 113{,}951

Maximum mortgage loan ≈\approx 113,951 (in the same currency as the 1000 monthly payment, which the question gives in dollars).

Answer checked by Claude Sonnet · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2082 · Differential Equation (MAT 202) (New Course) · 8 marks

Find the general solution of the given initial value problem y′′+4y=2sin⁡2t, y(0)=2, y′(0)=−1.y'' + 4y = 2\sin 2t,\ y(0) = 2,\ y'(0) = -1.
OR
A mass weighing 4 lb stretches 1.5 inches. The mass is displaced 2 inches in the positive direction from its equilibrium position and released with no initial velocity. Assuming that there is no damping and that the mass is acted on by external force of 2cos⁡3t2\cos 3t is formulated the IVP describing the motion and solve. [4+4]

Answer

Taking the first alternative:

Given [ y''+4y=2\sin 2t,\qquad y(0)=2,\quad y'(0)=-1. ]

Homogeneous equation: [ y''+4y=0. ] Characteristic equation: [ r^2+4=0 \implies r=\pm 2i. ] So [ y_c=c_1\cos 2t+c_2\sin 2t. ]

Since the right-hand side is (\sin 2t), which is already a solution of the homogeneous equation, we try [ y_p=t(A\cos 2t+B\sin 2t). ] Differentiating, [ y_p''+4y_p=-4A\sin 2t+4B\cos 2t. ] Comparing with (2\sin 2t), [ -4A=2 \implies A=-\frac12,\qquad 4B=0 \implies B=0. ] Hence [ y_p=-\frac{t}{2}\cos 2t. ]

Therefore the general solution is [ y=c_1\cos 2t+c_2\sin 2t-\frac{t}{2}\cos 2t. ]

Now apply the initial conditions.

[ y(0)=c_1=2. ]

Also, [ y'=-2c_1\sin 2t+2c_2\cos 2t-\frac12\cos 2t+t\sin 2t. ] At (t=0), [ y'(0)=2c_2-\frac12=-1. ] Thus [ 2c_2=-\frac12 \implies c_2=-\frac14. ]

So the solution of the initial value problem is [ \boxed{y(t)=2\cos 2t-\frac14\sin 2t-\frac{t}{2}\cos 2t.} ]

Check: [ y(0)=2,\qquad y'(0)=-\frac12-\frac12=-1. ] Hence the required solution is verified.

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2082 · Differential Equation (MAT 202) (New Course) · 8 marks

Find the general solution of the given initial value problem y′′+4y=t2+3et, y(0)=0, y′(0)=2y'' + 4y = t^2 + 3e^t,\ y(0) = 0,\ y'(0) = 2

Answer

For the differential equation [ y''+4y=t^2+3e^t, ] first solve the homogeneous equation: [ y''+4y=0. ] Its characteristic equation is [ r^2+4=0 \implies r=\pm 2i. ] So, [ y_h(t)=C_1\cos 2t+C_2\sin 2t. ]

Now find a particular solution.

For (t^2), take [ y_{p1}=At^2+Bt+C. ] Then [ y_{p1}''=2A. ] Substitute into (y''+4y=t^2): [ 2A+4(At^2+Bt+C)=t^2. ] Comparing coefficients: [ 4A=1 \implies A=\frac14, ] [ 4B=0 \implies B=0, ] [ 2A+4C=0 \implies \frac12+4C=0 \implies C=-\frac18. ] Thus, [ y_{p1}=\frac{t^2}{4}-\frac18. ]

For (3e^t), take [ y_{p2}=De^t. ] Then [ y_{p2}''+4y_{p2}=De^t+4De^t=5De^t. ] So, [ 5D=3 \implies D=\frac35. ] Hence, [ y_{p2}=\frac35e^t. ]

Therefore the general solution is [ y(t)=C_1\cos 2t+C_2\sin 2t+\frac{t^2}{4}-\frac18+\frac35e^t. ]

Now apply the initial conditions.

[ y(0)=0 ] gives [ C_1-\frac18+\frac35=0. ] So, [ C_1=\frac18-\frac35=-\frac{19}{40}. ]

Also, [ y'(t)=-2C_1\sin 2t+2C_2\cos 2t+\frac{t}{2}+\frac35e^t. ] Using (y'(0)=2): [ 2C_2+\frac35=2, ] so [ 2C_2=\frac75, ] [ C_2=\frac7{10}. ]

Hence the solution of the initial value problem is [ \boxed{ y(t)=-\frac{19}{40}\cos 2t+\frac7{10}\sin 2t+\frac{t^2}{4}-\frac18+\frac35e^t } ]

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2082 · Differential Equation (MAT 202) (New Course)

Group "B" 5x7=35
(a) How do you solve the homogenous differential equation of the form f(dydx,yx)=0f\left(\frac{dy}{dx}, \frac{y}{x}\right) = 0 ? [3] (b) Find complete and singular solution of the equation y=px+ap(1−p).y = px + ap(1 - p). [4]

Answer

(a) Solving a homogeneous differential equation of the form

[ f\left(\frac{dy}{dx},\frac{y}{x}\right)=0 ]

A first-order homogeneous differential equation can be written as
[ \frac{dy}{dx}=F\left(\frac{y}{x}\right) ] or, more generally, as a relation between (\frac{dy}{dx}) and (\frac{y}{x}).

Method:

  1. Put
    [ y=vx ] where (v) is a new variable.

  2. Differentiate with respect to (x): [ \frac{dy}{dx}=v+x\frac{dv}{dx}. ]

  3. Substitute (y=vx) and (\frac{dy}{dx}=v+x\frac{dv}{dx}) into the given equation. It becomes an equation involving (v) and (x) only: [ f\left(v+x\frac{dv}{dx},,v\right)=0. ]

  4. Solve this equation for (\frac{dv}{dx}). Usually it takes the separable form
    [ \frac{dv}{dx}=\phi(v). ]

  5. Separate the variables: [ \frac{dv}{\phi(v)}=\frac{dx}{x}. ]

  6. Integrate both sides: [ \int \frac{dv}{\phi(v)}=\int \frac{dx}{x}+C. ]

  7. Finally, replace (v) by (\frac{y}{x}) to get the required general solution.


(b) Complete and singular solution of

[ y=px+ap(1-p) ] where
[ p=\frac{dy}{dx}. ]

This is a Clairaut’s equation of the form
[ y=px+f(p), ] with
[ f(p)=ap(1-p)=ap-ap^2. ]

Complete solution

For Clairaut’s equation, the complete solution is obtained by replacing (p) by an arbitrary constant (c). Hence
[ \boxed{y=cx+ac(1-c)} ] or
[ \boxed{y=cx+ac-ac^2}. ]

Singular solution

Differentiate the given equation with respect to (p):
[ \frac{dy}{dp}=x+f'(p). ] But for the singular solution, we use
[ x+f'(p)=0. ]

Here,
[ f(p)=ap-ap^2 ] so
[ f'(p)=a-2ap=a(1-2p). ]

Thus,
[ x+a(1-2p)=0 ] [ x=-a(1-2p) ] [ x=a(2p-1). ]

Hence,
[ 2p-1=\frac{x}{a} ] [ p=\frac{x+a}{2a}. ]

Now substitute this value of (p) into the original equation
[ y=px+ap(1-p). ]

We have
[ p=\frac{x+a}{2a}. ]

Then
[ px=\frac{x(x+a)}{2a} ] and
[ 1-p=1-\frac{x+a}{2a}=\frac{a-x}{2a}. ]

Therefore,
[ ap(1-p)=a\cdot \frac{x+a}{2a}\cdot \frac{a-x}{2a} =\frac{(x+a)(a-x)}{4a} =\frac{a^2-x^2}{4a}. ]

So,
[ y=\frac{x(x+a)}{2a}+\frac{a^2-x^2}{4a}. ]

Taking L.C.M. (4a),
[ y=\frac{2x(x+a)+a^2-x^2}{4a} ] [ y=\frac{2x^2+2ax+a^2-x^2}{4a} ] [ y=\frac{x^2+2ax+a^2}{4a} ] [ \boxed{y=\frac{(x+a)^2}{4a}}. ]

Hence, the singular solution is
[ \boxed{y=\frac{(x+a)^2}{4a}}. ]


Final answers:

[ \boxed{\text{Complete solution: } y=cx+ac(1-c)} ]

[ \boxed{\text{Singular solution: } y=\frac{(x+a)^2}{4a}} ]

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2082 · Differential Equation (MAT 202) (New Course) · 7 marks

Find the equation of the integral surface satisfying 4yzp+4+2y=04yzp + 4 + 2y = 0 and passing through the circle y2+z2=1,x+z=2y^2 + z^2 = 1, x + z = 2.
P.T.O.

Answer

Write the equation as F(x,y,z,p,q)=4yzp+4+2y=0F(x,y,z,p,q) = 4yzp + 4 + 2y = 0, where p=∂z/∂xp = \partial z/\partial x and q=∂z/∂yq = \partial z/\partial y. It contains no qq, so Fq=0F_q = 0.

Characteristic equations (Lagrange–Charpit)

dxFp=dyFq=dzpFp+qFq=dx4yz=dy0=dz4yzp\frac{dx}{F_p} = \frac{dy}{F_q} = \frac{dz}{pF_p + qF_q} = \frac{dx}{4yz} = \frac{dy}{0} = \frac{dz}{4yzp}

From dy/0dy/0 we get dy=0dy = 0, so yy = constant along each characteristic. Then dz=p dxdz = p\,dx, and the equation itself gives

zp=−4+2y4y=−2+y2y.zp = -\frac{4+2y}{4y} = -\frac{2+y}{2y}.

Integration

With yy held constant, zp=12 ∂(z2)∂xzp = \tfrac12\,\dfrac{\partial (z^2)}{\partial x}, so

∂(z2)∂x=−2+yy⟹z2=−(2+y)y x+h(y)\frac{\partial (z^2)}{\partial x} = -\frac{2+y}{y} \quad\Longrightarrow\quad z^2 = -\frac{(2+y)}{y}\,x + h(y)

where h(y)h(y) is arbitrary (a complete integral in terms of the characteristic parameter). This is the general solution.

Applying the initial curve

The curve is y2+z2=1y^2 + z^2 = 1, x+z=2x + z = 2. Put z=w=1−y2z = w = \sqrt{1-y^2} (taking the branch of the circle used) and x=2−wx = 2 - w on it:

1−y2=−(2+y)(2−w)y+h(y)  ⟹  h(y)=1−y2+(2+y) (2−w)y.1 - y^2 = -\frac{(2+y)(2-w)}{y} + h(y) \;\Longrightarrow\; h(y) = 1 - y^2 + \frac{(2+y)\,(2-w)}{y}.

Substituting back:

z2=1−y2+(2+y)(2−x−1−y2)yz^2 = 1 - y^2 + \frac{(2+y)\left(2 - x - \sqrt{1-y^2}\right)}{y}

or

 yz2=y(1−y2)+(2+y)(2−x−1−y2) \boxed{\,y z^2 = y(1-y^2) + (2+y)\left(2 - x - \sqrt{1-y^2}\right)\,}

Check: on the curve x=2−zx = 2 - z with z=1−y2z = \sqrt{1-y^2}, the bracket 2−x−1−y2=z−z=02 - x - \sqrt{1-y^2} = z - z = 0, so yz2=y(1−y2)yz^2 = y(1-y^2), i.e. z2=1−y2z^2 = 1 - y^2, as required. This is the equation of the required integral surface.

Answer checked by Claude Sonnet · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2082 · Differential Equation (MAT 202) (New Course) · 7 marks

Reduce the equation ∂2z∂x2+2∂2z∂x∂y+∂2z∂y2=0\frac{\partial^2 z}{\partial x^2} + 2\frac{\partial^2 z}{\partial x \partial y} + \frac{\partial^2 z}{\partial y^2} = 0 to canonical form and solve it.

Answer

The given PDE is

∂2z∂x2+2∂2z∂x∂y+∂2z∂y2=0.\frac{\partial^2 z}{\partial x^2}+2\frac{\partial^2 z}{\partial x\partial y}+\frac{\partial^2 z}{\partial y^2}=0.

Comparing with

Azxx+2Bzxy+Czyy=0,A z_{xx}+2B z_{xy}+C z_{yy}=0,

we get

A=1,B=1,C=1.A=1,\quad B=1,\quad C=1.

Hence

B2−AC=1−1=0,B^2-AC=1-1=0,

so the equation is parabolic.

The characteristic equation is

A(dydx)2−2B(dydx)+C=0,A\left(\frac{dy}{dx}\right)^2-2B\left(\frac{dy}{dx}\right)+C=0,

i.e.

(dydx)2−2(dydx)+1=0.\left(\frac{dy}{dx}\right)^2-2\left(\frac{dy}{dx}\right)+1=0.

Thus

(dydx−1)2=0⇒dydx=1.\left(\frac{dy}{dx}-1\right)^2=0 \quad\Rightarrow\quad \frac{dy}{dx}=1.

So the characteristic curves are

y−x=c.y-x=c.

Let

ξ=y−x,η=x.\xi=y-x,\qquad \eta=x.

Then

zx=−zξ+zη,zy=zξ.z_x=-z_\xi+z_\eta,\qquad z_y=z_\xi.

Also,

zxx=zξξ−2zξη+zηη,z_{xx}=z_{\xi\xi}-2z_{\xi\eta}+z_{\eta\eta}, zxy=−zξξ+zξη,z_{xy}=-z_{\xi\xi}+z_{\xi\eta}, zyy=zξξ.z_{yy}=z_{\xi\xi}.

Substituting in the given PDE,

(zξξ−2zξη+zηη)+2(−zξξ+zξη)+zξξ=0.(z_{\xi\xi}-2z_{\xi\eta}+z_{\eta\eta}) +2(-z_{\xi\xi}+z_{\xi\eta})+z_{\xi\xi}=0.

The terms containing (z_{\xi\xi}) and (z_{\xi\eta}) cancel, leaving

zηη=0.z_{\eta\eta}=0.

This is the required canonical form.

Now integrate (z_{\eta\eta}=0) with respect to (\eta):

zη=ϕ(ξ),z_\eta=\phi(\xi),

where (\phi) is an arbitrary function of (\xi). Again integrating with respect to (\eta),

z=η ϕ(ξ)+ψ(ξ),z=\eta\,\phi(\xi)+\psi(\xi),

where (\psi) is another arbitrary function of (\xi).

Reverting to (x,y), since (\xi=y-x) and (\eta=x),

z=x ϕ(y−x)+ψ(y−x)\boxed{z=x\,\phi(y-x)+\psi(y-x)}

or equivalently,

z=y ϕ(y−x)+ψ(y−x).\boxed{z=y\,\phi(y-x)+\psi(y-x)}.

Hence the general solution is

z=xf(y−x)+g(y−x)\boxed{z=x f(y-x)+g(y-x)}

where (f) and (g) are arbitrary functions.

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2082 · Differential Equation (MAT 202) (New Course) · 7 marks

Sketch the graph of the given function for three periods and find the Fourier series for the given function. f(x)={x+1,−1≤x<01−x,0≤x<1f(x+2)=f(x)f(x) = \begin{cases} x + 1, & -1 \le x < 0 \\ 1 - x, & 0 \le x < 1 \end{cases} \quad f(x+2) = f(x)
OR
Solve r=a2tr = a^2 t by Monge's Method

Answer

Sketch

The function is f(x)=1−∣x∣f(x) = 1 - |x| on [−1,1][-1,1], repeated with period 2. Its graph is a zig-zag (triangular wave):

  • peaks of height 1 at x=0,±2,±4,…x = 0, \pm 2, \pm 4, \dots,
  • zeros at x=±1,±3,…x = \pm1, \pm3, \dots,
  • straight lines between them (slope +1+1 on [−1,0][-1,0] and −1-1 on [0,1][0,1]).

Over three periods, from x=−3x=-3 to x=3x=3, the graph rises from 0 at x=−3x=-3 to a peak of 1 at x=−2x=-2, falls to 0 at x=−1x=-1, rises to 1 at x=0x=0, falls to 0 at x=1x=1, rises to 1 at x=2x=2 and falls to 0 at x=3x=3.

Fourier series

Period 2L=22L = 2, so L=1L = 1. The function is even (f(−x)=f(x)f(-x) = f(x)), hence bn=0b_n = 0.

a0=2L∫01(1−x) dx=2[x−x22]01=1a_0 = \frac{2}{L}\int_0^1 (1-x)\,dx = 2\left[x - \frac{x^2}{2}\right]_0^1 = 1

an=2∫01(1−x)cos⁡nπx dx=2[(1−x)sin⁡nπxnπ]01+2nπ∫01sin⁡nπx dx=2 (1−cos⁡nπ)n2π2a_n = 2\int_0^1 (1-x)\cos n\pi x\,dx = 2\left[\frac{(1-x)\sin n\pi x}{n\pi}\right]_0^1 + \frac{2}{n\pi}\int_0^1 \sin n\pi x\,dx = \frac{2\,(1-\cos n\pi)}{n^2\pi^2}

So an=4n2π2a_n = \dfrac{4}{n^2\pi^2} for odd nn and an=0a_n = 0 for even nn.

f(x)=a02+∑ancos⁡nπx=12+4π2[cos⁡πx+cos⁡3πx32+cos⁡5πx52+⋯ ]f(x) = \frac{a_0}{2} + \sum a_n\cos n\pi x = \frac12 + \frac{4}{\pi^2}\left[\cos\pi x + \frac{\cos 3\pi x}{3^2} + \frac{\cos 5\pi x}{5^2} + \cdots\right]

Check at x=0x=0: 12+4π2⋅π28=1=f(0)\tfrac12 + \tfrac{4}{\pi^2}\cdot\tfrac{\pi^2}{8} = 1 = f(0).


OR: Solve r=a2tr = a^2 t by Monge's method

Here R=1, S=0, T=−a2R = 1,\ S = 0,\ T = -a^2 (with r=zxxr = z_{xx}, t=zyyt = z_{yy}). Monge's equations are

R dy2−S dx dy+T dx2=0andR dp dy+T dq dx=0.R\,dy^2 - S\,dx\,dy + T\,dx^2 = 0 \quad\text{and}\quad R\,dp\,dy + T\,dq\,dx = 0.

The first gives dy2−a2dx2=0dy^2 - a^2dx^2 = 0, i.e. dy=±a dxdy = \pm a\,dx, so y∓ax=y \mp ax = constant.

  • For dy=a dxdy = a\,dx: the second equation becomes a dp−a2 dq=0a\,dp - a^2\,dq = 0, so dp=a dqdp = a\,dq and p−aq=F(y−ax)p - aq = F(y - ax).
  • For dy=−a dxdy = -a\,dx: −a dp−a2 dq=0-a\,dp - a^2\,dq = 0, so p+aq=G(y+ax)p + aq = G(y + ax).

Solving for pp and qq gives p=12(F+G)p = \tfrac12(F+G) and q=12a(G−F)q = \tfrac1{2a}(G-F). Integrating dz=p dx+q dydz = p\,dx + q\,dy:

z=ϕ(y+ax)+ψ(y−ax)\boxed{z = \phi(y + ax) + \psi(y - ax)}

where ϕ\phi and ψ\psi are arbitrary functions. Check: zxx=a2(ϕ′′+ψ′′)z_{xx} = a^2(\phi'' + \psi'') and zyy=ϕ′′+ψ′′z_{yy} = \phi'' + \psi'', so r=a2tr = a^2 t.

Answer checked by Claude Sonnet · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2082 · Differential Equation (MAT 202) (New Course) · 7 marks

Find the solution of two dimensional Laplace's equation uxx+uyy=0u_{xx} + u_{yy} = 0 in a rectangle 0<x<a,0<y<b0 < x < a, 0 < y < b satisfying the following Dirichlet boundary conditions. u(x,0)=u(x,b)=0, 0<x<au(x, 0) = u(x, b) = 0,\ 0 < x < a u(0,y)=0,u(a,y)=f(y), 0≤y≤bu(0, y) = 0, u(a, y) = f(y),\ 0 \le y \le b Where f is a given function on 0≤y≤b0 \le y \le b
OR
Using the non-dimensional variable S=xLS = \frac{x}{L}, prove that the wave equation utt=a2uxxu_{tt} = a^2 u_{xx} can be changed into L2uss=a2uttL^2 u_{ss} = a^2 u_{tt}, so that it is non dimensionalized with respect to the spatial variable. [7]

Answer

Main part: Two-dimensional Laplace equation in a rectangle

We have to solve

uxx+uyy=0,0<x<a,0<y<bu_{xx}+u_{yy}=0,\qquad 0<x<a,\quad 0<y<b

with boundary conditions

u(x,0)=0,u(x,b)=0,0<x<a,u(x,0)=0,\qquad u(x,b)=0,\qquad 0<x<a, u(0,y)=0,u(a,y)=f(y),0≤y≤b.u(0,y)=0,\qquad u(a,y)=f(y),\qquad 0\le y\le b.

Use the method of separation of variables. Let

u(x,y)=X(x)Y(y).u(x,y)=X(x)Y(y).

Substituting into Laplace’s equation,

X′′Y+XY′′=0X''Y+XY''=0

so that

X′′X=−Y′′Y=−λ.\frac{X''}{X}=-\frac{Y''}{Y}=-\lambda.

Thus,

X′′+λX=0,Y′′−λY=0.X''+\lambda X=0,\qquad Y''-\lambda Y=0.

The boundary conditions in the (y)-direction are homogeneous:

Y(0)=0,Y(b)=0.Y(0)=0,\qquad Y(b)=0.

For nontrivial solutions we must take (\lambda<0). Put

λ=−k2.\lambda=-k^2.

Then

Y′′+k2Y=0.Y''+k^2Y=0.

Its general solution is

Y(y)=Acos⁡ky+Bsin⁡ky.Y(y)=A\cos ky+B\sin ky.

Using (Y(0)=0),

A=0.A=0.

Hence

Y(y)=Bsin⁡ky.Y(y)=B\sin ky.

Using (Y(b)=0),

Bsin⁡kb=0.B\sin kb=0.

For a nontrivial solution, (B\ne0), so

sin⁡kb=0.\sin kb=0.

Therefore,

kb=nπ,n=1,2,3,…kb=n\pi,\qquad n=1,2,3,\dots

so

k=nπb.k=\frac{n\pi}{b}.

Thus the (y)-eigenfunctions are

Yn(y)=sin⁡nπyb.Y_n(y)=\sin\frac{n\pi y}{b}.

Now the (x)-equation becomes

X′′−n2π2b2X=0.X''-\frac{n^2\pi^2}{b^2}X=0.

Its general solution is

X(x)=C1cosh⁡nπxb+C2sinh⁡nπxb.X(x)=C_1\cosh\frac{n\pi x}{b}+C_2\sinh\frac{n\pi x}{b}.

Using the boundary condition (u(0,y)=0), we get (X(0)=0). Hence

C1=0.C_1=0.

Therefore,

Xn(x)=C2sinh⁡nπxb.X_n(x)=C_2\sinh\frac{n\pi x}{b}.

So the elementary solutions are

un(x,y)=Ansinh⁡nπxbsin⁡nπyb.u_n(x,y)=A_n\sinh\frac{n\pi x}{b}\sin\frac{n\pi y}{b}.

By the principle of superposition, the required solution is

u(x,y)=∑n=1∞Ansinh⁡nπxbsin⁡nπyb.u(x,y)=\sum_{n=1}^{\infty}A_n\sinh\frac{n\pi x}{b}\sin\frac{n\pi y}{b}.

Now apply the non-homogeneous boundary condition

u(a,y)=f(y).u(a,y)=f(y).

Then

f(y)=∑n=1∞Ansinh⁡nπabsin⁡nπyb.f(y)=\sum_{n=1}^{\infty}A_n\sinh\frac{n\pi a}{b}\sin\frac{n\pi y}{b}.

This is the Fourier sine series of (f(y)) on (0<y<b). Hence

Ansinh⁡nπab=2b∫0bf(y)sin⁡nπyb dy.A_n\sinh\frac{n\pi a}{b} = \frac{2}{b}\int_0^b f(y)\sin\frac{n\pi y}{b}\,dy.

Therefore,

An=2bsinh⁡nπab∫0bf(y)sin⁡nπyb dy.A_n= \frac{2}{b\sinh\frac{n\pi a}{b}} \int_0^b f(y)\sin\frac{n\pi y}{b}\,dy.

Hence the solution of the given boundary value problem is

u(x,y)=∑n=1∞[2bsinh⁡nπab∫0bf(η)sin⁡nπηb dη]sinh⁡nπxbsin⁡nπyb\boxed{ u(x,y)= \sum_{n=1}^{\infty} \left[ \frac{2}{b\sinh\frac{n\pi a}{b}} \int_0^b f(\eta)\sin\frac{n\pi \eta}{b}\,d\eta \right] \sinh\frac{n\pi x}{b} \sin\frac{n\pi y}{b} }

where (\eta) is a dummy variable. This satisfies Laplace’s equation and all the given boundary conditions. For a classical solution, (f) should be piecewise smooth and satisfy the compatibility conditions (f(0)=f(b)=0).


OR part: Non-dimensionalization of the wave equation

Let

S=xL.S=\frac{x}{L}.

Then

x=LS.x=LS.

Write

u(x,t)=U(S,t).u(x,t)=U(S,t).

By the chain rule,

ux=∂U∂SdSdx=1LUS.u_x=\frac{\partial U}{\partial S}\frac{dS}{dx} =\frac{1}{L}U_S.

Differentiating again,

uxx=1L2USS.u_{xx}=\frac{1}{L^2}U_{SS}.

The wave equation is

utt=a2uxx.u_{tt}=a^2u_{xx}.

Substituting (u_{xx}=\dfrac{1}{L^2}U_{SS}),

Utt=a2L2USS.U_{tt}=\frac{a^2}{L^2}U_{SS}.

Multiplying by (L^2),

L2Utt=a2USS.L^2U_{tt}=a^2U_{SS}.

Rename (U) again as (u). Then

L2utt=a2uSS\boxed{L^2u_{tt}=a^2u_{SS}}

or, writing (S=s),

L2utt=a2uss.\boxed{L^2u_{tt}=a^2u_{ss}}.

Thus the wave equation becomes non-dimensional with respect to the spatial variable (S=x/L).
If the printed statement writes (L^2u_{ss}=a^2u_{tt}), then the subscripts (s) and (t) are interchanged; the correct transformed form is (L^2u_{tt}=a^2u_{ss}).

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2081 · Differential Equation (MAT 202) · 8 marks

Find the general solution of the given initial value problem. y′′+4y=2sin⁡t,y(0)=2,y′(0)=−1.y'' + 4y = 2\sin t, y(0) = 2, y'(0) = -1.
OR
A spring is stretched by 10 cm by a force of 2 N. A mass of 2 kg is hung from the spring and is also attached to a viscous damper that exerts a force of 4 N when the velocity is 1 m/s. The motion of the mass is driven by an external force of 4cos⁡2t4\cos 2t N. If the mass is initially at rest at equilibrium, find its position at any time t; identify the transient and steady state parts of the solution.

Answer

y′′+4y=2sin⁡t,y(0)=2,y′(0)=−1.y''+4y=2\sin t,\qquad y(0)=2,\quad y'(0)=-1.

Homogeneous equation: y′′+4y=0y''+4y=0 Auxiliary equation: r2+4=0  ⟹  r=±2ir^2+4=0 \implies r=\pm 2i So, yh=c1cos⁡2t+c2sin⁡2t.y_h=c_1\cos 2t+c_2\sin 2t.

Particular solution: Since the forcing term is (2\sin t), try yp=Asin⁡t+Bcos⁡t.y_p=A\sin t+B\cos t. Then yp′′=−Asin⁡t−Bcos⁡t.y_p''=-A\sin t-B\cos t. Substitute into (y''+4y=2\sin t): (−Asin⁡t−Bcos⁡t)+4(Asin⁡t+Bcos⁡t)=2sin⁡t(-A\sin t-B\cos t)+4(A\sin t+B\cos t)=2\sin t 3Asin⁡t+3Bcos⁡t=2sin⁡t.3A\sin t+3B\cos t=2\sin t. Comparing coefficients: 3A=2  ⟹  A=23,3B=0  ⟹  B=0.3A=2 \implies A=\frac23,\qquad 3B=0 \implies B=0. Hence, yp=23sin⁡t.y_p=\frac23\sin t.

General solution: y(t)=c1cos⁡2t+c2sin⁡2t+23sin⁡t.y(t)=c_1\cos 2t+c_2\sin 2t+\frac23\sin t.

Now apply the initial conditions.

y(0)=c1=2.y(0)=c_1=2.

Also, y′(t)=−2c1sin⁡2t+2c2cos⁡2t+23cos⁡t.y'(t)=-2c_1\sin 2t+2c_2\cos 2t+\frac23\cos t. Using (y'(0)=-1): 2c2+23=−12c_2+\frac23=-1 2c2=−532c_2=-\frac53 c2=−56.c_2=-\frac56.

Therefore the solution of the initial value problem is y(t)=2cos⁡2t−56sin⁡2t+23sin⁡t.\boxed{y(t)=2\cos 2t-\frac56\sin 2t+\frac23\sin t}.

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2081 · Differential Equation (MAT 202) · 8 marks

Find the general solution of the given IVP y′′+y′−2y=2t,y(0)=0,y′(0)=1y'' + y' - 2y = 2t, y(0) = 0, y'(0) = 1

Answer

The differential equation is

y′′+y′−2y=2t,y(0)=0,y′(0)=1.y''+y'-2y=2t,\qquad y(0)=0,\quad y'(0)=1.

Homogeneous solution

Auxiliary equation:

m2+m−2=0m^2+m-2=0

(m+2)(m−1)=0(m+2)(m-1)=0

m=1, −2.m=1,\ -2.

So,

yc=C1et+C2e−2t.y_c=C_1e^t+C_2e^{-2t}.

Particular solution

Since the right-hand side is (2t), take

yp=At+B.y_p=At+B.

Then

yp′=A,yp′′=0.y_p'=A,\qquad y_p''=0.

Substitute into the differential equation:

0+A−2(At+B)=2t0+A-2(At+B)=2t

−2At+(A−2B)=2t.-2At+(A-2B)=2t.

Comparing coefficients:

−2A=2  ⟹  A=−1,-2A=2 \implies A=-1,

A−2B=0  ⟹  −1−2B=0  ⟹  B=−12.A-2B=0 \implies -1-2B=0 \implies B=-\frac12.

Hence,

yp=−t−12.y_p=-t-\frac12.

Therefore the general solution is

y=C1et+C2e−2t−t−12.y=C_1e^t+C_2e^{-2t}-t-\frac12.

Using initial conditions

y(0)=C1+C2−12=0y(0)=C_1+C_2-\frac12=0

C1+C2=12.(1)C_1+C_2=\frac12. \quad (1)

Also,

y′=C1et−2C2e−2t−1.y'=C_1e^t-2C_2e^{-2t}-1.

y′(0)=C1−2C2−1=1y'(0)=C_1-2C_2-1=1

C1−2C2=2.(2)C_1-2C_2=2. \quad (2)

Subtracting (1) from (2):

3C2=−32  ⟹  C2=−12.3C_2=-\frac32 \implies C_2=-\frac12.

Then from (1):

C1−12=12  ⟹  C1=1.C_1-\frac12=\frac12 \implies C_1=1.

Final answer

y=et−12e−2t−t−12\boxed{y=e^t-\frac12e^{-2t}-t-\frac12}

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2081 · Differential Equation (MAT 202) · 3 marks

(a) How do you solve the Clairaut's differential equation y=px+f(p)y = px + f(p) where p=dydxp = \frac{dy}{dx}

Answer

For Clairaut's equation
y=px+f(p),p=dydx,y = px + f(p), \qquad p=\frac{dy}{dx}, differentiate both sides with respect to xx: dydx=p+xdpdx+f′(p)dpdx.\frac{dy}{dx}=p+x\frac{dp}{dx}+f'(p)\frac{dp}{dx}. Since dydx=p\frac{dy}{dx}=p, this gives p=p+[x+f′(p)]dpdxp=p+\bigl[x+f'(p)\bigr]\frac{dp}{dx} so that [x+f′(p)]dpdx=0.\bigl[x+f'(p)\bigr]\frac{dp}{dx}=0.

Hence two cases arise:

  1. General solution:
    dpdx=0  ⟹  p=c(arbitrary constant).\frac{dp}{dx}=0 \implies p=c \quad (\text{arbitrary constant}). Substituting p=cp=c in the original equation, y=cx+f(c).y=cx+f(c). This is the required general solution.

  2. Singular solution:
    x+f′(p)=0.x+f'(p)=0. Eliminate pp between this equation and the original equation y=px+f(p)y=px+f(p). The resulting relation between xx and yy, free from the arbitrary constant cc, is the singular solution.

Thus the complete solution consists of the general solution y=cx+f(c)y=cx+f(c) and, if it exists, the singular solution obtained by eliminating pp from x+f′(p)=0,y=px+f(p).x+f'(p)=0,\qquad y=px+f(p).

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2081 · Differential Equation (MAT 202) · 7 marks

(b) Find the complete and singular solution of the equation y=px+ap(1−p)y = px + ap(1-p)

Answer

Let
[ p=\frac{dy}{dx}. ] The given equation is [ y=px+ap(1-p). ] This is a Clairaut-type equation of the form [ y=px+f(p),\qquad f(p)=ap(1-p). ]

Complete solution

For Clairaut’s equation, replace (p) by an arbitrary constant (c). Hence the complete solution is [ \boxed{y=cx+ac(1-c)}. ]

Singular solution

Write [ F(x,y,p)=y-px-ap(1-p)=0. ] For the singular solution, eliminate (p) using [ \frac{\partial F}{\partial p}=0. ] Now, [ \frac{\partial F}{\partial p}=-x-a(1-2p)=0. ] So [ x=a(2p-1) ] or [ p=\frac{x+a}{2a}. ] Substitute this value of (p) into the complete solution: [ y=cx+ac(1-c), ] where (c=p=\dfrac{x+a}{2a}). Thus [ y=\frac{x+a}{2a}x+a\frac{x+a}{2a}\left(1-\frac{x+a}{2a}\right). ] Simplifying, [ y=\frac{(x+a)^2}{4a}. ]

Therefore the singular solution is [ \boxed{y=\frac{(x+a)^2}{4a}},\qquad a\ne 0. ]

If (a=0), the equation reduces to (y=px), whose complete solution is (y=cx), and there is no separate singular solution.

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2081 · Differential Equation (MAT 202) · 7 marks

Find the surface which is orthogonal to the one-parameter system. Z=cxy(x2+y2)Z = cxy(x^2 + y^2) and which passes through the hyperbola x2−y2=a2,Z=0.x^2 - y^2 = a^2, Z = 0.

Answer

Step 1: Differential equation of the given system. The one-parameter family is

z=c xy(x2+y2)  ⇒  c=zxy(x2+y2)z = c\,xy(x^2+y^2) \;\Rightarrow\; c = \frac{z}{xy(x^2+y^2)}

Differentiating z=c(x3y+xy3)z = c(x^3y + xy^3):

p1=c(3x2y+y3)=z(3x2+y2)x(x2+y2),q1=c(x3+3xy2)=z(x2+3y2)y(x2+y2)p_1 = c(3x^2y + y^3) = \frac{z(3x^2+y^2)}{x(x^2+y^2)},\qquad q_1 = c(x^3 + 3xy^2) = \frac{z(x^2+3y^2)}{y(x^2+y^2)}

The normal to the family is (p1,q1,−1)(p_1, q_1, -1).

Step 2: Orthogonal surfaces. A surface z=f(x,y)z=f(x,y) is orthogonal to the family if its normal (p,q,−1)(p,q,-1) is perpendicular to (p1,q1,−1)(p_1,q_1,-1):

p1p+q1q+1=0 ⇒ p1p+q1q=−1p_1 p + q_1 q + 1 = 0\ \Rightarrow\ p_1 p + q_1 q = -1

This is a Lagrange equation, with auxiliary equations

dxp1=dyq1=dz−1\frac{dx}{p_1} = \frac{dy}{q_1} = \frac{dz}{-1}

Step 3: First integral. From dxp1=dyq1\dfrac{dx}{p_1} = \dfrac{dy}{q_1}:

x dx3x2+y2=y dyx2+3y2 ⇒ (x3+3xy2) dx=(3x2y+y3) dy\frac{x\,dx}{3x^2+y^2} = \frac{y\,dy}{x^2+3y^2}\ \Rightarrow\ (x^3+3xy^2)\,dx = (3x^2y+y^3)\,dy

Put u=x2, v=y2u=x^2,\ v=y^2: (u+3v) du=(3u+v) dv(u+3v)\,du = (3u+v)\,dv. With s=u+vs=u+v and d=u−vd=u-v, this is 2s dd=d ds2s\,dd = d\,ds, giving d2=c1sd^2 = c_1 s:

(x2−y2)2x2+y2=c1\frac{(x^2-y^2)^2}{x^2+y^2} = c_1

Step 4: Second integral. Note that xp1+yq1=z[(3x2+y2)+(x2+3y2)]x2+y2=4zx p_1 + y q_1 = \dfrac{z[(3x^2+y^2)+(x^2+3y^2)]}{x^2+y^2} = 4z. Using multipliers x,yx, y:

x dx+y dy4z=dz−1 ⇒ x dx+y dy+4z dz=0 ⇒ x2+y2+4z2=c2\frac{x\,dx + y\,dy}{4z} = \frac{dz}{-1}\ \Rightarrow\ x\,dx+y\,dy+4z\,dz = 0\ \Rightarrow\ x^2+y^2+4z^2 = c_2

Step 5: Apply the condition. The surface must pass through x2−y2=a2, z=0x^2-y^2=a^2,\ z=0. On this curve, with x2+y2=sx^2+y^2=s:

c1=a4s,c2=s⇒c1c2=a4c_1 = \frac{a^4}{s},\qquad c_2 = s \quad\Rightarrow\quad c_1c_2 = a^4

Hence the required orthogonal surface is

(x2−y2)2 (x2+y2+4z2)=a4 (x2+y2)\boxed{(x^2-y^2)^2\,(x^2+y^2+4z^2) = a^4\,(x^2+y^2)}

Answer checked by Claude Sonnet · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2081 · Differential Equation (MAT 202) · 7 marks

Solve p2x+q2y=zp^2x + q^2y = z using Charpit's method.

Answer

Let F(x,y,z,p,q)=p2x+q2y−z=0F(x,y,z,p,q) = p^2x + q^2y - z = 0, where p=∂z/∂xp=\partial z/\partial x, q=∂z/∂yq=\partial z/\partial y.

Charpit's auxiliary equations:

dpFx+pFz=dqFy+qFz=dz−pFp−qFq=dx−Fp=dy−Fq\frac{dp}{F_x+pF_z} = \frac{dq}{F_y+qF_z} = \frac{dz}{-pF_p-qF_q} = \frac{dx}{-F_p} = \frac{dy}{-F_q}

Here Fx=p2, Fy=q2, Fz=−1, Fp=2px, Fq=2qyF_x = p^2,\ F_y = q^2,\ F_z = -1,\ F_p = 2px,\ F_q = 2qy, so

dpp2−p=dqq2−q=dz−2(p2x+q2y)=dz−2z=dx−2px=dy−2qy\frac{dp}{p^2-p} = \frac{dq}{q^2-q} = \frac{dz}{-2(p^2x+q^2y)} = \frac{dz}{-2z} = \frac{dx}{-2px} = \frac{dy}{-2qy}

Finding an integral. Take dpp(p−1)=dx−2px\dfrac{dp}{p(p-1)} = \dfrac{dx}{-2px}, which gives dpp−1=−dx2x\dfrac{dp}{p-1} = -\dfrac{dx}{2x}, and take dpp(p−1)=dz−2z\dfrac{dp}{p(p-1)} = \dfrac{dz}{-2z}. Using these two,

2 dpp+dxx−dzz=0 ⇒ p2x=a2z(a constant)\frac{2\,dp}{p} + \frac{dx}{x} - \frac{dz}{z} = 0\ \Rightarrow\ p^2x = a^2 z \qquad (a\ \text{constant})

Verification: along the auxiliary equations dp:dx:dz=(p2−p):(−2px):(−2z)dp:dx:dz = (p^2-p) : (-2px) : (-2z), so 2(p−1)−2p+2=02(p-1) - 2p + 2 = 0, as required.

Finding pp and qq.

p=azxp = a\sqrt{\frac{z}{x}}

Substituting in p2x+q2y=zp^2x + q^2y = z: a2z+q2y=za^2z + q^2y = z, so

q=1−a2 zyq = \sqrt{1-a^2}\,\sqrt{\frac{z}{y}}

Solving dz=p dx+q dydz = p\,dx + q\,dy:

dzz=a dxx+1−a2 dyy\frac{dz}{\sqrt z} = a\,\frac{dx}{\sqrt x} + \sqrt{1-a^2}\,\frac{dy}{\sqrt y}

Integrating,

2z=2ax+21−a2y+2b2\sqrt z = 2a\sqrt x + 2\sqrt{1-a^2}\sqrt y + 2b

The complete integral is therefore

z=ax+1−a2 y+borz=(ax+1−a2y+b)2\boxed{\sqrt z = a\sqrt x + \sqrt{1-a^2}\,\sqrt y + b}\quad\text{or}\quad z = \big(a\sqrt x + \sqrt{1-a^2}\sqrt y + b\big)^2

where aa and bb are arbitrary constants.

Answer checked by Claude Sonnet · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2081 · Differential Equation (MAT 202) · 7 marks

Sketch the graph of the function for three periods. Find the Fourier series for the given function f(x)={x+L,−L≤x≤0,L,0<x<L,f(x+2L)=f(x)f(x) = \begin{cases} x + L, & -L \le x \le 0, \\ L, & 0 < x < L, \end{cases} \quad f(x+2L) = f(x)
OR Find the eigenvalues and eigenfunctions of the given boundary value problem. Assume that all eigenvalues are real y′′+λy=0,y′(0)=0,y'' + \lambda y = 0, y'(0) = 0, y(π)=0.y(\pi) = 0.

Answer

Fourier series of f(x)={x+L,−L≤x≤0L,0<x<Lf(x) = \begin{cases} x+L, & -L\le x\le 0 \\ L, & 0<x<L\end{cases},  f(x+2L)=f(x)\ f(x+2L)=f(x).

Graph (three periods). On [−L,0][-L,0] the graph is a straight line of slope 1 rising from 00 at x=−Lx=-L to LL at x=0x=0. On (0,L)(0,L) it is the horizontal line y=Ly=L. At x=Lx=L the function drops from LL back to 00 (a jump discontinuity), and the same pattern repeats on every interval of length 2L2L, for example on [−3L,−L][-3L,-L], [−L,L][-L,L] and [L,3L][L,3L]. The curve is continuous at x=0x=0, ±2L\pm2L, ... and has jumps at x=±L,±3L,…x=\pm L,\pm3L,\dots

Coefficients.

a0=1L[∫−L0(x+L) dx+∫0LL dx]=1L[L22+L2]=3L2a_0 = \frac1L\left[\int_{-L}^{0}(x+L)\,dx + \int_0^L L\,dx\right] = \frac1L\left[\frac{L^2}{2}+L^2\right] = \frac{3L}{2}

an=1L[∫−L0(x+L)cos⁡nπxL dx+∫0LLcos⁡nπxL dx]a_n = \frac1L\left[\int_{-L}^{0}(x+L)\cos\frac{n\pi x}{L}\,dx + \int_0^L L\cos\frac{n\pi x}{L}\,dx\right]

The second integral is zero. In the first, put u=x+Lu=x+L, so cos⁡nπxL=(−1)ncos⁡nπuL\cos\frac{n\pi x}{L} = (-1)^n\cos\frac{n\pi u}{L}:

an=(−1)nL∫0Lucos⁡nπuL du=(−1)nL⋅L2n2π2[(−1)n−1]=L [1−(−1)n]n2π2a_n = \frac{(-1)^n}{L}\int_0^L u\cos\frac{n\pi u}{L}\,du = \frac{(-1)^n}{L}\cdot\frac{L^2}{n^2\pi^2}\big[(-1)^n-1\big] = \frac{L\,[1-(-1)^n]}{n^2\pi^2}

so an=2Ln2π2a_n = \dfrac{2L}{n^2\pi^2} for odd nn and 00 for even nn.

bn=1L[∫−L0(x+L)sin⁡nπxL dx+L∫0Lsin⁡nπxL dx]b_n = \frac1L\left[\int_{-L}^{0}(x+L)\sin\frac{n\pi x}{L}\,dx + L\int_0^L\sin\frac{n\pi x}{L}\,dx\right]

∫−L0(x+L)sin⁡nπxL dx=(−1)n∫0Lusin⁡nπuL du=(−1)n⋅(−L2(−1)nnπ)=−L2nπ\int_{-L}^{0}(x+L)\sin\frac{n\pi x}{L}\,dx = (-1)^n\int_0^L u\sin\frac{n\pi u}{L}\,du = (-1)^n\cdot\left(-\frac{L^2(-1)^n}{n\pi}\right) = -\frac{L^2}{n\pi}

L∫0Lsin⁡nπxL dx=L2nπ[1−(−1)n]L\int_0^L\sin\frac{n\pi x}{L}\,dx = \frac{L^2}{n\pi}\big[1-(-1)^n\big]

bn=1L[−L2nπ+L2nπ−L2(−1)nnπ]=(−1)n+1Lnπb_n = \frac1L\left[-\frac{L^2}{n\pi}+\frac{L^2}{n\pi}-\frac{L^2(-1)^n}{n\pi}\right] = \frac{(-1)^{n+1}L}{n\pi}

Fourier series:

f(x)=3L4+∑n=1∞[L (1−(−1)n)n2π2cos⁡nπxL+(−1)n+1Lnπsin⁡nπxL]f(x) = \frac{3L}{4} + \sum_{n=1}^\infty\left[\frac{L\,(1-(-1)^n)}{n^2\pi^2}\cos\frac{n\pi x}{L} + \frac{(-1)^{n+1}L}{n\pi}\sin\frac{n\pi x}{L}\right]

Check at x=0x=0: 3L4+2Lπ2(1+19+125+⋯ )=3L4+2Lπ2⋅π28=L=f(0)\dfrac{3L}{4}+\dfrac{2L}{\pi^2}\left(1+\dfrac1{9}+\dfrac1{25}+\cdots\right) = \dfrac{3L}{4}+\dfrac{2L}{\pi^2}\cdot\dfrac{\pi^2}{8} = L = f(0).


OR: Eigenvalues and eigenfunctions of y′′+λy=0, y′(0)=0, y(π)=0y''+\lambda y=0,\ y'(0)=0,\ y(\pi)=0.

  • λ=0\lambda = 0: y=Ax+By = Ax+B. y′(0)=A=0y'(0)=A=0 and y(π)=B=0y(\pi)=B=0, so only the trivial solution.
  • λ<0\lambda<0 (let λ=−k2\lambda=-k^2): y=Acosh⁡kx+Bsinh⁡kxy = A\cosh kx + B\sinh kx. y′(0)=Bk=0⇒B=0y'(0)=Bk=0\Rightarrow B=0, and y(π)=Acosh⁡kπ=0⇒A=0y(\pi)=A\cosh k\pi=0\Rightarrow A=0: trivial.
  • λ>0\lambda>0 (let λ=k2\lambda=k^2, k>0k>0): y=Acos⁡kx+Bsin⁡kxy = A\cos kx + B\sin kx. y′(0)=Bk=0⇒B=0y'(0)=Bk=0\Rightarrow B=0. Then y(π)=Acos⁡kπ=0y(\pi)=A\cos k\pi=0 with A≠0A\ne0 requires kπ=(2n−1)π2k\pi = \dfrac{(2n-1)\pi}{2}.

Eigenvalues:

λn=(2n−12)2,n=1,2,3,…\lambda_n = \left(\frac{2n-1}{2}\right)^2,\quad n=1,2,3,\dots

Eigenfunctions:

yn(x)=cos⁡(2n−1)x2,n=1,2,3,…y_n(x) = \cos\frac{(2n-1)x}{2},\quad n=1,2,3,\dots

Answer checked by Claude Sonnet · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2081 · Differential Equation (MAT 202) · 7 marks

Find the solution u(x,y)u(x, y) of two dimensional Laplace's equation uxx+uyy=0u_{xx} + u_{yy} = 0 in the rectangle 0<x<a,0<y<b0 < x < a, 0 < y < b that satisfies the boundary conditions u(0,y)=0,u(a,y)=0,0<y<bu(0, y) = 0, u(a, y) = 0, 0 < y < b u(x,0)=h(x),u(x,b)=0,0≤x≤a.u(x, 0) = h(x), u(x, b) = 0, 0 \le x \le a.
OR Find the solution u(r,θ)u(r, \theta) of Laplace equation in semi-circular region r<a,0<θ<πr < a, 0 < \theta < \pi that satisfies the boundary conditions u(r,0)=0,u(r,π)=0,0≤r<au(r, 0) = 0, u(r, \pi) = 0, 0 \le r < a u(a,θ)=f(θ),0≤θ≤πu(a, \theta) = f(\theta), 0 \le \theta \le \pi by assuming that u is single valued and bounded in the region.

Answer

Solution of the first alternative

We have to solve

[ u_{xx}+u_{yy}=0,\qquad 0<x<a,\quad 0<y<b ]

with

[ u(0,y)=0,\quad u(a,y)=0,\qquad 0<y<b, ] [ u(x,0)=h(x),\quad u(x,b)=0,\qquad 0\le x\le a. ]

Use separation of variables. Let

[ u(x,y)=X(x)Y(y). ]

Substituting into Laplace’s equation,

[ X''Y+XY''=0 ]

so that

[ \frac{X''}{X}=-\frac{Y''}{Y}=-\lambda. ]

Thus,

[ X''+\lambda X=0, ] [ Y''-\lambda Y=0. ]

The boundary conditions in (x) are

[ X(0)=0,\qquad X(a)=0. ]

Hence we have the Sturm–Liouville problem

[ X''+\lambda X=0,\qquad X(0)=X(a)=0. ]

Its eigenvalues and eigenfunctions are

[ \lambda_n=\left(\frac{n\pi}{a}\right)^2,\qquad X_n(x)=\sin\left(\frac{n\pi x}{a}\right),\qquad n=1,2,3,\dots ]

For (Y), we get

[ Y_n''-\left(\frac{n\pi}{a}\right)^2Y_n=0. ]

Let

[ \mu_n=\frac{n\pi}{a}. ]

Then

[ Y_n(y)=A_n\cosh(\mu_n y)+B_n\sinh(\mu_n y). ]

Using the boundary condition

[ u(x,b)=0, ]

we require

[ Y_n(b)=0. ]

A convenient solution satisfying this is

[ Y_n(y)=\sinh\bigl(\mu_n(b-y)\bigr). ]

Therefore, by superposition,

[ u(x,y)=\sum_{n=1}^{\infty} C_n \sin\left(\frac{n\pi x}{a}\right) \sinh\left(\frac{n\pi}{a}(b-y)\right). ]

Now apply the remaining boundary condition

[ u(x,0)=h(x). ]

This gives

[ h(x)=\sum_{n=1}^{\infty} C_n \sinh\left(\frac{n\pi b}{a}\right) \sin\left(\frac{n\pi x}{a}\right). ]

This is the Fourier sine series of (h(x)) on (0<x<a). Hence,

[ C_n\sinh\left(\frac{n\pi b}{a}\right)

\frac{2}{a}\int_0^a h(x) \sin\left(\frac{n\pi x}{a}\right),dx. ]

So,

[ C_n= \frac{2}{a\sinh\left(\frac{n\pi b}{a}\right)} \int_0^a h(x) \sin\left(\frac{n\pi x}{a}\right),dx. ]

Therefore the required solution is

[ \boxed{ u(x,y)= \sum_{n=1}^{\infty} C_n \sin\left(\frac{n\pi x}{a}\right) \sinh\left(\frac{n\pi}{a}(b-y)\right) } ]

where

[ \boxed{ C_n= \frac{2}{a\sinh\left(\frac{n\pi b}{a}\right)} \int_0^a h(x) \sin\left(\frac{n\pi x}{a}\right),dx } ]


OR alternative

For the semi-circular region (r<a,\ 0<\theta<\pi), with

[ u(r,0)=0,\qquad u(r,\pi)=0, ] [ u(a,\theta)=f(\theta), ]

the solution is

[ \boxed{ u(r,\theta)=\sum_{n=1}^{\infty} A_n r^n \sin(n\theta) } ]

where

[ \boxed{ A_n= \frac{2}{\pi a^n} \int_0^\pi f(\theta)\sin(n\theta),d\theta } ]

This satisfies Laplace’s equation, the boundary conditions on (\theta=0,\pi), and remains bounded at (r=0).

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

Questions from TU BSc 2nd Year 2082 Applied Statistics (2nd Year · Applied Statistics · 2082), TU 4 Yrs BSc 2nd Yr 2083 Regular (2nd Year · Probability and Inference-I / Stat.201 · 2083), TU 4 Yrs BSc 2nd Yr 2083 Regular Basic (2nd Year · Basic Chemistry II / CHE 201 · 2083), TU BSc 2nd Yr 2083 Regular Zoology (2nd Year · Zoology (Zool. 201) · 2083), TU BSc 2nd Yr 2083 Regular Geology (Bachelor Level (4 Yrs.) / Science & Tech. / II Year · Geology Geo.201 · 2083), TU BSc 2nd Yr 2083 Regular ENV 201 (Bachelor Level (4 yrs.)/ II Year /Sc. & Tech. · ENV.201 · 2083), TU 4 Yrs BSc 2nd Yr 2083 Regular (Bachelor Level (4 yrs.) II Year /Sc. & Tech. · Biochemistry and Microbial Biotechnology (MB.201) · 2083), TU 4 Yrs BSc 2nd Yr 2083 Regular Physics (Bachelor Level (4 Yrs.) Science & Tech. / II Year · Optics, Modern Physics and Electronics; PHY. 201 · 2083), TU BSc 2nd Yrs 2083 Regular Applied (2nd Year · Applied Statistics · 2083), TU 4 Yrs BSc 2nd Yr Differential (2nd Year · Differential Equation (MAT 202) (New Course) · 2082), TU 4 Yrs BSc 2nd Yr Probability and (2nd year · Probability and Inference-I · 2082), TU 4 Years BSc 2nd Yr PHY 201 New (2nd year · PHY. 201 · 2082), TU 4 Years BSc 2nd Yr Applied (2nd Year · Applied Statistics (New Course) · 2082), TU 4 Years BSc 2nd Yr Probability and (2nd year · Probability and Inference - I (Stat. 201) · 2081), TU 4 Years BSc 2nd Yr Linear Algebra (BSc 2nd Year · Linear Algebra (MAT 201) · 2081), TU 4 Years BSc 2nd Yr Differential (BSc 2nd Year · Differential Equation (MAT 202) · 2081), TU 4 Years BSc 2nd Yr Applied (2nd Year · Applied Statistics · 2081), TU 4 Years BSc 2nd Yr Basic Chemistry (2nd Year · Basic Chemistry II (CHE 201) · 2081) and TU 4 Years BSc 2nd Yr Basic Chemistry (2nd Year · General Chemistry II (CHE 201) · 2078). Answers are written for this site; check them against your class notes.