%Initial values einzelner Stufen Xv0 = [0.3*10^6 ,0.3*10^6 ,0.3*10^6 ,0.3*10^6 ,0.4*10^6 ,0.4*10^6 ,0.25*10^6 ,0.3*10^6 ,0.38*10^6 ,0.38*10^6]; Xd0 = [6*10^3 ,6*10^3 ,6*10^3 ,6*10^3 ,6*10^3 ,6*10^3 ,6*10^3 ,6*10^3 ,6*10^3 ,6*10^3];% ,0.3*10^6]; cGlc0 = [6.8 ,6.8 ,6.8 ,6.8 ,6.8 ,6.8 ,6.8 ,6.8 ,6.8 ,6.8]; cLac0 = [0.3 ,0.3 ,0.3 ,0.3 ,0.3 ,0.3 ,0.3 ,0.3 ,0.3 ,0.3]; cGln0 = [0.9 ,0.9 ,0.9 ,0.9 ,0.9 ,0.9 ,0.9 ,0.9 ,0.9 ,0.9]; cAmm0 = [1 ,1 ,1 ,1 ,1 ,1 ,1 ,1 ,1 ,1]; v0 = [0.25 ,0.5 ,1 ,2 ,2 ,10 ,50 ,80 ,400 ,2000]; %% %Parameters() Parameters(1) = 0.055; % my_max [h^-1] Maximum Specific Growth Rate Parameters(2) = 0.03; % my_d_max [h^-1] Maximum Specific Death Rate Parameters(3) = 0.002; % my_d_min [h^-1] Minimum Specific Death Rate Parameters(4) = 1.8*10^-8; % qGlc_max [mmol/cells *h] Maximum Cell-specific Glucose Uptake Rate Parameters(5) = 0.8*10^-8; % qGln_max [mmol/cells *h] Maximum Cell-specific Glutamine Uptake Rate Parameters(6) = 1.2*10^-12;% qLac_max [mmol/cells *h] Maximum Cell-specific Lactate Uptake Rate Parameters(7) = 0.19; % kGlc [mmol/L] Monod Kinetic Constant for Glucose Uptake Parameters(8) = 0.3; % kGln [mmol/L] Monod Kinetic Constant for Glutamine Uptake Parameters(9) = 0.03; % KsGlc [mmol/L] Monod Kinetic Constant for Glucose Parameters(10) = 0.03; % KsGln [mmol/L] Monod Kinetic Constant for Glutamine %% %constants %Constants(1) = 0.01; % KLys [h^-1] Constant for Cell Lysing Constants(2) = 4*10^-12; % qAmm_uptake [mmol/cells *h] Maximum Cell-specific Ammonia Uptake Rate Constants(3) = 0.4; % YAmm/Gln [mmol/mmol] Kinetic production Constant for Ammonia by consuming Glutamine Constants(4) = 1.6; %YLac/Glc [mmol/mmol] Kinetic production Constant for Lactate by consuming Glucose Constants(5) = 0; % [L] Sample Volume Constants(6) = 0.5; % Correction Factor for Ammonia uptake %% n = 10; %Anzahl an Stufen %Laufdauer t=0; y=[0,0,0,0,0,0,0]; Laufdauer =[0 72;72 144 ;144 216;216 288;288 360;360 432;432 504;504 576;576 648;648 1059]; %Laufdauer der einzelnen Stufen for k = 1:n tRange= Laufdauer(k,:); c0 = [ Xv0(k); Xd0(k); cGlc0(k); cLac0(k); cGln0(k); cAmm0(k); v0(k)]; %Vector that contains all the Initial values %% %ODE45 [tSpeicher,ySpeicher] = ode15s(@(t,c) Model_BI130_1(t,c,Parameters,Constants), tRange, c0); t = vertcat(t,tSpeicher); y= vertcat(y,ySpeicher); end %Total Cell Density Xv = y (:,1); Xd = y (:,2); Xt = Xv + Xd; %Viability Viability = (Xv./Xt)*100; [x,fval,exitflag,output] = particleswarm(@(c)Model_BI130_1([],c,Parameters,Constants),10,[],[])