Can I get help on building this code? Thank you!

Provided is the BM code for the injection and orally administered drug Theophylline. Theophylline is a component in tea, related to caffeine, and used in high concentrations as a bronchodilator to treat Asthma and COPD. The volume of distribution is 35 L and the clearance kc is 2.88 L/hr.

1a) Convert the Berkeley Madonna code for the injection to Matlab (include your m.file). Constants are provided in mg and L. Runtime is from 0 to 96 h, with a dt of 0.02. The initial dose is 100 mg, administered every 24 hrs. Provide a plot of the concentration using the data provided (include graphs in a word-file).

1b) Modify the dose and the number of injections to establish plasma concentrations (C) in a range of 10-20 mg/L (in this range the drug is effective and has manageable toxic side effects). It is understood that some injections maybe needed to establish a more steady behavior within the desired therapeutic range. Provide the dose and interval and plot the concentration.

1c) Implement the model for oral delivery which includes the GI tract absorbing the compound (include your m.file). Provide a plot of the concentration C using the data provided. Runtime is again from 0 to 96 h, with a dt of 0.02, and an initial dose of 100 mg.

1d) Modify the dose and the number of injections to establish plasma concentrations (C) in a range of 10-20 mg/L. It is understood that multiple pills maybe needed to establish a more steady behavior within the therapeutic range. Provide the dose and interval and plot the concentration C.

1e) Summarize and compare your findings.

{Injection Model}

{Reservoirs} d/dt (m) = + Inj - Ex INIT m = 0

{Flows} Inj = PULSE(Dose,1,Interval) Ex = kc*C

{Functions} Vd = 35 C = m/Vd kc = 2.88 Dose = 100 Interval = 24

{Globals} {End Globals}

{Pill-GI-Model}

{Reservoirs} d/dt (m) = - Ex + Egm INIT m = 0 d/dt (GI) = - Egm + Inj INIT GI = 0

{Flows} Ex = kc*C Egm = kgm*GI Inj = Pulse(Dose,1,Interval)

{Functions} Vd = 35 C = m/Vd kc = 2.88 Dose = 100 Interval = 24 kgm = 0.1

{Globals} {End Globals}

Provided is the BM code for the injection and orally administered drug Theophylline. Theophylline is a component in tea, related to caffeine, and used in high concentrations as a bronchodilator to treat Asthma and COPD. The volume of distribution is 35L and the clearance kc is 2.88L/hr. 1a) Convert the Berkeley Madonna code for the injection to Matlab (include your m.file). Constants are provided in mg and L. Runtime is from 0 to 96h, with a dt of 0.02. The initial dose is 100mg, administered every 24hrs. Provide a plot of the concentration using the data provided (include graphs in a word-file). 1b) Modify the dose and the number of injections to establish plasma concentrations (C ) in a range of 1020mg/L (in this range the drug is effective and has manageable toxic side effects). It is understood that some injections maybe needed to establish a more steady behavior within the desired therapeutic range. Provide the dose and interval and plot the concentration. 1c) Implement the model for oral delivery which includes the Gl tract absorbing the compound (include your m.file). Provide a plot of the concentration C using the data provided. Runtime is again from 0 to 96h, with a dt of 0.02 , and an initial dose of 100mg. 1d) Modify the dose and the number of injections to establish plasma concentrations (C) in a range of 1020mg/L. It is understood that multiple pills maybe needed to establish a more steady behavior within the therapeutic range. Provide the dose and interval and plot the concentration C. 1e) Summarize and compare your findings. Inj=PULSE(Dose,1,Interval)Ex=kcC{Functions}Vd=35C=m/Vdkc=2.88Dose=100 \{Reservoirs } d/dt(m)= Ex + Egm INIT m=0 d/dt(GI)= Egm + Inj INIT GI =0 { Flows } Ex = kc*C Egm = kgm*GI Inj = Pulse(Dose,1,Interval) \{Functions } Vd =35 C=m/Vd kc =2.88 Dose =100 Interval =24 kgm =0.1 \{Globals } \{End Globals } \{Reservoirs d/dt(m)=Ex+ Egm INIT m=0 d/dt(GI)=Egm+Inj INIT GI =0 \{Flows Ex=kcC Egm=kgmGI Inj = Pulse(Dose,1,Interval) \{Functions Vd=35 C=m/Vd kc=2.88 Dose =100 Interval =24 kgm=0.1 \{Globals\} \{End Globals }