(* 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.9241, aBHBO -> -0.2199, aBHBT -> -0.2641, aBHBU -> -0.2284, aBHBV -> -0.131, aBHCA -> -0.312, aBHCH -> -0.3284, aBHDP -> -0.9182, aBHEL -> 0.1821, aBHER -> -0.5488, aBHFP -> -0.2367, aBHPC -> -0.5294, aBOBH -> 0.0, aBOBO -> -0.7134, aBOBT -> -0.5769, aBOBU -> -0.5064, aBOBV -> -0.5089, aBOCA -> 0.0, aBOCH -> -0.4668, aBODP -> 0.0, aBOEL -> 1.8096, aBOER -> -0.0391, aBOFP -> -0.1758, aBOPC -> -0.0001, aBTBH -> 0.0582, aBTBO -> -0.8237, aBTBT -> -0.9597, aBTBU -> -0.972, aBTBV -> -0.6819, aBTCA -> 0.0, aBTCH -> -0.0727, aBTDP -> 0.0, aBTEL -> 1.9668, aBTER -> 0.0, aBTFP -> -0.4498, aBTPC -> 0.0, aBUBH -> 0.0, aBUBO -> -0.7487, aBUBT -> -0.9377, aBUBU -> -0.9059, aBUBV -> -0.5581, aBUCA -> 0.0, aBUCH -> 0.0564, aBUDP -> 0.0, aBUEL -> 3.0816, aBUER -> 0.0, aBUFP -> -0.8309, aBUPC -> 0.0, aBVBH -> 0.1337, aBVBO -> -0.5768, aBVBT -> -0.6133, aBVBU -> -0.5773, aBVBV -> -0.6578, aBVCA -> -0.654, aBVCH -> -0.0464, aBVDP -> -0.108, aBVEL -> 1.2848, aBVER -> 0.0, aBVFP -> -0.6349, aBVPC -> 0.0, aCABH -> 0.4301, aCABO -> -0.283, aCABT -> -0.2737, aCABU -> -0.2612, aCABV -> -0.1701, aCACA -> -0.818, aCACH -> 0.3033, aCADP -> 0.0, aCAEL -> -0.4766, aCAER -> 0.0, aCAFP -> -1.1212, aCAPC -> -0.685, aCHBH -> -0.3846, aCHBO -> -0.4679, aCHBT -> -0.5906, aCHBU -> -0.2163, aCHBV -> 0.0, aCHCA -> -0.3186, aCHCH -> -1.242, aCHDP -> -2.3006, aCHEL -> 0.0, aCHER -> -0.5077, aCHFP -> -0.6943, aCHPC -> -0.7845, aDPBH -> 0.0, aDPBO -> -0.187, aDPBT -> -0.1915, aDPBU -> -0.1628, aDPBV -> -0.0172, aDPCA -> 0.0532, aDPCH -> -0.1472, aDPDP -> -1.3228, aDPEL -> 0.1989, aDPER -> 0.0, aDPFP -> 0.0, aDPPC -> -0.4407, aELBH -> 0.8285, aELBO -> -0.0932, aELBT -> -0.099, aELBU -> -0.0776, aELBV -> -0.0229, aELCA -> -0.9925, aELCH -> 0.0, aELDP -> -0.0003, aELEL -> -2.6202, aELER -> 0.0, aELFP -> 0.0, aELPC -> -1.091, aERBH -> 1.4712, aERBO -> 0.0, aERBT -> -0.0315, aERBU -> -0.0285, aERBV -> -0.0178, aERCA -> 0.0, aERCH -> 1.1905, aERDP -> 0.0, aEREL -> 0.0, aERER -> -1.3219, aERFP -> -0.0319, aERPC -> 0.0, aFPBH -> 0.9091, aFPBO -> -0.114, aFPBT -> -0.0742, aFPBU -> 0.1988, aFPBV -> 0.7253, aFPCA -> 0.0, aFPCH -> 0.4306, aFPDP -> 0.9762, aFPEL -> -0.8005, aFPER -> 0.0, aFPFP -> -0.9946, aFPPC -> -0.423, aPCBH -> 0.0, aPCBO -> -0.2718, aPCBT -> -0.3041, aPCBU -> -0.3241, aPCBV -> -0.2039, aPCCA -> -0.5827, aPCCH -> 0.2406, aPCDP -> -1.0009, aPCEL -> -0.9852, aPCER -> -0.9156, aPCFP -> -0.4513, aPCPC -> -0.6154, uBH -> 0.2472, uBO -> 0.4624, uBT -> 0.6259, uBU -> 0.5991, uBV -> 0.4562, uCA -> 0.2464, uCH -> 0.4682, uDP -> 0.2361, uEL -> 0.4052, uER -> 0.1507, uFP -> 0.2164, uPC -> 0.2381, 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 = {
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]), 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]), 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]), 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]}]