Skip to content

STEP2 2000 -- Pure Mathematics

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

All questions are pure mathematics. Model solutions 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

(a−N)(b−N)=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+1≠N(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 a≠ba \neq b, multiply through by NabNab:

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

ab−Na−Nb=0ab - Na - Nb = 0

Add N2N^2 to both sides:

ab−Na−Nb+N2=N2ab - Na - Nb + N^2 = N^2

(a−N)(b−N)=N2(a - N)(b - N) = N^2

Since 1a<1N\frac{1}{a} < \frac{1}{N} (because 1a=1N−1b\frac{1}{a} = \frac{1}{N} - \frac{1}{b} and b>0b > 0), we have a>Na > N, so a−N>0a - N > 0. Similarly b>Nb > N, so b−N>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 d1⩽d2d_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 a≠ba \neq b. Multiply through by NabNab:

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

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

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

Add N2N^2 to both sides:

4ab−2Na−2Nb+N2=N24ab - 2Na - 2Nb + N^2 = N^2

(2a−N)(2b−N)=N2(2a - N)(2b - N) = N^2

Since 1a<2N\frac{1}{a} < \frac{2}{N}, we have a>N/2a > N/2, so 2a−N>02a - N > 0. Similarly 2b−N>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): 2a−N=N2a - N = N and 2b−N=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 a≠ba \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 (x−a)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 (x−a)4(x - a)^4 is a factor of p(x)p(x).

Given that the polynomial

x6+4x5−5x4−40x3−40x2+32x+kx^6 + 4x^5 - 5x^4 - 40x^3 - 40x^2 + 32x + k

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

Model Solution

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

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

Evaluating at x=ax = a:

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


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

p′(x)=4(x−a)3 q(x)+(x−a)4 q′(x)=(x−a)3[4q(x)+(x−a)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)=(x−a)3s(x)p'(x) = (x - a)^3 s(x), which means (x−a)2(x - a)^2 is a factor of p′(x)p'(x). By the result just proved (applied to p′p' instead of pp):

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

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

p′′(x)=3(x−a)2 s(x)+(x−a)3 s′(x)=(x−a)2[3s(x)+(x−a)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)=(x−a)2t(x)p''(x) = (x - a)^2 t(x), meaning (x−a)2(x - a)^2 is a factor of p′′(x)p''(x). Applying the same result to p′′p'':

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

Conclusion: if (x−a)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+4x5−5x4−40x3−40x2+32x+kp(x) = x^6 + 4x^5 - 5x^4 - 40x^3 - 40x^2 + 32x + k. Compute the first three derivatives:

p′(x)=6x5+20x4−20x3−120x2−80x+32p'(x) = 6x^5 + 20x^4 - 20x^3 - 120x^2 - 80x + 32

p′′(x)=30x4+80x3−60x2−240x−80p''(x) = 30x^4 + 80x^3 - 60x^2 - 240x - 80

p′′′(x)=120x3+240x2−120x−240p'''(x) = 120x^3 + 240x^2 - 120x - 240

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

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

Factor by grouping:

a2(a+2)−1(a+2)=0  ⟹  (a+2)(a2−1)=0  ⟹  (a+2)(a−1)(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=480−640−240+480−80=0✓p''(-2) = 30(16) + 80(-8) - 60(4) - 240(-2) - 80 = 480 - 640 - 240 + 480 - 80 = 0 \checkmark

p′′(−1)=30−80−60+240−80=50≠0p''(-1) = 30 - 80 - 60 + 240 - 80 = 50 \neq 0

p′′(1)=30+80−60−240−80=−270≠0p''(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+160−480+160+32=0✓= -192 + 320 + 160 - 480 + 160 + 32 = 0 \checkmark

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

p(−2)=64−128−80+320−160−64+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(x2−4x+3)=x6+4x5−5x4−40x3−40x2+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 x2−4x+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^6−4x5-4x^53x43x^4
8x38x^38x58x^5−32x4-32x^424x324x^3
24x224x^224x424x^4−96x3-96x^372x272x^2
32x32x32x332x^3−128x2-128x^296x96x
161616x216x^2−64x-64x4848

Collecting:

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

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

p(x)=(x+2)4(x−1)(x−3)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 a≈7+ηa \approx 7 + \eta, where η=(13 ϵ1−2 ϵ2+243 ϵ3)/14\eta = (13 \, \epsilon_1 - 2 \, \epsilon_2 + 24\sqrt{3} \, \epsilon_3)/14.

Given now that

∣ϵ1∣⩽2×10−3,∣ϵ2∣⩽4.9×10−2,∣ϵ3∣⩽3×10−3,|\epsilon_1| \leqslant 2 \times 10^{-3}, \quad |\epsilon_2| \leqslant 4.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 a≈7+ηa \approx 7 + \eta

By the cosine rule,

a2=b2+c2−2bccos⁡A.(*)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+9−48cos⁡(π/3)=73−24=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:

2a da=2b db+2c dc−2cos⁡A (c db+b dc)−2bc(−sin⁡A) dA.2a\,da = 2b\,db + 2c\,dc - 2\cos A\,(c\,db + b\,dc) - 2bc(-\sin A)\,dA.

Note the last term: ∂∂A(−2bccos⁡A)=2bcsin⁡A\frac{\partial}{\partial A}(-2bc\cos A) = 2bc\sin A. Dividing by 2:

a da=b db+c dc−cos⁡A (c db+b dc)+bcsin⁡A dA.a\,da = b\,db + c\,dc - \cos A\,(c\,db + b\,dc) + bc\sin A\,dA.

Collecting terms in dbdb, dcdc, dAdA:

a da=(b−ccos⁡A) db+(c−bcos⁡A) dc+bcsin⁡A dA.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:

  • b−ccos⁡A=8−3⋅12=132b - c\cos A = 8 - 3 \cdot \tfrac{1}{2} = \tfrac{13}{2}
  • c−bcos⁡A=3−8⋅12=−1c - b\cos A = 3 - 8 \cdot \tfrac{1}{2} = -1
  • bcsin⁡A=24⋅32=123bc\sin A = 24 \cdot \tfrac{\sqrt{3}}{2} = 12\sqrt{3}

Therefore

7 dη=132 dϵ1−dϵ2+123 dϵ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ϵ1−2ϵ2+243 ϵ314.\eta = \frac{13\epsilon_1 - 2\epsilon_2 + 24\sqrt{3}\,\epsilon_3}{14}.

Hence a≈7+ηa \approx 7 + \eta where η=(13ϵ1−2ϵ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ϵ1−2ϵ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 ∣ϵ3∣14.|\eta| \leqslant \frac{13|\epsilon_1| + 2|\epsilon_2| + 24\sqrt{3}\,|\epsilon_3|}{14}.

Substituting the given bounds:

∣η∣⩽13×2×10−3+2×4.9×10−2+243×3×10−314.|\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×10−3=26×10−313 \times 2 \times 10^{-3} = 26 \times 10^{-3}
  • 2×4.9×10−2=98×10−32 \times 4.9 \times 10^{-2} = 98 \times 10^{-3}
  • 243×3×10−3=24×3×10−3=72×10−324\sqrt{3} \times \sqrt{3} \times 10^{-3} = 24 \times 3 \times 10^{-3} = 72 \times 10^{-3}

Therefore

∣η∣⩽(26+98+72)×10−314=196×10−314=14×10−3=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×10−3\epsilon_1 = 2 \times 10^{-3}, ϵ2=−4.9×10−2\epsilon_2 = -4.9 \times 10^{-2}, ϵ3=3×10−3\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=cos⁡nθ+isin⁡nθ.(\cos \theta + i \sin \theta)^n = \cos n\theta + i \sin n\theta. By considering (5−i)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 tan⁡−1θ\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=cos⁡nθ+isin⁡nθ(\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=cos⁡kθ+isin⁡kθ(\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) =(cos⁡kθ+isin⁡kθ)(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 (5−i)2(1+i)(5 - i)^2(1 + i).

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

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

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

(24−10i)(1+i)=24+24i−10i−10i2=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 5−i5 - i has positive real part and negative imaginary part, so arg⁡(5−i)=−arctan⁡(1/5)\arg(5 - i) = -\arctan(1/5). By de Moivre’s theorem:

arg⁡((5−i)2(1+i))=2arg⁡(5−i)+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×1985 i+992i2= 993 \times 1985 + 993i + 992 \times 1985\,i + 992i^2 =(993×1985−992)+(993+992×1985)i.= (993 \times 1985 - 992) + (993 + 992 \times 1985)i.

Evaluating:

  • Real part: 993×1985−992=1971105−992=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 0⩽x⩽10 \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 λ=3∫01xf(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))2dx−13λ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) lim⁡n→∞R=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)2 dxI(\lambda) = \displaystyle\int_{0}^{1}(f(x) - \lambda x)^2\,dx. Expand the integrand:

I(λ)=∫01f(x)2 dx−2λ∫01xf(x) dx+λ2∫01x2 dxI(\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:

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

So

I(λ)=∫01f(x)2 dx−2λ∫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λ=−2∫01xf(x) dx+2λ3=0\frac{dI}{d\lambda} = -2\int_{0}^{1}xf(x)\,dx + \frac{2\lambda}{3} = 0

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

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

From Step 1, we have λ=3∫01xf(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)2 dx−2λ⋅λ3+λ23R^2 = I(\lambda) = \int_{0}^{1}f(x)^2\,dx - 2\lambda\cdot\frac{\lambda}{3} + \frac{\lambda^2}{3}

=∫01f(x)2 dx−2λ23+λ23=∫01f(x)2 dx−λ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.

λ=3∫01xsin⁡ ⁣(π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:

λ≈3∫01x(πxn−π3x36n3)dx=3∫01(π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 lim⁡n→∞R=0\lim_{n\to\infty}R = 0.

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

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

∫01sin⁡2 ⁣(πxn)dx=n2π∫0π/n(1−cos⁡2u) du=n2π[u−sin⁡2u2]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π(πn−sin⁡(2π/n)2)=12−nsin⁡(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πn−8π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)2 dx=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 n→∞n \to \infty, so lim⁡n→∞R=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⁡θ=1−t21+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 cos⁡2(θ/2)\cos^2(\theta/2):

sin⁡θ=2sin⁡(θ/2)cos⁡(θ/2)1=2sin⁡(θ/2)cos⁡(θ/2)cos⁡2(θ/2)+sin⁡2(θ/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 cos⁡2(θ/2)\cos^2(\theta/2):

=2tan⁡(θ/2)1+tan⁡2(θ/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⁡θ=cos⁡2(θ/2)−sin⁡2(θ/2)\cos\theta = \cos^2(\theta/2) - \sin^2(\theta/2):

cos⁡θ=cos⁡2(θ/2)−sin⁡2(θ/2)cos⁡2(θ/2)+sin⁡2(θ/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 cos⁡2(θ/2)\cos^2(\theta/2):

=1−tan⁡2(θ/2)1+tan⁡2(θ/2)=1−t21+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⁡θ=2cos⁡2(θ/2)−1\cos\theta = 2\cos^2(\theta/2) - 1, so 1+cos⁡θ=2cos⁡2(θ/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⁡θ=2cos⁡2(θ/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+t2 dtd\theta = \dfrac{2}{1+t^2}\,dt. The limits transform as: θ=0⇒t=0\theta = 0 \Rightarrow t = 0; θ=π/2⇒t=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+t2⋅21+t2 dt\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⁡α+1 dt= \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+1−cos⁡2α=(t+cos⁡α)2+sin⁡2α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+sin⁡2α 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=0⇒u=cot⁡αt = 0 \Rightarrow u = \cot\alpha; t=1⇒u=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⁡αsin⁡2α (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⁡α[arctan⁡u]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⁡(cot⁡x)=arctan⁡(tan⁡(π/2−x))=π/2−x\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+(2 sin⁡θ) 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=μ(a i+b j+c k)\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 x−z=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=a i+b j+c k\mathbf{t} = a\,\mathbf{i} + b\,\mathbf{j} + c\,\mathbf{k}.

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

∣s∣2=(cos⁡θ+3)2+2sin⁡2θ+(cos⁡θ−3)2|\mathbf{s}|^2 = (\cos\theta + \sqrt{3})^2 + 2\sin^2\theta + (\cos\theta - \sqrt{3})^2

=cos⁡2θ+23cos⁡θ+3+2sin⁡2θ+cos⁡2θ−23cos⁡θ+3= \cos^2\theta + 2\sqrt{3}\cos\theta + 3 + 2\sin^2\theta + \cos^2\theta - 2\sqrt{3}\cos\theta + 3

=2cos⁡2θ+2sin⁡2θ+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⁡ϕ=∣s⋅t∣∣s∣∣t∣=∣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⁡θ+2 bsin⁡θ+3(a−c)∣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=p i+q j+r k\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:

s⋅m=(p+r)cos⁡θ+2 qsin⁡θ+3(p−r)\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=p i−p k=p(i−k)\mathbf{m} = p\,\mathbf{i} - p\,\mathbf{k} = p(\mathbf{i} - \mathbf{k}).

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

cos⁡α=∣23 p∣22⋅∣p∣2=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 i−k\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 x−z=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 x−z=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|:

∣CP∣2=(2cos⁡θ+23−23)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

=4cos⁡2θ+8sin⁡2θ+4cos⁡2θ=8cos⁡2θ+8sin⁡2θ=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+4x e−x2(y+3)12=0(x⩾0),\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 y→1y \to 1 as x→∞x \to \infty.

(ii) Let yy be any solution of the differential equation dydx−x e6x2(y+3)1−k=0(x⩾0).\frac{dy}{dx} - x \, e^{6x^2}(y + 3)^{1-k} = 0 \quad (x \geqslant 0). Find a value of kk such that, as x→∞x \to \infty, e−3x2ye^{-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=−4x e−x2(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/2 dydx=−4x e−x2(y + 3)^{-1/2}\,\frac{dy}{dx} = -4x\,e^{-x^2}

Integrate both sides with respect to xx:

∫(y+3)−1/2 dy=∫−4x e−x2 dx\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:

∫u−1/2 du=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=−2x dxdv = -2x\,dx, hence −4x dx=2 dv-4x\,dx = 2\,dv:

∫2ev dv=2ev=2e−x2\int 2e^v\,dv = 2e^v = 2e^{-x^2}

Combining:

2(y+3)1/2=2e−x2+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  ⟹  2⋅3=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=2e−x2+42(y + 3)^{1/2} = 2e^{-x^2} + 4

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

Squaring both sides:

y+3=(e−x2+2)2=e−2x2+4e−x2+4y + 3 = (e^{-x^2} + 2)^2 = e^{-2x^2} + 4e^{-x^2} + 4

y=e−2x2+4e−x2+1y = e^{-2x^2} + 4e^{-x^2} + 1

As x→∞x \to \infty, both e−2x2→0e^{-2x^2} \to 0 and e−x2→0e^{-x^2} \to 0, so y→0+0+1=1y \to 0 + 0 + 1 = 1.

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

The equation is:

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

Separate variables:

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

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

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

Case k≠0k \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 C≥0C \geq 0:

(y+3)k≈k12e6x2(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 e−3x2ye^{-3x^2}y:

e−3x2y≈(k12)1/ke(6/k−3)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 x→∞x \to \infty, the exponent of ee must be exactly zero:

6k−3=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 x→∞x \to \infty:

y+3=e3x261+6Ce−6x2→e3x26y + 3 = \frac{e^{3x^2}}{\sqrt{6}}\sqrt{1 + 6Ce^{-6x^2}} \to \frac{e^{3x^2}}{\sqrt{6}}

Therefore:

e−3x2(y+3)→16  ⟹  e−3x2y=e−3x2(y+3)−3e−3x2→16−0=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}}.