Skip to content

STEP2 2000 -- Pure Mathematics

STEP2 2000 — Section A (Pure Mathematics)

Section titled “STEP2 2000 — Section A (Pure Mathematics)”

Exam: STEP2  |  Year: 2000  |  Questions: Q1—Q8  |  Total marks per question: 20

All questions are pure mathematics. Solutions and examiner commentary are included below.

QTopicDifficultyKey Techniques
1数论 Number TheoryStandard因数分解,代数变换,唯一性论证
2代数 AlgebraStandard多项式导数,重根条件,方程组求解
3微积分 CalculusStandard偏微分,线性近似,误差传播
4复数 Complex NumbersChallengingDe Moivre定理,复数乘法,辐角原理
5微积分 CalculusStandard最小二乘法,积分求极值,极限分析
6积分 IntegrationChallenging万能代换,三角恒等式,定积分计算
7向量几何 Vector GeometryChallenging向量方程,直线夹角公式,平面方程,轨迹分析
8微分方程 Differential EquationsChallenging分离变量法,极限分析,指数函数

Topic: 数论 Number Theory  |  Difficulty: Standard  |  Marks: 20

1 A number of the form 1/N1/N, where NN is an integer greater than 1, is called a unit fraction.

Noting that

12=13+16and13=14+112,\frac{1}{2} = \frac{1}{3} + \frac{1}{6} \quad \text{and} \quad \frac{1}{3} = \frac{1}{4} + \frac{1}{12},

guess a general result of the form

1N=1a+1b(*)\frac{1}{N} = \frac{1}{a} + \frac{1}{b} \qquad \text{(*)}

and hence prove that any unit fraction can be expressed as the sum of two distinct unit fractions.

By writing (*) in the form

(aN)(bN)=N2(a - N)(b - N) = N^2

and by considering the factors of N2N^2, show that if NN is prime, then there is only one way of expressing 1/N1/N as the sum of two distinct unit fractions.

Prove similarly that any fraction of the form 2/N2/N, where NN is prime number greater than 2, can be expressed uniquely as the sum of two distinct unit fractions.

Model Solution

Guessing the general result. From the two examples:

N=2:a=3=2+1,b=6=2×3N = 2: \quad a = 3 = 2+1, \quad b = 6 = 2 \times 3

N=3:a=4=3+1,b=12=3×4N = 3: \quad a = 4 = 3+1, \quad b = 12 = 3 \times 4

the pattern suggests:

1N=1N+1+1N(N+1)(**)\frac{1}{N} = \frac{1}{N+1} + \frac{1}{N(N+1)} \qquad \text{(**)}

Proof. Compute the right-hand side with a common denominator:

1N+1+1N(N+1)=NN(N+1)+1N(N+1)=N+1N(N+1)=1N\frac{1}{N+1} + \frac{1}{N(N+1)} = \frac{N}{N(N+1)} + \frac{1}{N(N+1)} = \frac{N+1}{N(N+1)} = \frac{1}{N}

Since N>1N > 1, we have N+1N(N+1)N + 1 \neq N(N+1), so the two unit fractions are distinct. This proves that any unit fraction can be expressed as the sum of two distinct unit fractions. \square


Uniqueness when NN is prime. Starting from 1N=1a+1b\frac{1}{N} = \frac{1}{a} + \frac{1}{b} with aba \neq b, multiply through by NabNab:

ab=N(a+b)ab = N(a + b)

abNaNb=0ab - Na - Nb = 0

Add N2N^2 to both sides:

abNaNb+N2=N2ab - Na - Nb + N^2 = N^2

(aN)(bN)=N2(a - N)(b - N) = N^2

Since 1a<1N\frac{1}{a} < \frac{1}{N} (because 1a=1N1b\frac{1}{a} = \frac{1}{N} - \frac{1}{b} and b>0b > 0), we have a>Na > N, so aN>0a - N > 0. Similarly b>Nb > N, so bN>0b - N > 0.

If NN is prime, the positive divisors of N2N^2 are 1,N,N21, N, N^2. The factor pairs (d1,d2)(d_1, d_2) with d1d2=N2d_1 d_2 = N^2 and d1d2d_1 \leqslant d_2 are:

  • (1,N2)(1, N^2): gives a=N+1a = N + 1, b=N2+N=N(N+1)b = N^2 + N = N(N+1). These are distinct since N>1N > 1.
  • (N,N)(N, N): gives a=b=2Na = b = 2N. This is excluded since we require distinct fractions.

(The pair (N2,1)(N^2, 1) gives the same solution with aa and bb swapped.)

Therefore, when NN is prime, the only expression of 1/N1/N as the sum of two distinct unit fractions (up to order) is:

1N=1N+1+1N(N+1)\frac{1}{N} = \frac{1}{N+1} + \frac{1}{N(N+1)}

\square


The case 2/N2/N with NN an odd prime. Suppose 2N=1a+1b\frac{2}{N} = \frac{1}{a} + \frac{1}{b} with aba \neq b. Multiply through by NabNab:

2ab=N(a+b)2ab = N(a + b)

4ab=2N(a+b)4ab = 2N(a + b)

4ab2Na2Nb=04ab - 2Na - 2Nb = 0

Add N2N^2 to both sides:

4ab2Na2Nb+N2=N24ab - 2Na - 2Nb + N^2 = N^2

(2aN)(2bN)=N2(2a - N)(2b - N) = N^2

Since 1a<2N\frac{1}{a} < \frac{2}{N}, we have a>N/2a > N/2, so 2aN>02a - N > 0. Similarly 2bN>02b - N > 0.

Since NN is an odd prime (N>2N > 2), the positive factor pairs of N2N^2 are the same as before: (1,N2)(1, N^2), (N,N)(N, N), (N2,1)(N^2, 1).

  • (N,N)(N, N): 2aN=N2a - N = N and 2bN=N2b - N = N, giving a=b=Na = b = N. Not distinct, so excluded.
  • (1,N2)(1, N^2): 2a=N+12a = N + 1 and 2b=N2+N=N(N+1)2b = N^2 + N = N(N+1). Since NN is odd, N+1N + 1 is even, so a=(N+1)/2a = (N+1)/2 is a positive integer. Since N+1N + 1 is even, N(N+1)N(N+1) is also even, so b=N(N+1)/2b = N(N+1)/2 is a positive integer. And aba \neq b since N>1N > 1.

(The pair (N2,1)(N^2, 1) gives the same solution with aa and bb swapped.)

Therefore 2/N2/N can be expressed uniquely (up to order) as:

2N=1(N+1)/2+1N(N+1)/2\frac{2}{N} = \frac{1}{(N+1)/2} + \frac{1}{N(N+1)/2}

Verification:

1(N+1)/2+1N(N+1)/2=2N+1+2N(N+1)=2N+2N(N+1)=2(N+1)N(N+1)=2N\frac{1}{(N+1)/2} + \frac{1}{N(N+1)/2} = \frac{2}{N+1} + \frac{2}{N(N+1)} = \frac{2N + 2}{N(N+1)} = \frac{2(N+1)}{N(N+1)} = \frac{2}{N} \checkmark

\square


Topic: 代数 Algebra  |  Difficulty: Standard  |  Marks: 20

2 Prove that if (xa)2(x - a)^2 is a factor of the polynomial p(x)p(x), then p(a)=0p'(a) = 0. Prove a corresponding result if (xa)4(x - a)^4 is a factor of p(x)p(x).

Given that the polynomial

x6+4x55x440x340x2+32x+kx^6 + 4x^5 - 5x^4 - 40x^3 - 40x^2 + 32x + k

has a factor of the form (xa)4(x - a)^4, find kk.

Model Solution

Proof for (xa)2(x - a)^2. If (xa)2(x - a)^2 is a factor of p(x)p(x), then p(x)=(xa)2q(x)p(x) = (x - a)^2 q(x) for some polynomial q(x)q(x). Differentiating by the product rule:

p(x)=2(xa)q(x)+(xa)2q(x)p'(x) = 2(x - a) \, q(x) + (x - a)^2 \, q'(x)

Evaluating at x=ax = a:

p(a)=2(0)q(a)+02q(a)=0p'(a) = 2(0) \cdot q(a) + 0^2 \cdot q'(a) = 0 \qquad \square


Corresponding result for (xa)4(x - a)^4. If (xa)4(x - a)^4 is a factor of p(x)p(x), then p(x)=(xa)4q(x)p(x) = (x - a)^4 q(x) for some polynomial q(x)q(x). Differentiate:

p(x)=4(xa)3q(x)+(xa)4q(x)=(xa)3[4q(x)+(xa)q(x)]s(x)p'(x) = 4(x - a)^3 \, q(x) + (x - a)^4 \, q'(x) = (x - a)^3 \underbrace{\left[ 4q(x) + (x - a) q'(x) \right]}_{s(x)}

So p(x)=(xa)3s(x)p'(x) = (x - a)^3 s(x), which means (xa)2(x - a)^2 is a factor of p(x)p'(x). By the result just proved (applied to pp' instead of pp):

p(a)=0p''(a) = 0

Now differentiate p(x)=(xa)3s(x)p'(x) = (x - a)^3 s(x):

p(x)=3(xa)2s(x)+(xa)3s(x)=(xa)2[3s(x)+(xa)s(x)]t(x)p''(x) = 3(x - a)^2 \, s(x) + (x - a)^3 \, s'(x) = (x - a)^2 \underbrace{\left[ 3s(x) + (x - a) s'(x) \right]}_{t(x)}

So p(x)=(xa)2t(x)p''(x) = (x - a)^2 t(x), meaning (xa)2(x - a)^2 is a factor of p(x)p''(x). Applying the same result to pp'':

p(a)=0p'''(a) = 0

Conclusion: if (xa)4(x - a)^4 is a factor of p(x)p(x), then p(a)=p(a)=p(a)=0p'(a) = p''(a) = p'''(a) = 0. \square


Finding kk. Let p(x)=x6+4x55x440x340x2+32x+kp(x) = x^6 + 4x^5 - 5x^4 - 40x^3 - 40x^2 + 32x + k. Compute the first three derivatives:

p(x)=6x5+20x420x3120x280x+32p'(x) = 6x^5 + 20x^4 - 20x^3 - 120x^2 - 80x + 32

p(x)=30x4+80x360x2240x80p''(x) = 30x^4 + 80x^3 - 60x^2 - 240x - 80

p(x)=120x3+240x2120x240p'''(x) = 120x^3 + 240x^2 - 120x - 240

Step 1. Solve p(a)=0p'''(a) = 0:

120(a3+2a2a2)=0120(a^3 + 2a^2 - a - 2) = 0

Factor by grouping:

a2(a+2)1(a+2)=0    (a+2)(a21)=0    (a+2)(a1)(a+1)=0a^2(a + 2) - 1(a + 2) = 0 \implies (a + 2)(a^2 - 1) = 0 \implies (a + 2)(a - 1)(a + 1) = 0

So the candidates are a=2,1,1a = -2, \, -1, \, 1.

Step 2. Check which candidates satisfy p(a)=0p''(a) = 0:

p(2)=30(16)+80(8)60(4)240(2)80=480640240+48080=0p''(-2) = 30(16) + 80(-8) - 60(4) - 240(-2) - 80 = 480 - 640 - 240 + 480 - 80 = 0 \checkmark

p(1)=308060+24080=500p''(-1) = 30 - 80 - 60 + 240 - 80 = 50 \neq 0

p(1)=30+806024080=2700p''(1) = 30 + 80 - 60 - 240 - 80 = -270 \neq 0

So a=2a = -2.

Step 3. Verify p(2)=0p'(-2) = 0:

p(2)=6(32)+20(16)20(8)120(4)80(2)+32p'(-2) = 6(-32) + 20(16) - 20(-8) - 120(4) - 80(-2) + 32 =192+320+160480+160+32=0= -192 + 320 + 160 - 480 + 160 + 32 = 0 \checkmark

Step 4. Find kk from p(2)=0p(-2) = 0:

p(2)=6412880+32016064+k=48+k=0p(-2) = 64 - 128 - 80 + 320 - 160 - 64 + k = -48 + k = 0

k=48\boxed{k = 48}

Verification. We check that (x+2)4(x24x+3)=x6+4x55x440x340x2+32x+48(x + 2)^4 (x^2 - 4x + 3) = x^6 + 4x^5 - 5x^4 - 40x^3 - 40x^2 + 32x + 48.

(x+2)4=x4+8x3+24x2+32x+16(x + 2)^4 = x^4 + 8x^3 + 24x^2 + 32x + 16 and x24x+3x^2 - 4x + 3. Multiplying term by term:

Term from (x+2)4(x+2)^4×x2\times \, x^2×(4x)\times \, (-4x)×3\times \, 3
x4x^4x6x^64x5-4x^53x43x^4
8x38x^38x58x^532x4-32x^424x324x^3
24x224x^224x424x^496x3-96x^372x272x^2
32x32x32x332x^3128x2-128x^296x96x
161616x216x^264x-64x4848

Collecting:

  • x6x^6: 11
  • x5x^5: 4+8=4-4 + 8 = 4
  • x4x^4: 332+24=53 - 32 + 24 = -5
  • x3x^3: 2496+32=4024 - 96 + 32 = -40
  • x2x^2: 72128+16=4072 - 128 + 16 = -40
  • xx: 9664=3296 - 64 = 32
  • constant: 4848

This matches the original polynomial with k=48k = 48. Note also that x24x+3=(x1)(x3)x^2 - 4x + 3 = (x-1)(x-3), so the full factorisation is:

p(x)=(x+2)4(x1)(x3)p(x) = (x + 2)^4 (x - 1)(x - 3)


Topic: 微积分 Calculus  |  Difficulty: Standard  |  Marks: 20

3 The lengths of the sides BC,CA,ABBC, CA, AB of the triangle ABCABC are denoted by a,b,ca, b, c, respectively. Given that

b=8+ϵ1,c=3+ϵ2,A=π/3+ϵ3,b = 8 + \epsilon_1, c = 3 + \epsilon_2, A = \pi/3 + \epsilon_3,

where ϵ1,ϵ2\epsilon_1, \epsilon_2, and ϵ3\epsilon_3 are small, show that a7+ηa \approx 7 + \eta, where η=(13ϵ12ϵ2+243ϵ3)/14\eta = (13 \, \epsilon_1 - 2 \, \epsilon_2 + 24\sqrt{3} \, \epsilon_3)/14.

Given now that

ϵ12×103,ϵ249×102,ϵ33×103,|\epsilon_1| \leqslant 2 \times 10^{-3}, \quad |\epsilon_2| \leqslant 4 \cdot 9 \times 10^{-2}, \quad |\epsilon_3| \leqslant \sqrt{3} \times 10^{-3},

find the range of possible values of η\eta.

Model Solution

Part (a): Showing a7+ηa \approx 7 + \eta

By the cosine rule,

a2=b2+c22bccosA.(*)a^2 = b^2 + c^2 - 2bc\cos A. \qquad \text{(*)}

When ϵ1=ϵ2=ϵ3=0\epsilon_1 = \epsilon_2 = \epsilon_3 = 0, we have b=8b = 8, c=3c = 3, A=π/3A = \pi/3, so

a2=64+948cos(π/3)=7324=49,a^2 = 64 + 9 - 48\cos(\pi/3) = 73 - 24 = 49,

giving a=7a = 7.

To find the first-order perturbation, we differentiate (*). Since aa, bb, cc, and AA all vary:

2ada=2bdb+2cdc2cosA(cdb+bdc)2bc(sinA)dA.2a\,da = 2b\,db + 2c\,dc - 2\cos A\,(c\,db + b\,dc) - 2bc(-\sin A)\,dA.

Note the last term: A(2bccosA)=2bcsinA\frac{\partial}{\partial A}(-2bc\cos A) = 2bc\sin A. Dividing by 2:

ada=bdb+cdccosA(cdb+bdc)+bcsinAdA.a\,da = b\,db + c\,dc - \cos A\,(c\,db + b\,dc) + bc\sin A\,dA.

Collecting terms in dbdb, dcdc, dAdA:

ada=(bccosA)db+(cbcosA)dc+bcsinAdA.a\,da = (b - c\cos A)\,db + (c - b\cos A)\,dc + bc\sin A\,dA.

Now substitute db=dϵ1db = d\epsilon_1, dc=dϵ2dc = d\epsilon_2, dA=dϵ3dA = d\epsilon_3, and evaluate at b=8b = 8, c=3c = 3, A=π/3A = \pi/3, a=7a = 7:

  • bccosA=8312=132b - c\cos A = 8 - 3 \cdot \tfrac{1}{2} = \tfrac{13}{2}
  • cbcosA=3812=1c - b\cos A = 3 - 8 \cdot \tfrac{1}{2} = -1
  • bcsinA=2432=123bc\sin A = 24 \cdot \tfrac{\sqrt{3}}{2} = 12\sqrt{3}

Therefore

7dη=132dϵ1dϵ2+123dϵ3,7\,d\eta = \tfrac{13}{2}\,d\epsilon_1 - d\epsilon_2 + 12\sqrt{3}\,d\epsilon_3,

which integrates (to first order) to

η=13ϵ12ϵ2+243ϵ314.\eta = \frac{13\epsilon_1 - 2\epsilon_2 + 24\sqrt{3}\,\epsilon_3}{14}.

Hence a7+ηa \approx 7 + \eta where η=(13ϵ12ϵ2+243ϵ3)/14\eta = (13\epsilon_1 - 2\epsilon_2 + 24\sqrt{3}\,\epsilon_3)/14.

Part (b): Finding the range of η\eta

Since η=(13ϵ12ϵ2+243ϵ3)/14\eta = (13\epsilon_1 - 2\epsilon_2 + 24\sqrt{3}\,\epsilon_3)/14 is linear in the ϵi\epsilon_i, its maximum and minimum occur at the extreme values of each ϵi\epsilon_i.

The maximum of η|\eta| is achieved when each term 13ϵ113\epsilon_1, 2ϵ2-2\epsilon_2, and 243ϵ324\sqrt{3}\,\epsilon_3 has the same sign. By the triangle inequality:

η13ϵ1+2ϵ2+243ϵ314.|\eta| \leqslant \frac{13|\epsilon_1| + 2|\epsilon_2| + 24\sqrt{3}\,|\epsilon_3|}{14}.

Substituting the given bounds:

η13×2×103+2×4.9×102+243×3×10314.|\eta| \leqslant \frac{13 \times 2 \times 10^{-3} + 2 \times 4.9 \times 10^{-2} + 24\sqrt{3} \times \sqrt{3} \times 10^{-3}}{14}.

Evaluating each term:

  • 13×2×103=26×10313 \times 2 \times 10^{-3} = 26 \times 10^{-3}
  • 2×4.9×102=98×1032 \times 4.9 \times 10^{-2} = 98 \times 10^{-3}
  • 243×3×103=24×3×103=72×10324\sqrt{3} \times \sqrt{3} \times 10^{-3} = 24 \times 3 \times 10^{-3} = 72 \times 10^{-3}

Therefore

η(26+98+72)×10314=196×10314=14×103=0.014.|\eta| \leqslant \frac{(26 + 98 + 72) \times 10^{-3}}{14} = \frac{196 \times 10^{-3}}{14} = 14 \times 10^{-3} = 0.014.

This bound is attained (for example, when ϵ1=2×103\epsilon_1 = 2 \times 10^{-3}, ϵ2=4.9×102\epsilon_2 = -4.9 \times 10^{-2}, ϵ3=3×103\epsilon_3 = \sqrt{3} \times 10^{-3}), so

0.014η0.014.-0.014 \leqslant \eta \leqslant 0.014.


Topic: 复数 Complex Numbers  |  Difficulty: Challenging  |  Marks: 20

4 Prove that (cosθ+isinθ)(cosϕ+isinϕ)=cos(θ+ϕ)+isin(θ+ϕ)(\cos \theta + i \sin \theta)(\cos \phi + i \sin \phi) = \cos(\theta + \phi) + i \sin(\theta + \phi) and that, for every positive integer nn, (cosθ+isinθ)n=cosnθ+isinnθ.(\cos \theta + i \sin \theta)^n = \cos n\theta + i \sin n\theta. By considering (5i)2(1+i)(5 - i)^2(1 + i), or otherwise, prove that arctan(7/17)+2arctan(1/5)=π/4.\arctan(7/17) + 2 \arctan(1/5) = \pi/4. Prove also that 3arctan(1/4)+arctan(1/20)+arctan(1/1985)=π/4.3 \arctan(1/4) + \arctan(1/20) + \arctan(1/1985) = \pi/4. [Note that arctanθ\arctan \theta is another notation for tan1θ\tan^{-1} \theta.]

Model Solution

Part (a): Proving the multiplication formula

(cosθ+isinθ)(cosϕ+isinϕ)(\cos\theta + i\sin\theta)(\cos\phi + i\sin\phi) =cosθcosϕ+icosθsinϕ+isinθcosϕ+i2sinθsinϕ= \cos\theta\cos\phi + i\cos\theta\sin\phi + i\sin\theta\cos\phi + i^2\sin\theta\sin\phi =(cosθcosϕsinθsinϕ)+i(sinθcosϕ+cosθsinϕ).= (\cos\theta\cos\phi - \sin\theta\sin\phi) + i(\sin\theta\cos\phi + \cos\theta\sin\phi).

By the angle addition formulae cos(θ+ϕ)=cosθcosϕsinθsinϕ\cos(\theta + \phi) = \cos\theta\cos\phi - \sin\theta\sin\phi and sin(θ+ϕ)=sinθcosϕ+cosθsinϕ\sin(\theta + \phi) = \sin\theta\cos\phi + \cos\theta\sin\phi, this equals

cos(θ+ϕ)+isin(θ+ϕ).(*)\cos(\theta + \phi) + i\sin(\theta + \phi). \qquad \text{(*)}

Part (b): Proving de Moivre’s theorem by induction

We prove (cosθ+isinθ)n=cosnθ+isinnθ(\cos\theta + i\sin\theta)^n = \cos n\theta + i\sin n\theta for all positive integers nn.

Base case (n=1n = 1): (cosθ+isinθ)1=cosθ+isinθ(\cos\theta + i\sin\theta)^1 = \cos\theta + i\sin\theta. True.

Inductive step: Assume the result holds for n=kn = k, i.e., (cosθ+isinθ)k=coskθ+isinkθ(\cos\theta + i\sin\theta)^k = \cos k\theta + i\sin k\theta. Then

(cosθ+isinθ)k+1=(cosθ+isinθ)k(cosθ+isinθ)(\cos\theta + i\sin\theta)^{k+1} = (\cos\theta + i\sin\theta)^k \cdot (\cos\theta + i\sin\theta) =(coskθ+isinkθ)(cosθ+isinθ).= (\cos k\theta + i\sin k\theta)(\cos\theta + i\sin\theta).

Applying (*) with ϕ=θ\phi = \theta:

=cos(kθ+θ)+isin(kθ+θ)=cos(k+1)θ+isin(k+1)θ.= \cos(k\theta + \theta) + i\sin(k\theta + \theta) = \cos(k+1)\theta + i\sin(k+1)\theta.

By induction, the result holds for all positive integers nn.

Part (c): Proving arctan(7/17)+2arctan(1/5)=π/4\arctan(7/17) + 2\arctan(1/5) = \pi/4

Consider the complex number (5i)2(1+i)(5 - i)^2(1 + i).

First, compute (5i)2(5 - i)^2:

(5i)2=2510i+i2=2410i.(5 - i)^2 = 25 - 10i + i^2 = 24 - 10i.

Then multiply by (1+i)(1 + i):

(2410i)(1+i)=24+24i10i10i2=24+14i+10=34+14i.(24 - 10i)(1 + i) = 24 + 24i - 10i - 10i^2 = 24 + 14i + 10 = 34 + 14i.

Now we find the argument of this product in two ways.

Method 1 (direct): Since 34+14i34 + 14i has positive real and imaginary parts,

arg(34+14i)=arctan ⁣(1434)=arctan ⁣(717).\arg(34 + 14i) = \arctan\!\left(\frac{14}{34}\right) = \arctan\!\left(\frac{7}{17}\right).

Method 2 (via arguments of factors): We have arg(1+i)=π/4\arg(1 + i) = \pi/4. Also 5i5 - i has positive real part and negative imaginary part, so arg(5i)=arctan(1/5)\arg(5 - i) = -\arctan(1/5). By de Moivre’s theorem:

arg((5i)2(1+i))=2arg(5i)+arg(1+i)=2arctan(1/5)+π/4.\arg\bigl((5-i)^2(1+i)\bigr) = 2\arg(5-i) + \arg(1+i) = -2\arctan(1/5) + \pi/4.

Equating the two expressions:

arctan(7/17)=2arctan(1/5)+π/4,\arctan(7/17) = -2\arctan(1/5) + \pi/4,

hence

arctan(7/17)+2arctan(1/5)=π/4.\arctan(7/17) + 2\arctan(1/5) = \pi/4. \qquad \blacksquare

Part (d): Proving 3arctan(1/4)+arctan(1/20)+arctan(1/1985)=π/43\arctan(1/4) + \arctan(1/20) + \arctan(1/1985) = \pi/4

Consider the product (4+i)3(20+i)(1985+i)(4 + i)^3(20 + i)(1985 + i).

Step 1: Compute (4+i)3(4 + i)^3.

(4+i)2=16+8i+i2=15+8i.(4 + i)^2 = 16 + 8i + i^2 = 15 + 8i.

(4+i)3=(15+8i)(4+i)=60+15i+32i+8i2=52+47i.(4 + i)^3 = (15 + 8i)(4 + i) = 60 + 15i + 32i + 8i^2 = 52 + 47i.

Step 2: Multiply by (20+i)(20 + i).

(52+47i)(20+i)=1040+52i+940i+47i2=993+992i.(52 + 47i)(20 + i) = 1040 + 52i + 940i + 47i^2 = 993 + 992i.

Step 3: Multiply by (1985+i)(1985 + i).

(993+992i)(1985+i)(993 + 992i)(1985 + i) =993×1985+993i+992×1985i+992i2= 993 \times 1985 + 993i + 992 \times 1985\,i + 992i^2 =(993×1985992)+(993+992×1985)i.= (993 \times 1985 - 992) + (993 + 992 \times 1985)i.

Evaluating:

  • Real part: 993×1985992=1971105992=1970113993 \times 1985 - 992 = 1971105 - 992 = 1970113.
  • Imaginary part: 993+992×1985=993+1969120=1970113993 + 992 \times 1985 = 993 + 1969120 = 1970113.

So

(4+i)3(20+i)(1985+i)=1970113(1+i).(4 + i)^3(20 + i)(1985 + i) = 1970113(1 + i).

Since arg(1+i)=π/4\arg(1 + i) = \pi/4, we have

arg((4+i)3(20+i)(1985+i))=π/4.\arg\bigl((4+i)^3(20+i)(1985+i)\bigr) = \pi/4.

On the other hand, since each factor 4+i4 + i, 20+i20 + i, 1985+i1985 + i has positive real and imaginary parts, their arguments lie in (0,π/2)(0, \pi/2), and

arg(4+i)=arctan(1/4),arg(20+i)=arctan(1/20),arg(1985+i)=arctan(1/1985).\arg(4+i) = \arctan(1/4), \quad \arg(20+i) = \arctan(1/20), \quad \arg(1985+i) = \arctan(1/1985).

The sum 3arctan(1/4)+arctan(1/20)+arctan(1/1985)<3π/4+π/2<2π3\arctan(1/4) + \arctan(1/20) + \arctan(1/1985) < 3 \cdot \pi/4 + \pi/2 < 2\pi, so there is no ambiguity in adding the arguments. Therefore

3arctan(1/4)+arctan(1/20)+arctan(1/1985)=π/4.3\arctan(1/4) + \arctan(1/20) + \arctan(1/1985) = \pi/4. \qquad \blacksquare


Topic: 微积分 Calculus  |  Difficulty: Standard  |  Marks: 20

5 It is required to approximate a given function f(x)f(x), over the interval 0x10 \leqslant x \leqslant 1, by the linear function λx\lambda x, where λ\lambda is chosen to minimise 01(f(x)λx)2dx.\int_{0}^{1} \left( f(x) - \lambda x \right)^2 dx. Show that λ=301xf(x)dx.\lambda = 3 \int_{0}^{1} xf(x) dx. The residual error, RR, of this approximation process is such that R2=01(f(x)λx)2dx.R^2 = \int_{0}^{1} \left( f(x) - \lambda x \right)^2 dx. Show that R2=01(f(x))2dx13λ2.R^2 = \int_{0}^{1} \left( f(x) \right)^2 dx - \frac{1}{3} \lambda^2. Given now that f(x)=sin(πx/n)f(x) = \sin(\pi x/n), show that (i) for large nn, λπ/n\lambda \approx \pi/n and (ii) limnR=0\lim_{n \to \infty} R = 0. Explain why, prior to any calculation, these results are to be expected. [You may assume that, when θ\theta is small, sinθθθ3/6\sin \theta \approx \theta - \theta^3/6 and cosθ1θ2/2\cos \theta \approx 1 - \theta^2/2.]

Model Solution

Step 1: Find λ\lambda.

Define I(λ)=01(f(x)λx)2dxI(\lambda) = \displaystyle\int_{0}^{1}(f(x) - \lambda x)^2\,dx. Expand the integrand:

I(λ)=01f(x)2dx2λ01xf(x)dx+λ201x2dxI(\lambda) = \int_{0}^{1}f(x)^2\,dx - 2\lambda\int_{0}^{1}xf(x)\,dx + \lambda^2\int_{0}^{1}x^2\,dx

Evaluating the last integral:

01x2dx=[x33]01=13\int_{0}^{1}x^2\,dx = \left[\frac{x^3}{3}\right]_0^1 = \frac{1}{3}

So

I(λ)=01f(x)2dx2λ01xf(x)dx+λ23I(\lambda) = \int_{0}^{1}f(x)^2\,dx - 2\lambda\int_{0}^{1}xf(x)\,dx + \frac{\lambda^2}{3}

This is a quadratic in λ\lambda with positive leading coefficient, so it has a unique minimum. Differentiate and set to zero:

dIdλ=201xf(x)dx+2λ3=0\frac{dI}{d\lambda} = -2\int_{0}^{1}xf(x)\,dx + \frac{2\lambda}{3} = 0

λ=301xf(x)dx(as required)\therefore\quad \lambda = 3\int_{0}^{1}xf(x)\,dx \qquad \text{(as required)}

Step 2: Show R2=01f(x)2dx13λ2R^2 = \displaystyle\int_0^1 f(x)^2\,dx - \tfrac{1}{3}\lambda^2.

From Step 1, we have λ=301xf(x)dx\lambda = 3\displaystyle\int_0^1 xf(x)\,dx, so 01xf(x)dx=λ3\displaystyle\int_0^1 xf(x)\,dx = \frac{\lambda}{3}.

Substituting back:

R2=I(λ)=01f(x)2dx2λλ3+λ23R^2 = I(\lambda) = \int_{0}^{1}f(x)^2\,dx - 2\lambda\cdot\frac{\lambda}{3} + \frac{\lambda^2}{3}

=01f(x)2dx2λ23+λ23=01f(x)2dxλ23(as required)= \int_{0}^{1}f(x)^2\,dx - \frac{2\lambda^2}{3} + \frac{\lambda^2}{3} = \int_{0}^{1}f(x)^2\,dx - \frac{\lambda^2}{3} \qquad \text{(as required)}

Step 3: Now let f(x)=sin(πx/n)f(x) = \sin(\pi x/n). Part (i) — show λπ/n\lambda \approx \pi/n for large nn.

λ=301xsin ⁣(πxn)dx\lambda = 3\int_{0}^{1}x\sin\!\left(\frac{\pi x}{n}\right)dx

For large nn, use sinθθθ3/6\sin\theta \approx \theta - \theta^3/6 with θ=πx/n\theta = \pi x/n:

λ301x(πxnπ3x36n3)dx=301(πx2nπ3x46n3)dx\lambda \approx 3\int_{0}^{1}x\left(\frac{\pi x}{n} - \frac{\pi^3 x^3}{6n^3}\right)dx = 3\int_{0}^{1}\left(\frac{\pi x^2}{n} - \frac{\pi^3 x^4}{6n^3}\right)dx

=3[πx33nπ3x530n3]01=3(π3nπ330n3)=πnπ310n3= 3\left[\frac{\pi x^3}{3n} - \frac{\pi^3 x^5}{30n^3}\right]_0^1 = 3\left(\frac{\pi}{3n} - \frac{\pi^3}{30n^3}\right) = \frac{\pi}{n} - \frac{\pi^3}{10n^3}

For large nn, the second term is negligible, so λπn\lambda \approx \dfrac{\pi}{n}.

Step 4: Part (ii) — show limnR=0\lim_{n\to\infty}R = 0.

Compute 01f(x)2dx=01sin2 ⁣(πxn)dx\displaystyle\int_0^1 f(x)^2\,dx = \int_0^1 \sin^2\!\left(\frac{\pi x}{n}\right)dx.

Using sin2u=12(1cos2u)\sin^2 u = \tfrac{1}{2}(1 - \cos 2u) and substituting u=πx/nu = \pi x/n:

01sin2 ⁣(πxn)dx=n2π0π/n(1cos2u)du=n2π[usin2u2]0π/n\int_0^1 \sin^2\!\left(\frac{\pi x}{n}\right)dx = \frac{n}{2\pi}\int_0^{\pi/n}(1 - \cos 2u)\,du = \frac{n}{2\pi}\left[u - \frac{\sin 2u}{2}\right]_0^{\pi/n}

=n2π(πnsin(2π/n)2)=12nsin(2π/n)4π= \frac{n}{2\pi}\left(\frac{\pi}{n} - \frac{\sin(2\pi/n)}{2}\right) = \frac{1}{2} - \frac{n\sin(2\pi/n)}{4\pi}

Expand sin(2π/n)=2πn(2π)36n3+\sin(2\pi/n) = \dfrac{2\pi}{n} - \dfrac{(2\pi)^3}{6n^3} + \cdots:

nsin(2π/n)4π=n4π(2πn8π36n3+)=12π23n2+\frac{n\sin(2\pi/n)}{4\pi} = \frac{n}{4\pi}\left(\frac{2\pi}{n} - \frac{8\pi^3}{6n^3} + \cdots\right) = \frac{1}{2} - \frac{\pi^2}{3n^2} + \cdots

So

01f(x)2dx=12(12π23n2+)=π23n2+O ⁣(1n4)\int_0^1 f(x)^2\,dx = \frac{1}{2} - \left(\frac{1}{2} - \frac{\pi^2}{3n^2} + \cdots\right) = \frac{\pi^2}{3n^2} + O\!\left(\frac{1}{n^4}\right)

And 13λ2=13(πn+O ⁣(1n3))2=π23n2+O ⁣(1n4)\frac{1}{3}\lambda^2 = \frac{1}{3}\left(\frac{\pi}{n} + O\!\left(\frac{1}{n^3}\right)\right)^2 = \frac{\pi^2}{3n^2} + O\!\left(\frac{1}{n^4}\right).

Therefore

R2=π23n2π23n2+O ⁣(1n4)=O ⁣(1n4)0R^2 = \frac{\pi^2}{3n^2} - \frac{\pi^2}{3n^2} + O\!\left(\frac{1}{n^4}\right) = O\!\left(\frac{1}{n^4}\right) \to 0

as nn \to \infty, so limnR=0\lim_{n\to\infty}R = 0.

Explanation. For large nn, sin(πx/n)πx/n\sin(\pi x/n) \approx \pi x/n over [0,1][0,1], since πx/n\pi x/n is small. But πx/n\pi x/n is already a linear function of xx, so the best linear approximation λx\lambda x fits f(x)f(x) almost perfectly and the residual error tends to zero.


Topic: 积分 Integration  |  Difficulty: Challenging  |  Marks: 20

6 Show that sinθ=2t1+t2,cosθ=1t21+t2,1+cosθsinθ=tan(π/2θ/2),\sin \theta = \frac{2t}{1 + t^2}, \quad \cos \theta = \frac{1 - t^2}{1 + t^2}, \quad \frac{1 + \cos \theta}{\sin \theta} = \tan(\pi/2 - \theta/2), where t=tan(θ/2)t = \tan(\theta/2).

Use the substitution t=tan(θ/2)t = \tan(\theta/2) to show that, for 0<α<π/20 < \alpha < \pi/2, 0π211+cosαsinθdθ=αsinα,\int_{0}^{\frac{\pi}{2}} \frac{1}{1 + \cos \alpha \sin \theta} \, d\theta = \frac{\alpha}{\sin \alpha}, and deduce a similar result for 0π211+sinαcosθdθ.\int_{0}^{\frac{\pi}{2}} \frac{1}{1 + \sin \alpha \cos \theta} \, d\theta.

Model Solution

Part (a): Proving the three identities with t=tan(θ/2)t = \tan(\theta/2).

For the first two identities, start from sinθ=2sin(θ/2)cos(θ/2)\sin\theta = 2\sin(\theta/2)\cos(\theta/2). Divide numerator and denominator by cos2(θ/2)\cos^2(\theta/2):

sinθ=2sin(θ/2)cos(θ/2)1=2sin(θ/2)cos(θ/2)cos2(θ/2)+sin2(θ/2)\sin\theta = \frac{2\sin(\theta/2)\cos(\theta/2)}{1} = \frac{2\sin(\theta/2)\cos(\theta/2)}{\cos^2(\theta/2) + \sin^2(\theta/2)}

Dividing top and bottom by cos2(θ/2)\cos^2(\theta/2):

=2tan(θ/2)1+tan2(θ/2)=2t1+t2.(i)= \frac{2\tan(\theta/2)}{1 + \tan^2(\theta/2)} = \frac{2t}{1 + t^2}. \qquad \text{(i)}

Similarly for cosine, start from cosθ=cos2(θ/2)sin2(θ/2)\cos\theta = \cos^2(\theta/2) - \sin^2(\theta/2):

cosθ=cos2(θ/2)sin2(θ/2)cos2(θ/2)+sin2(θ/2)\cos\theta = \frac{\cos^2(\theta/2) - \sin^2(\theta/2)}{\cos^2(\theta/2) + \sin^2(\theta/2)}

Dividing top and bottom by cos2(θ/2)\cos^2(\theta/2):

=1tan2(θ/2)1+tan2(θ/2)=1t21+t2.(ii)= \frac{1 - \tan^2(\theta/2)}{1 + \tan^2(\theta/2)} = \frac{1 - t^2}{1 + t^2}. \qquad \text{(ii)}

For the third identity, use cosθ=2cos2(θ/2)1\cos\theta = 2\cos^2(\theta/2) - 1, so 1+cosθ=2cos2(θ/2)1 + \cos\theta = 2\cos^2(\theta/2). Also sinθ=2sin(θ/2)cos(θ/2)\sin\theta = 2\sin(\theta/2)\cos(\theta/2).

1+cosθsinθ=2cos2(θ/2)2sin(θ/2)cos(θ/2)=cos(θ/2)sin(θ/2)=cot(θ/2)\frac{1 + \cos\theta}{\sin\theta} = \frac{2\cos^2(\theta/2)}{2\sin(\theta/2)\cos(\theta/2)} = \frac{\cos(\theta/2)}{\sin(\theta/2)} = \cot(\theta/2)

Since cot(θ/2)=tan(π/2θ/2)\cot(\theta/2) = \tan(\pi/2 - \theta/2), we have

1+cosθsinθ=tan(π/2θ/2).(iii)\frac{1 + \cos\theta}{\sin\theta} = \tan(\pi/2 - \theta/2). \qquad \text{(iii)}

Part (b): Evaluating 0π/2dθ1+cosαsinθ\displaystyle\int_0^{\pi/2}\frac{d\theta}{1 + \cos\alpha\,\sin\theta}.

Use the substitution t=tan(θ/2)t = \tan(\theta/2), so dθ=21+t2dtd\theta = \dfrac{2}{1+t^2}\,dt. The limits transform as: θ=0t=0\theta = 0 \Rightarrow t = 0; θ=π/2t=tan(π/4)=1\theta = \pi/2 \Rightarrow t = \tan(\pi/4) = 1.

Substituting sinθ=2t1+t2\sin\theta = \dfrac{2t}{1+t^2} from (i):

0π/2dθ1+cosαsinθ=0111+cosα2t1+t221+t2dt\int_0^{\pi/2}\frac{d\theta}{1 + \cos\alpha\,\sin\theta} = \int_0^1 \frac{1}{1 + \cos\alpha\cdot\dfrac{2t}{1+t^2}}\cdot\frac{2}{1+t^2}\,dt

=012(1+t2)+2tcosαdt=012t2+2tcosα+1dt= \int_0^1 \frac{2}{(1+t^2) + 2t\cos\alpha}\,dt = \int_0^1 \frac{2}{t^2 + 2t\cos\alpha + 1}\,dt

Complete the square in the denominator:

t2+2tcosα+1=(t+cosα)2+1cos2α=(t+cosα)2+sin2αt^2 + 2t\cos\alpha + 1 = (t + \cos\alpha)^2 + 1 - \cos^2\alpha = (t + \cos\alpha)^2 + \sin^2\alpha

So the integral becomes

012(t+cosα)2+sin2αdt\int_0^1 \frac{2}{(t+\cos\alpha)^2 + \sin^2\alpha}\,dt

Substitute u=t+cosαsinαu = \dfrac{t + \cos\alpha}{\sin\alpha}, so du=dtsinαdu = \dfrac{dt}{\sin\alpha}, i.e., dt=sinαdudt = \sin\alpha\,du.

Limits: t=0u=cotαt = 0 \Rightarrow u = \cot\alpha; t=1u=1+cosαsinα=cot(α/2)t = 1 \Rightarrow u = \dfrac{1 + \cos\alpha}{\sin\alpha} = \cot(\alpha/2) by identity (iii).

=cotαcot(α/2)2sinαsin2α(u2+1)du=2sinαcotαcot(α/2)duu2+1= \int_{\cot\alpha}^{\cot(\alpha/2)} \frac{2\sin\alpha}{\sin^2\alpha\,(u^2 + 1)}\,du = \frac{2}{\sin\alpha}\int_{\cot\alpha}^{\cot(\alpha/2)}\frac{du}{u^2 + 1}

=2sinα[arctanu]cotαcot(α/2)=2sinα(arctan(cot(α/2))arctan(cotα))= \frac{2}{\sin\alpha}\Big[\arctan u\Big]_{\cot\alpha}^{\cot(\alpha/2)} = \frac{2}{\sin\alpha}\left(\arctan(\cot(\alpha/2)) - \arctan(\cot\alpha)\right)

Now use the identity: for 0<x<π/20 < x < \pi/2, arctan(cotx)=arctan(tan(π/2x))=π/2x\arctan(\cot x) = \arctan(\tan(\pi/2 - x)) = \pi/2 - x.

=2sinα(π2α2π2+α)=2sinαα2=αsinα.= \frac{2}{\sin\alpha}\left(\frac{\pi}{2} - \frac{\alpha}{2} - \frac{\pi}{2} + \alpha\right) = \frac{2}{\sin\alpha}\cdot\frac{\alpha}{2} = \frac{\alpha}{\sin\alpha}. \qquad \blacksquare

Part (c): Deducing 0π/2dθ1+sinαcosθ\displaystyle\int_0^{\pi/2}\frac{d\theta}{1 + \sin\alpha\,\cos\theta}.

Substitute ϕ=π/2θ\phi = \pi/2 - \theta in the integral. Then cosθ=sinϕ\cos\theta = \sin\phi, dθ=dϕd\theta = -d\phi, and the limits reverse: θ=0ϕ=π/2\theta = 0 \to \phi = \pi/2, θ=π/2ϕ=0\theta = \pi/2 \to \phi = 0.

0π/2dθ1+sinαcosθ=0π/2dϕ1+sinαsinϕ\int_0^{\pi/2}\frac{d\theta}{1 + \sin\alpha\,\cos\theta} = \int_0^{\pi/2}\frac{d\phi}{1 + \sin\alpha\,\sin\phi}

This has exactly the form 0π/2dϕ1+cosβsinϕ\displaystyle\int_0^{\pi/2}\frac{d\phi}{1 + \cos\beta\,\sin\phi} with cosβ=sinα\cos\beta = \sin\alpha.

Since 0<α<π/20 < \alpha < \pi/2, we have cosβ=sinα=cos(π/2α)\cos\beta = \sin\alpha = \cos(\pi/2 - \alpha), so β=π/2α\beta = \pi/2 - \alpha.

By the result of Part (b):

0π/2dϕ1+sinαsinϕ=βsinβ=π/2αsin(π/2α)=π/2αcosα\int_0^{\pi/2}\frac{d\phi}{1 + \sin\alpha\,\sin\phi} = \frac{\beta}{\sin\beta} = \frac{\pi/2 - \alpha}{\sin(\pi/2 - \alpha)} = \frac{\pi/2 - \alpha}{\cos\alpha}

Therefore

0π/2dθ1+sinαcosθ=π/2αcosα.\int_0^{\pi/2}\frac{d\theta}{1 + \sin\alpha\,\cos\theta} = \frac{\pi/2 - \alpha}{\cos\alpha}. \qquad \blacksquare


Topic: 向量几何 Vector Geometry  |  Difficulty: Challenging  |  Marks: 20

7 The line ll has vector equation r=λs\mathbf{r} = \lambda\mathbf{s}, where s=(cosθ+3)i+(2sinθ)j+(cosθ3)k\mathbf{s} = (\cos \theta + \sqrt{3} \, ) \, \mathbf{i} + (\sqrt{2} \, \sin \theta) \, \mathbf{j} + (\cos \theta - \sqrt{3} \, ) \, \mathbf{k} and λ\lambda is a scalar parameter. Find an expression for the angle between ll and the line r=μ(ai+bj+ck)\mathbf{r} = \mu(a \, \mathbf{i} + b \, \mathbf{j} + c \, \mathbf{k}). Show that there is a line mm through the origin such that, whatever the value of θ\theta, the acute angle between ll and mm is π/6\pi/6.

A plane has equation xz=43x - z = 4\sqrt{3}. The line ll meets this plane at PP. Show that, as θ\theta varies, PP describes a circle, with its centre on mm. Find the radius of this circle.

Model Solution

Part 1: Angle between ll and a general line through the origin.

The direction vector of ll is s=(cosθ+3)i+2sinθj+(cosθ3)k\mathbf{s} = (\cos\theta + \sqrt{3})\,\mathbf{i} + \sqrt{2}\sin\theta\,\mathbf{j} + (\cos\theta - \sqrt{3})\,\mathbf{k}. Let the second line have direction t=ai+bj+ck\mathbf{t} = a\,\mathbf{i} + b\,\mathbf{j} + c\,\mathbf{k}.

First, compute s|\mathbf{s}|:

s2=(cosθ+3)2+2sin2θ+(cosθ3)2|\mathbf{s}|^2 = (\cos\theta + \sqrt{3})^2 + 2\sin^2\theta + (\cos\theta - \sqrt{3})^2

=cos2θ+23cosθ+3+2sin2θ+cos2θ23cosθ+3= \cos^2\theta + 2\sqrt{3}\cos\theta + 3 + 2\sin^2\theta + \cos^2\theta - 2\sqrt{3}\cos\theta + 3

=2cos2θ+2sin2θ+6=2+6=8= 2\cos^2\theta + 2\sin^2\theta + 6 = 2 + 6 = 8

So s=22|\mathbf{s}| = 2\sqrt{2}.

The acute angle ϕ\phi between the two lines is given by

cosϕ=stst=a(cosθ+3)+b2sinθ+c(cosθ3)22a2+b2+c2\cos\phi = \frac{|\mathbf{s} \cdot \mathbf{t}|}{|\mathbf{s}||\mathbf{t}|} = \frac{|a(\cos\theta + \sqrt{3}) + b\sqrt{2}\sin\theta + c(\cos\theta - \sqrt{3})|}{2\sqrt{2}\sqrt{a^2 + b^2 + c^2}}

=(a+c)cosθ+2bsinθ+3(ac)22a2+b2+c2(*)= \frac{|(a + c)\cos\theta + \sqrt{2}\,b\sin\theta + \sqrt{3}(a - c)|}{2\sqrt{2}\sqrt{a^2 + b^2 + c^2}} \qquad \text{(*)}

Part 2: Showing line mm exists with constant angle π/6\pi/6.

We need a fixed direction m=pi+qj+rk\mathbf{m} = p\,\mathbf{i} + q\,\mathbf{j} + r\,\mathbf{k} such that the angle between ll and mm is π/6\pi/6 for all θ\theta. From (*), setting a=p,b=q,c=ra = p, b = q, c = r:

sm=(p+r)cosθ+2qsinθ+3(pr)\mathbf{s} \cdot \mathbf{m} = (p + r)\cos\theta + \sqrt{2}\,q\sin\theta + \sqrt{3}(p - r)

For this dot product to be independent of θ\theta, the coefficients of cosθ\cos\theta and sinθ\sin\theta must vanish:

p+r=0andq=0p + r = 0 \quad \text{and} \quad q = 0

So r=pr = -p and q=0q = 0, giving m=pipk=p(ik)\mathbf{m} = p\,\mathbf{i} - p\,\mathbf{k} = p(\mathbf{i} - \mathbf{k}).

With this choice, sm=3(p(p))=23p\mathbf{s} \cdot \mathbf{m} = \sqrt{3}(p - (-p)) = 2\sqrt{3}\,p, and m=p2|\mathbf{m}| = |p|\sqrt{2}. The angle satisfies:

cosα=23p22p2=234=32\cos\alpha = \frac{|2\sqrt{3}\,p|}{2\sqrt{2} \cdot |p|\sqrt{2}} = \frac{2\sqrt{3}}{4} = \frac{\sqrt{3}}{2}

Since cosα=3/2\cos\alpha = \sqrt{3}/2, we get α=π/6\alpha = \pi/6 for every value of θ\theta. The line mm through the origin has direction ik\mathbf{i} - \mathbf{k}.

Part 3: Locus of PP and the circle.

Points on ll have coordinates x=λ(cosθ+3)x = \lambda(\cos\theta + \sqrt{3}), y=λ2sinθy = \lambda\sqrt{2}\sin\theta, z=λ(cosθ3)z = \lambda(\cos\theta - \sqrt{3}). Substituting into the plane equation xz=43x - z = 4\sqrt{3}:

λ(cosθ+3)λ(cosθ3)=43\lambda(\cos\theta + \sqrt{3}) - \lambda(\cos\theta - \sqrt{3}) = 4\sqrt{3}

23λ=43    λ=22\sqrt{3}\,\lambda = 4\sqrt{3} \implies \lambda = 2

So PP has coordinates:

P=(2cosθ+23,  22sinθ,  2cosθ23)P = \big(2\cos\theta + 2\sqrt{3},\; 2\sqrt{2}\sin\theta,\; 2\cos\theta - 2\sqrt{3}\big)

The centre of the circle lies on both mm and the plane. A general point on mm is (t,0,t)(t, 0, -t); substituting into xz=43x - z = 4\sqrt{3} gives t(t)=43t - (-t) = 4\sqrt{3}, so t=23t = 2\sqrt{3}. The centre is C=(23,0,23)C = (2\sqrt{3}, 0, -2\sqrt{3}).

The radius is CP|CP|:

CP2=(2cosθ+2323)2+(22sinθ0)2+(2cosθ23(23))2|CP|^2 = (2\cos\theta + 2\sqrt{3} - 2\sqrt{3})^2 + (2\sqrt{2}\sin\theta - 0)^2 + (2\cos\theta - 2\sqrt{3} - (-2\sqrt{3}))^2

=(2cosθ)2+(22sinθ)2+(2cosθ)2= (2\cos\theta)^2 + (2\sqrt{2}\sin\theta)^2 + (2\cos\theta)^2

=4cos2θ+8sin2θ+4cos2θ=8cos2θ+8sin2θ=8= 4\cos^2\theta + 8\sin^2\theta + 4\cos^2\theta = 8\cos^2\theta + 8\sin^2\theta = 8

So CP=22|CP| = 2\sqrt{2}, which is independent of θ\theta. As θ\theta varies, PP traces out a circle of radius 222\sqrt{2} with centre on mm.


Topic: 微分方程 Differential Equations  |  Difficulty: Challenging  |  Marks: 20

8 (i) Let yy be the solution of the differential equation dydx+4xex2(y+3)12=0(x0),\frac{dy}{dx} + 4x \, e^{-x^2}(y + 3)^{\frac{1}{2}} = 0 \quad (x \geqslant 0), that satisfies the condition y=6y = 6 when x=0x = 0. Find yy in terms of xx and show that y1y \to 1 as xx \to \infty.

(ii) Let yy be any solution of the differential equation dydxxe6x2(y+3)1k=0(x0).\frac{dy}{dx} - x \, e^{6x^2}(y + 3)^{1-k} = 0 \quad (x \geqslant 0). Find a value of kk such that, as xx \to \infty, e3x2ye^{-3x^2}y tends to a finite non-zero limit, which you should determine.

Model Solution

Part (i): Solving the ODE and finding the limit.

The equation is:

dydx=4xex2(y+3)1/2\frac{dy}{dx} = -4x\,e^{-x^2}(y + 3)^{1/2}

This is separable. Divide both sides by (y+3)1/2(y + 3)^{1/2}:

(y+3)1/2dydx=4xex2(y + 3)^{-1/2}\,\frac{dy}{dx} = -4x\,e^{-x^2}

Integrate both sides with respect to xx:

(y+3)1/2dy=4xex2dx\int (y + 3)^{-1/2}\,dy = \int -4x\,e^{-x^2}\,dx

The left side: let u=y+3u = y + 3, so du=dydu = dy:

u1/2du=2u1/2=2(y+3)1/2\int u^{-1/2}\,du = 2u^{1/2} = 2(y + 3)^{1/2}

The right side: let v=x2v = -x^2, so dv=2xdxdv = -2x\,dx, hence 4xdx=2dv-4x\,dx = 2\,dv:

2evdv=2ev=2ex2\int 2e^v\,dv = 2e^v = 2e^{-x^2}

Combining:

2(y+3)1/2=2ex2+C2(y + 3)^{1/2} = 2e^{-x^2} + C

Apply the initial condition y=6y = 6 when x=0x = 0:

2(6+3)1/2=2e0+C    23=2+C    C=42(6 + 3)^{1/2} = 2e^0 + C \implies 2 \cdot 3 = 2 + C \implies C = 4

So:

2(y+3)1/2=2ex2+42(y + 3)^{1/2} = 2e^{-x^2} + 4

(y+3)1/2=ex2+2(y + 3)^{1/2} = e^{-x^2} + 2

Squaring both sides:

y+3=(ex2+2)2=e2x2+4ex2+4y + 3 = (e^{-x^2} + 2)^2 = e^{-2x^2} + 4e^{-x^2} + 4

y=e2x2+4ex2+1y = e^{-2x^2} + 4e^{-x^2} + 1

As xx \to \infty, both e2x20e^{-2x^2} \to 0 and ex20e^{-x^2} \to 0, so y0+0+1=1y \to 0 + 0 + 1 = 1.

Part (ii): Finding kk for a finite non-zero limit of e3x2ye^{-3x^2}y.

The equation is:

dydx=xe6x2(y+3)1k\frac{dy}{dx} = x\,e^{6x^2}(y + 3)^{1-k}

Separate variables:

(y+3)k1dy=xe6x2dx(y + 3)^{k-1}\,dy = x\,e^{6x^2}\,dx

Integrate both sides. The right side: let w=6x2w = 6x^2, so dw=12xdxdw = 12x\,dx:

xe6x2dx=112ewdw=112e6x2\int x\,e^{6x^2}\,dx = \frac{1}{12}\int e^w\,dw = \frac{1}{12}e^{6x^2}

Case k0k \neq 0: The left side integrates to (y+3)kk\frac{(y+3)^k}{k}, giving:

(y+3)kk=112e6x2+C1\frac{(y + 3)^k}{k} = \frac{1}{12}e^{6x^2} + C_1

    (y+3)k=k12e6x2+Cwhere C=kC1\implies (y + 3)^k = \frac{k}{12}e^{6x^2} + C \qquad \text{where } C = kC_1

For large xx, if C0C \geq 0:

(y+3)kk12e6x2(y + 3)^k \approx \frac{k}{12}e^{6x^2}

y+3(k12)1/ke6x2/ky + 3 \approx \left(\frac{k}{12}\right)^{1/k} e^{6x^2/k}

y(k12)1/ke6x2/ky \approx \left(\frac{k}{12}\right)^{1/k} e^{6x^2/k}

Now examine e3x2ye^{-3x^2}y:

e3x2y(k12)1/ke(6/k3)x2e^{-3x^2}y \approx \left(\frac{k}{12}\right)^{1/k} e^{(6/k - 3)x^2}

For this to tend to a finite non-zero limit as xx \to \infty, the exponent of ee must be exactly zero:

6k3=0    k=2\frac{6}{k} - 3 = 0 \implies k = 2

Verification and determination of the limit. With k=2k = 2:

(y+3)2=16e6x2+C(y + 3)^2 = \frac{1}{6}e^{6x^2} + C

y+3=16e6x2+C(taking positive root since y+3>0)y + 3 = \sqrt{\frac{1}{6}e^{6x^2} + C} \qquad \text{(taking positive root since } y + 3 > 0\text{)}

As xx \to \infty:

y+3=e3x261+6Ce6x2e3x26y + 3 = \frac{e^{3x^2}}{\sqrt{6}}\sqrt{1 + 6Ce^{-6x^2}} \to \frac{e^{3x^2}}{\sqrt{6}}

Therefore:

e3x2(y+3)16    e3x2y=e3x2(y+3)3e3x2160=16e^{-3x^2}(y + 3) \to \frac{1}{\sqrt{6}} \implies e^{-3x^2}y = e^{-3x^2}(y+3) - 3e^{-3x^2} \to \frac{1}{\sqrt{6}} - 0 = \frac{1}{\sqrt{6}}

The value of kk is 22, and the limit is 16\dfrac{1}{\sqrt{6}}.