(* 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 = {
ATIII[t], II[t], IIa[t], IIa\[LetterSpace]ATIII[t], IX[t], IXa[t], IXa\[LetterSpace]ATIII[t], IXa\[LetterSpace]VIIIa[t], IXa\[LetterSpace]VIIIa\[LetterSpace]X[t], TF[t], TFPI[t], TF\[LetterSpace]VII[t], TF\[LetterSpace]VIIa[t], TF\[LetterSpace]VIIa\[LetterSpace]ATIII[t], TF\[LetterSpace]VIIa\[LetterSpace]IX[t], TF\[LetterSpace]VIIa\[LetterSpace]X[t], TF\[LetterSpace]VIIa\[LetterSpace]Xa[t], TF\[LetterSpace]VIIa\[LetterSpace]Xa\[LetterSpace]TFPI[t], V[t], VII[t], VIII[t], VIIIa[t], VIIIa1\[LetterSpace]L[t], VIIIa2[t], VIIa[t], Va[t], X[t], Xa[t], Xa\[LetterSpace]ATIII[t], Xa\[LetterSpace]TFPI[t], Xa\[LetterSpace]Va[t], Xa\[LetterSpace]Va\[LetterSpace]II[t], mIIa[t], mIIa\[LetterSpace]ATIII[t]
};

initialValues = {
ATIII[0] == 3.4*^-06, II[0] == 1.4*^-06, IIa[0] == 1*^-09, IIa\[LetterSpace]ATIII[0] == 0.0, IX[0] == 9*^-08, IXa[0] == 2.1*^-10, IXa\[LetterSpace]ATIII[0] == 0.0, IXa\[LetterSpace]VIIIa[0] == 0.0, IXa\[LetterSpace]VIIIa\[LetterSpace]X[0] == 0.0, TF[0] == 0.0, TFPI[0] == 2.5*^-09, TF\[LetterSpace]VII[0] == 0.0, TF\[LetterSpace]VIIa[0] == 0.0, TF\[LetterSpace]VIIa\[LetterSpace]ATIII[0] == 0.0, TF\[LetterSpace]VIIa\[LetterSpace]IX[0] == 0.0, TF\[LetterSpace]VIIa\[LetterSpace]X[0] == 0.0, TF\[LetterSpace]VIIa\[LetterSpace]Xa[0] == 0.0, TF\[LetterSpace]VIIa\[LetterSpace]Xa\[LetterSpace]TFPI[0] == 0.0, V[0] == 2*^-08, VII[0] == 1*^-08, VIII[0] == 7*^-10, VIIIa[0] == 0.0, VIIIa1\[LetterSpace]L[0] == 0.0, VIIIa2[0] == 0.0, VIIa[0] == 1*^-10, Va[0] == 0.0, X[0] == 1.6*^-07, Xa[0] == 9.4*^-11, Xa\[LetterSpace]ATIII[0] == 0.0, Xa\[LetterSpace]TFPI[0] == 0.0, Xa\[LetterSpace]Va[0] == 0.0, Xa\[LetterSpace]Va\[LetterSpace]II[0] == 0.0, mIIa[0] == 0.0, mIIa\[LetterSpace]ATIII[0] == 0.0
};

rates = {
R1, R10, R11, R12, R12b, R13, R14, R15, R16, R17, R18, R18b, R19, R2, R20, R21, R22, R23, R24, R25, R26, R27, R28, R29, R3, R4, R5, R6, R6b, R7, R8, R8b, R9
};

rateEquations = {
R1 -> compartment\[LetterSpace]1*(-(k1*TF\[LetterSpace]VII[t]) + k2*TF[t]*VII[t]), R10 -> compartment\[LetterSpace]1*k17*IIa[t]*VIII[t], R11 -> compartment\[LetterSpace]1*(-(k18*IXa\[LetterSpace]VIIIa[t]) + k19*IXa[t]*VIIIa[t]), R12 -> compartment\[LetterSpace]1*(-(k20*IXa\[LetterSpace]VIIIa\[LetterSpace]X[t]) + k21*IXa\[LetterSpace]VIIIa[t]*X[t]), R12b -> compartment\[LetterSpace]1*k22*IXa\[LetterSpace]VIIIa\[LetterSpace]X[t], R13 -> compartment\[LetterSpace]1*(k24*VIIIa[t] - k23*VIIIa1\[LetterSpace]L[t]*VIIIa2[t]), R14 -> compartment\[LetterSpace]1*k25*IXa\[LetterSpace]VIIIa\[LetterSpace]X[t], R15 -> compartment\[LetterSpace]1*k25*IXa\[LetterSpace]VIIIa[t], R16 -> compartment\[LetterSpace]1*k26*IIa[t]*V[t], R17 -> compartment\[LetterSpace]1*(k28*Va[t]*Xa[t] - k27*Xa\[LetterSpace]Va[t]), R18 -> compartment\[LetterSpace]1*(k30*II[t]*Xa\[LetterSpace]Va[t] - k29*Xa\[LetterSpace]Va\[LetterSpace]II[t]), R18b -> compartment\[LetterSpace]1*k31*Xa\[LetterSpace]Va\[LetterSpace]II[t], R19 -> compartment\[LetterSpace]1*k32*mIIa[t]*Xa\[LetterSpace]Va[t], R2 -> compartment\[LetterSpace]1*(-(k3*TF\[LetterSpace]VIIa[t]) + k4*TF[t]*VIIa[t]), R20 -> compartment\[LetterSpace]1*(k34*TFPI[t]*Xa[t] - k33*Xa\[LetterSpace]TFPI[t]), R21 -> compartment\[LetterSpace]1*(k36*TFPI[t]*TF\[LetterSpace]VIIa\[LetterSpace]Xa[t] - k35*TF\[LetterSpace]VIIa\[LetterSpace]Xa\[LetterSpace]TFPI[t]), R22 -> compartment\[LetterSpace]1*k37*TF\[LetterSpace]VIIa[t]*Xa\[LetterSpace]TFPI[t], R23 -> compartment\[LetterSpace]1*k38*ATIII[t]*Xa[t], R24 -> compartment\[LetterSpace]1*k39*ATIII[t]*mIIa[t], R25 -> compartment\[LetterSpace]1*k40*ATIII[t]*IXa[t], R26 -> compartment\[LetterSpace]1*k41*ATIII[t]*IIa[t], R27 -> compartment\[LetterSpace]1*k42*ATIII[t]*TF\[LetterSpace]VIIa[t], R28 -> compartment\[LetterSpace]1*k43*IXa[t]*X[t], R29 -> compartment\[LetterSpace]1*k44*mIIa[t]*V[t], R3 -> compartment\[LetterSpace]1*k5*TF\[LetterSpace]VIIa[t]*VII[t], R4 -> compartment\[LetterSpace]1*k6*VII[t]*Xa[t], R5 -> compartment\[LetterSpace]1*k7*IIa[t]*VII[t], R6 -> compartment\[LetterSpace]1*(-(k8*TF\[LetterSpace]VIIa\[LetterSpace]X[t]) + k9*TF\[LetterSpace]VIIa[t]*X[t]), R6b -> compartment\[LetterSpace]1*k10*TF\[LetterSpace]VIIa\[LetterSpace]X[t], R7 -> compartment\[LetterSpace]1*(-(k11*TF\[LetterSpace]VIIa\[LetterSpace]Xa[t]) + k12*TF\[LetterSpace]VIIa[t]*Xa[t]), R8 -> compartment\[LetterSpace]1*(k14*IX[t]*TF\[LetterSpace]VIIa[t] - k13*TF\[LetterSpace]VIIa\[LetterSpace]IX[t]), R8b -> compartment\[LetterSpace]1*k15*TF\[LetterSpace]VIIa\[LetterSpace]IX[t], R9 -> compartment\[LetterSpace]1*k16*II[t]*Xa[t]
};

parameters = {
k1 -> 0.0031, k10 -> 6.0, k11 -> 19.0, k12 -> 22000000.0, k13 -> 2.4, k14 -> 10000000.0, k15 -> 1.8, k16 -> 7500.0, k17 -> 20000000.0, k18 -> 0.005, k19 -> 10000000.0, k2 -> 3200000.0, k20 -> 0.001, k21 -> 100000000.0, k22 -> 8.2, k23 -> 22000.0, k24 -> 0.006, k25 -> 0.001, k26 -> 20000000.0, k27 -> 0.2, k28 -> 400000000.0, k29 -> 103.0, k3 -> 0.0031, k30 -> 100000000.0, k31 -> 63.5, k32 -> 15000000.0, k33 -> 0.00036, k34 -> 900000.0, k35 -> 0.00011, k36 -> 320000000.0, k37 -> 50000000.0, k38 -> 1500.0, k39 -> 7100.0, k4 -> 23000000.0, k40 -> 490.0, k41 -> 7100.0, k42 -> 230.0, k43 -> 5700.0, k44 -> 3000000.0, k5 -> 440000.0, k6 -> 13000000.0, k7 -> 23000.0, k8 -> 1.05, k9 -> 25000000.0, compartment\[LetterSpace]1 -> 1.0
};

assignments = {
IIa\[LetterSpace]plus\[LetterSpace]1\[LetterSpace]2mIIa -> IIa[t] + 1.2*mIIa[t]
};

events = {

};

speciesAnnotations = {
ATIII[t]->"http://identifiers.org/uniprot/P01008", II[t]->"http://identifiers.org/uniprot/P00734", IIa[t]->"http://identifiers.org/uniprot/P00734", IX[t]->"http://identifiers.org/uniprot/P00740", IXa[t]->"http://identifiers.org/uniprot/P00740", TF[t]->"http://identifiers.org/uniprot/P13726", TFPI[t]->"http://identifiers.org/uniprot/P10646", V[t]->"http://identifiers.org/uniprot/P12259", VII[t]->"http://identifiers.org/uniprot/P08709", VIII[t]->"http://identifiers.org/uniprot/P00451", VIIIa[t]->"http://identifiers.org/uniprot/P00451", VIIIa1\[LetterSpace]L[t]->"http://identifiers.org/uniprot/P00451", VIIIa2[t]->"http://identifiers.org/uniprot/P00451", VIIa[t]->"http://identifiers.org/uniprot/P08709", Va[t]->"http://identifiers.org/uniprot/P12259", X[t]->"http://identifiers.org/uniprot/P00742", Xa[t]->"http://identifiers.org/uniprot/P00742", mIIa[t]->"http://identifiers.org/uniprot/P00734"
};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
ATIII'[t] ==  -1.0*R23 -1.0*R24 -1.0*R25 -1.0*R26 -1.0*R27, II'[t] ==  -1.0*R9 -1.0*R18, IIa'[t] == 1.0*R5 +1.0*R9 +1.0*R10 +1.0*R16 +1.0*R19 -1.0*R5 -1.0*R10 -1.0*R16 -1.0*R26, IIa\[LetterSpace]ATIII'[t] == 1.0*R26 , IX'[t] ==  -1.0*R8, IXa'[t] == 1.0*R14 +1.0*R15 +1.0*R8b +1.0*R28 -1.0*R11 -1.0*R25 -1.0*R28, IXa\[LetterSpace]ATIII'[t] == 1.0*R25 , IXa\[LetterSpace]VIIIa'[t] == 1.0*R11 +1.0*R12b -1.0*R12 -1.0*R15, IXa\[LetterSpace]VIIIa\[LetterSpace]X'[t] == 1.0*R12 -1.0*R14 -1.0*R12b, TF'[t] ==  -1.0*R1 -1.0*R2, TFPI'[t] ==  -1.0*R20 -1.0*R21, TF\[LetterSpace]VII'[t] == 1.0*R1 , TF\[LetterSpace]VIIa'[t] == 1.0*R2 +1.0*R3 +1.0*R8b -1.0*R3 -1.0*R6 -1.0*R7 -1.0*R8 -1.0*R22 -1.0*R27, TF\[LetterSpace]VIIa\[LetterSpace]ATIII'[t] == 1.0*R27 , TF\[LetterSpace]VIIa\[LetterSpace]IX'[t] == 1.0*R8 -1.0*R8b, TF\[LetterSpace]VIIa\[LetterSpace]X'[t] == 1.0*R6 -1.0*R6b, TF\[LetterSpace]VIIa\[LetterSpace]Xa'[t] == 1.0*R7 +1.0*R6b -1.0*R21, TF\[LetterSpace]VIIa\[LetterSpace]Xa\[LetterSpace]TFPI'[t] == 1.0*R21 +1.0*R22 , V'[t] ==  -1.0*R16 -1.0*R29, VII'[t] ==  -1.0*R1 -1.0*R3 -1.0*R4 -1.0*R5, VIII'[t] ==  -1.0*R10, VIIIa'[t] == 1.0*R10 -1.0*R11 -1.0*R13, VIIIa1\[LetterSpace]L'[t] == 1.0*R13 +1.0*R14 +1.0*R15 , VIIIa2'[t] == 1.0*R13 +1.0*R14 +1.0*R15 , VIIa'[t] == 1.0*R3 +1.0*R4 +1.0*R5 -1.0*R2, Va'[t] == 1.0*R16 +1.0*R29 -1.0*R17, X'[t] == 1.0*R14 -1.0*R6 -1.0*R12 -1.0*R28, Xa'[t] == 1.0*R4 +1.0*R9 +1.0*R12b +1.0*R28 -1.0*R4 -1.0*R7 -1.0*R9 -1.0*R17 -1.0*R20 -1.0*R23, Xa\[LetterSpace]ATIII'[t] == 1.0*R23 , Xa\[LetterSpace]TFPI'[t] == 1.0*R20 -1.0*R22, Xa\[LetterSpace]Va'[t] == 1.0*R17 +1.0*R19 +1.0*R18b -1.0*R18 -1.0*R19, Xa\[LetterSpace]Va\[LetterSpace]II'[t] == 1.0*R18 -1.0*R18b, mIIa'[t] == 1.0*R18b +1.0*R29 -1.0*R19 -1.0*R24 -1.0*R29, mIIa\[LetterSpace]ATIII'[t] == 1.0*R24 
};
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]}]