(* 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.802, aBHBO -> -0.1407, aBHBT -> -0.1829, aBHBU -> -0.1492, aBHBV -> -0.0542, aBHCA -> -0.1616, aBHCH -> -0.0674, aBHDP -> -0.5281, aBHEL -> 1.1317, aBHER -> 0.0, aBHFP -> -0.0714, aBHPC -> -0.0016, aBOBH -> 0.0002, aBOBO -> -0.8101, aBOBT -> -0.6393, aBOBU -> -0.5679, aBOBV -> -0.5467, aBOCA -> 0.0, aBOCH -> -0.5602, aBODP -> 0.0, aBOEL -> 1.9238, aBOER -> 0.0, aBOFP -> -0.172, aBOPC -> -0.0022, aBTBH -> 0.0, aBTBO -> -0.894, aBTBT -> -1.0301, aBTBU -> -1.0671, aBTBV -> -0.7607, aBTCA -> 0.0, aBTCH -> 0.0, aBTDP -> 0.0, aBTEL -> 1.7918, aBTER -> 0.0, aBTFP -> -0.5519, aBTPC -> 0.0, aBUBH -> 0.0, aBUBO -> -0.7768, aBUBT -> -0.9798, aBUBU -> -0.9673, aBUBV -> -0.5769, aBUCA -> 0.0, aBUCH -> 0.0509, aBUDP -> 0.0, aBUEL -> 2.9487, aBUER -> 0.0, aBUFP -> -0.8667, aBUPC -> 0.0, aBVBH -> 0.2644, aBVBO -> -0.6707, aBVBT -> -0.7173, aBVBU -> -0.6842, aBVBV -> -0.7864, aBVCA -> -0.0638, aBVCH -> -0.0365, aBVDP -> -0.0024, aBVEL -> 1.5582, aBVER -> 0.0, aBVFP -> -0.7502, aBVPC -> 0.0, aCABH -> 0.3093, aCABO -> -0.3117, aCABT -> -0.2976, aCABU -> -0.2897, aCABV -> -0.1381, aCACA -> -0.8541, aCACH -> 0.2757, aCADP -> 0.0, aCAEL -> -0.3389, aCAER -> 0.0, aCAFP -> -0.9018, aCAPC -> -0.552, aCHBH -> 0.0, aCHBO -> 0.2353, aCHBT -> -0.7164, aCHBU -> -0.2685, aCHBV -> 0.0, aCHCA -> -0.3833, aCHCH -> -1.7395, aCHDP -> 0.0, aCHEL -> 0.0, aCHER -> -0.5151, aCHFP -> -0.8546, aCHPC -> -0.498, aDPBH -> 0.0, aDPBO -> -0.1243, aDPBT -> -0.132, aDPBU -> -0.1096, aDPBV -> -0.0218, aDPCA -> 0.0902, aDPCH -> -0.124, aDPDP -> -0.8923, aDPEL -> 0.2009, aDPER -> 0.0, aDPFP -> 0.0, aDPPC -> 0.0, aELBH -> 0.7903, aELBO -> -0.0508, aELBT -> -0.0554, aELBU -> -0.0282, aELBV -> 0.0, aELCA -> -0.0021, aELCH -> 0.0, aELDP -> -0.0002, aELEL -> -3.4228, aELER -> 0.0, aELFP -> 0.0, aELPC -> 0.0, aERBH -> 0.0003, aERBO -> 0.0008, aERBT -> 0.0, aERBU -> 0.0, aERBV -> 0.0009, aERCA -> 0.0, aERCH -> 0.8619, aERDP -> 0.0, aEREL -> 0.0, aERER -> -0.7175, aERFP -> 0.0, aERPC -> 0.0, aFPBH -> 0.8821, aFPBO -> -0.1144, aFPBT -> -0.0707, aFPBU -> 0.1825, aFPBV -> 0.6556, aFPCA -> -0.0001, aFPCH -> 0.3854, aFPDP -> 0.9454, aFPEL -> -0.7612, aFPER -> 0.0, aFPFP -> -0.9455, aFPPC -> -0.504, aPCBH -> 0.1861, aPCBO -> -0.1822, aPCBT -> -0.2322, aPCBU -> -0.2527, aPCBV -> -0.119, aPCCA -> 0.4026, aPCCH -> 0.3598, aPCDP -> 0.0, aPCEL -> -0.1694, aPCER -> 0.0, aPCFP -> 0.0, aPCPC -> -0.8721, uBH -> 0.1995, uBO -> 0.5128, uBT -> 0.6794, uBU -> 0.6371, uBV -> 0.5279, uCA -> 0.2639, uCH -> 0.6287, uDP -> 0.1919, uEL -> 0.4773, uER -> 0.1066, uFP -> 0.2144, uPC -> 0.2144, 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 = {
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]), 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]), 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]), 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]), 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]), 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]), 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]), 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]), 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]}]