You may use Gamma or Beta function to express the final answer. Please show intermediate steps, otherwise there are no points.
1. Compute the integral (5 points) ∫011−x3x4dx
Replace variable u=x3, then
I=∫01(1−u)1/2u4/3du1/3=(1/3)∫01(1−u)−1/2u2/3du=(1/3)∫01(1−u)−1/2u2/3du=(1/3)B(5/3,1/2).
2. Compute the integral (5 points) ∫01x2(−lnx)2dx
Replace variable u=−lnx, then
∫0∞e−2uu2e−udu=(1/3)3Γ(3).
3. Compute the special value of Gamma functions (the result should not be expressed using Gamma function).
Γ(−1/2)=Γ(1/2)/(1/2)=2π (2 points)
∣Γ(1/2+i)∣=? (3 points) [Hint: Γ(zˉ)=Γ(z), and Γ(p)Γ(1−p)=π/sin(pπ) ]
Take p=1/2+i, then
∣Γ(1/2+i)∣2=Γ(1/2+i)Γ(1/2−i)=π/sin(π(1/2+i))=π/cos(iπ)=π/cosh(pi)
2. Legendre Function (25 points)
1. Compute P3(x) using Rodrigue formula. (5 points)
By Rodrigue formula
P3(x)=233!1∂x3(x2−1)3=233!1∂x3(x2−1)3=233!1∂x3(x6−3x4+3x2−1)=481(120x3−72x)=(5x3−3x)/2.
2. Prove the recursion relation 5.8( c) (10 points)
Pl′(x)−xPl−1′(x)=lPl−1(x)
using the generating function
Φ(x,h)=n=0∑∞hnPn(x)=1−2xh+h21
If we multiply both sides of the equation by hl and sum over l, we get
LHS=∂xΦ(x,h)−xh∂xΦ(x,h)RHS=h∂h(hΦ(x,h))
We plug in Φ(x,h) to test whether this is true.
LHS=(1−xh)(1−2xh+h2)3/2h
and
RHS=hΦ+h2(−1/2)(−2x+2h)(1−2xh+h2)−3/2=[h(1−2xh+h2)+h2(x−h)](1−2xh+h2)−3/2=h(1−hx)(1−2xh+h2)−3/2
OK, the same, done.
3. Compute ∫−11xnPn(x)dx in the following steps: ( 10 points)
Find the constant c, such that Pn(x)=cxn+ lower order terms. (try Rodrigue formula to get the leading term)
Show that ∫−11Pn(x)2dx=∫−11cxnPn(x)dx
Look up ∫−11Pn(x)Pn(x)dx.
The leading term of Pn(x) can be obtained as
2nn!1(2n)(2n−1)⋯(n+1)=2−n(n2n)
Because lower order term is orthogonal to Pn(x).
∫−11Pn(x)Pn(x)dx=2n+12
So finally, we get (2n+1)(2n)!2n+1(n!)2
3. Bessel Function (30 points)
1. Problem 12.1. Show by ratio test that the series for Jp(x) converges for all x. (10 points)
2. Problem 15.6 (7 points)
3. Problem 19.1 (7 points)
4. Problem 20.3, 6, 7 (6 points)
4. Solving PDE with separation of variables (30 points)
1. Solve the steady state heat equation on 2D square [0,1]2. (10 point)
Δu(x,y)=0,0≤x,y≤1
with boundary condition that
u(0,y)=0,u(1,y)=1,u(x,0)=0,u(x,1)=1.