(* Generated by JWS Online *)

(* This is an experimental feature of JWS Online. Please report any mistakes.*)

(* Note that the following notable SBML entities or features are not supported in notebook outputyet: *)
    (* Events *)
    (* Constraints *)
    (* Units and UnitDefinitions *)
    (* AlgebraicRules *)
    (* conversionFactors *)

variables = {
species\[LetterSpace]10[t], species\[LetterSpace]11[t], species\[LetterSpace]13[t], species\[LetterSpace]14[t], species\[LetterSpace]15[t], species\[LetterSpace]17[t], species\[LetterSpace]18[t], species\[LetterSpace]2[t], species\[LetterSpace]20[t], species\[LetterSpace]21[t], species\[LetterSpace]22[t], species\[LetterSpace]23[t], species\[LetterSpace]24[t], species\[LetterSpace]25[t], species\[LetterSpace]27[t], species\[LetterSpace]28[t], species\[LetterSpace]29[t], species\[LetterSpace]30[t], species\[LetterSpace]31[t], species\[LetterSpace]33[t], species\[LetterSpace]34[t], species\[LetterSpace]4[t], species\[LetterSpace]5[t], species\[LetterSpace]6[t], species\[LetterSpace]7[t], species\[LetterSpace]8[t]
};

initialValues = {
species\[LetterSpace]10[0] == 100.0, species\[LetterSpace]11[0] == 57516.0, species\[LetterSpace]13[0] == 10000.0, species\[LetterSpace]14[0] == 100.0, species\[LetterSpace]15[0] == 100.0, species\[LetterSpace]17[0] == 20.0, species\[LetterSpace]18[0] == 100.0, species\[LetterSpace]2[0] == 3150.0, species\[LetterSpace]20[0] == 0.0, species\[LetterSpace]21[0] == 20.0, species\[LetterSpace]22[0] == 100.0, species\[LetterSpace]23[0] == 100.0, species\[LetterSpace]24[0] == 100.0, species\[LetterSpace]25[0] == 100.0, species\[LetterSpace]27[0] == 0.0, species\[LetterSpace]28[0] == 9363.0, species\[LetterSpace]29[0] == 0.0, species\[LetterSpace]30[0] == 45.0, species\[LetterSpace]31[0] == 100.0, species\[LetterSpace]33[0] == 100.0, species\[LetterSpace]34[0] == 100.0, species\[LetterSpace]4[0] == 400.0, species\[LetterSpace]5[0] == 467.0, species\[LetterSpace]6[0] == 0.0, species\[LetterSpace]7[0] == 60000.0, species\[LetterSpace]8[0] == 0.0
};

rates = {
reaction\[LetterSpace]1, reaction\[LetterSpace]10, reaction\[LetterSpace]11, reaction\[LetterSpace]12, reaction\[LetterSpace]13, reaction\[LetterSpace]14, reaction\[LetterSpace]15, reaction\[LetterSpace]16, reaction\[LetterSpace]17, reaction\[LetterSpace]18, reaction\[LetterSpace]19, reaction\[LetterSpace]2, reaction\[LetterSpace]20, reaction\[LetterSpace]21, reaction\[LetterSpace]22, reaction\[LetterSpace]23, reaction\[LetterSpace]24, reaction\[LetterSpace]25, reaction\[LetterSpace]26, reaction\[LetterSpace]27, reaction\[LetterSpace]28, reaction\[LetterSpace]29, reaction\[LetterSpace]3, reaction\[LetterSpace]30, reaction\[LetterSpace]31, reaction\[LetterSpace]32, reaction\[LetterSpace]33, reaction\[LetterSpace]34, reaction\[LetterSpace]35, reaction\[LetterSpace]4, reaction\[LetterSpace]5, reaction\[LetterSpace]6, reaction\[LetterSpace]7, reaction\[LetterSpace]8, reaction\[LetterSpace]9
};

rateEquations = {
reaction\[LetterSpace]1 -> multiplier*reaction\[LetterSpace]1\[LetterSpace]k1*species\[LetterSpace]1, reaction\[LetterSpace]10 -> function\[LetterSpace]6[reaction\[LetterSpace]10\[LetterSpace]BCRmax, reaction\[LetterSpace]10\[LetterSpace]BCRt, species\[LetterSpace]7[t], reaction\[LetterSpace]10\[LetterSpace]BS], reaction\[LetterSpace]11 -> function\[LetterSpace]7[reaction\[LetterSpace]11\[LetterSpace]HCSmax, species\[LetterSpace]7[t], reaction\[LetterSpace]11\[LetterSpace]HCSt, reaction\[LetterSpace]11\[LetterSpace]HS], reaction\[LetterSpace]12 -> function\[LetterSpace]8[reaction\[LetterSpace]12\[LetterSpace]k9, species\[LetterSpace]14[t], species\[LetterSpace]7[t]], reaction\[LetterSpace]13 -> function\[LetterSpace]9[reaction\[LetterSpace]13\[LetterSpace]k10, species\[LetterSpace]15[t], species\[LetterSpace]13[t]], reaction\[LetterSpace]14 -> function\[LetterSpace]10[reaction\[LetterSpace]14\[LetterSpace]k11, species\[LetterSpace]11[t]], reaction\[LetterSpace]15 -> reaction\[LetterSpace]15\[LetterSpace]k1*species\[LetterSpace]7[t], reaction\[LetterSpace]16 -> function\[LetterSpace]11[reaction\[LetterSpace]16\[LetterSpace]khrs, species\[LetterSpace]19, species\[LetterSpace]7[t]], reaction\[LetterSpace]17 -> reaction\[LetterSpace]17\[LetterSpace]k1*species\[LetterSpace]18[t], reaction\[LetterSpace]18 -> reaction\[LetterSpace]18\[LetterSpace]k1*species\[LetterSpace]17[t], reaction\[LetterSpace]19 -> function\[LetterSpace]12[reaction\[LetterSpace]19\[LetterSpace]k15, species\[LetterSpace]17[t], species\[LetterSpace]22[t]], reaction\[LetterSpace]2 -> function\[LetterSpace]1[reaction\[LetterSpace]2\[LetterSpace]ICSmax, species\[LetterSpace]2[t], reaction\[LetterSpace]2\[LetterSpace]ICt, reaction\[LetterSpace]2\[LetterSpace]IS], reaction\[LetterSpace]20 -> reaction\[LetterSpace]20\[LetterSpace]k1*species\[LetterSpace]21[t], reaction\[LetterSpace]21 -> function\[LetterSpace]13[reaction\[LetterSpace]21\[LetterSpace]k17, species\[LetterSpace]21[t], species\[LetterSpace]24[t]], reaction\[LetterSpace]22 -> function\[LetterSpace]14[reaction\[LetterSpace]22\[LetterSpace]k18, species\[LetterSpace]23[t], species\[LetterSpace]18[t]], reaction\[LetterSpace]23 -> reaction\[LetterSpace]23\[LetterSpace]k1*species\[LetterSpace]23[t], reaction\[LetterSpace]24 -> function\[LetterSpace]15[reaction\[LetterSpace]24\[LetterSpace]k20, species\[LetterSpace]25[t], species\[LetterSpace]23[t]], reaction\[LetterSpace]25 -> reaction\[LetterSpace]25\[LetterSpace]k1*species\[LetterSpace]23[t], reaction\[LetterSpace]26 -> function\[LetterSpace]16[reaction\[LetterSpace]26\[LetterSpace]kprs, species\[LetterSpace]26, species\[LetterSpace]11[t]], reaction\[LetterSpace]27 -> reaction\[LetterSpace]27\[LetterSpace]k1*species\[LetterSpace]25[t], reaction\[LetterSpace]28 -> function\[LetterSpace]17[reaction\[LetterSpace]28\[LetterSpace]k23, species\[LetterSpace]14[t], species\[LetterSpace]11[t]], reaction\[LetterSpace]29 -> function\[LetterSpace]18[reaction\[LetterSpace]29\[LetterSpace]k24, species\[LetterSpace]15[t], species\[LetterSpace]28[t]], reaction\[LetterSpace]3 -> reaction\[LetterSpace]3\[LetterSpace]k1*species\[LetterSpace]4[t], reaction\[LetterSpace]30 -> reaction\[LetterSpace]30\[LetterSpace]k1*species\[LetterSpace]11[t], reaction\[LetterSpace]31 -> function\[LetterSpace]19[reaction\[LetterSpace]31\[LetterSpace]k26, species\[LetterSpace]11[t], species\[LetterSpace]10[t], species\[LetterSpace]31[t]], reaction\[LetterSpace]32 -> function\[LetterSpace]20[reaction\[LetterSpace]32\[LetterSpace]PCSmax, species\[LetterSpace]11[t], reaction\[LetterSpace]32\[LetterSpace]PPCt, reaction\[LetterSpace]32\[LetterSpace]PCSS], reaction\[LetterSpace]33 -> function\[LetterSpace]21[reaction\[LetterSpace]33\[LetterSpace]k27, species\[LetterSpace]30[t], species\[LetterSpace]33[t]], reaction\[LetterSpace]34 -> function\[LetterSpace]22[reaction\[LetterSpace]34\[LetterSpace]k28, species\[LetterSpace]30[t], species\[LetterSpace]33[t]], reaction\[LetterSpace]35 -> function\[LetterSpace]23[reaction\[LetterSpace]35\[LetterSpace]k29, species\[LetterSpace]30[t], species\[LetterSpace]34[t]], reaction\[LetterSpace]4 -> reaction\[LetterSpace]4\[LetterSpace]k1*species\[LetterSpace]5[t], reaction\[LetterSpace]5 -> reaction\[LetterSpace]5\[LetterSpace]k1*species\[LetterSpace]5[t], reaction\[LetterSpace]6 -> function\[LetterSpace]2[reaction\[LetterSpace]6\[LetterSpace]k5, species\[LetterSpace]7[t], species\[LetterSpace]4[t]], reaction\[LetterSpace]7 -> function\[LetterSpace]3[reaction\[LetterSpace]7\[LetterSpace]k6, species\[LetterSpace]2[t], species\[LetterSpace]5[t]], reaction\[LetterSpace]8 -> function\[LetterSpace]4[reaction\[LetterSpace]8\[LetterSpace]k7, species\[LetterSpace]2[t], species\[LetterSpace]5[t]], reaction\[LetterSpace]9 -> function\[LetterSpace]5[reaction\[LetterSpace]9\[LetterSpace]k8, species\[LetterSpace]11[t]]
};

parameters = {
multiplier -> 1.0, species\[LetterSpace]1 -> 304.0, species\[LetterSpace]12 -> 0.0, species\[LetterSpace]16 -> 0.0, species\[LetterSpace]19 -> 600.0, species\[LetterSpace]26 -> 575.16, species\[LetterSpace]3 -> 0.0, species\[LetterSpace]32 -> 0.0, species\[LetterSpace]9 -> 0.0, reaction\[LetterSpace]10\[LetterSpace]BCRmax -> 2000.0, reaction\[LetterSpace]1\[LetterSpace]k1 -> 1.0, reaction\[LetterSpace]2\[LetterSpace]ICSmax -> 100.0, reaction\[LetterSpace]2\[LetterSpace]ICt -> 3120.0, reaction\[LetterSpace]2\[LetterSpace]IS -> 5.0, reaction\[LetterSpace]3\[LetterSpace]k1 -> 6.0, reaction\[LetterSpace]4\[LetterSpace]k1 -> 4.29, reaction\[LetterSpace]5\[LetterSpace]k1 -> 0.856, reaction\[LetterSpace]6\[LetterSpace]k5 -> 2.66, reaction\[LetterSpace]7\[LetterSpace]k6 -> 0.0005286, reaction\[LetterSpace]8\[LetterSpace]k7 -> 0.0005286, reaction\[LetterSpace]9\[LetterSpace]k8 -> 0.0005, reaction\[LetterSpace]10\[LetterSpace]BCRt -> 55326.0, reaction\[LetterSpace]10\[LetterSpace]BS -> 5.0, reaction\[LetterSpace]11\[LetterSpace]HCSmax -> 500.0, reaction\[LetterSpace]11\[LetterSpace]HCSt -> 93925.0, reaction\[LetterSpace]11\[LetterSpace]HS -> 5.0, reaction\[LetterSpace]12\[LetterSpace]k9 -> 1.0, reaction\[LetterSpace]13\[LetterSpace]k10 -> 5.998, reaction\[LetterSpace]14\[LetterSpace]k11 -> 0.005, reaction\[LetterSpace]15\[LetterSpace]k1 -> 0.016, reaction\[LetterSpace]16\[LetterSpace]khrs -> 100.0, reaction\[LetterSpace]17\[LetterSpace]k1 -> 0.01, reaction\[LetterSpace]18\[LetterSpace]k1 -> 0.0496, reaction\[LetterSpace]19\[LetterSpace]k15 -> 0.43, reaction\[LetterSpace]20\[LetterSpace]k1 -> 0.054, reaction\[LetterSpace]21\[LetterSpace]k17 -> 0.38, reaction\[LetterSpace]22\[LetterSpace]k18 -> 0.068, reaction\[LetterSpace]23\[LetterSpace]k1 -> 0.005, reaction\[LetterSpace]24\[LetterSpace]k20 -> 0.00675, reaction\[LetterSpace]25\[LetterSpace]k1 -> 5*^-06, reaction\[LetterSpace]26\[LetterSpace]kprs -> 100.0, reaction\[LetterSpace]27\[LetterSpace]k1 -> 0.01, reaction\[LetterSpace]28\[LetterSpace]k23 -> 0.017386, reaction\[LetterSpace]29\[LetterSpace]k24 -> 0.1068, reaction\[LetterSpace]30\[LetterSpace]k1 -> 0.0005, reaction\[LetterSpace]31\[LetterSpace]k26 -> 1.5*^-05, reaction\[LetterSpace]32\[LetterSpace]PCSmax -> 500.0, reaction\[LetterSpace]32\[LetterSpace]PPCt -> 80342.0, reaction\[LetterSpace]32\[LetterSpace]PCSS -> 5.0, reaction\[LetterSpace]33\[LetterSpace]k27 -> 0.01, reaction\[LetterSpace]34\[LetterSpace]k28 -> 0.001, reaction\[LetterSpace]35\[LetterSpace]k29 -> 0.05, Excreted -> 1.0, HepaticTissue -> 1.0, Intake -> 1.0, Intestine -> 1.0, PeripheralTissue -> 1.0, Plasma -> 1.0
};

assignments = {
function\[LetterSpace]4[k7_,IC_,IBS_] -> IBS*IC*k7, function\[LetterSpace]5[k8_,PFC_] -> k8*PFC, function\[LetterSpace]6[BCRmax_,BCRt_,HFC_,BS_] -> BCRmax/(1 + (BCRt/HFC)^BS), function\[LetterSpace]10[k11_,PFC_] -> k11*PFC, function\[LetterSpace]8[k9_,ACAT_,HFC_] -> ACAT*HFC*k9, function\[LetterSpace]16[kprs_,PLDLRsS_,PFC_] -> (kprs*PLDLRsS)/PFC, function\[LetterSpace]15[k20_,PLDLRs_,LDLC_] -> k20*LDLC*PLDLRs, function\[LetterSpace]14[k18_,LDLC_,HLDLRs_] -> HLDLRs*k18*LDLC, function\[LetterSpace]13[k17_,IDLC_,HSL_] -> HSL*IDLC*k17, function\[LetterSpace]12[k15_,VLDLC_,LPL_] -> k15*LPL*VLDLC, function\[LetterSpace]11[khrs_,HLDLRsS_,HFC_] -> (HLDLRsS*khrs)/HFC, function\[LetterSpace]1[ICSmax_,IC_,ICt_,IS_] -> ICSmax/(1 + (IC/ICt)^IS), function\[LetterSpace]2[k5_,HFC_,HBS_] -> (HFC*k5)/HBS, function\[LetterSpace]3[k6_,IC_,IBS_] -> IBS*IC*k6, function\[LetterSpace]7[HCSmax_,HFC_,HCSt_,HS_] -> HCSmax/(1 + (HFC/HCSt)^HS), function\[LetterSpace]9[k10_,CEH_,HCE_] -> CEH*HCE*k10, function\[LetterSpace]23[k29_,HDLC_,SRB1_] -> HDLC*k29*SRB1, function\[LetterSpace]22[k28_,HDLC_,CETP_] -> CETP*HDLC*k28, function\[LetterSpace]21[k27_,HDLC_,CETP_] -> CETP*HDLC*k27, function\[LetterSpace]20[PCSmax_,PFC_,PPCt_,PCSS_] -> PCSmax/(1 + (PFC/PPCt)^PCSS), function\[LetterSpace]19[k26_,PFC_,NHDL_,LCAT_] -> k26*LCAT*NHDL*PFC, function\[LetterSpace]18[k24_,CEH_,PCE_] -> CEH*k24*PCE, function\[LetterSpace]17[k23_,ACAT_,PFC_] -> ACAT*k23*PFC
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

units = {
{"time" -> "", "metabolite" -> "<substance_units>/volume", "extent" -> "<substance_units>/volume"}
};

(* Time evolution *)
odes = {
species\[LetterSpace]10'[t] == 1.0*reaction\[LetterSpace]9 +1.0*reaction\[LetterSpace]14 -1.0*reaction\[LetterSpace]31, species\[LetterSpace]11'[t] == 1.0*reaction\[LetterSpace]24 +1.0*reaction\[LetterSpace]25 +1.0*reaction\[LetterSpace]29 +1.0*reaction\[LetterSpace]32 -1.0*reaction\[LetterSpace]28 -1.0*reaction\[LetterSpace]30 -1.0*reaction\[LetterSpace]31, species\[LetterSpace]13'[t] == 1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]13, species\[LetterSpace]14'[t] == 0.0 , species\[LetterSpace]15'[t] == 0.0 , species\[LetterSpace]17'[t] == 1.0*reaction\[LetterSpace]15 +1.0*reaction\[LetterSpace]33 -1.0*reaction\[LetterSpace]18 -1.0*reaction\[LetterSpace]19, species\[LetterSpace]18'[t] == 1.0*reaction\[LetterSpace]16 -1.0*reaction\[LetterSpace]17, species\[LetterSpace]2'[t] == 1.0*reaction\[LetterSpace]1 +1.0*reaction\[LetterSpace]2 +1.0*reaction\[LetterSpace]10 -1.0*reaction\[LetterSpace]7 -1.0*reaction\[LetterSpace]8, species\[LetterSpace]20'[t] == 1.0*reaction\[LetterSpace]17 , species\[LetterSpace]21'[t] == 1.0*reaction\[LetterSpace]19 -1.0*reaction\[LetterSpace]20 -1.0*reaction\[LetterSpace]21, species\[LetterSpace]22'[t] == 0.0 , species\[LetterSpace]23'[t] == 1.0*reaction\[LetterSpace]21 +1.0*reaction\[LetterSpace]34 -1.0*reaction\[LetterSpace]22 -1.0*reaction\[LetterSpace]23 -1.0*reaction\[LetterSpace]24 -1.0*reaction\[LetterSpace]25, species\[LetterSpace]24'[t] == 0.0 , species\[LetterSpace]25'[t] == 1.0*reaction\[LetterSpace]26 -1.0*reaction\[LetterSpace]27, species\[LetterSpace]27'[t] == 1.0*reaction\[LetterSpace]27 , species\[LetterSpace]28'[t] == 1.0*reaction\[LetterSpace]28 -1.0*reaction\[LetterSpace]29, species\[LetterSpace]29'[t] == 1.0*reaction\[LetterSpace]30 , species\[LetterSpace]30'[t] == 1.0*reaction\[LetterSpace]31 -1.0*reaction\[LetterSpace]33 -1.0*reaction\[LetterSpace]34 -1.0*reaction\[LetterSpace]35, species\[LetterSpace]31'[t] == 0.0 , species\[LetterSpace]33'[t] == 0.0 , species\[LetterSpace]34'[t] == 0.0 , species\[LetterSpace]4'[t] == 1.0*reaction\[LetterSpace]4 +1.0*reaction\[LetterSpace]6 -1.0*reaction\[LetterSpace]3, species\[LetterSpace]5'[t] == 1.0*reaction\[LetterSpace]3 -1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]5, species\[LetterSpace]6'[t] == 1.0*reaction\[LetterSpace]5 , species\[LetterSpace]7'[t] == 1.0*reaction\[LetterSpace]7 +1.0*reaction\[LetterSpace]11 +1.0*reaction\[LetterSpace]13 +1.0*reaction\[LetterSpace]18 +1.0*reaction\[LetterSpace]20 +1.0*reaction\[LetterSpace]22 +1.0*reaction\[LetterSpace]23 +1.0*reaction\[LetterSpace]35 -1.0*reaction\[LetterSpace]6 -1.0*reaction\[LetterSpace]10 -1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]15, species\[LetterSpace]8'[t] == 1.0*reaction\[LetterSpace]8 
};
timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}];

(* Steady-state solution initialized with result of time evolution *)
findRootEquations = odes /.D[_[t],t]->0;
findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2];
steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100]
fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations

(* Plot the time evolution of the variables *)
plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]