(* 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.7905, aBHBO -> 0.0, aBHBT -> 0.0, aBHBU -> 0.0, aBHBV -> 0.0, aBHCA -> 0.0, aBHCH -> 0.0, aBHDP -> 0.0, aBHEL -> 0.0, aBHER -> 0.0, aBHFP -> 0.0, aBHPC -> 0.0, aBOBH -> 0.0, aBOBO -> -1.082, aBOBT -> 0.0, aBOBU -> 0.0, aBOBV -> 0.0, aBOCA -> 0.0, aBOCH -> 0.0, aBODP -> 0.0, aBOEL -> 0.0, aBOER -> 0.0, aBOFP -> 0.0, aBOPC -> 0.0, aBTBH -> 0.0, aBTBO -> 0.0, aBTBT -> -1.1639, aBTBU -> 0.0, aBTBV -> 0.0, aBTCA -> 0.0, aBTCH -> 0.0, aBTDP -> 0.0, aBTEL -> 0.0, aBTER -> 0.0, aBTFP -> 0.0, aBTPC -> 0.0, aBUBH -> 0.0, aBUBO -> 0.0, aBUBT -> 0.0, aBUBU -> -1.1205, aBUBV -> 0.0, aBUCA -> 0.0, aBUCH -> 0.0, aBUDP -> 0.0, aBUEL -> 0.0, aBUER -> 0.0, aBUFP -> 0.0, aBUPC -> 0.0, aBVBH -> 0.0, aBVBO -> 0.0, aBVBT -> 0.0, aBVBU -> 0.0, aBVBV -> -0.7535, aBVCA -> 0.0, aBVCH -> 0.0, aBVDP -> 0.0, aBVEL -> 0.0, aBVER -> 0.0, aBVFP -> 0.0, aBVPC -> 0.0, aCABH -> 0.0, aCABO -> 0.0, aCABT -> 0.0, aCABU -> 0.0, aCABV -> 0.0, aCACA -> -0.7839, aCACH -> 0.0, aCADP -> 0.0, aCAEL -> 0.0, aCAER -> 0.0, aCAFP -> 0.0, aCAPC -> 0.0, aCHBH -> 0.0, aCHBO -> 0.0, aCHBT -> 0.0, aCHBU -> 0.0, aCHBV -> 0.0, aCHCA -> 0.0, aCHCH -> -1.4921, aCHDP -> 0.0, aCHEL -> 0.0, aCHER -> 0.0, aCHFP -> 0.0, aCHPC -> 0.0, aDPBH -> 0.0, aDPBO -> 0.0, aDPBT -> 0.0, aDPBU -> 0.0, aDPBV -> 0.0, aDPCA -> 0.0, aDPCH -> 0.0, aDPDP -> -0.6553, aDPEL -> 0.0, aDPER -> 0.0, aDPFP -> 0.0, aDPPC -> 0.0, aELBH -> 0.0, aELBO -> 0.0, aELBT -> 0.0, aELBU -> 0.0, aELBV -> 0.0, aELCA -> 0.0, aELCH -> 0.0, aELDP -> 0.0, aELEL -> -0.6457, aELER -> 0.0, aELFP -> 0.0, aELPC -> 0.0, aERBH -> 0.0, aERBO -> 0.0, aERBT -> 0.0, aERBU -> 0.0, aERBV -> 0.0, aERCA -> 0.0, aERCH -> 0.0, aERDP -> 0.0, aEREL -> 0.0, aERER -> -0.5581, aERFP -> 0.0, aERPC -> 0.0, aFPBH -> 0.0, aFPBO -> 0.0, aFPBT -> 0.0, aFPBU -> 0.0, aFPBV -> 0.0, aFPCA -> 0.0, aFPCH -> 0.0, aFPDP -> 0.0, aFPEL -> 0.0, aFPER -> 0.0, aFPFP -> -0.5218, aFPPC -> 0.0, aPCBH -> 0.0, aPCBO -> 0.0, aPCBT -> 0.0, aPCBU -> 0.0, aPCBV -> 0.0, aPCCA -> 0.0, aPCCH -> 0.0, aPCDP -> 0.0, aPCEL -> 0.0, aPCER -> 0.0, aPCFP -> 0.0, aPCPC -> -0.5146, uBH -> 0.2972, uBO -> 0.7332, uBT -> 0.8777, uBU -> 0.7689, uBV -> 0.5834, uCA -> 0.3454, uCH -> 0.7103, uDP -> 0.2082, uEL -> 0.1715, uER -> 0.0904, uFP -> 0.1619, uPC -> 0.1262, 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 = {
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]), 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]), 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]), 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]), 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]}]