1 The coordinates of a particle at time t are x and y. For t⩾0, they satisfy the pair of coupled differential equations
x˙=−x−kyy˙=x−y
where k is a constant. When t=0, x=1 and y=0.
(i) Let k=1. Find x and y in terms of t and sketch y as a function of t.
Sketch the path of the particle in the x-y plane, giving the coordinates of the point at which y is greatest and the coordinates of the point at which x is least.
(ii) Instead, let k=0. Find x and y in terms of t and sketch the path of the particle in the x-y plane.
Hint
(i) x˙=−x−y
So y=−x˙−x
As y˙=x−y, −x¨−x˙=x+x˙+x
x¨+2x˙+2x=0
[M1]
AQE λ2+2λ+2=0, λ=−1±i so x=e−t(Acost+Bsint) [M1 A1 cao]
[Alternatively, x=y˙+y leading to y¨+2y˙+2y=0 and y=e−t(Ccost+Dsint) etc]
So y=e−t(Acost+Bsint)−e−t(−Asint+Bcost)−e−t(Acost+Bsint)
Applying initial conditions: x(0)=1 gives A=1, and y(0)=0 gives −B=0, so B=0.
x=e−tcost,y=e−tsint
Sketch of y as a function of t:
y(0)=0, and y=0 whenever sint=0, i.e. at t=0,π,2π,… The envelope ±e−t decays exponentially. y reaches its first positive maximum at t=4π, its first negative minimum at t=45π, and the oscillations die away as t→∞.
Path in the x-y plane:
Note that x2+y2=e−2t(cos2t+sin2t)=e−2t, so the distance from the origin is r=e−t, which decreases monotonically. The particle spirals inward toward the origin.
Maximum y: We need y˙=0, i.e. x−y=0, so x=y:
e−tcost=e−tsint⟹tant=1⟹t=4π+nπ
For the first (and largest) positive y, t=4π:
(2e−π/4,2e−π/4)
Minimum x: We need x˙=0, i.e. −x−y=0, so x=−y:
e−tcost=−e−tsint⟹tant=−1⟹t=43π+nπ
For the first such value, t=43π:
(−2e−3π/4,2e−3π/4)
Part (ii) (k=0)
With k=0, x˙=−x and y˙=x−y.
From the first equation: x˙=−x gives x=Ae−t. With x(0)=1: A=1, so x=e−t.
Substituting into the second equation: y˙+y=e−t.
The integrating factor is et:
dtd(ety)=et⋅e−t=1
ety=t+c
With y(0)=0: c=0, so:
x=e−t,y=te−t
Path in the x-y plane:
The path starts at (1,0) when t=0. As t increases, x=e−t decreases monotonically while y=te−t increases from 0 to a maximum and then decreases back to 0.
To find the maximum y on the curve, we compute dxdy=x˙y˙=−e−t(1−t)e−t=t−1.
Setting dxdy=0: t=1, giving the point (e−1,e−1).
As t→∞, dxdy=t−1→∞, meaning the tangent becomes vertical (approaching the origin along the positive x-axis from above).
The path starts at (1,0), rises to (e−1,e−1), then curves back toward the origin.
As h→0: g(h)→g(0)=0, so 1+g(x)g(h)→1, and hg(h)→g′(0)=k.
g′(x)=k(1−g(x)2)
Finding g(x):
1−g(x)2g′(x)=k
Integrating both sides. Since ∣g(x)∣<1, we use the standard result ∫1−u2du=tanh−1u:
tanh−1(g(x))=kx+c
With g(0)=0: tanh−1(0)=0, so c=0.
g(x)=tanh(kx)
We can verify: tanh satisfies the addition formula tanh(x+y)=1+tanhxtanhytanhx+tanhy, ∣tanh(kx)∣<1 for all x, and dxdtanh(kx)x=0=ksech2(0)=k. All conditions are satisfied.
(i) You are given that the transformation represented by A has a line L1 of invariant points (so that each point on L1 is transformed to itself). Let (x,y) be a point on L1. Show that ((a−1)(d−1)−bc)xy=0.
Show further that (a−1)(d−1)=bc.
What can be said about A if L1 does not pass through the origin?
(ii) By considering the cases b=0 and b=0 separately, show that if (a−1)(d−1)=bc then the transformation represented by A has a line of invariant points. You should identify the line in the different cases that arise.
(iii) You are given instead that the transformation represented by A has an invariant line L2 (so that each point on L2 is transformed to a point on L2) and that L2 does not pass through the origin. If L2 has the form y=mx+k, show that (a−1)(d−1)=bc.
Hint
(i) (acbd)(xy)=(xy)⇒ax+by=xcx+dy=yM1
(a−1)x=−by
−cx=(d−1)y
Thus (a−1)x(d−1)y=bycx, that is (a−1)(d−1)xy−bcxy=0
So ((a−1)(d−1)−bc)xy=0M1
Thus ((a−1)(d−1)−bc)=0 or x=0 or y=0M1
If L1 is x=0, then both by=0 and dy=y,∀y
Thus, b=0 and d=1, meaning that ((a−1)(d−1)−bc)=0E1
Similarly, if L1 is y=0, then both cx=0 and ax=x,∀x
Then c=0 and a=1, meaning that ((a−1)(d−1)−bc)=0
In all three cases, (a−1)(d−1)=bcE1
Alternatively,
(acbd)(xy)=(xy)⇒(acbd)(xy)=(1001)(xy)M1
So (acbd)(xy)−(1001)(xy)=(00)
That is (a−1cbd−1)(xy)=(00)M1
As this is true for a line of invariant points, it does not have a unique solution and so E1
det(a−1cbd−1)=0M1
and so ((a−1)(d−1)−bc)=0 which implies both ((a−1)(d−1)−bc)xy=0 and (a−1)(d−1)=bcE1
If L1 does not pass through the origin then either L1 is a) y=mx+k with k=0
or b) x=k with k=0E1
For a) (acbd)(xmx+k)=(xmx+k)∀x
Thus ax+b(mx+k)=x and cx+d(mx+k)=mx+k
As these apply for all x, and k=0, bk=0 which implies b=0 and a+bm=1 and thus a=1
Also dk=k implying d=1 and c+dm=m which thus gives c=0M1
Let (x,y) be a point on L1. Since it is an invariant point, A maps (x,y) to itself:
(acbd)(xy)=(xy)
This gives the two equations:
ax+by=x⇒(a−1)x=−by(1)
cx+dy=y⇒cx=−(d−1)y(2)
From (1): (a−1)x=−by. From (2): (d−1)y=−cx.
Multiplying these two equations:
(a−1)(d−1)xy=(−by)(−cx)=bcxy
((a−1)(d−1)−bc)xy=0
((a−1)(d−1)−bc)xy=0
as required. ■
Showing (a−1)(d−1)=bc:
The equation ((a−1)(d−1)−bc)xy=0 holds for every point on L1. Since L1 is a line (not just the origin), it contains infinitely many points. We consider three cases:
Case 1: If there exists a point on L1 with x=0 and y=0, then we can divide by xy to get (a−1)(d−1)−bc=0, i.e.\ (a−1)(d−1)=bc.
Case 2:L1 is the line x=0 (the y-axis). Then every point (0,y) is invariant. From (1): by=0 for all y, giving b=0. From (2): dy=y for all y, giving d=1. Hence (a−1)(d−1)=(a−1)(0)=0 and bc=0, so (a−1)(d−1)=bc.
Case 3:L1 is the line y=0 (the x-axis). Then every point (x,0) is invariant. From (1): (a−1)x=0 for all x, giving a=1. From (2): cx=0 for all x, giving c=0. Hence (a−1)(d−1)=0 and bc=0, so again (a−1)(d−1)=bc.
In all cases, (a−1)(d−1)=bc. ■
If L1 does not pass through the origin: Suppose L1 is the line y=mx+k with k=0. Every point (x,mx+k) satisfies A(mx+kx)=(mx+kx), so:
(a−1)x+b(mx+k)=0for all x
The constant term gives bk=0. Since k=0, we get b=0. Then the x-coefficient gives a−1+bm=0, so a=1.
(d−1)(mx+k)+cx=0for all x
The constant term gives (d−1)k=0, so d=1. The x-coefficient gives c+(d−1)m=0, so c=0.
Therefore A=I. The only matrix with an invariant line not through the origin is the identity matrix.
Part (ii)
Case b=0: We are given (a−1)(d−1)=bc. Consider the line (a−1)x+by=0, i.e.\ y=−b(a−1)x.
For any point (x,y) on this line:
First component: ax+by=ax+b⋅(−ba−1x)=ax−(a−1)x=x. ✓
Second component: cx+dy=cx−bd(a−1)x. We need this to equal y=−ba−1x, so:
cx−bd(a−1)x=−ba−1x
Dividing by x (for x=0, and the case x=0 gives y=0 which is trivially invariant):
c−bd(a−1)=−ba−1
cb−d(a−1)=−(a−1)
bc=(a−1)(d−1)
This is precisely our hypothesis. So the line of invariant points is (a−1)x+by=0, equivalently cx+(d−1)y=0 (these represent the same line since (a−1)(d−1)=bc).
Case b=0: The condition (a−1)(d−1)=bc=0 means a=1 or d=1.
Subcase a=1:A=(1c0d). The invariant point equations become 0=0 and cx+(d−1)y=0. The line of invariant points is cx+(d−1)y=0.
Subcase d=1, a=1:A=(ac01). The invariant point equations become (a−1)x=0 and cx=0. Since a=1, we need x=0 (and cx=0 is then automatic). The line of invariant points is x=0.
In every case, A has a line of invariant points. ■
Part (iii)
L2:y=mx+k with k=0 is an invariant line, meaning A maps every point on L2 to a (possibly different) point on L2. For any x, the image of (x,mx+k) must be of the form (x′,mx′+k):
ax+b(mx+k)=x′(3)
cx+d(mx+k)=mx′+k(4)
Case m=0: From (3): x′=(a+bm)x+bk. Substitute into (4):
cx+dmx+dk=m(a+bm)x+mbk+k
(c+dm)x+dk=(ma+m2b)x+mbk+k
Equating constant terms (x=0):
dk=mbk+k
Since k=0, divide by k:
d=mb+1⇒d−1=mb(5)
Equating coefficients of x:
c+dm=ma+m2b
Using (5), dm=mb2+m:
c+mb2+m=ma+m2b
c=ma−m=m(a−1)(6)
Multiplying (5) and (6):
(d−1)⋅m(a−1)=mb⋅c
m(a−1)(d−1)=mbc
Since m=0, divide by m:
(a−1)(d−1)=bc■
Case m=0: The line is y=k with k=0. For any (x,k), the image is (x′,k):
ax+bk=x′(7)
cx+dk=k(8)
From (8), since this must hold for all x: the coefficient of x gives c=0, and the constant term gives dk=k, so d=1 (since k=0).
Therefore (a−1)(d−1)=(a−1)⋅0=0 and bc=b⋅0=0, giving (a−1)(d−1)=bc. ■
4 The nth degree polynomial P(x) is said to be reflexive if:
(a) P(x) is of the form xn−a1xn−1+a2xn−2−⋯+(−1)nan where n⩾1;
(b) a1,a2,…,an are real;
(c) the n (not necessarily distinct) roots of the equation P(x)=0 are a1,a2,…,an.
(i) Find all reflexive polynomials of degree less than or equal to 3.
(ii) For a reflexive polynomial with n>3, show that
2a2=−a22−a32−⋯−an2.
Deduce that, if all the coefficients of a reflexive polynomial of degree n are integers and an=0, then n⩽3.
(iii) Determine all reflexive polynomials with integer coefficients.
Hint
(i) Degree 1, x−a1=0 has root x=a1 and so x−a1 is reflexive. [B1]
Degree 2, x2−a1x+a2=0 has to have roots a1 and a2 to be reflexive.
Thus a12−a1a1+a2=0 (giving a2=0) and a22−a1a2+a2=0 [M1] which is consistent. Thus, x2−a1x (+ 0) [B1]
a1=a1+a2 and a2=a1a2 giving a2=0 and consistent for any a1 [M1]
Thus, x2−a1x (+ 0) [B1]
Degree 3, x3−a1x2+a2x−a3=0
a1=a1+a2+a3
a2=a1a2+a2a3+a3a1
a3=a1a2a3 [M1]
The first equation implies that a2+a3=0 [M1] or in other words a2=−a3
This result substituted into the second equation implies that a2=a2a3 [M1]
Continuing a2a3−a2=0, a2(a3−1)=0 [M1] so either a2=0 and thus a3=0 in which case the equations are consistent for any a1 [M1] or a3=1 and thus a2=−1 and a1=−1 from the third equation. [M1]
Alternatively, from the third equation, a3=a1a2a3, a1a2a3−a3=0, a3(a1a2−1)=0
so a3=0, a2=0 and consistent for any a1 or a1a2=1 In the latter case it is simpler to use equation two route again.
Yielding x3−a1x2 [A1] or x3+x2−x−1 [A1]
Alternative approach a13−a1a12+a2a1−a3=0 (A), a23−a1a22+a2a2−a3=0 (B), and
a33−a1a32+a2a3−a3=0 (C) [M1]
So (A) implies a3=a1a2 and thus (B) becomes a23−a1a22+a22−a1a2=0
x3−a1x2−x+a1=(x−a1)(x2−1) for which the equation has roots a1 and ±1. Thus a3=−a1=1 so the polynomial is x3+x2−x−1 for which the equation has roots −1, −1 and 1 and so is reflexive M1 A1
or a2=a1 and so A gives a3=a12 But C is a33−a1a32+a2a3−a3=0 and so
which gives a root of zero plus the roots of the bracketed expression. Thus we require the bracketed expression to be itself a reflexive polynomial. This can only happen if either the bracketed expression is of degree 3 or, in turn, an−1=0, and so on. E1
Hence, we have x−a1, x2−a1x (+ 0), x3−a1x2, (with a1 integer), x3+x2−x−1, or these multiplied by xr.
That is (x−a1)xr or (x+1)2(x−1)xr with r=0,1,2,…A1[4]
Degree 1:P(x)=x−a1. The equation x−a1=0 has root x=a1. The single coefficient is a1, which equals the single root. So x−a1 is reflexive for every real a1.
Degree 2:P(x)=x2−a1x+a2. By Vieta’s formulas, the roots a1,a2 satisfy:
a1+a2=a1⇒a2=0(1)
a1⋅a2=a2⇒a2(a1−1)=0(2)
From (1), a2=0, which satisfies (2) for any a1. So the reflexive polynomials of degree 2 are:
x2−a1xfor any real a1
Verification:x2−a1x=x(x−a1) has roots 0 and a1, which are indeed a2=0 and a1. ✓
Degree 3:P(x)=x3−a1x2+a2x−a3. By Vieta’s formulas:
a1+a2+a3=a1⇒a2+a3=0⇒a3=−a2(3)
a1a2+a2a3+a3a1=a2(4)
a1a2a3=a3⇒a3(a1a2−1)=0(5)
From (5), either a3=0 or a1a2=1.
Subcase a3=0: From (3), a2=0. Equation (4) becomes 0=0. So a1 is free, giving the polynomial x3−a1x2.
Verification:x3−a1x2=x2(x−a1) has roots 0,0,a1, which are a2=0, a3=0, a1. ✓
Subcase a1a2=1, a3=0: From (3), a3=−a2. Substitute into (4):
a1a2+a2(−a2)+(−a2)a1=a2
a1a2−a22−a1a2=a2
−a22=a2
a2(a2+1)=0
Since a3=−a2=0, we have a2=0, so a2=−1. Then a3=−(−1)=1 and a1=a21=−1.
Verification:P(x)=x3+x2−x−1=(x+1)2(x−1), with roots −1,−1,1. These are a1=−1, a2=−1, a3=1. ✓
So the reflexive polynomials of degree ⩽3 are:
x−a1,x2−a1x,x3−a1x2(any real a1),x3+x2−x−1
Part (ii)
For a reflexive polynomial of degree n, the roots are a1,a2,…,an. By Vieta’s formulas:
a1=∑r=1nar⇒∑r=2nar=0(6)
a2=∑1⩽i<j⩽naiaj(7)
Consider the square of ∑r=2nar:
(∑r=2nar)2=∑r=2nar2+2∑2⩽i<j⩽naiaj
From (6), the left side is 0, so:
0=∑r=2nar2+2∑2⩽i<j⩽naiaj(8)
Now we relate ∑2⩽i<j⩽naiaj to a2. From (7):
a2=∑1⩽i<j⩽naiaj=a1∑r=2nar+∑2⩽i<j⩽naiaj
By (6), a1∑r=2nar=a1⋅0=0, so:
a2=∑2⩽i<j⩽naiaj
Substituting into (8):
0=∑r=2nar2+2a2
2a2=−∑r=2nar2=−a22−a32−⋯−an2
2a2=−a22−a32−⋯−an2■
Deduction: Rearranging:
a22+2a2+a32+⋯+an2=0
(a2+1)2+a32+⋯+an2=1
Since squares are non-negative:
a32+a42+⋯+an2=1−(a2+1)2⩽1
If all coefficients are integers and an=0, then an2⩾1. Combined with a32+⋯+an2⩽1, we must have an2=1 and a3=a4=⋯=an−1=0.
From (6): a2+an=0 (since a3=⋯=an−1=0), so a2=−an=∓1.
From a1a2⋯an=an (the product of roots equals an), dividing by an=0:
a1a2⋯an−1=1
Since a3=⋯=an−1=0 (and n>3 means at least one such coefficient exists), the left side is 0. But the right side is 1. Contradiction.
Therefore, if all coefficients are integers and an=0, then n⩽3. ■
Part (iii)
From part (ii), any reflexive polynomial with integer coefficients and n>3 must have an=0. Then:
P(x)=xn−a1xn−1+⋯+(−1)n−1an−1x=x⋅Q(x)
where Q(x)=xn−1−a1xn−2+⋯+(−1)n−1an−1.
Since P(x)=0 has root x=0=an, and the remaining roots are a1,…,an−1, the polynomial Q(x) must have roots a1,…,an−1. Moreover, Q(x) has the correct sign alternation (the coefficient of xn−1−k is (−1)kak). So Q(x) is itself reflexive of degree n−1.
By induction, we keep factoring out x until we reach a reflexive polynomial of degree ⩽3 with integer coefficients. From part (i), these are:
Degree 1:x−a1 with a1∈Z
Degree 2:x2−a1x with a1∈Z
Degree 3:x3−a1x2 with a1∈Z, or x3+x2−x−1=(x+1)2(x−1)
But x2−a1x=x(x−a1) is just x times a degree-1 reflexive, and x3−a1x2=x2(x−a1) is x2 times a degree-1 reflexive. So the complete list of reflexive polynomials with integer coefficients is:
(x−a1)⋅xror(x+1)2(x−1)⋅xrwith a1∈Z,r=0,1,2,…
Equivalently: xr(x−a1) for any integer a1 and non-negative integer r, or xr(x+1)2(x−1) for any non-negative integer r. ■
5 (i) Let
f(x)=x2+px,
where p is a non-zero constant. Sketch the curve y=f(x) for x⩾0 in the case p>0.
(ii) Let
I=∫(b2−y2)c2−y21dy,
where b and c are positive constants. Use the substitution y=x2+pbx, where p is a suitably chosen constant, to show that
I=∫b2+(b2−c2)x21dx.
Evaluate
∫12(3−y2)2−y21dy.
[ Note:∫a2+x21dx=a1tan−1ax+ constant. ]
Hence evaluate
∫211(3y2−1)2y2−1ydy.
(iii) By means of a suitable substitution, evaluate
∫211(3y2−1)2y2−11dy.
Since p>0, f′(x)>0 for all x⩾0, so f is strictly increasing.
We can also note that for x>0:
f′′(x)=(x2+p)5/2−3px<0
so the curve is concave (curving towards the asymptote from below).
The sketch shows the curve passing through the origin, increasing monotonically, and approaching the horizontal asymptote y=1 from below. The curve is concave throughout x>0.
Part (ii)
Note: The question paper contains a misprint. The substitution should be y=x2+pcx (with c not b). We use the corrected substitution.
So:
∫1/21(3y2−1)2y2−1ydy=∫21(3−x2)2−x2x2⋅(−x21)dx=∫12(3−x2)2−x21dx
This is identical to the integral evaluated above, so:
∫1/21(3y2−1)2y2−1ydy=33π
Part (iii)
We evaluate ∫1/21(3y2−1)2y2−11dy.
The integrand involves y2 throughout (no y in the numerator), so the substitution y=1/x from part (ii) does not reduce it to the same form. Instead we use y=x2+pbx with suitably chosen b and p.
With y=x2+pbx:
y2=x2+pb2x2
3y2−1=x2+p3b2x2−(x2+p)=x2+p(3b2−1)x2−p
2y2−1=x2+p2b2x2−(x2+p)=x2+p(2b2−1)x2−p
dxdy=(x2+p)3/2bp
The integral becomes:
∫(3b2−1)x2−px2+p⋅(2b2−1)x2−px2+p⋅(x2+p)3/2bpdx=∫[(3b2−1)x2−p](2b2−1)x2−pbpdx
We choose 2b2−1=0, i.e., b2=21, b=21, and p=−1.
Then:
bp=21⋅(−1)=−21
(3b2−1)x2−p=21x2+1=2x2+2
(2b2−1)x2−p=0+1=1
The integral simplifies to:
∫2x2+2⋅1−1/2dx=∫x2+2−2dx
6 The point P in the Argand diagram is represented by the the complex number z, which satisfies
zz∗−az∗−a∗z+aa∗−r2=0.
Here, r is a positive real number and r2=a∗a. By writing ∣z−a∣2 as (z−a)(z−a)∗, show that the locus of P is a circle, C, the radius and the centre of which you should give.
(i) The point Q is represented by ω, and is related to P by ω=z1. Let C′ be the locus of Q. Show that C′ is also a circle, and give its radius and centre.
If C and C′ are the same circle, show that
(∣a∣2−r2)2=1
and that either a is real or a is imaginary. Give sketches to indicate the position of C in these two cases.
(ii) Suppose instead that the point Q is represented by ω, where ω=z∗1. If the locus of Q is C, is it the case that either a is real or a is imaginary?
Hint
∣z−a∣2=(z−a)(z−a)∗=(z−a)(z∗−a∗)=zz∗−az∗−a∗z+aa∗
M1
So z satisfies ∣z−a∣2=r2 which means that the locus of P is a circle (C) centre a and radius r (which does not pass through the origin as r2=aa∗.) A1 [2]
(i) As w=z1, z=w1 so
w1w∗1−aw∗1−a∗w1+aa∗−r2=0
Recall that ∣z−a∣2=(z−a)(z−a)∗=(z−a)(z∗−a∗)=zz∗−az∗−a∗z+aa∗.
The given equation zz∗−az∗−a∗z+aa∗−r2=0 is therefore equivalent to:
∣z−a∣2=r2
This is the equation of a circle C with centre a and radius r.
Since r2=aa∗=∣a∣2, the circle does not pass through the origin. ■
Part (i)
Since ω=1/z, we have z=1/ω and z∗=1/ω∗. Substituting into the equation of C:
ω1⋅ω∗1−a⋅ω∗1−a∗⋅ω1+aa∗−r2=0
Multiplying through by ωω∗:
1−aω−a∗ω∗+(aa∗−r2)ωω∗=0
Since r2=aa∗, we may divide by (aa∗−r2):
ωω∗−aa∗−r2aω−aa∗−r2a∗ω∗+aa∗−r21=0
Note that aa∗−r2a=(aa∗−r2a∗)∗ since aa∗−r2 is real.
Completing the square:
ω−aa∗−r2a∗2=aa∗−r2a∗2−aa∗−r21=(aa∗−r2)2aa∗−aa∗−r21
=(aa∗−r2)2aa∗−(aa∗−r2)=(aa∗−r2)2r2
So C′ is a circle with centre aa∗−r2a∗ and radius ∣aa∗−r2∣r. ■
If C and C′ are the same circle:
The centres must be equal:
a=aa∗−r2a∗...(1)
The radii must be equal:
r2=(aa∗−r2)2r2...(2)
From (2), since r=0, we divide by r2:
(aa∗−r2)2=1⟹(∣a∣2−r2)2=1
So ∣a∣2−r2=±1.
Case 1:∣a∣2−r2=1
From (1): a=1a∗=a∗.
Writing a=c+di: c+di=c−di, so d=0 and a is real.
Case 2:∣a∣2−r2=−1
From (1): a=−1a∗=−a∗.
Writing a=c+di: c+di=−c+di, so c=0 and a is imaginary.
Sketches:
When a is real (Case 1): ∣a∣2=r2+1, so ∣a∣>r. The circle C is centred on the real axis at a with radius r<∣a∣, so it does not enclose the origin.
When a is imaginary (Case 2): ∣a∣2=r2−1, so r>∣a∣. The circle C is centred on the imaginary axis at a with radius r>∣a∣, so it encloses the origin.
Part (ii)
Now ω=1/z∗, so z∗=1/ω and z=1/ω∗. Substituting:
ω∗1⋅ω1−a⋅ω1−a∗⋅ω∗1+aa∗−r2=0
Multiplying by ωω∗:
1−aω∗−a∗ω+(aa∗−r2)ωω∗=0
Dividing by (aa∗−r2):
ωω∗−aa∗−r2a∗ω−aa∗−r2aω∗+aa∗−r21=0
Completing the square:
ω−aa∗−r2a2=(aa∗−r2)2∣a∣2−aa∗−r21=(aa∗−r2)2∣a∣2−(aa∗−r2)=(aa∗−r2)2r2
So the locus of Q is a circle with centre ∣a∣2−r2a and radius ∣a∣2−r2r.
If this is to be the circle C:
Centre: a=∣a∣2−r2a, which gives ∣a∣2−r2=1 (since a=0).
Radius: r2=(∣a∣2−r2)2r2, which gives (∣a∣2−r2)2=1.
These are consistent: ∣a∣2−r2=1 satisfies (∣a∣2−r2)2=1.
But the centre condition a=∣a∣2−r2a reduces to ∣a∣2−r2=1, which places no constraint on the argument of a. Any complex number with ∣a∣2=r2+1 works.
For example, a=1+i with r=1 gives ∣a∣2=2=r2+1. The circle C has centre 1+i and radius 1. Since ∣a∣2−r2=1, the transformed circle C′ also has centre 11+i=1+i and radius 11=1, so C′=C. But a=1+i is neither real nor imaginary.
Therefore, it is not the case that either a is real or a is imaginary. ■
7 The Devil’s Curve is given by
y2(y2−b2)=x2(x2−a2),
where a and b are positive constants.
(i) In the case a=b, sketch the Devil’s Curve.
(ii) Now consider the case a=2 and b=5, and x⩾0,y⩾0.
(a) Show by considering a quadratic equation in x2 that either 0⩽y⩽1 or y⩾2.
(b) Describe the curve very close to and very far from the origin.
(c) Find the points at which the tangent to the curve is parallel to the x-axis and the point at which the tangent to the curve is parallel to the y-axis.
Sketch the Devil’s Curve in this case.
(iii) Sketch the Devil’s Curve in the case a=2 and b=5 again, but with −∞<x<∞ and −∞<y<∞.
Factorising the left side as a difference of squares:
(y2−x2)(y2+x2)=a2(y2−x2).
(y2−x2)(y2+x2−a2)=0.
This gives two factors:
y2−x2=0, i.e.\ y=±x (a pair of straight lines through the origin), or
x2+y2=a2 (a circle of radius a centred at the origin).
So the Devil’s Curve with a=b consists of the circle x2+y2=a2 together with the lines y=x and y=−x.
Part (ii)(a)
With a=2 and b=5, the equation is
y2(y2−5)=x2(x2−4).
Expanding:
y4−5y2=x4−4x2.
Rearranging as a quadratic in x2:
(x2)2−4x2−y2(y2−5)=0.
For real solutions in x, the discriminant must be non-negative:
16+4y2(y2−5)⩾0.
Dividing by 4:
4+y4−5y2⩾0,
y4−5y2+4⩾0.
Factorising:
(y2−1)(y2−4)⩾0.
Since y⩾0, we can write this as
(y−1)(y+1)(y−2)(y+2)⩾0.
As y+1>0 and y+2>0 for y⩾0, the sign is determined by (y−1)(y−2). This product is non-negative when y⩽1 or y⩾2. Therefore either 0⩽y⩽1 or y⩾2.
Part (ii)(b)
Very close to the origin: When x and y are small, the higher-order terms x4 and y4 are negligible compared to x2 and y2. The equation
y4−5y2=x4−4x2
becomes approximately
−5y2≈−4x2,
so y2≈54x2, giving y≈52x (taking the positive root since y⩾0). The curve passes through the origin with slope 52.
Very far from the origin: When x and y are large, the x4 and y4 terms dominate. The equation becomes approximately
y4≈x4,
so y≈x. The curve approaches the line y=x asymptotically.
Part (ii)(c)
Differentiating y4−5y2=x4−4x2 implicitly with respect to x:
(4y3−10y)dxdy=4x3−8x.
Tangent parallel to the x-axis (dxdy=0):
4x3−8x=0,
4x(x2−2)=0,
so x=0 or x=2 (taking x⩾0).
x=0: y4−5y2=0 gives y2(y2−5)=0, so y=0 or y=5.
x=2: y4−5y2=4−8=−4, so y4−5y2+4=0, giving (y2−1)(y2−4)=0, so y=1 or y=2.
At (0,0) the tangent behaviour is degenerate (the curve passes through the origin with a definite slope 52, but dxdy takes the indeterminate form 0/0). Excluding (0,0), the points with horizontal tangent are:
(0,5),(2,1),(2,2).
Tangent parallel to the y-axis (dydx=0, equivalently 4y3−10y=0):
2y(2y2−5)=0,
so y=0 or y=25.
y=0: x4−4x2=0 gives x2(x2−4)=0, so x=0 or x=2.
y=25: y2=25, so y2(y2−5)=25⋅(−25)=−425. Then x4−4x2+425=0, with discriminant 16−25=−9<0, so no real x exists.
Excluding (0,0), the point with vertical tangent is (2,0).
Sketch (first quadrant): The curve has two branches. The first branch starts at the origin (slope 52), passes through (2,1) (horizontal tangent), reaches (2,0) (vertical tangent), then continues to (0,5) (horizontal tangent). The second branch starts near the line y=x for large values and comes down to (2,2) (horizontal tangent), then connects to the asymptotic behaviour.
Part (iii)
The full curve for all x,y is obtained by exploiting the symmetry of the equation. Since the equation involves only x2 and y2 (even powers), the curve is symmetric under x→−x and y→−y independently. That is, it is symmetric in both the x-axis and the y-axis (and hence also under 180∘ rotation about the origin).
Reflecting the first-quadrant sketch into all four quadrants, the full Devil’s Curve consists of:
A closed loop surrounding the origin, passing through (±2,0) and (0,±5), with four “indentations” at (±2,±1).
A separate component in each quadrant for y⩾2, approaching the lines y=±x asymptotically. These four outer branches form a shape resembling a four-pointed star extending to infinity.
8 A pyramid has a horizontal rectangular base ABCD and its vertex V is vertically above the centre of the base. The acute angle between the face AVB and the base is α, the acute angle between the face BVC and the base is β and the obtuse angle between the faces AVB and BVC is π−θ.
(i) The edges AB and BC are parallel to the unit vectors i and j, respectively, and the unit vector k is vertical. Find a unit vector that is perpendicular to the face AVB.
Show that
cosθ=cosαcosβ.
(ii) The edge BV makes an angle ϕ with the base. Show that
Set up coordinates with the centre of the base at the origin W=(0,0,0). Let the half-lengths of the rectangle be a along i and b along j, and let the height of the pyramid be h. Then:
A=(−a,−b,0),B=(a,−b,0),C=(a,b,0),V=(0,0,h).
Normal to face AVB: The midpoint of AB is M=(0,−b,0). The direction from M to V is (0,b,h), which makes angle α with the base (the horizontal), so tanα=bh.
We need a vector perpendicular to face AVB that lies in the plane perpendicular to AB. Since AB is along i, the normal has no i-component and lies in the yz-plane.
Taking the direction from M to V as (0,b,h), a perpendicular direction in the yz-plane is (0,−h,b). Normalising:
n^1=b2+h21(0,−h,b).
This normal points inward (toward the interior of the pyramid). We can verify the angle: n^1⋅k=b2+h2b=cosα, confirming α is the angle between the face and the base.
Normal to face BVC: Similarly, the midpoint of BC is N=(a,0,0), and the direction from N to V is (−a,0,h), making angle β with the base, so tanβ=ah.
A perpendicular to BC (which is along j) in the xz-plane is (h,0,a). Normalising:
n^2=a2+h21(h,0,a).
Angle between the faces: The obtuse dihedral angle between faces AVB and BVC is π−θ, so the acute angle between the inward normals is θ. Therefore: