The left hand column displays the number of hours since the immune system came into effect, and the right hand column displays the particle count. On top the formula input by me is displayed. 122 hours after the immune system responds the particle count is 1.0071 × 1012 which I feel is pretty close to 1 × 1012, so I keep my answer as 122 hours. However, this is the amount of time it would take for the victim to die from the time the immune system responds. To determine the time of death from first infected I need to add 26.58 and 122, which gives us 148.6 (4 s.f.) hours.
At this point I decide to see the effect of an antiviral medication. I postulate that the medicine will not have an effect on the growth of the particles but it, along with the immune system, will eliminate 1,200,000 particles every hour. Despite the huge number of particles being expelled from the victim’s system I believe that there is a point where the growth of the particles outnumber expulsion of particles leading to inevitable death. I assume the point when the number of particles are so large that the medicine is useless, is between 9 million and 10 million. Therefore I decide to plug both numbers into the following formula:
un(1.60.25) – 1,200,000
Plugging 9,000,000 I find that:
9,000,000(1.60.25) – 1,200,000 = 8922000 (4 s.f.)
This shows that the number of particles in the system have actually decreased implying that the victim is recovering.
However when I plug 10,000,000 into the equation I find that:
10,000,000(1.60.25) – 1,200,000 = 10,050,000 (4 s.f.)
This, on the other hand, displays that despite the medicine being in the victim’s system (seen by the 1,200,000 particle expulsion) the victim’s particle count has still increased deeming the medicine useless. The following graphs can illustrate this further:
The graph above shows how if the medicine is administered at 9,000,000 particles, the victim begins to recover. However, the graph below shows how the victim’s particle count still rises if the medicine is administered at 10,000,000 particles.
These two graphs display that the medicine MUST be administered before 9,000,000 particles are reached if the victim wants to make a full recovery.
Now I presume that realistically a person requires time to adapt to medicine so therefore I decided to give the victim a 4 hour time period to adapt to the medicine during which he is on continuous intravenous dosing, whereby he gets a drip of medicine every minute. However, I presume that during this time the kidney is removing 2.5% of the medicine per hour. I also assume that the victim requires at least 90 micrograms of medicine at the end of the 4 hours to begin and maintain the elimination of 1,200,000 particles. Using this information I decide to determine how much of the medicine the patient receives every minute and every hour to reach 90 micrograms by the end of 4 hours.
The victim loses 2.5% of the medicine every hour, but since takes place over the course of the hour I cannot simply subtract 2.5% from total amount of medicine at the end of every hour. However, finding the exact amount lost is very complicated so in order to find a decently accurate I presume that the decay is linear, therefore (2.5 ÷ 60 = 0.0416)% of the medicine is eliminated every minute. This enables to create the following equation:
D(0.9996240) + D(0.9996239) + … + D(0.99961) + D
This formula helps determine the dosage because D represents the dosage per minute, and multiplying with 0.999583 is the equivalent of removing 0.0416%. This number has a power of 240 because 0.0416% is removed 240 times by the end. For the next dose injected a minute later, this takes place 239 times because it is in the body for 1 minute less. Similarly, the dose injected 1 minute before the end of the 4 hours only has a power of 1 because 0.0416% is only removed one time before the end.
This equation can then be rearranged to put the smallest powers at the beginning and the largest powers at the end but either way I have a geometric sequence. I then use the formula to find the sum of a geometric sequence:
a(rn – 1) ÷ (r – 1) where a is the initial term (D, if the equation is rearranged) and r is the ratio of multiplication (0.9996, if the equation is rearranged) and n is the number of times the ratio is multiplied (240).
This equation must multiply out to 90 so that 90 micrograms are ensured in the system at the end of the four hours. Plugging the values in I get the following:
D(0.9996240 – 1) ÷ (0.9996 – 1) = 90
D(-0.095) ÷ (-0.000416) = 90
D(-0.095) = 90 × (-0.000416) = -0.038
D = (=0.038) ÷ (-0.095) = 0.39 micrograms
I presume that once the level of medication has reached 90 micrograms, at the end of the 4 hours, the body will have adapted to the medicine therefore the victim will be taken off the intravenous phase and he will be given injections every 4 hours. Nevertheless, the kidneys will still be removing 2.5% of the medication every hour, therefore I need to calculate the extra dosage which must be given every 4 hours so that the amount of medicine in the victims bloodstream is never lower than 90 micrograms. The following equation can help me come up with my answer:
(90 + D) × 0.9754 ≥ 90
(90 + D) × 0.904 ≥ 90
90 + D ≥ 90 ÷ 0.904 ≥ 99.59
D ≥ 99.59 – 90
D ≥ 9.59
I used this equation because I presumed that a dose, D, was added before the 4 hour period begun, then 2.5% was removed every hour, for 4 hours.
Now I decide to go into more detail regarding the last point when the victim can take medicine and survive. To do this first I have to find out the number of particles the person must have in his body at this point. To do this I looked at each aspect of the following formula:
un(1.60.25) – 1,200,000
This formula shows the number of particles in the next hour. If I have un as the starting amount of particles and the finishing amount of particles is 1 less than un, then un is the number of particles the person would have in his body at the last possible time to take medicine and recover. In other words, the gain in particles would be 1,199,999 and the loss would be 1,200,000 so 1 particle is lost starting the recovery process. Using this I can equate the formula above to:
un(1.60.25) – 1,200,000 = un – 1
un(1.60.25) – 1,200,000 - un = -1
un(1.60.25) – un = 1,199,999
un(1.60.25 – 1) = 1,199,999
un(0.1246) = 1,199,999
un = 9624426
Using this knowledge I can utilize the general formula to calculate the last possible time from the onset of the infection to start the medication.
General formula: a(r(n)) – e[(r(n) – 1) ÷ (r – 1)]
In this formula, a is the initial term, r is the ratio of multiplication, n is the number of hours, and e is the amount being eliminated by the body every hour.
I insert values in this formula to determine the time it takes to reach 9624426 from the moment the immune system responds to the virus:
a = 1,000,000
r = 1.60.25
e = 50,000
rn = 1.60.25n because of indices rules
1,000,000(1.6(0.25n)) – 50,000[(1.6(0.25n) – 1) ÷ (1.60.25 – 1)] = 9624426
Factorizing and clearing this up, I get:
20(1.6(0.25n)) – [(1.6(0.25n) – 1) ÷ (0.1246)] = 192.5 (4 s.f.)
(1.6(0.25n) – 1) ÷ (0.1246) = 20(1.6(0.25n)) – 192.5
(1.6(0.25n) – 1) = 2.49(1.6(0.25n)) – 23.99
(1.6(0.25n) – 1) - 2.49(1.6(0.25n)) = - 23.99
1.6(0.25n) - 2.49(1.6(0.25n)) = -22.99
1.6(0.25n)(1 – 2.49) = -22.99
1.6(0.25n)(-1.49) = -22.99
1.6(0.25n) = (-22.99) ÷ (-1.49) = 15.43
Using log:
(0.25n)log1.6 = log15.43
0.25n = log15.43 ÷ log1.6 = 5.82
n = 23.28 hours (last possible time from response of immune system to start medication)
However, I need to find the last possible time from the onset of infection to start the medication therefore I need to add 26.58 to 23.28 which gives us 49.87 hours.
In order to complete the analysis of this victim I also need to calculate the time it will take him to recover from this point in time. In order to do this I will manipulate the general formula again and equate it to 0 (giving the time he has 0 particles left in his system):
9624426(1.60.25n) – 1,200,000[(1.60.25n – 1) ÷ (1.60.25 – 1)] = 0
1,200,000[8.0203(1.60.25n) – [(1.60.25n – 1) ÷ (1.60.25 – 1)] = 0
8.02(1.60.25n) – [(1.60.25n – 1) ÷ (0.1246)] = 0
8.02(1.60.25n) = (1.60.25n – 1) ÷ (0.1246)
0.9999(1.60.25n) = (1.60.25n) – 1
0.9999(1.60.25n) - (1.60.25n) = – 1
(1.60.25n)(0.9999 – 1) = – 1
(1.60.25n)(-0.0000881) = - 1
(1.60.25n) = (- 1) ÷ (-0.00000881) = 1135073
Using log:
(0.25n)log1.6 = log1135073
0.25n = log1135073 ÷ log 1.6 = 29.66 hours
n = 118.65 hours
This displays that it will take the victim 118.65 hours from the start of the regimen of medicine to recover completely. From this I can determine that the victim is infected for a total time of 118.65 + 23.28 + 26.58 = 168.51 hours.
In my investigation I have made many assumptions which I have listed before utilizing them, throughout the course of the investigation. I think the strength of the model is that it gives each answer very precisely as the general formula doesn’t give an approximation, however this is also a weakness because it is highly unlikely that any syringe can drop 0.39 micrograms of medicine every minute. There are a lot of variables in this investigation that would make the answers very unrealistic – deeming the investigation very unreliable.
I have displayed how this affects a young adult, but this model can be modified if the victim was a child instead. This can be done by changing many of the variables – the initial amount could be 25000 instead of 10000 since a child has lower resistance, the rate of multiplication could be 3 instead of 2 at the beginning for the same reason. The immune system could kick in at 3,000,000 instead of 1,000,000 because it would take longer for the body to acknowledge the presence of the virus. The amount of particles expelled could be 20000 instead of 50000 as the child’s immune system isn’t as strong. Furthermore, the child – having a smaller body – would take it a lot less medicine dropping all those figures even further. Finally, the amount of particles needed for death would be a lot lower, 109 instead of 1012.
I will now explain how I derived the general formula. First I presume that I are just finding a general formula for the time period after the immune system kicks in so I get the following:
{[1,000,000(1.60.25) – 50,000](1.60.25) – 50,000}(1.60.25) – 50,000 …
Opening this up I get the following:
1,000,000(1.60.25)(1.60.25)(1.60.25)–50,000(1.60.25)(1.60.25)–50,000(1.60.25)–50,000 …
For the first part I can immediately make out that the number of (1.60.25) attached to the 1,000,000 are equal to the number of hours passed making it, where n is the number of hours passed:
1,000,000(1.60.25n)
For the second half we can factorize 50,000 out leaving us with:
X – 50,000{(1.60.25)(1.60.25) + (1.60.25) + 1}
X – 50,000{(1.60.25(2)) + (1.60.25(1)) + (1.60.25(0))}
I then see a geometric equation forming with a being 1 because (1.60.25(0)) = 1 and r being (1.60.25). Using the sum of a geometric formula:
a(rn – 1) ÷ (r – 1)
1[(1.60.25n) – 1] ÷ [(1.60.25) – 1]
From this we can derive the general formula after the immune system kicks in to be:
1,000,000(1.60.25n) – 50,000[(1.60.25n) – 1] ÷ [(1.60.25) – 1])
We can then use this format to derive a general formula that would work for any numbers:
P(r0.25n) – e[(r0.25n – 1) ÷ (r – 1)]
In this formula P is the starting number of particles (so 10,000 if we start from when he was first affected, or 1,000,000 if we start from when the immune system kicks in). r is the rate of multiplication every 4 hours, and e is the number of particles eliminated from the body every hour. From this formula we can plug in any values for any of the factors and determine the number of particles.