A seminal breakthrough in HIV treatment was the use of multiple-drug combination therapy. The idea behind...

Transcribed Image Text:

A seminal breakthrough in HIV treatment was the use of multiple-drug combination therapy. The idea behind the therapy is that multiple drugs, introduced to the body at specific times, fight the virus off better than one drug alone. A key step in the development of multiple-drug therapy is determining the best time to take the next drug. One way this is done is to take the next drug when the previous drug is at its maximum in the body (and beginning to decrease from there). Drug absorption and elimination in the body is modeled using compartmental models. Consider a model whose solution is given by 11 x(t) = —— (e-¹/¹² - e-t), -t/12 (1) where r(t) is the amount of the drug in the body at time t. Assume that we want to administer the next drug when r(t) is at its maximum. We call this time tmax in what follows. Recall that we have only talked about minimizing functions so you will have to restate the problem as a minimization problem. 1 (a) Find tmax and the maximum of the function r(t) by taking the derivative (by hand) and using scipy. optimize. fsolve to solve the resulting equation (see an example of how to use this in the Homework 4 template), with an initial guess of 20 = 1.5, to find the root of x' (t). First save a' (1.5) to the variable A1 to make sure you have r' (t) defined correctly. Then save tmax to the variable A2 and the maximum of r(t) to the variable A3. Note: scipy. optimize.fsolve returns an array, make sure you get the single number out of that array for your answers here. (b) Use scipy. optimize.fminbound with the interval [0, 10] to find tmax and the maximum of r(t). Create an array with 2 elements: tmax in the first component and the maximum in the second component. Save this array to the variable A4. A seminal breakthrough in HIV treatment was the use of multiple-drug combination therapy. The idea behind the therapy is that multiple drugs, introduced to the body at specific times, fight the virus off better than one drug alone. A key step in the development of multiple-drug therapy is determining the best time to take the next drug. One way this is done is to take the next drug when the previous drug is at its maximum in the body (and beginning to decrease from there). Drug absorption and elimination in the body is modeled using compartmental models. Consider a model whose solution is given by 11 x(t) = —— (e-¹/¹² - e-t), -t/12 (1) where r(t) is the amount of the drug in the body at time t. Assume that we want to administer the next drug when r(t) is at its maximum. We call this time tmax in what follows. Recall that we have only talked about minimizing functions so you will have to restate the problem as a minimization problem. 1 (a) Find tmax and the maximum of the function r(t) by taking the derivative (by hand) and using scipy. optimize. fsolve to solve the resulting equation (see an example of how to use this in the Homework 4 template), with an initial guess of 20 = 1.5, to find the root of x' (t). First save a' (1.5) to the variable A1 to make sure you have r' (t) defined correctly. Then save tmax to the variable A2 and the maximum of r(t) to the variable A3. Note: scipy. optimize.fsolve returns an array, make sure you get the single number out of that array for your answers here. (b) Use scipy. optimize.fminbound with the interval [0, 10] to find tmax and the maximum of r(t). Create an array with 2 elements: tmax in the first component and the maximum in the second component. Save this array to the variable A4.

Answer rating: 100% (QA)

a We have the function xt 3 Differentiating this function we find x t et24 e12 1210 1271... View the full answer