Skip to main content

Chapter 9

Differential Equations

TU BBS past exam questions

Past questions and answers

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

  • Asked 2 times
  • 2082 · Differential Equation (MAT 202) (New Course) · 4+4 marks
  • 2081 · Differential Equation (MAT 202) · 4+4 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.

Answer

(a) Transform into 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 into the second equation: [ \frac{x_1'-x_1''}{2}=3x_1-4\left(\frac{x_1-x_1'}{2}\right) ] [ \frac{x_1'-x_1''}{2}=3x_1-2x_1+2x_1' ] [ \frac{x_1'-x_1''}{2}=x_1+2x_1' ]

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

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

Initial conditions: [ x_1(0)=1 ] and [ x_1'(0)=x_1(0)-2x_2(0)=1-2(2)=-3 ]

So, [ \boxed{x_1(0)=1,\qquad x_1'(0)=-3} ]


(b) Find (x_1) and (x_2)

The auxiliary equation is [ m^2+3m+2=0 ] [ (m+1)(m+2)=0 ] [ m=-1,\ -2 ]

Hence, [ x_1=c_1e^{-t}+c_2e^{-2t} ]

Using (x_1(0)=1): [ c_1+c_2=1 ]

Differentiate: [ x_1'=-c_1e^{-t}-2c_2e^{-2t} ]

Using (x_1'(0)=-3): [ -c_1-2c_2=-3 ] [ c_1+2c_2=3 ]

Solving: [ c_1+c_2=1 ] [ c_1+2c_2=3 ]

Subtracting gives [ c_2=2 ] and then [ c_1=-1 ]

Therefore, [ \boxed{x_1=-e^{-t}+2e^{-2t}} ]

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

Since [ x_1'=e^{-t}-4e^{-2t} ] we get [ x_2=\frac{(-e^{-t}+2e^{-2t})-(e^{-t}-4e^{-2t})}{2} ] [ x_2=\frac{-2e^{-t}+6e^{-2t}}{2} ] [ \boxed{x_2=-e^{-t}+3e^{-2t}} ]

Hence the required solution is [ \boxed{x_1(t)=-e^{-t}+2e^{-2t}} ] [ \boxed{x_2(t)=-e^{-t}+3e^{-2t}} ]

Check: [ x_1(0)=-1+2=1,\qquad x_2(0)=-1+3=2 ] which satisfy the given initial conditions.

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=2sin⁡2t, y(0)=2, y′(0)=−1y'' + 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.

Answer

First alternative

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

Homogeneous equation: [ y''+4y=0. ] Its solution is [ y_h=C_1\cos 2t+C_2\sin 2t. ]

Since the forcing term (2\sin 2t) is resonant, take [ y_p=t(A\cos 2t+B\sin 2t). ] Then [ y_p''+4y_p=-4A\sin 2t+4B\cos 2t. ] Comparing with (2\sin 2t), [ -4A=2,\qquad 4B=0. ] So [ A=-\frac12,\qquad B=0. ] Thus [ y_p=-\frac12 t\cos 2t. ]

General solution: [ y=C_1\cos 2t+C_2\sin 2t-\frac12 t\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. ] So [ 2C_2=-\frac12 \quad\Rightarrow\quad C_2=-\frac14. ]

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


OR alternative

Mass weighing (4) lb: [ m=\frac{W}{g}=\frac{4}{32}=\frac18\ \text{slug}. ] Spring stretch: [ s=1.5\text{ in}=\frac{1.5}{12}=\frac18\ \text{ft}. ] Spring constant: [ k=\frac{W}{s}=\frac{4}{1/8}=32\ \text{lb/ft}. ]

Let (u(t)) be displacement from equilibrium in feet. No damping, external force (2\cos 3t), so [ m u''+k u=2\cos 3t. ] Thus [ \frac18 u''+32u=2\cos 3t. ] Multiplying by (8), [ u''+256u=16\cos 3t. ]

Initial conditions: [ u(0)=\frac{2}{12}=\frac16\ \text{ft},\qquad u'(0)=0. ]

Homogeneous solution: [ u_h=C_1\cos 16t+C_2\sin 16t. ]

For the particular solution, try [ u_p=A\cos 3t+B\sin 3t. ] Substituting, [ (256-9)(A\cos 3t+B\sin 3t)=16\cos 3t. ] So [ 247A=16,\qquad B=0, ] giving [ A=\frac{16}{247}. ] Hence [ u_p=\frac{16}{247}\cos 3t. ]

General solution: [ u=C_1\cos 16t+C_2\sin 16t+\frac{16}{247}\cos 3t. ]

Using (u(0)=\frac16), [ C_1+\frac{16}{247}=\frac16 ] [ C_1=\frac16-\frac{16}{247}=\frac{151}{1482}. ]

Also, [ u'=-16C_1\sin 16t+16C_2\cos 16t-\frac{48}{247}\sin 3t. ] At (t=0), [ u'(0)=16C_2=0 \quad\Rightarrow\quad C_2=0. ]

Therefore, [ \boxed{u(t)=\frac{151}{1482}\cos 16t+\frac{16}{247}\cos 3t.} ]

Equivalently, [ \boxed{u(t)=\left(\frac16-\frac{16}{247}\right)\cos 16t+\frac{16}{247}\cos 3t.} ]

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

The differential equation is y′′+4y=t2+3et.y''+4y=t^2+3e^t.

1. Homogeneous solution

Solve y′′+4y=0.y''+4y=0. Auxiliary equation: r2+4=0  ⟹  r=±2i.r^2+4=0 \implies r=\pm 2i. So, yh=C1cos⁡2t+C2sin⁡2t.y_h=C_1\cos 2t+C_2\sin 2t.

2. Particular solution

For (t^2), try yp1=At2+Bt+C.y_{p1}=At^2+Bt+C. Then yp1′′=2A.y_{p1}''=2A. Substitute: 2A+4(At2+Bt+C)=t22A+4(At^2+Bt+C)=t^2 4At2+4Bt+(2A+4C)=t2.4At^2+4Bt+(2A+4C)=t^2. Comparing coefficients: 4A=1  ⟹  A=14,4A=1 \implies A=\frac14, 4B=0  ⟹  B=0,4B=0 \implies B=0, 2A+4C=0  ⟹  12+4C=0  ⟹  C=−18.2A+4C=0 \implies \frac12+4C=0 \implies C=-\frac18. Thus, yp1=14t2−18.y_{p1}=\frac14t^2-\frac18.

For (3e^t), try yp2=Det.y_{p2}=De^t. Then yp2′′+4yp2=Det+4Det=5Det.y_{p2}''+4y_{p2}=De^t+4De^t=5De^t. So, 5D=3  ⟹  D=35.5D=3 \implies D=\frac35. Thus, yp2=35et.y_{p2}=\frac35e^t.

Hence the general solution is y=C1cos⁡2t+C2sin⁡2t+14t2−18+35et.y=C_1\cos 2t+C_2\sin 2t+\frac14t^2-\frac18+\frac35e^t.

3. Apply initial conditions

Given (y(0)=0): C1−18+35=0C_1-\frac18+\frac35=0 C1=18−35=5−2440=−1940.C_1=\frac18-\frac35=\frac{5-24}{40}=-\frac{19}{40}.

Now differentiate: y′=−2C1sin⁡2t+2C2cos⁡2t+12t+35et.y'=-2C_1\sin 2t+2C_2\cos 2t+\frac12t+\frac35e^t. Given (y'(0)=2): 2C2+35=22C_2+\frac35=2 2C2=752C_2=\frac75 C2=710.C_2=\frac7{10}.

Therefore, the solution of the initial value problem is y=−1940cos⁡2t+710sin⁡2t+14t2−18+35et.\boxed{y=-\frac{19}{40}\cos 2t+\frac7{10}\sin 2t+\frac14t^2-\frac18+\frac35e^t.}

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

  • 2082 · Differential Equation (MAT 202) (New Course) · 3+4 marks

(a) How do you solve the homogeneous differential equation of the form f(dydx,yx)=0f\left(\frac{dy}{dx}, \frac{y}{x}\right) = 0? (b) Find complete and singular solution of the equation y=px+ap(1−p)y = px + ap(1-p).

Answer

(a) Equation of the form f(dydx,yx)=0f\left(\dfrac{dy}{dx}, \dfrac{y}{x}\right)=0

Write p=dydxp = \dfrac{dy}{dx}. The equation relates pp and y/xy/x only, so it is homogeneous. Two cases:

Method 1 (solve for pp). Solve f(p,y/x)=0f(p, y/x)=0 for pp to get p=ϕ ⁣(yx)p = \phi\!\left(\dfrac{y}{x}\right). Put y=vxy = vx, so dydx=v+xdvdx\dfrac{dy}{dx} = v + x\dfrac{dv}{dx}:

v+xdvdx=ϕ(v)  ⇒  dvϕ(v)−v=dxxv + x\frac{dv}{dx} = \phi(v)\;\Rightarrow\; \frac{dv}{\phi(v)-v} = \frac{dx}{x}

The variables are separated. Integrate:

∫dvϕ(v)−v=ln⁡x+c\int\frac{dv}{\phi(v)-v} = \ln x + c

then replace vv by y/xy/x.

Method 2 (solve for y/xy/x). If it is easier, solve for yx=ψ(p)\dfrac{y}{x} = \psi(p), i.e. y=xψ(p)y = x\psi(p). Differentiating with respect to xx:

p=ψ(p)+xψ′(p)dpdx  ⇒  dxx=ψ′(p) dpp−ψ(p)p = \psi(p) + x\psi'(p)\frac{dp}{dx}\;\Rightarrow\; \frac{dx}{x} = \frac{\psi'(p)\,dp}{p-\psi(p)}

Integrating gives xx in terms of pp (with a constant cc); together with y=xψ(p)y = x\psi(p) this gives the solution in parametric form (eliminate pp if possible).

(b) y=px+ap(1−p)y = px + ap(1-p)

This is Clairaut's equation y=px+F(p)y = px + F(p) with F(p)=ap(1−p)F(p) = ap(1-p).

Complete solution: replace pp by an arbitrary constant cc:

y=cx+ac(1−c)y = cx + ac(1-c)

Singular solution: differentiate the complete solution with respect to cc:

0=x+a(1−2c)  ⇒  c=x+a2a0 = x + a(1-2c)\;\Rightarrow\; c = \frac{x+a}{2a}

Substituting back, with y=c(x+a)−ac2y = c(x+a) - ac^2:

y=(x+a)22a−(x+a)24a=(x+a)24ay = \frac{(x+a)^2}{2a} - \frac{(x+a)^2}{4a} = \frac{(x+a)^2}{4a}

(x+a)2=4ay\boxed{(x+a)^2 = 4ay}

This parabola is the envelope of the family of straight lines and is the singular solution.

Answer checked by Claude Sonnet · 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.

Answer

Let p=∂z/∂xp = \partial z/\partial x. The equation is

4yz p=−(4+2y)  ⇒  2z ∂z∂x=−2+yy(i.e. ∂(z2)∂x=−2+yy)4yz\,p = -(4+2y)\;\Rightarrow\; 2z\,\frac{\partial z}{\partial x} = -\frac{2+y}{y}\quad\Big(\text{i.e. } \frac{\partial (z^2)}{\partial x} = -\frac{2+y}{y}\Big)

It contains no qq, so yy acts as a constant and we integrate with respect to xx:

z2+2+yy x=f(y)z^2 + \frac{2+y}{y}\,x = f(y)

where ff is arbitrary (Lagrange's auxiliary equations dx4yz=dy0=dz−(4+2y)\dfrac{dx}{4yz} = \dfrac{dy}{0} = \dfrac{dz}{-(4+2y)} give y=c1y=c_1 and z2+2+yyx=c2z^2 + \frac{2+y}{y}x = c_2). The integral surfaces are therefore ϕ ⁣(y, z2+2+yyx)=0\phi\!\left(y,\ z^2 + \dfrac{2+y}{y}x\right)=0.

Fitting the curve y2+z2=1, x+z=2y^2+z^2=1,\ x+z=2.

On the curve, x=2−zx = 2 - z and z2=1−y2z^2 = 1-y^2. Put c=2+yyc = \dfrac{2+y}{y}. Then

c2=z2+c x=1−y2+c(2−z)  ⇒  c2−(1−y2)−2c=−c zc_2 = z^2 + c\,x = 1 - y^2 + c(2-z)\;\Rightarrow\; c_2 - (1-y^2) - 2c = -c\,z

Squaring removes zz (using z2=1−y2z^2 = 1-y^2 on the curve):

[c2−(1−y2)−2c]2=c2(1−y2)\big[c_2 - (1-y^2) - 2c\big]^2 = c^2(1-y^2)

Now replace c2c_2 by z2+cxz^2 + cx and cc by 2+yy\dfrac{2+y}{y}:

[z2+y2−1+c(x−2)]2=c2(1−y2)\big[z^2 + y^2 - 1 + c(x-2)\big]^2 = c^2(1-y^2)

Multiplying by y2y^2:

[y(y2+z2−1)+(2+y)(x−2)]2=(2+y)2(1−y2)\boxed{\big[y(y^2+z^2-1) + (2+y)(x-2)\big]^2 = (2+y)^2(1-y^2)}

Check: on the curve y2+z2−1=0y^2+z^2-1=0 and x−2=−zx-2=-z, so the left side is (2+y)2z2=(2+y)2(1−y2)(2+y)^2z^2 = (2+y)^2(1-y^2), which equals the right side. The surface therefore contains the given circle and satisfies the differential equation (it is a function of yy and z2+2+yyxz^2+\frac{2+y}{y}x only).

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

zxx+2zxy+zyy=0.z_{xx}+2z_{xy}+z_{yy}=0.

Here,

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

So,

B2−4AC=4−4=0.B^2-4AC=4-4=0.

Hence the equation is parabolic.

The characteristic equation is

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

So,

dydx=1\frac{dy}{dx}=1

and the characteristic curves are

y−x=corx−y=c.y-x=c \quad \text{or} \quad x-y=c.

Choose the new variables

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

Then

∂∂x=∂∂ξ+∂∂η,∂∂y=∂∂ξ−∂∂η.\frac{\partial}{\partial x}=\frac{\partial}{\partial \xi}+\frac{\partial}{\partial \eta}, \qquad \frac{\partial}{\partial y}=\frac{\partial}{\partial \xi}-\frac{\partial}{\partial \eta}.

Therefore,

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

Substituting in the given equation,

zxx+2zxy+zyy=4zξξ=0.z_{xx}+2z_{xy}+z_{yy} =4z_{\xi\xi}=0.

Thus the canonical form is

zξξ=0.\boxed{z_{\xi\xi}=0}.

Now integrate with respect to (\xi):

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

where (\phi) is an arbitrary function of (\eta). Again integrating,

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

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

Returning to the original variables (\xi=x+y,\ \eta=x-y),

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

Hence the general solution of the given PDE is

z=(x+y)F(x−y)+G(x−y)\boxed{z=(x+y)F(x-y)+G(x-y)}

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

Fourier Series of the Given Function

Given
[ f(x)= \begin{cases} x+1, & -1\le x<0,\ 1-x, & 0\le x<1, \end{cases} \qquad f(x+2)=f(x). ]

The function is periodic with period (2), so (2L=2\Rightarrow L=1).
Also, [ f(-x)=f(x), ] so (f(x)) is an even function. Hence all sine coefficients (b_n=0), and the Fourier series contains only cosine terms: [ f(x)=\frac{a_0}{2}+\sum_{n=1}^{\infty}a_n\cos(n\pi x). ]

1. Find (a_0): [ a_0=\frac1L\int_{-L}^{L}f(x),dx =\int_{-1}^{1}f(x),dx. ] Since (f(x)) is even, [ a_0=2\int_0^1(1-x),dx =2\left[x-\frac{x^2}{2}\right]_0^1 =2\left(1-\frac12\right)=1. ] Thus, [ \frac{a_0}{2}=\frac12. ]

2. Find (a_n): [ a_n=\frac1L\int_{-1}^{1}f(x)\cos(n\pi x),dx =2\int_0^1(1-x)\cos(n\pi x),dx. ] Let [ I=\int_0^1(1-x)\cos(n\pi x),dx. ] Using integration by parts, [ I=\left[(1-x)\frac{\sin(n\pi x)}{n\pi}\right]_0^1 +\frac1{n\pi}\int_0^1\sin(n\pi x),dx. ] The boundary term is zero. Hence, [ I=\frac1{n\pi}\left[-\frac{\cos(n\pi x)}{n\pi}\right]_0^1 =\frac{1-\cos(n\pi)}{n^2\pi^2} =\frac{1-(-1)^n}{n^2\pi^2}. ] Therefore, [ a_n=2I=\frac{2(1-(-1)^n)}{n^2\pi^2}. ] For even (n), (a_n=0).
For odd (n), (a_n=\dfrac{4}{n^2\pi^2}).

3. Fourier series: [ \boxed{ f(x)=\frac12+\frac{4}{\pi^2}\sum_{k=0}^{\infty} \frac{\cos\big((2k+1)\pi x\big)}{(2k+1)^2} } ] or [ \boxed{ f(x)=\frac12+\frac{4}{\pi^2} \left( \frac{\cos\pi x}{1^2} +\frac{\cos3\pi x}{3^2} +\frac{\cos5\pi x}{5^2} +\cdots \right) } ]

Sketch of the graph for three periods

Over one period ([-1,1]), the function is [ f(x)=1-|x|. ] It has a peak of height (1) at (x=0) and zeros at (x=-1,1).
Because the period is (2), the graph repeats every (2) units. Thus:

  • peaks occur at even integers: (x=\ldots,-2,0,2,\ldots)
  • zeros occur at odd integers: (x=\ldots,-3,-1,1,3,\ldots)

For three periods, say from (x=-3) to (x=3), the graph is a triangular wave:

      /\        /\        /\
     /  \      /  \      /  \
____/    \____/    \____/    \____
-3   -2   -1    0    1    2    3

The straight-line segments join the points [ (-3,0),(-2,1),(-1,0),(0,1),(1,0),(2,1),(3,0). ]

Hence the required Fourier series is [ \boxed{ f(x)=\frac12+\frac{4}{\pi^2} \left( \cos\pi x+\frac{\cos3\pi x}{3^2} +\frac{\cos5\pi x}{5^2} +\cdots \right) } ] and the graph is the triangular wave shown above.

Answer checked by Claude Haiku · 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 L2a2uss=utt\frac{L^2}{a^2} u_{ss} = u_{tt}, so that it is non dimensionalized with respect to the spatial variable.

Answer

Let
[ u(x,y)=X(x)Y(y) ] be a separated solution of Laplace’s equation [ u_{xx}+u_{yy}=0. ] Substituting gives [ X''Y+XY''=0 ] or [ \frac{X''}{X}=-\frac{Y''}{Y}=\mu^2 \quad \text{(separation constant)}. ] Hence [ X''-\mu^2 X=0,\qquad Y''+\mu^2 Y=0. ]

The boundary conditions in the (y)-direction are [ u(x,0)=u(x,b)=0, ] so [ Y(0)=Y(b)=0. ] The non-trivial solutions of [ Y''+\mu^2Y=0 ] with these boundary conditions are [ Y_n(y)=\sin\left(\frac{n\pi y}{b}\right),\qquad \mu=\frac{n\pi}{b},\quad n=1,2,3,\dots ]

For each (n), the equation for (X) becomes [ X_n''-\left(\frac{n\pi}{b}\right)^2X_n=0. ] Its general solution is [ X_n(x)=A_n e^{n\pi x/b}+B_n e^{-n\pi x/b}. ] Using the boundary condition [ u(0,y)=0 \implies X_n(0)=0, ] we get (A_n+B_n=0), so [ X_n(x)=C_n\sinh\left(\frac{n\pi x}{b}\right). ]

Therefore the separated solutions are [ u_n(x,y)=C_n\sinh\left(\frac{n\pi x}{b}\right) \sin\left(\frac{n\pi y}{b}\right). ] By superposition, [ u(x,y)=\sum_{n=1}^{\infty}C_n\sinh\left(\frac{n\pi x}{b}\right) \sin\left(\frac{n\pi y}{b}\right). ]

Now apply the remaining boundary condition [ u(a,y)=f(y). ] Then [ f(y)=\sum_{n=1}^{\infty}C_n\sinh\left(\frac{n\pi a}{b}\right) \sin\left(\frac{n\pi y}{b}\right). ] This is the Fourier sine series of (f(y)) on (0<y<b). Hence [ C_n\sinh\left(\frac{n\pi a}{b}\right)

\frac{2}{b}\int_0^b f(y)\sin\left(\frac{n\pi y}{b}\right),dy. ] Thus [ \boxed{ C_n= \frac{2}{b\sinh\left(\frac{n\pi a}{b}\right)} \int_0^b f(y)\sin\left(\frac{n\pi y}{b}\right),dy } ] and the required solution is [ \boxed{ u(x,y)=\sum_{n=1}^{\infty}C_n \sinh\left(\frac{n\pi x}{b}\right) \sin\left(\frac{n\pi y}{b}\right) } ] with (C_n) as above. This solution satisfies all the given Dirichlet boundary conditions, provided (f) is piecewise smooth and (f(0)=f(b)=0).

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⁡2t,y(0)=2,y′(0)=−1.y'' + 4y = 2\sin 2t, 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

First alternative

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

Homogeneous equation: y′′+4y=0.y''+4y=0. Its characteristic equation is r2+4=0  ⟹  r=±2i.r^2+4=0 \implies r=\pm 2i. So yh=C1cos⁡2t+C2sin⁡2t.y_h=C_1\cos 2t+C_2\sin 2t.

Since the forcing term (2\sin 2t) is already in (y_h), take yp=t(Acos⁡2t+Bsin⁡2t).y_p=t(A\cos 2t+B\sin 2t).

Differentiate: yp′=Acos⁡2t+Bsin⁡2t+t(−2Asin⁡2t+2Bcos⁡2t),y_p'=A\cos 2t+B\sin 2t+t(-2A\sin 2t+2B\cos 2t), yp′′=−4Asin⁡2t+4Bcos⁡2t−4t(Acos⁡2t+Bsin⁡2t).y_p''=-4A\sin 2t+4B\cos 2t-4t(A\cos 2t+B\sin 2t).

Substitute into (y''+4y=2\sin 2t): −4Asin⁡2t+4Bcos⁡2t=2sin⁡2t.-4A\sin 2t+4B\cos 2t=2\sin 2t. Hence −4A=2  ⟹  A=−12,4B=0  ⟹  B=0.-4A=2 \implies A=-\frac12,\qquad 4B=0 \implies B=0. Thus yp=−t2cos⁡2t.y_p=-\frac{t}{2}\cos 2t.

General solution: y(t)=C1cos⁡2t+C2sin⁡2t−t2cos⁡2t.y(t)=C_1\cos 2t+C_2\sin 2t-\frac{t}{2}\cos 2t.

Use (y(0)=2): C1=2.C_1=2.

Now y′(t)=−2C1sin⁡2t+2C2cos⁡2t−12cos⁡2t+tsin⁡2t.y'(t)=-2C_1\sin 2t+2C_2\cos 2t-\frac12\cos 2t+t\sin 2t. Using (y'(0)=-1): 2C2−12=−1  ⟹  2C2=−12  ⟹  C2=−14.2C_2-\frac12=-1 \implies 2C_2=-\frac12 \implies C_2=-\frac14.

Therefore, y(t)=2cos⁡2t−14sin⁡2t−t2cos⁡2t.\boxed{y(t)=2\cos 2t-\frac14\sin 2t-\frac{t}{2}\cos 2t.}


Second alternative

Spring constant: k=Fx=20.10=20 N/m.k=\frac{F}{x}=\frac{2}{0.10}=20\ \text{N/m}.

Mass: m=2 kg.m=2\ \text{kg}.

Damping coefficient: c=Fv=41=4 N s/m.c=\frac{F}{v}=\frac{4}{1}=4\ \text{N s/m}.

External force: F(t)=4cos⁡2t.F(t)=4\cos 2t.

Equation of motion: mx′′+cx′+kx=F(t),m x''+c x'+k x=F(t), so 2x′′+4x′+20x=4cos⁡2t.2x''+4x'+20x=4\cos 2t. Dividing by (2): x′′+2x′+10x=2cos⁡2t.x''+2x'+10x=2\cos 2t.

Homogeneous equation: x′′+2x′+10x=0.x''+2x'+10x=0. Characteristic equation: r2+2r+10=0  ⟹  r=−1±3i.r^2+2r+10=0 \implies r=-1\pm 3i. Hence xh=e−t(C1cos⁡3t+C2sin⁡3t).x_h=e^{-t}(C_1\cos 3t+C_2\sin 3t).

For the particular solution, take xp=Acos⁡2t+Bsin⁡2t.x_p=A\cos 2t+B\sin 2t. Then xp′=−2Asin⁡2t+2Bcos⁡2t,x_p'=-2A\sin 2t+2B\cos 2t, xp′′=−4Acos⁡2t−4Bsin⁡2t.x_p''=-4A\cos 2t-4B\sin 2t.

Substitute into (x''+2x'+10x=2\cos 2t): (6A+4B)cos⁡2t+(−4A+6B)sin⁡2t=2cos⁡2t.(6A+4B)\cos 2t+(-4A+6B)\sin 2t=2\cos 2t. Thus 6A+4B=2,−4A+6B=0.6A+4B=2,\qquad -4A+6B=0. Solving, A=313,B=213.A=\frac{3}{13},\qquad B=\frac{2}{13}. So xp=313cos⁡2t+213sin⁡2t.x_p=\frac{3}{13}\cos 2t+\frac{2}{13}\sin 2t.

General solution: x(t)=e−t(C1cos⁡3t+C2sin⁡3t)+313cos⁡2t+213sin⁡2t.x(t)=e^{-t}(C_1\cos 3t+C_2\sin 3t)+\frac{3}{13}\cos 2t+\frac{2}{13}\sin 2t.

Initial conditions: initially at rest at equilibrium means x(0)=0,x′(0)=0.x(0)=0,\qquad x'(0)=0.

From (x(0)=0): C1+313=0  ⟹  C1=−313.C_1+\frac{3}{13}=0 \implies C_1=-\frac{3}{13}.

Now x′(t)=e−t[(−C1+3C2)cos⁡3t+(−C2−3C1)sin⁡3t]−613sin⁡2t+413cos⁡2t.x'(t)=e^{-t}\left[(-C_1+3C_2)\cos 3t+(-C_2-3C_1)\sin 3t\right]-\frac{6}{13}\sin 2t+\frac{4}{13}\cos 2t. Using (x'(0)=0): −C1+3C2+413=0.-C_1+3C_2+\frac{4}{13}=0. Since (C_1=-\frac{3}{13}), 313+3C2+413=0  ⟹  3C2=−713  ⟹  C2=−739.\frac{3}{13}+3C_2+\frac{4}{13}=0 \implies 3C_2=-\frac{7}{13} \implies C_2=-\frac{7}{39}.

Therefore, x(t)=e−t(−313cos⁡3t−739sin⁡3t)+313cos⁡2t+213sin⁡2t.\boxed{x(t)=e^{-t}\left(-\frac{3}{13}\cos 3t-\frac{7}{39}\sin 3t\right)+\frac{3}{13}\cos 2t+\frac{2}{13}\sin 2t.}

The position is in metres.

Transient part: xtr(t)=e−t(−313cos⁡3t−739sin⁡3t).x_{\text{tr}}(t)=e^{-t}\left(-\frac{3}{13}\cos 3t-\frac{7}{39}\sin 3t\right).

Steady-state part: xss(t)=313cos⁡2t+213sin⁡2t.x_{\text{ss}}(t)=\frac{3}{13}\cos 2t+\frac{2}{13}\sin 2t.

Equivalently, xss(t)=1313cos⁡(2t−δ),x_{\text{ss}}(t)=\frac{\sqrt{13}}{13}\cos(2t-\delta), where δ=tan⁡−1(23).\delta=\tan^{-1}\left(\frac{2}{3}\right).

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.

1. Complementary function

Solve the homogeneous equation: y′′+y′−2y=0.y''+y'-2y=0.

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.

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

2. Particular integral

Since the RHS 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.

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

3. General solution

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

4. Apply initial conditions

Given (y(0)=0): C1+C2−12=0C_1+C_2-\frac12=0 C1+C2=12.(1)C_1+C_2=\frac12. \quad (1)

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

Given (y'(0)=1): C1−2C2−1=1C_1-2C_2-1=1 C1−2C2=2.(2)C_1-2C_2=2. \quad (2)

Subtracting (2) from (1): (C1+C2)−(C1−2C2)=12−2(C_1+C_2)-(C_1-2C_2)=\frac12-2 3C2=−323C_2=-\frac32 C2=−12.C_2=-\frac12.

From (1): C1−12=12  ⟹  C1=1.C_1-\frac12=\frac12 \implies C_1=1.

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

This is the required solution of the given initial value problem.

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

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

(a) How do you solve the Clairaut's differential equation y=px+f(p)y = p x + f(p) where p=dydxp = \frac{dy}{dx}
(b) Find the complete and singular solution of the equation y=px+ap(1−p)y = px + ap(1 - p)

Answer

(a) Method for solving Clairaut’s equation

The Clairaut equation is
[ y=px+f(p),\qquad p=\frac{dy}{dx}. ]

Step 1: Differentiate both sides with respect to (x):
[ p=p+x\frac{dp}{dx}+f'(p)\frac{dp}{dx}. ]

Step 2: Cancel (p) from both sides:
[ \left[x+f'(p)\right]\frac{dp}{dx}=0. ]

So either
[ \frac{dp}{dx}=0 ] or
[ x+f'(p)=0. ]

Step 3: If (\frac{dp}{dx}=0), then (p=c), where (c) is an arbitrary constant. Substituting (p=c) in the original equation gives the complete solution:
[ y=cx+f(c). ]

Step 4: If (x+f'(p)=0), eliminate (p) between this equation and the original equation. The resulting relation, free from arbitrary constants, is the singular solution. It is the envelope of the family (y=cx+f(c)).


(b) Complete and singular solution of

[ y=px+ap(1-p) ]

This is a Clairaut equation with
[ f(p)=ap(1-p). ]

Complete solution

Replace (p) by the arbitrary constant (c):
[ y=cx+ac(1-c). ]

Hence the complete solution is
[ \boxed{y=cx+ac-ac^2} ] or
[ \boxed{y=cx+ac(1-c)}. ]

Singular solution

For the singular solution, use
[ x+f'(p)=0. ]

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

Thus
[ x+a(1-2p)=0 ] [ x+a-2ap=0 ] [ 2ap=x+a ] [ p=\frac{x+a}{2a}. ]

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

So
[ y=x\left(\frac{x+a}{2a}\right)+a\left(\frac{x+a}{2a}\right)\left(1-\frac{x+a}{2a}\right). ]

Now
[ 1-\frac{x+a}{2a}=\frac{2a-x-a}{2a}=\frac{a-x}{2a}. ]

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

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

Simplify:
[ 2x(x+a)+(x+a)(a-x) =(x+a)\left[2x+a-x\right] =(x+a)(x+a) =(x+a)^2. ]

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

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

So, for (a\neq 0),
[ \boxed{\text{Complete solution: } y=cx+ac(1-c)} ] [ \boxed{\text{Singular solution: } y=\frac{(x+a)^2}{4a}}. ]

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

Setting up. The given one-parameter family is

z=cxy(x2+y2)⟹u(x,y,z)=zxy(x2+y2)=c.z=cxy(x^2+y^2)\quad\Longrightarrow\quad u(x,y,z)=\frac{z}{xy(x^2+y^2)}=c.

Its normal direction is (ux,uy,uz)(u_x,u_y,u_z). With D=xy(x2+y2)=x3y+xy3D=xy(x^2+y^2)=x^3y+xy^3:

ux=−z y(3x2+y2)D2,uy=−z x(x2+3y2)D2,uz=1D,u_x=-\frac{z\,y(3x^2+y^2)}{D^2},\qquad u_y=-\frac{z\,x(x^2+3y^2)}{D^2},\qquad u_z=\frac1D,

so the normal is proportional to

P=−zy(3x2+y2),Q=−zx(x2+3y2),R=xy(x2+y2).P=-zy(3x^2+y^2),\qquad Q=-zx(x^2+3y^2),\qquad R=xy(x^2+y^2).

A surface z=ϕ(x,y)z=\phi(x,y) cuts every member of the family orthogonally when its normal (p,q,−1)(p,q,-1) is perpendicular to (P,Q,R)(P,Q,R), i.e. when it satisfies the Lagrange equation

Pp+Qq=R.Pp+Qq=R.

Auxiliary equations.

dx−zy(3x2+y2)=dy−zx(x2+3y2)=dzxy(x2+y2).\frac{dx}{-zy(3x^2+y^2)}=\frac{dy}{-zx(x^2+3y^2)}=\frac{dz}{xy(x^2+y^2)}.

First integral. Taking the first two ratios:

x(x2+3y2) dx=y(3x2+y2) dy.x(x^2+3y^2)\,dx=y(3x^2+y^2)\,dy .

Put u=x+y, v=x−yu=x+y,\ v=x-y. Then x3+3xy2=12(u3+v3)x^3+3xy^2=\tfrac12(u^3+v^3) and 3x2y+y3=12(u3−v3)3x^2y+y^3=\tfrac12(u^3-v^3), and the equation becomes

u3 dv+v3 du=0 ⇒ duu3+dvv3=0 ⇒ 1u2+1v2=const.u^3\,dv+v^3\,du=0\ \Rightarrow\ \frac{du}{u^3}+\frac{dv}{v^3}=0\ \Rightarrow\ \frac1{u^2}+\frac1{v^2}=\text{const}.

Since u2+v2=2(x2+y2)u^2+v^2=2(x^2+y^2) and uv=x2−y2uv=x^2-y^2,

x2+y2(x2−y2)2=C1.\frac{x^2+y^2}{(x^2-y^2)^2}=C_1 .

Second integral. Use multipliers x, y, 4zx,\ y,\ 4z. Then

x P+y Q=−4zxy(x2+y2)=−4z R,x\,P+y\,Q=-4zxy(x^2+y^2)=-4z\,R,

so

x dx+y dy−4z xy(x2+y2)=dzxy(x2+y2) ⇒ x dx+y dy+4z dz=0,\frac{x\,dx+y\,dy}{-4z\,xy(x^2+y^2)}=\frac{dz}{xy(x^2+y^2)}\ \Rightarrow\ x\,dx+y\,dy+4z\,dz=0,

x2+y2+4z2=C2.x^2+y^2+4z^2=C_2 .

General orthogonal surface.

x2+y2+4z2=φ ⁣[x2+y2(x2−y2)2].x^2+y^2+4z^2=\varphi\!\left[\frac{x^2+y^2}{(x^2-y^2)^2}\right].

Condition. The surface must pass through x2−y2=a2, z=0x^2-y^2=a^2,\ z=0. On this curve C2=x2+y2C_2=x^2+y^2 and C1=x2+y2a4C_1=\dfrac{x^2+y^2}{a^4}, so C2=a4C1C_2=a^4C_1, i.e. φ(t)=a4t\varphi(t)=a^4t.

Required surface:

(x2+y2+4z2)(x2−y2)2=a4(x2+y2)\boxed{(x^2+y^2+4z^2)(x^2-y^2)^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/∂x, q=∂z/∂yp=\partial z/\partial x,\ q=\partial z/\partial y.

Charpit's auxiliary equations

dxFp=dyFq=dzpFp+qFq=−dpFx+pFz=−dqFy+qFz.\frac{dx}{F_p}=\frac{dy}{F_q}=\frac{dz}{pF_p+qF_q}=\frac{-dp}{F_x+pF_z}=\frac{-dq}{F_y+qF_z}.

Here Fp=2px, Fq=2qy, Fx=p2, Fy=q2, Fz=−1F_p=2px,\ F_q=2qy,\ F_x=p^2,\ F_y=q^2,\ F_z=-1, so pFp+qFq=2(p2x+q2y)=2zpF_p+qF_q=2(p^2x+q^2y)=2z, Fx+pFz=p2−pF_x+pF_z=p^2-p and Fy+qFz=q2−qF_y+qF_z=q^2-q:

dx2px=dy2qy=dz2z=−dpp(p−1)=−dqq(q−1).\frac{dx}{2px}=\frac{dy}{2qy}=\frac{dz}{2z}=\frac{-dp}{p(p-1)}=\frac{-dq}{q(q-1)}.

Finding pp and qq. From dx2px=−dpp(p−1)\dfrac{dx}{2px}=\dfrac{-dp}{p(p-1)}:

dx2x=−dpp−1 ⇒ 12ln⁡x=−ln⁡(p−1)+ln⁡a ⇒ p=1+ax.\frac{dx}{2x}=-\frac{dp}{p-1}\ \Rightarrow\ \tfrac12\ln x=-\ln(p-1)+\ln a\ \Rightarrow\ p=1+\frac{a}{\sqrt x}.

In the same way, from dy2qy=−dqq(q−1)\dfrac{dy}{2qy}=\dfrac{-dq}{q(q-1)}:

q=1+by.q=1+\frac{b}{\sqrt y}.

Check against the PDE. With these values, p2x=(x+a)2p^2x=(\sqrt x+a)^2 and q2y=(y+b)2q^2y=(\sqrt y+b)^2, so p2x+q2y=zp^2x+q^2y=z gives z=(x+a)2+(y+b)2z=(\sqrt x+a)^2+(\sqrt y+b)^2. Here aa and bb are independent constants, so this is a two-parameter solution. It can be confirmed by integrating dz=p dx+q dydz=p\,dx+q\,dy.

Integrating dz=p dx+q dydz=p\,dx+q\,dy

dz=(1+ax)dx+(1+by)dydz=\left(1+\frac{a}{\sqrt x}\right)dx+\left(1+\frac{b}{\sqrt y}\right)dy

z=x+y+2ax+2by+C.z=x+y+2a\sqrt x+2b\sqrt y+C.

Substituting p,qp,q in p2x+q2y=zp^2x+q^2y=z:

z=(x+2ax+a2)+(y+2by+b2)=x+y+2ax+2by+a2+b2,z=(x+2a\sqrt x+a^2)+(y+2b\sqrt y+b^2)=x+y+2a\sqrt x+2b\sqrt y+a^2+b^2,

so C=a2+b2C=a^2+b^2.

Complete integral:

z=x+y+2ax+2by+a2+b2=(x+a)2+(y+b)2\boxed{z=x+y+2a\sqrt x+2b\sqrt y+a^2+b^2=(\sqrt x+a)^2+(\sqrt y+b)^2}

where a,ba,b 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(π)=0.y'' + \lambda y = 0, y'(0) = 0, y(\pi) = 0.

Answer

Sketch of f(x)f(x) over three periods (period 2L2L, shown on −3L≤x≤3L-3L\le x\le 3L)

On −L≤x≤0-L\le x\le 0 the graph is the straight line y=x+Ly=x+L, rising from 00 at x=−Lx=-L to LL at x=0x=0. On 0<x<L0<x<L it is the horizontal line y=Ly=L. This pattern repeats every 2L2L:

 f
 L |      ____       ____       ____
   |     /    |     /    |     /    |
   |    /     |    /     |    /     |
 0 |___/      |___/      |___/      |___
   -3L  -2L  -L   0    L   2L   3L

Each period is a ramp from 00 to LL followed by a flat level at LL. The function is continuous at x=0,±2L,…x=0,\pm2L,\dots and has a jump from LL down to 00 at x=±L,±3L,…x=\pm L,\pm3L,\dots

Fourier series. Period 2L2L:

f(x)=a02+∑n=1∞(ancos⁡nπxL+bnsin⁡nπxL).f(x)=\frac{a_0}{2}+\sum_{n=1}^{\infty}\left(a_n\cos\frac{n\pi x}{L}+b_n\sin\frac{n\pi x}{L}\right).

Constant term.

a0=1L[∫−L0(x+L) dx+∫0LL dx]=1L[L22+L2]=3L2,a02=3L4.a_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},\qquad \frac{a_0}{2}=\frac{3L}{4}.

Cosine coefficients. Write k=nπ/Lk=n\pi/L. The integral over (0,L)(0,L) of Lcos⁡kxL\cos kx is 00. Integrating by parts,

an=1L∫−L0(x+L)cos⁡kx dx=1L[(x+L)sin⁡kxk]−L0−1Lk∫−L0sin⁡kx dx=1−(−1)nL k2a_n=\frac1L\int_{-L}^{0}(x+L)\cos kx\,dx=\frac1L\left[\frac{(x+L)\sin kx}{k}\right]_{-L}^{0}-\frac1{Lk}\int_{-L}^{0}\sin kx\,dx =\frac{1-(-1)^n}{L\,k^2}

an=L [1−(−1)n]n2π2(=2Ln2π2 for odd n, 0 for even n).a_n=\frac{L\,[1-(-1)^n]}{n^2\pi^2}\quad\left(=\frac{2L}{n^2\pi^2}\text{ for odd }n,\ 0\text{ for even }n\right).

Sine coefficients.

bn=1L[∫−L0(x+L)sin⁡kx dx+∫0LLsin⁡kx dx].b_n=\frac1L\left[\int_{-L}^{0}(x+L)\sin kx\,dx+\int_0^L L\sin kx\,dx\right].

First integral: [−(x+L)cos⁡kxk]−L0+1k∫−L0cos⁡kx dx=−Lk\left[-\dfrac{(x+L)\cos kx}{k}\right]_{-L}^{0}+\dfrac1k\int_{-L}^0\cos kx\,dx=-\dfrac Lk. Second integral: L[1−(−1)n]k\dfrac{L[1-(-1)^n]}{k}. Hence

bn=1L⋅Lk(−1+1−(−1)n)=(−1)n+1Lnπ.b_n=\frac1L\cdot\frac{L}{k}\left(-1+1-(-1)^n\right)=\frac{(-1)^{n+1}L}{n\pi}.

Result

f(x)=3L4+∑n=1∞[L (1−(−1)n)n2π2cos⁡nπxL+(−1)n+1Lnπsin⁡nπxL]\boxed{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 (continuity point): 3L4+2Lπ2∑n odd1n2=3L4+2Lπ2⋅π28=L=f(0)\frac{3L}{4}+\frac{2L}{\pi^2}\sum_{n\ \text{odd}}\frac1{n^2}=\frac{3L}{4}+\frac{2L}{\pi^2}\cdot\frac{\pi^2}{8}=L=f(0). At the jump x=±Lx=\pm L the series converges to L/2L/2, the mean of LL and 00.


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

Case λ=0\lambda=0: y=A+Bxy=A+Bx. y′(0)=B=0y'(0)=B=0 and y(π)=A=0y(\pi)=A=0, so only y≡0y\equiv0. Not an eigenvalue.

Case λ=−μ2<0\lambda=-\mu^2<0: y=Acosh⁡μx+Bsinh⁡μxy=A\cosh\mu x+B\sinh\mu x. y′(0)=Bμ=0⇒B=0y'(0)=B\mu=0\Rightarrow B=0, and y(π)=Acosh⁡μπ=0⇒A=0y(\pi)=A\cosh\mu\pi=0\Rightarrow A=0. No eigenvalue.

Case λ=μ2>0\lambda=\mu^2>0: y=Acos⁡μx+Bsin⁡μxy=A\cos\mu x+B\sin\mu x, y′=−Aμsin⁡μx+Bμcos⁡μxy'=-A\mu\sin\mu x+B\mu\cos\mu x.

y′(0)=Bμ=0⇒B=0y'(0)=B\mu=0\Rightarrow B=0. Then y(π)=Acos⁡μπ=0y(\pi)=A\cos\mu\pi=0 with A≠0A\ne0 requires cos⁡μπ=0\cos\mu\pi=0:

μπ=(2n−1)π2 ⇒ μ=2n−12,n=1,2,3,…\mu\pi=\frac{(2n-1)\pi}{2}\ \Rightarrow\ \mu=\frac{2n-1}{2},\quad n=1,2,3,\dots

λn=(2n−12)2=(2n−1)24,yn(x)=cos⁡(2n−1)x2(n=1,2,… )\boxed{\lambda_n=\left(\frac{2n-1}{2}\right)^2=\frac{(2n-1)^2}{4},\qquad y_n(x)=\cos\frac{(2n-1)x}{2}}\quad(n=1,2,\dots)

The eigenvalues are 14,94,254,…\tfrac14,\tfrac94,\tfrac{25}{4},\dots, and each eigenfunction is determined up to a constant multiple.

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

Choosing the first alternative.

Let
[ u(x,y)=X(x)Y(y) ] be a separated solution of Laplace’s equation [ u_{xx}+u_{yy}=0. ] Substituting gives [ \frac{X''}{X}=-\frac{Y''}{Y}=-\lambda. ] Hence [ X''+\lambda X=0,\qquad Y''-\lambda Y=0. ]

The boundary conditions [ u(0,y)=0,\qquad u(a,y)=0 ] give [ X(0)=0,\qquad X(a)=0. ] Therefore the eigenfunctions and eigenvalues are [ X_n(x)=\sin\frac{n\pi x}{a},\qquad \lambda_n=\left(\frac{n\pi}{a}\right)^2,\qquad n=1,2,3,\dots ]

For (Y(y)), we have [ Y_n(y)=A_n\cosh\frac{n\pi y}{a}+B_n\sinh\frac{n\pi y}{a}. ] Using the condition [ u(x,b)=0 ] gives [ Y_n(b)=0. ] Thus we may write [ Y_n(y)=C_n\sinh\frac{n\pi(b-y)}{a}. ]

Hence the required solution has the form [ u(x,y)=\sum_{n=1}^{\infty}C_n\sin\frac{n\pi x}{a}\sinh\frac{n\pi(b-y)}{a}. ]

Now apply the remaining boundary condition [ u(x,0)=h(x). ] Then [ h(x)=\sum_{n=1}^{\infty}C_n\sinh\frac{n\pi b}{a}\sin\frac{n\pi x}{a}. ] This is the Fourier sine series of (h(x)) on (0<x<a). Therefore [ C_n\sinh\frac{n\pi b}{a}

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

Therefore the solution is [ \boxed{ u(x,y)= \sum_{n=1}^{\infty} \left[ \frac{2}{a\sinh\frac{n\pi b}{a}} \int_0^a h(\xi)\sin\frac{n\pi \xi}{a},d\xi \right] \sin\frac{n\pi x}{a} \sinh\frac{n\pi(b-y)}{a} } ] for (0<x<a,\ 0<y<b).


OR alternative: semi-circular region

In polar coordinates, Laplace’s equation is [ u_{rr}+\frac1r u_r+\frac1{r^2}u_{\theta\theta}=0. ] Let [ u(r,\theta)=R(r)\Theta(\theta). ] Then [ \Theta''+\lambda \Theta=0. ] Using [ u(r,0)=0,\qquad u(r,\pi)=0, ] we get [ \Theta(0)=0,\qquad \Theta(\pi)=0. ] Hence [ \Theta_n(\theta)=\sin n\theta,\qquad \lambda_n=n^2,\qquad n=1,2,3,\dots ]

The radial equation is [ r^2R''+rR'-n^2R=0, ] whose general solution is [ R_n(r)=A_nr^n+B_nr^{-n}. ] Since (u) is bounded at (r=0), we must have (B_n=0). Thus [ R_n(r)=A_nr^n. ]

Therefore [ u(r,\theta)=\sum_{n=1}^{\infty}A_nr^n\sin n\theta. ] Using the boundary condition [ u(a,\theta)=f(\theta), ] we get [ f(\theta)=\sum_{n=1}^{\infty}A_na^n\sin n\theta. ] This is the Fourier sine series of (f(\theta)) on (0<\theta<\pi). Hence [ A_na^n=\frac{2}{\pi}\int_0^\pi f(\theta)\sin n\theta,d\theta. ] So [ A_n=\frac{2}{\pi a^n}\int_0^\pi f(\theta)\sin n\theta,d\theta. ]

Thus the required solution is [ \boxed{ u(r,\theta)= \sum_{n=1}^{\infty} \left[ \frac{2}{\pi a^n} \int_0^\pi f(\phi)\sin n\phi,d\phi \right] r^n\sin n\theta } ] for (0\le r<a,\ 0\le \theta\le \pi).

Answer checked by Claude Haiku · Oct 21, 2026

Test yourself · 2 MCQs

Loading…

Questions from TU BSc 2nd Year 2082 Applied Statistics (BSc 2nd Year · Applied Statistics · 2082), TU 4 Yrs BSc 2nd Yr 2083 Regular (BSc 2nd year · Probability and Inference-I / Stat.201 · 2083), TU 4 Yrs BSc 2nd Yr 2083 Regular Basic (Bachelor Level (4 Yrs.)/Science & Tech./II Year · Basic Chemistry II (CHE 201) · 2083), TU BSc 2nd Yr 2083 Regular Zoology (BSc 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 (BSc 2nd year · 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 (BSc 2nd Year · Applied Statistics · 2083), TU 4 Yrs BSc 2nd Yr Differential (BSc 2nd year · Differential Equation (MAT 202) (New Course) · 2082), TU 4 Yrs BSc 2nd Yr Probability and (BSc 2nd year · Probability and Inference-I · 2082), TU 4 Yrs BSc 2nd Yr Linear Algebra 2082 (Bachelor Level (4 Yrs.) Science & Tech. II Year · Linear Algebra (MAT 201) (New Course) · 2082), TU 4 Years BSc 2nd Yr PHY 201 New (BSc 2nd year · PHY. 201 (New Course) · 2082), TU 4 Years BSc 2nd Yr Applied (BSc 2nd Year · Applied Statistics (New Course) · 2082), TU 4 Years BSc 2nd Yr Probability and (BSc 2nd Year · Probability and Inference - I (Stat. 201) · 2081), TU 4 Years BSc 2nd Yr PHY 201 2081 (Bachelor Level (4 Yrs.)/ Science & Tech. / II Year · PHY. 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 (BSc 2nd year · Applied Statistics · 2081), TU 4 Years BSc 2nd Yr Basic Chemistry (BSc 2nd year · Basic Chemistry II (CHE 201) · 2081) and TU 4 Years BSc 2nd Yr Basic Chemistry (BSc 2nd Year · General Chemistry II (CHE 201) · 2078). Answers are written for this site; check them against your class notes.