1b) Use code from Berkeley Madonna listed below as a starting point to develop your own Matlab implementation. Use ode45 (RK4 solver). Of note, the code is simplified and is tuned towards a specific ethanol dose (5g) but should deliver a qualitative similar behavior to what is published in the paper. Submit the following items (you can combine your items b-e): a) Submit your Matlab code b) Include your Forrester Diagram (integrated in doc-file or pdf) c) Plot the outcome of your simulation for Cg, Cl, Cc and Cm (doc/image or pdf) d) Discuss the behavior of these parameters with respect to the model structure and parameters (doc/pdf, approximately 5 sentences). e) Change one parameter to improve clearance by the liver compartment (simulating a drug). If you do, how does this improve the overall ethanol clearance? Explain what you did and provide a simulation. Here is the Berkeley Madonna code: METHOD RK4 STARTTIME = 0 STOPTIME = 300 {minutes} DT = 0.02 Km = 0.0184 Vg = 2.4 Vl = 1.08 Vc = 11.56 Vm = 25.76 FL = 1.35 Fm = 0.95 D=0.2 {D is a function of 5g ethanol dose and 100 kg body weight} ks = -0.049*(D) + 0.0545 INIT Vs = 1 d/dt(Vs) = -ks*Vs INIT Cg = 0 INIT Cl = 0 INIT Cc = 0 INIT Cm = 0 Vmax = 0.1012 d/dt(Cg) =((((2/3)*FL)*(Cc - Cg) +(ks*Vs))/Vg) rel = (Vmax*Cl)/(Km + Cl) d/dt(Cl) = ((FL*(0.33*Cc + .66*Cg - Cl) - rel)/Vl) d/dt(Cc) = ((-FL*(Cc - Cl) - Fm*(Cc - Cm))/Vc) d/dt(Cm) = (Fm*(Cc-Cm)/Vm)

Please help!! Thank you!

1b) Use code from Berkeley Madonna listed below as a starting point to develop your own Matlab implementation. Use ode45 (RK4 solver). Of note, the code is simplified and is tuned towards a specific ethanol dose (5g) but should deliver a qualitative similar behavior to what is published in the paper. Submit the following items (you can combine your items b-e): a) Submit your Matlab code b) Include your Forrester Diagram (integrated in doc-file or pdf) c) Plot the outcome of your simulation for Cg,Cl,Cc and Cm (doc/image or pdf) d) Discuss the behavior of these parameters with respect to the model structure and parameters (doc/pdf, approximately 5 sentences). e) Change one parameter to improve clearance by the liver compartment (simulating a drug). If you do, how does this improve the overall ethanol clearance? Explain what you did and provide a simulation. STOPTIME = 300 \{minutes } DT=0.02Km=0.0184Vg=2.4Vl=1.08VC=11.56Vm=25.76FL=1.35Fm=0.95D=0.2 \{D is a function of 5g ethand ks=0.049(D)+0.0545 INIT Vs =1 d/dt(Vs)=ksVs INIT Cg =0 INIT Cl =0 INIT CC =0 INIT Cm =0 INIT Cg =0 INIT Cl =0 INIT CC =0 INIT Cm =0 Vmax =0.1012 d/dt(Cg)=((((2/3)FL)(CcCg)+(ksVs))/Vg)rel=(VmaxCl)/(Km+Cl) d/dt(Cl)=((FL(0.33Cc+.66CgCl)rel)/Vl)d/dt(Cc)=((FL(CcCl)Fm(CcCm))/Vc)d/dt(Cm)=(Fm(CcCm)/Vm)