Page 394 - Control Theory in Biomedical Engineering

360 Appendices for Chapter 2

k12=3.6; % the rate constant of cholesterol exchange between
compartments (from compartment I to compartment II)
k21=1.0; % the rate constant of cholesterol exchange between
compartments (from compartment II to compartment I)
mdiet=0; % dietary cholesterol, 0-3g daily intake, but only up to
0.33g/day is absorbed
min=0.85; % cholesterol obtained from intestine
mout=1.2; % cholesterol used as a precursor of bile acid and being
a part of bile
mtis=0.243; % cholesterol transported to muscles and peripheral
tissue; 350mg/day
v1=6.5; % volume of blood plasma in the liver compartment in dl
v2=23.5; % volume of blood plasma in the peripheral blood com-
partment in dl
% in doi:10.1152/advan.00141.2014 v1 and v2 referred to whole blood
inf = msgbox('Before you run the calculations you have to input: the
number of time intervals, tmax for each interval, values of model
parameters.','Info','warn');
pause(2)
disp('Physiological Parameters:')
disp(strcat('1. k=',num2str(k)))
disp(strcat('2. k12=',num2str(k12)))
disp(strcat('3. k21=',num2str(k21)))
disp(strcat('4. mdiet=',num2str(mdiet)))
disp(strcat('5. min=',num2str(min)))
disp(strcat('6. mout=',num2str(mout)))
disp(strcat('7. mtis=',num2str(mtis)))
prompty = 'Do you want to disturb this system? Y/N: ';
stra1 = input(prompty, 's');
ada=0;
if stra1=='Y'j stra1=='y'
ada=1;
end
while ada>0
ada=0;
str2=input('enter the number of parameter you want to change:');
switch str2
case 1
str3=input('Enter a new value:');
k=str3;