|
|||||
Jodium-124 is een radioactieve
stof die onder het uitzenden van straling langzaam vervalt naar
Tellurium-124. Het wordt bijvoorbeeld gebruikt om schildklierkanker op
te sporen. Een patiënt krijgt een hoeveelheid Jodium-124 toegediend en
die verzamelt zich na een poosje vooral in de schildklier. Met een zgn. PET-CT scan kan dan de straling van dat Jodium gemeten worden en daarmee kunnen tumoren worden opgespoord. De halfwaardetijd van Jodium-124 is 4,18 dagen. Dat betekent dat elke dag ongeveer 15,3% van de hoeveelheid Jodium verandert in Tellurium (reken dat zelf maar na). |
|||||
Stel dat een patiënt op t = 0 een hoeveelheid van 12 milligram I-124 krijgt toegediend, en dat we daarna elke dag meten hoeveel gram Tellurium in zijn lichaam aanwezig is. Dat begint op 0, en al uiteindelijk naar de 12 toelopen, als langzamerhand alle Jodium vervalt. De grafiek zal er daarom uitzien als hieronder. | |||||
|
|||||
De recursievergelijking is
Tn = Tn - 1 + 0,153 • (12 - Tn
-1) Die laatste term geeft aan dat 15,3% van de Jodium (12 - T) wordt omgezet in Tellurium, dus er bijkomt als Tellurium. Deze soort groei, waarbij de waarde naar een grenswaarde toeloopt heet asymptotische groei (of ook wel begrensde groei). Een recursievergelijking daarvan ziet er altijd zó uit: |
|||||
|
|||||
G is de grenswaarde waar de
hoeveelheid uiteindelijk zal eindigen. c is één of andere
constante, die per geval verschillend is en heet ook wel de
groeivoet. Er geldt wel altijd 0 < c < 1. In plaats van toename kan er natuurlijk ook afname naar een bepaalde grenswaarde G zijn. Al we bijvoorbeeld in bovenstaand voorbeeld de hoeveelheid Jodium-124 hadden bekeken, dan was die begonnen bij 12 mg en afgenomen naar de grenswaarde G = 0. |
|||||
Een directe formule. | |||||
Nou, die is makkelijk te vinden,
als je je maar realiseert dat dit soort groei een speciaal geval
is van de al eerder besproken lineaire lineaire
differentievergelijkingen. Kijk maar: un = un-1 + c • (G - un - 1) ⇒ un = un - 1 + cG - c • un - 1 ⇒ un = (1 - c)• un - 1 + cG Daar staat gewoon un = a • un-1 + b met a = 1 - c en b = cG. Daarbij moet wel gelden dat 0 < a < 1. De directe formule hebben we bij lineaire differentievergelijkingen in opgave1 al afgeleid. In dit geval wordt dat: |
|||||
|
|||||
Voor bovenstaand voorbeeld geeft dat:
T(n) = 12(1 - 0,847n) en hiernaast zie je hoe mooi die grafiek inderdaad past bij de berekende punten. |
|||||
OPGAVEN | |||||