(* 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 = {
xBH[t], xBO[t], xBT[t], xBU[t], xBV[t], xCA[t], xCH[t], xDP[t], xEL[t], xER[t], xFP[t], xPC[t]
};

initialValues = {
xBH[0] == 0.01, xBO[0] == 0.01, xBT[0] == 0.01, xBU[0] == 0.01, xBV[0] == 0.01, xCA[0] == 0.01, xCH[0] == 0.01, xDP[0] == 0.01, xEL[0] == 0.01, xER[0] == 0.01, xFP[0] == 0.01, xPC[0] == 0.01
};

rates = {
v1, v10, v11, v12, v2, v3, v4, v5, v6, v7, v8, v9
};

rateEquations = {
v1 -> xBH[t]*(uBH + aBHBH*xBH[t] + aBHBO*xBO[t] + aBHBT*xBT[t] + aBHBU*xBU[t] + aBHBV*xBV[t] + aBHCA*xCA[t] + aBHCH*xCH[t] + aBHDP*xDP[t] + aBHEL*xEL[t] + aBHER*xER[t] + aBHFP*xFP[t] + aBHPC*xPC[t]), v10 -> xCH[t]*(uCH + aCHBH*xBH[t] + aCHBO*xBO[t] + aCHBT*xBT[t] + aCHBU*xBU[t] + aCHBV*xBV[t] + aCHCA*xCA[t] + aCHCH*xCH[t] + aCHDP*xDP[t] + aCHEL*xEL[t] + aCHER*xER[t] + aCHFP*xFP[t] + aCHPC*xPC[t]), v11 -> xDP[t]*(uDP + aDPBH*xBH[t] + aDPBO*xBO[t] + aDPBT*xBT[t] + aDPBU*xBU[t] + aDPBV*xBV[t] + aDPCA*xCA[t] + aDPCH*xCH[t] + aDPDP*xDP[t] + aDPEL*xEL[t] + aDPER*xER[t] + aDPFP*xFP[t] + aDPPC*xPC[t]), v12 -> xER[t]*(uER + aERBH*xBH[t] + aERBO*xBO[t] + aERBT*xBT[t] + aERBU*xBU[t] + aERBV*xBV[t] + aERCA*xCA[t] + aERCH*xCH[t] + aERDP*xDP[t] + aEREL*xEL[t] + aERER*xER[t] + aERFP*xFP[t] + aERPC*xPC[t]), v2 -> xCA[t]*(uCA + aCABH*xBH[t] + aCABO*xBO[t] + aCABT*xBT[t] + aCABU*xBU[t] + aCABV*xBV[t] + aCACA*xCA[t] + aCACH*xCH[t] + aCADP*xDP[t] + aCAEL*xEL[t] + aCAER*xER[t] + aCAFP*xFP[t] + aCAPC*xPC[t]), v3 -> xBU[t]*(uBU + aBUBH*xBH[t] + aBUBO*xBO[t] + aBUBT*xBT[t] + aBUBU*xBU[t] + aBUBV*xBV[t] + aBUCA*xCA[t] + aBUCH*xCH[t] + aBUDP*xDP[t] + aBUEL*xEL[t] + aBUER*xER[t] + aBUFP*xFP[t] + aBUPC*xPC[t]), v4 -> xPC[t]*(uPC + aPCBH*xBH[t] + aPCBO*xBO[t] + aPCBT*xBT[t] + aPCBU*xBU[t] + aPCBV*xBV[t] + aPCCA*xCA[t] + aPCCH*xCH[t] + aPCDP*xDP[t] + aPCEL*xEL[t] + aPCER*xER[t] + aPCFP*xFP[t] + aPCPC*xPC[t]), v5 -> xBO[t]*(uBO + aBOBH*xBH[t] + aBOBO*xBO[t] + aBOBT*xBT[t] + aBOBU*xBU[t] + aBOBV*xBV[t] + aBOCA*xCA[t] + aBOCH*xCH[t] + aBODP*xDP[t] + aBOEL*xEL[t] + aBOER*xER[t] + aBOFP*xFP[t] + aBOPC*xPC[t]), v6 -> xBV[t]*(uBV + aBVBH*xBH[t] + aBVBO*xBO[t] + aBVBT*xBT[t] + aBVBU*xBU[t] + aBVBV*xBV[t] + aBVCA*xCA[t] + aBVCH*xCH[t] + aBVDP*xDP[t] + aBVEL*xEL[t] + aBVER*xER[t] + aBVFP*xFP[t] + aBVPC*xPC[t]), v7 -> xBT[t]*(uBT + aBTBH*xBH[t] + aBTBO*xBO[t] + aBTBT*xBT[t] + aBTBU*xBU[t] + aBTBV*xBV[t] + aBTCA*xCA[t] + aBTCH*xCH[t] + aBTDP*xDP[t] + aBTEL*xEL[t] + aBTER*xER[t] + aBTFP*xFP[t] + aBTPC*xPC[t]), v8 -> xEL[t]*(uEL + aELBH*xBH[t] + aELBO*xBO[t] + aELBT*xBT[t] + aELBU*xBU[t] + aELBV*xBV[t] + aELCA*xCA[t] + aELCH*xCH[t] + aELDP*xDP[t] + aELEL*xEL[t] + aELER*xER[t] + aELFP*xFP[t] + aELPC*xPC[t]), v9 -> xFP[t]*(uFP + aFPBH*xBH[t] + aFPBO*xBO[t] + aFPBT*xBT[t] + aFPBU*xBU[t] + aFPBV*xBV[t] + aFPCA*xCA[t] + aFPCH*xCH[t] + aFPDP*xDP[t] + aFPEL*xEL[t] + aFPER*xER[t] + aFPFP*xFP[t] + aFPPC*xPC[t])
};

parameters = {
aBHBH -> -0.9118, aBHBO -> -0.2145, aBHBT -> -0.2718, aBHBU -> -0.2275, aBHBV -> -0.1294, aBHCA -> -0.3058, aBHCH -> -0.3478, aBHDP -> -0.9002, aBHEL -> 0.1764, aBHER -> -0.5455, aBHFP -> -0.2307, aBHPC -> -0.5286, aBOBH -> 0.0, aBOBO -> -0.7339, aBOBT -> -0.6235, aBOBU -> -0.6317, aBOBV -> -0.515, aBOCA -> 0.0, aBOCH -> -0.5069, aBODP -> 0.0, aBOEL -> 1.7566, aBOER -> 0.0, aBOFP -> -0.2087, aBOPC -> 0.0, aBTBH -> 0.0, aBTBO -> -0.8193, aBTBT -> -0.9067, aBTBU -> -0.7538, aBTBV -> -0.7552, aBTCA -> 0.0, aBTCH -> -0.0864, aBTDP -> 0.0, aBTEL -> 2.2719, aBTER -> -0.7375, aBTFP -> -0.7033, aBTPC -> 0.0, aBUBH -> 0.0, aBUBO -> -0.9208, aBUBT -> -0.8155, aBUBU -> -0.8804, aBUBV -> -0.5498, aBUCA -> 0.0, aBUCH -> 0.0656, aBUDP -> 0.0, aBUEL -> 3.3782, aBUER -> 0.0, aBUFP -> -0.7822, aBUPC -> 0.0, aBVBH -> 0.1366, aBVBO -> -0.5556, aBVBT -> -0.6416, aBVBU -> -0.5837, aBVBV -> -0.6597, aBVCA -> -0.6566, aBVCH -> -0.0468, aBVDP -> -0.1062, aBVEL -> 1.3026, aBVER -> 0.0, aBVFP -> -0.6379, aBVPC -> 0.0, aCABH -> 0.4526, aCABO -> -0.2776, aCABT -> -0.2736, aCABU -> -0.2614, aCABV -> -0.1679, aCACA -> -0.829, aCACH -> 0.3118, aCADP -> 0.0, aCAEL -> -0.4475, aCAER -> 0.0, aCAFP -> -1.121, aCAPC -> -0.6709, aCHBH -> 0.0, aCHBO -> -0.4645, aCHBT -> -0.6319, aCHBU -> -0.1511, aCHBV -> 0.039, aCHCA -> -0.2413, aCHCH -> -1.4543, aCHDP -> -2.157, aCHEL -> 0.0239, aCHER -> -0.4394, aCHFP -> -0.5069, aCHPC -> -0.771, aDPBH -> 0.0, aDPBO -> -0.2028, aDPBT -> -0.1999, aDPBU -> -0.1763, aDPBV -> -0.026, aDPCA -> 0.042, aDPCH -> -0.151, aDPDP -> -1.2535, aDPEL -> 0.1756, aDPER -> 0.0, aDPFP -> 0.0, aDPPC -> -0.4333, aELBH -> 0.6924, aELBO -> -0.1141, aELBT -> -0.1688, aELBU -> -0.1241, aELBV -> -0.0493, aELCA -> -1.098, aELCH -> 0.0, aELDP -> -0.4084, aELEL -> -2.4418, aELER -> 0.0, aELFP -> -0.1531, aELPC -> -1.0774, aERBH -> 1.3433, aERBO -> -0.0203, aERBT -> 0.0, aERBU -> -0.0613, aERBV -> -0.0296, aERCA -> 0.0, aERCH -> 1.0831, aERDP -> 0.0, aEREL -> -0.1385, aERER -> -1.2705, aERFP -> -0.1676, aERPC -> 0.0, aFPBH -> 0.9613, aFPBO -> -0.0993, aFPBT -> -0.0685, aFPBU -> 0.2313, aFPBV -> 0.7585, aFPCA -> 0.0, aFPCH -> 0.4481, aFPDP -> 1.0147, aFPEL -> -0.768, aFPER -> 0.0, aFPFP -> -1.0382, aFPPC -> -0.402, aPCBH -> 0.0, aPCBO -> -0.2646, aPCBT -> -0.3033, aPCBU -> -0.3237, aPCBV -> -0.2017, aPCCA -> -0.5596, aPCCH -> 0.265, aPCDP -> -0.9771, aPCEL -> -0.9041, aPCER -> -0.8171, aPCFP -> -0.4053, aPCPC -> -0.6217, uBH -> 0.2453, uBO -> 0.4777, uBT -> 0.5979, uBU -> 0.5841, uBV -> 0.4573, uCA -> 0.2464, uCH -> 0.5025, uDP -> 0.2321, uEL -> 0.4021, uER -> 0.1558, uFP -> 0.2192, uPC -> 0.2375, x1 -> 1.0, x10 -> 1.0, x11 -> 1.0, x12 -> 1.0, x2 -> 1.0, x3 -> 1.0, x4 -> 1.0, x5 -> 1.0, x6 -> 1.0, x7 -> 1.0, x8 -> 1.0, x9 -> 1.0, default -> 1.0
};

assignments = {
fBO -> xBO[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fCH -> xCH[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fBT -> xBT[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fEL -> xEL[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fFP -> xFP[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fBH -> xBH[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fBU -> xBU[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fBV -> xBV[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fPC -> xPC[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fER -> xER[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fDP -> xDP[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t]), fCA -> xCA[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

units = {
{"time" -> "", "metabolite" -> "", "extent" -> ""}
};

(* Time evolution *)
odes = {
xBH'[t] == 1.0*v1 , xBO'[t] == 1.0*v5 , xBT'[t] == 1.0*v7 , xBU'[t] == 1.0*v3 , xBV'[t] == 1.0*v6 , xCA'[t] == 1.0*v2 , xCH'[t] == 1.0*v10 , xDP'[t] == 1.0*v11 , xEL'[t] == 1.0*v8 , xER'[t] == 1.0*v12 , xFP'[t] == 1.0*v9 , xPC'[t] == 1.0*v4 
};
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]}]