Porisma

Original question & working

Logistic proportions and the first whole week · source reference

[2021/HCI/Prelim/I/10] · Question pp. 7 · Working pp. 29, 30, 31, 32, 33

Open the interactive lesson

The working and teaching notes are included in the interactive lesson.

View source review notes
Full original question, working & notes

Printed question

A flu virus is spreading in a community which has a fixed population of NN people. Scientists discover that at time tt weeks from the beginning, the rate of change of the number of people who are infected is proportional to the product of the number of people infected, II, and the number of people who are not infected.

(i) Write down a differential equation relating I,tI,t. [1]

It is noted that at the beginning, 1% of the people in the community is infected with the virus. 12 weeks later, the proportion of people infected in the community has risen to 25%.

(ii) By solving the differential equation obtained in part (i), show that I=N99(33−t/12)+1I=\frac N{99(33^{-t/12})+1}. [6]

(iii) Find the least number of weeks needed for the proportion of people infected to exceed 50%, giving your answer to the nearest integer. [2]

(iv) Sketch the graph of II against tt. Explain, with justification, what happens to the number of people infected if this situation continues indefinitely. [2]

(v) State a possible limitation of the model in part (i). [1]

[2021/HCI/Prelim/I/10]

Printed working

I′=kI(N−I)I'=kI(N-I).

Method 1: partial fractions

dt/dI=1/[kI(N−I)]=[1/I+1/(N−I)]/(kN)dt/dI=1/[kI(N-I)]=[1/I+1/(N-I)]/(kN).

t=∫[1/I+1/(N−I)]/(kN)dI=[ln⁡I−ln⁡(N−I)]/(kN)+Ct=\int[1/I+1/(N-I)]/(kN)dI=[\ln I-\ln(N-I)]/(kN)+C (0<I<N0<I<N).

kN(t−C)=ln⁡[I/(N−I)]kN(t-C)=\ln[I/(N-I)]; ekNt−kNC=I/(N−I)e^{kNt-kNC}=I/(N-I); AekNt=I/(N−I)Ae^{kNt}=I/(N-I), A=e−kNCA=e^{-kNC}.

AekNt(N−I)=IAe^{kNt}(N-I)=I; AekNtN=I(1+AekNt)Ae^{kNt}N=I(1+Ae^{kNt}); I=AekNtN/(1+AekNt)=AN/(e−kNt+A)I=Ae^{kNt}N/(1+Ae^{kNt})=AN/(e^{-kNt}+A).

At 0, I=0.01NI=0.01N: 0.01(1+A)=A⇒A=1/990.01(1+A)=A\Rightarrow A=1/99.

At 12, I=0.25NI=0.25N: 0.25(e−12Nk+1/99)=1/99⇒e−12Nk=1/33⇒Nk=ln⁡33/120.25(e^{-12Nk}+1/99)=1/99\Rightarrow e^{-12Nk}=1/33\Rightarrow Nk=\ln33/12.

I=N/99e−tln⁡33/12+1/99=N99e−tln⁡33/12+1=N99(33−t/12)+1I=\frac{N/99}{e^{-t\ln33/12}+1/99}=\frac N{99e^{-t\ln33/12}+1}=\frac N{99(33^{-t/12})+1}.

Method 2: MF27 / complete the square

dt/dI=1/[kI(N−I)]=1/{k[(N/2)2−(I−N/2)2]}dt/dI=1/[kI(N-I)]=1/\{k[(N/2)^2-(I-N/2)^2]\}.

t=1kNln⁡∣N/2+I−N/2N/2−I+N/2∣+C=1kNln⁡∣I/(N−I)∣+Ct=\frac1{kN}\ln\left|\frac{N/2+I-N/2}{N/2-I+N/2}\right|+C=\frac1{kN}\ln|I/(N-I)|+C.

The printed solution repeats the exponentiation, initial conditions A=1/99A=1/99, calibration Nk=ln⁡33/12Nk=\ln33/12, and final formula above, giving the same result.

Threshold method 1: algebra

I>0.50N⇒N/[99(33−t/12)+1]>0.50NI>0.50N\Rightarrow N/[99(33^{-t/12})+1]>0.50N.

99(33−t/12)+1<2⇒99(33−t/12)<1⇒33−t/12<1/9999(33^{-t/12})+1<2\Rightarrow99(33^{-t/12})<1\Rightarrow33^{-t/12}<1/99.

−tln⁡33/12<−ln⁡99⇒t>15.8-t\ln33/12<-\ln99\Rightarrow t>15.8. Least whole number of weeks: 16.

Threshold method 2: graphical calculator

The original calculator figures display mode NORMAL FLOAT AUTO REAL RADIAN MP. Diagram 1 displays CALC INTERSECT, and Diagram 2 displays CALC ZERO; their intersection/zero readouts are recreated in the native figures.

Plot y=1/[99(33−t/12)+1]y=1/[99(33^{-t/12})+1] with y=0.5y=0.5, or plot y=1/[99(33−t/12)+1]−0.5y=1/[99(33^{-t/12})+1]-0.5 and locate its zero. Both give threshold t=15.770433…t=15.770433\ldots, so t>15.8t>15.8 rounded and whole week 16.

As t→∞t\to\infty, 33−t/12→0⇒I→N33^{-t/12}\to0\Rightarrow I\to N. Almost everyone in the community will be infected if this continues indefinitely.

Limitations printed: the model assumes everyone will get the virus, but some may have natural immunity; or, deaths may occur during some weeks so NN is no longer constant.

Printed answer key

(i) dI/dt=kI(N−I)dI/dt=kI(N-I). (iii) 16. (iv) Native sketch below. If this situation continues indefinitely, almost everyone in the community will be infected by the virus. (v) One possible limitation is that the model assumes everyone will get the virus which in actual fact some may have natural immunity against the flu virus.

Teacher note 1

Given dI/dt=f(I)dI/dt=f(I), t=∫dI/f(I)t=\int dI/f(I). Let 1/[kI(N−I)]=(1/k)[A/I+B/(N−I)]1/[kI(N-I)]=(1/k)[A/I+B/(N-I)]; 1=A(N−I)+BI1=A(N-I)+BI. At I=0I=0, A=1/NA=1/N; at I=NI=N, B=1/NB=1/N. Hence 1/[kI(N−I)]=(1/(kN))[1/I+1/(N−I)]1/[kI(N-I)]=(1/(kN))[1/I+1/(N-I)].

Teacher note 2

Complete the square for kI(N−I)=−k(I2−NI)kI(N-I)=-k(I^2-NI), where k,Nk,N are constants.

Teacher note 3

For Diagram 1, sketch y=1/[99(33−t/12)+1]y=1/[99(33^{-t/12})+1] and y=0.5y=0.5. On GC go to 2nd [calc] → 5:intersect to find intersection between the 2 curves.

For Diagram 2, sketch y=1/[99(33−t/12)+1]−0.5y=1/[99(33^{-t/12})+1]-0.5. On GC go to 2nd [calc] → 2:zero to find intersection of curve and horizontal axis.

Diagram 10-18-population

Flu infection population

Diagram 10-18-GC-intersection

GC Diagram 1: intersection with 0.5

Diagram 10-18-GC-zero

GC Diagram 2: zero of shifted fraction

Mathematical review

(i) I′=kI(N−I)I'=kI(N-I), k>0.

(ii) Separating gives ln⁡[I/(N−I)]=kNt+C\ln[I/(N-I)]=kNt+C. At t=0 the odds are 1/99. At t=12 they are 1/3, so e12kN=33e^{12kN}=33 and kN=ln⁡33/12kN=\ln33/12. Thus I=N/[1+99e−tln⁡33/12]=N/[1+99(33)−t/12]I=N/[1+99e^{-t\ln33/12}]=N/[1+99(33)^{-t/12}].

(iii) I>N/2I>N/2 requires 99(33)−t/12<199(33)^{-t/12}<1, hence t>12ln⁡99/ln⁡33≈15.7704t>12\ln99/\ln33\approx15.7704. The nearest integer and first complete integer week are both 16.

(iv) I begins at 0.01N, increases sigmoidally, has an inflection at I=N/2 and approaches the asymptote I=N. (v) The model assumes no recovery or immunity and a fixed total population. Either assumption can fail in an actual epidemic; its mathematical limiting behavior is conditional on these assumptions.

Provenance

TJC JC2 2025 H2 Maths (9758) MSM. Native transcription checked against original visible question, working and sidebars. Question page(s): 7. Working page(s): 29, 30, 31, 32, 33. Original records remain private for QA; this native article has no PDF or raster worksheet dependency.