Page 395 - Control Theory in Biomedical Engineering
P. 395

Appendices for Chapter 2  361


                      case 2
                        str3=input('Enter a new value:');
                        k12=str3;
                      case 3
                        str3=input('Enter a new value:');
                        k21=str3;
                      case 4
                        str3=input('Enter a new value:');
                        mdiet=str3;
                      case 5
                        str3=input('Enter a new value:');
                        min=str3;
                      case 6
                        str3=input('Enter a new value:');
                        mout=str3;
                      case 7
                        str3=input('Enter a new value:');
                        mtis=str3;
                      otherwise
                      disp('bad command')
                   end
              disp('New 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)))
                      str4=input(prompty,'s');
                      if str4=='Y'j str4=='y'
                        ada=1;
                      end
                   end
                   m1=k/(mout+mtis-mdiet-min);

                   m2=(mdiet 2+mdiet*min-mdiet*mout-k*k12-2*mdiet*mtis-

              min*mtis+mout*mtis+mtis 2)/(k21*(mdiet+min-mout-mtis));
                   c1=m1/v1;
                   c2=m2/v2;
   390   391   392   393   394   395   396   397   398   399   400