Page 396 - Control Theory in Biomedical Engineering

362 Appendices for Chapter 2

disp(strcat('m1=',num2str(m1)))
disp(strcat('c1=',num2str(c1)))
disp(strcat('m2=',num2str(m2)))
disp(strcat('c2=',num2str(c2)))

n=input('Enter the number of time intervals= ');
if n<1
hg = msgbox('Wrong value. Please try again', 'Error', 'warn');
n=input('Enter the time intervals= ');
end
tmax=[];
parameter=zeros(7,n);
for i=1:n
format compact
interval=i
t=input('tmax[minute]= ');
if t<0
hg = msgbox('Wrong value. Please try again', 'Error',
'warn');
t=input('tmax[minute]= ');
end
tmax=[tmax;t];
parameter(1,i)=k;
parameter(2,i)=k12;
parameter(3,i)=k21;
parameter(4,i)=mdiet;
parameter(5,i)=min;
parameter(6,i)=mout;
parameter(7,i)=mtis;
disp('Parameters:')
disp(strcat('1. k=',num2str(parameter(1,i))))
disp(strcat('2. k12=',num2str(parameter(2,i))))
disp(strcat('3. k21=',num2str(parameter(3,i))))
disp(strcat('4. mdiet=',num2str(parameter(4,i))))
disp(strcat('5. min=',num2str(parameter(5,i))))
disp(strcat('6. mout=',num2str(parameter(6,i))))
disp(strcat('7. mtis=',num2str(parameter(7,i))))
prompt = 'Do you want to change parameters in the current time inter-
val? Y/N: ';