(* 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 = {
ERKMEKPP[t], ERKP[t], ERKPMEKPP[t], ERKPP[t], ERKPPPase3[t], ERKPPase3[t], GS[t], MEKP[t], MEKPP[t], MEKPPPase2[t], MEKPPase2[t], MEKPRafa[t], MEKRafa[t], PLCgP[t], PLCgPI[t], R2[t], R2P[t], RG[t], RGAP[t], RGS[t], RGSr[t], RGSrD[t], RGSrT[t], RGrT[t], RPL[t], RPLP[t], RSh[t], RShGSr[t], RShGSrD[t], RShGSrT[t], RShP[t], RShPG[t], RShPGS[t], Ra[t], Rafa[t], RafaPase1[t], RafrT[t], ShP[t], ShPG[t], ShPGS[t], rasD[t], rasT[t]
};

initialValues = {
ERKMEKPP[0] == 0.0, ERKP[0] == 0.0, ERKPMEKPP[0] == 0.0, ERKPP[0] == 0.0, ERKPPPase3[0] == 0.0, ERKPPase3[0] == 0.0, GS[0] == 27.0125, MEKP[0] == 0.0, MEKPP[0] == 0.0, MEKPPPase2[0] == 0.0, MEKPPase2[0] == 0.0, MEKPRafa[0] == 0.0, MEKRafa[0] == 0.0, PLCgP[0] == 0.0, PLCgPI[0] == 0.0, R2[0] == 0.0, R2P[0] == 0.0, RG[0] == 0.0, RGAP[0] == 0.0, RGS[0] == 0.0, RGSr[0] == 0.0, RGSrD[0] == 0.0, RGSrT[0] == 0.0, RGrT[0] == 0.0, RPL[0] == 0.0, RPLP[0] == 0.0, RSh[0] == 0.0, RShGSr[0] == 0.0, RShGSrD[0] == 0.0, RShGSrT[0] == 0.0, RShP[0] == 0.0, RShPG[0] == 0.0, RShPGS[0] == 0.0, Ra[0] == 0.0, Rafa[0] == 0.0, RafaPase1[0] == 0.0, RafrT[0] == 0.0, ShP[0] == 0.0, ShPG[0] == 0.0, ShPGS[0] == 0.0, rasD[0] == 77.3499, rasT[0] == 7.65
};

rates = {
v\[LetterSpace]1, v\[LetterSpace]10, v\[LetterSpace]11, v\[LetterSpace]12, v\[LetterSpace]13, v\[LetterSpace]14, v\[LetterSpace]15, v\[LetterSpace]16, v\[LetterSpace]17, v\[LetterSpace]18, v\[LetterSpace]19, v\[LetterSpace]2, v\[LetterSpace]20, v\[LetterSpace]21, v\[LetterSpace]22, v\[LetterSpace]23, v\[LetterSpace]24, v\[LetterSpace]25, v\[LetterSpace]26, v\[LetterSpace]27, v\[LetterSpace]28, v\[LetterSpace]29, v\[LetterSpace]3, v\[LetterSpace]30, v\[LetterSpace]31, v\[LetterSpace]32, v\[LetterSpace]33, v\[LetterSpace]34, v\[LetterSpace]35, v\[LetterSpace]36, v\[LetterSpace]37, v\[LetterSpace]38, v\[LetterSpace]39, v\[LetterSpace]4, v\[LetterSpace]40, v\[LetterSpace]41, v\[LetterSpace]42, v\[LetterSpace]43, v\[LetterSpace]44, v\[LetterSpace]45, v\[LetterSpace]46, v\[LetterSpace]47, v\[LetterSpace]48, v\[LetterSpace]49, v\[LetterSpace]5, v\[LetterSpace]50, v\[LetterSpace]51, v\[LetterSpace]52, v\[LetterSpace]53, v\[LetterSpace]54, v\[LetterSpace]55, v\[LetterSpace]56, v\[LetterSpace]57, v\[LetterSpace]58, v\[LetterSpace]59, v\[LetterSpace]6, v\[LetterSpace]60, v\[LetterSpace]61, v\[LetterSpace]7, v\[LetterSpace]8, v\[LetterSpace]9
};

rateEquations = {
v\[LetterSpace]1 -> EGF*kon1*R - koff1*Ra[t], v\[LetterSpace]10 -> Grb*kon10*R2P[t] - koff10*RG[t], v\[LetterSpace]11 -> kon11*SOS*RG[t] - koff11*RGS[t], v\[LetterSpace]12 -> -(koff12*GS[t]*R2P[t]) + kon12*RGS[t], v\[LetterSpace]13 -> -(Grb*koff13*SOS) + kon13*GS[t], v\[LetterSpace]14 -> kon14*Shc*R2P[t] - koff14*RSh[t], v\[LetterSpace]15 -> kon15*RSh[t] - koff15*RShP[t], v\[LetterSpace]16 -> kon16*RShP[t] - koff16*R2P[t]*ShP[t], v\[LetterSpace]17 -> Grb*kon17*RShP[t] - koff17*RShPG[t], v\[LetterSpace]18 -> kon18*RShPG[t] - koff18*R2P[t]*ShPG[t], v\[LetterSpace]19 -> kon19*SOS*RShPG[t] - koff19*RShPGS[t], v\[LetterSpace]2 -> -(koff2*R2[t]) + kon2*Ra[t]^2, v\[LetterSpace]20 -> kon20*RShPGS[t] - koff20*R2P[t]*ShPGS[t], v\[LetterSpace]21 -> kon21*SOS*ShPG[t] - koff21*ShPGS[t], v\[LetterSpace]22 -> Grb*kon22*ShP[t] - koff22*ShPG[t], v\[LetterSpace]23 -> (vmax23*ShP[t])/(kmax23 + ShP[t]), v\[LetterSpace]24 -> -(koff24*GS[t]*ShP[t]) + kon24*ShPGS[t], v\[LetterSpace]25 -> kon25*GS[t]*RShP[t] - koff25*RShPGS[t], v\[LetterSpace]26 -> -(GDP*koff26*ras) + kon26*rasD[t], v\[LetterSpace]27 -> GTP*kon27*ras - koff27*rasT[t], v\[LetterSpace]28 -> kon28*rasT[t], v\[LetterSpace]29 -> -(koff29*vd*RShGSrD[t]) + (kon29*va*rasD[t]*RShPGS[t])/(1 + (kii*ERKPP[t])^kn), v\[LetterSpace]3 -> kon3*R2[t] - koff3*R2P[t], v\[LetterSpace]30 -> -(GDP*koff30*RShGSr[t]) + kon30*RShGSrD[t], v\[LetterSpace]31 -> GTP*kon31*RShGSr[t] - koff31*RShGSrT[t], v\[LetterSpace]32 -> kon32*vd*RShGSrT[t] - koff32*va*rasT[t]*RShPGS[t], v\[LetterSpace]33 -> -(koff33*vd*RShGSr[t]) + (kon33*ras*va*RShPGS[t])/(1 + (kii*ERKPP[t])^kn), v\[LetterSpace]34 -> (kon34*va*rasD[t]*RGS[t])/(1 + (kii*ERKPP[t])^kn) - koff34*vd*RGSrD[t], v\[LetterSpace]35 -> -(GDP*koff35*RGSr[t]) + kon35*RGSrD[t], v\[LetterSpace]36 -> GTP*kon36*RGSr[t] - koff36*RGSrT[t], v\[LetterSpace]37 -> -(koff37*va*rasT[t]*RGS[t]) + kon37*vd*RGSrT[t], v\[LetterSpace]38 -> (kon38*ras*va*RGS[t])/(1 + (kii*ERKPP[t])^kn) - koff38*vd*RGSr[t], v\[LetterSpace]39 -> -(koff39*RafrT[t]) + kon39*Raf*rasT[t], v\[LetterSpace]4 -> (vmax4*R2P[t])/(kmax4 + R2P[t]), v\[LetterSpace]40 -> kon40*RafrT[t] - koff40*Rafa[t]*rasT[t], v\[LetterSpace]41 -> kon41*Pase1*Rafa[t] - koff41*RafaPase1[t], v\[LetterSpace]42 -> kon42*RafaPase1[t], v\[LetterSpace]43 -> -(koff43*MEKRafa[t]) + kon43*MEK*Rafa[t], v\[LetterSpace]44 -> kon44*MEKRafa[t], v\[LetterSpace]45 -> -(koff45*MEKPRafa[t]) + kon45*MEKP[t]*Rafa[t], v\[LetterSpace]46 -> kon46*MEKPRafa[t], v\[LetterSpace]47 -> kon47*Pase2*MEKPP[t] - koff47*MEKPPPase2[t], v\[LetterSpace]48 -> kon48*MEKPPPase2[t], v\[LetterSpace]49 -> kon49*Pase2*MEKP[t] - koff49*MEKPPase2[t], v\[LetterSpace]5 -> kon5*PLCg*R2P[t] - koff5*RPL[t], v\[LetterSpace]50 -> kon50*MEKPPase2[t], v\[LetterSpace]51 -> -(koff51*ERKMEKPP[t]) + ERK*kon51*MEKPP[t], v\[LetterSpace]52 -> kon52*ERKMEKPP[t], v\[LetterSpace]53 -> -(koff53*ERKPMEKPP[t]) + kon53*ERKP[t]*MEKPP[t], v\[LetterSpace]54 -> kon54*ERKPMEKPP[t], v\[LetterSpace]55 -> kon55*Pase3*ERKPP[t] - koff55*ERKPPPase3[t], v\[LetterSpace]56 -> kon56*ERKPPPase3[t], v\[LetterSpace]57 -> kon57*Pase3*ERKP[t] - koff57*ERKPPase3[t], v\[LetterSpace]58 -> kon58*ERKPPase3[t], v\[LetterSpace]59 -> GAP*ka*kon59*R2P[t] - koff59*RGAP[t], v\[LetterSpace]6 -> kon6*RPL[t] - koff6*RPLP[t], v\[LetterSpace]60 -> ka*kon60*va*rasT[t]*RGAP[t] - koff60*vd*RGrT[t], v\[LetterSpace]61 -> ka*kon61*RGrT[t], v\[LetterSpace]7 -> -(koff7*PLCgP[t]*R2P[t]) + kon7*RPLP[t], v\[LetterSpace]8 -> (vmax8*PLCgP[t])/(kmax8 + PLCgP[t]), v\[LetterSpace]9 -> kon9*PLCgP[t] - koff9*PLCgPI[t]
};

parameters = {
Pase3 -> 0.0, Pool1 -> 680.0, Pool10 -> 200.0, Pool11 -> 200.0, Pool12 -> 50.0, Pool13 -> 50.0, Pool14 -> 100.0, Pool2 -> 100.0, Pool3 -> 105.0, Pool4 -> 34.0, Pool5 -> 150.0, Pool6 -> 85.0, Pool7 -> 85.0, Pool8 -> 12.0, Pool9 -> 50.0, ka -> 0.83, kii -> 3.0, kmax23 -> 340.0, kmax4 -> 50.0, kmax8 -> 100.0, kn -> 1.0, koff1 -> 0.06, koff10 -> 0.2, koff11 -> 0.06, koff12 -> 0.0028, koff13 -> 0.0001, koff14 -> 0.6, koff15 -> 0.06, koff16 -> 0.0009, koff17 -> 0.1, koff18 -> 0.0009, koff19 -> 0.0214, koff2 -> 0.1, koff20 -> 0.00024, koff21 -> 0.064, koff22 -> 0.1, koff24 -> 0.021, koff25 -> 0.0429, koff26 -> 0.00027, koff27 -> 0.00078, koff29 -> 0.76, koff3 -> 0.01, koff30 -> 0.093, koff31 -> 2.4, koff32 -> 0.0063, koff33 -> 0.001, koff34 -> 1.2, koff35 -> 0.1, koff36 -> 80.0, koff37 -> 0.005, koff38 -> 0.0016, koff39 -> 0.0053, koff40 -> 0.0007, koff41 -> 0.2, koff43 -> 0.01833, koff45 -> 0.01833, koff47 -> 0.8, koff49 -> 0.5, koff5 -> 0.2, koff51 -> 0.033, koff53 -> 0.033, koff55 -> 0.6, koff57 -> 0.5, koff59 -> 0.0, koff6 -> 0.05, koff60 -> 0.0, koff7 -> 0.006, koff9 -> 0.03, kon1 -> 0.003, kon10 -> 0.0015, kon11 -> 0.01, kon12 -> 0.15, kon13 -> 0.0015, kon14 -> 0.09, kon15 -> 6.0, kon16 -> 0.3, kon17 -> 0.003, kon18 -> 0.3, kon19 -> 0.01, kon2 -> 0.01, kon20 -> 0.12, kon21 -> 0.03, kon22 -> 0.003, kon24 -> 0.1, kon25 -> 0.009, kon26 -> 5.4*^-06, kon27 -> 0.078, kon28 -> 0.015, kon29 -> 1.9*^-05, kon3 -> 1.0, kon30 -> 46.5, kon31 -> 0.003, kon32 -> 806.4, kon33 -> 0.000625, kon34 -> 3*^-05, kon35 -> 50.0, kon36 -> 0.1, kon37 -> 640.0, kon38 -> 0.001, kon39 -> 0.01, kon40 -> 1.0, kon41 -> 0.0717, kon42 -> 1.0, kon43 -> 0.0111, kon44 -> 3.5, kon45 -> 0.0111, kon46 -> 2.9, kon47 -> 0.0143, kon48 -> 0.058, kon49 -> 0.00025, kon5 -> 0.06, kon50 -> 0.058, kon51 -> 0.01, kon52 -> 16.0, kon53 -> 0.01, kon54 -> 5.7, kon55 -> 0.0145, kon56 -> 0.27, kon57 -> 0.05, kon58 -> 0.3, kon59 -> 0.0, kon6 -> 1.0, kon60 -> 0.0, kon61 -> 0.0, kon7 -> 0.3, kon9 -> 1.0, va -> 250.0, vd -> 1.0, vmax23 -> 1.7, vmax4 -> 450.0, vmax8 -> 1.0, GDP -> 500.0, GTP -> 10000.0, sink -> 1.0, source -> 1.0, default\[LetterSpace]compartment -> 1.0
};

assignments = {
Shc -> Pool5 - RSh[t] - RShGSr[t] - RShGSrD[t] - RShGSrT[t] - RShP[t] - RShPG[t] - RShPGS[t] - ShP[t] - ShPG[t] - ShPGS[t], SOS -> Pool4 - GS[t] - RGS[t] - RGSr[t] - RGSrD[t] - RGSrT[t] - RShGSr[t] - RShGSrD[t] - RShGSrT[t] - RShPGS[t] - ShPGS[t], PLCg -> Pool3 - PLCgP[t] - PLCgPI[t] - RPL[t] - RPLP[t], Pase1 -> Pool12 - RafaPase1[t], Pase2 -> Pool13 - MEKPPase2[t] - MEKPPPase2[t], ERK -> Pool11 - ERKMEKPP[t] - ERKP[t] - ERKPMEKPP[t] - ERKPP[t] - ERKPPase3[t] - ERKPPPase3[t], GAP -> Pool8 - RGAP[t] - RGrT[t], R -> Pool2 - Ra[t] - 2*(R2[t] + R2P[t] + RG[t] + RGAP[t] + RGrT[t] + RGS[t] + RGSr[t] + RGSrD[t] + RGSrT[t] + RPL[t] + RPLP[t] + RSh[t] + RShGSr[t] + RShGSrD[t] + RShGSrT[t] + RShP[t] + RShPG[t] + RShPGS[t]), EGF -> Pool1 - Ra[t] - 2*(R2[t] + R2P[t] + RG[t] + RGAP[t] + RGrT[t] + RGS[t] + RGSr[t] + RGSrD[t] + RGSrT[t] + RPL[t] + RPLP[t] + RSh[t] + RShGSr[t] + RShGSrD[t] + RShGSrT[t] + RShP[t] + RShPG[t] + RShPGS[t]), Grb -> Pool6 - GS[t] - RG[t] - RGS[t] - RGSr[t] - RGSrD[t] - RGSrT[t] - RShGSr[t] - RShGSrD[t] - RShGSrT[t] - RShPG[t] - RShPGS[t] - ShPG[t] - ShPGS[t], ras -> Pool7 - RafrT[t] - rasD[t] - rasT[t] - RGrT[t] - RGSr[t] - RGSrD[t] - RGSrT[t] - RShGSr[t] - RShGSrD[t] - RShGSrT[t], Raf -> Pool9 - MEKPRafa[t] - MEKRafa[t] - Rafa[t] - RafaPase1[t] - RafrT[t], MEK -> Pool10 - ERKMEKPP[t] - ERKPMEKPP[t] - MEKP[t] - MEKPP[t] - MEKPPase2[t] - MEKPPPase2[t] - MEKPRafa[t] - MEKRafa[t]
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
ERKMEKPP'[t] == 1.0*v\[LetterSpace]51 -1.0*v\[LetterSpace]52, ERKP'[t] == 1.0*v\[LetterSpace]56 +1.0*v\[LetterSpace]52 -1.0*v\[LetterSpace]53 -1.0*v\[LetterSpace]57, ERKPMEKPP'[t] == 1.0*v\[LetterSpace]53 -1.0*v\[LetterSpace]54, ERKPP'[t] == 1.0*v\[LetterSpace]54 -1.0*v\[LetterSpace]55, ERKPPPase3'[t] == 1.0*v\[LetterSpace]55 -1.0*v\[LetterSpace]56, ERKPPase3'[t] == 1.0*v\[LetterSpace]57 -1.0*v\[LetterSpace]58, GS'[t] == 1.0*v\[LetterSpace]12 +1.0*v\[LetterSpace]24 -1.0*v\[LetterSpace]13 -1.0*v\[LetterSpace]25, MEKP'[t] == 1.0*v\[LetterSpace]44 +1.0*v\[LetterSpace]48 -1.0*v\[LetterSpace]45 -1.0*v\[LetterSpace]49, MEKPP'[t] == 1.0*v\[LetterSpace]46 +1.0*v\[LetterSpace]54 +1.0*v\[LetterSpace]52 -1.0*v\[LetterSpace]51 -1.0*v\[LetterSpace]53 -1.0*v\[LetterSpace]47, MEKPPPase2'[t] == 1.0*v\[LetterSpace]47 -1.0*v\[LetterSpace]48, MEKPPase2'[t] == 1.0*v\[LetterSpace]49 -1.0*v\[LetterSpace]50, MEKPRafa'[t] == 1.0*v\[LetterSpace]45 -1.0*v\[LetterSpace]46, MEKRafa'[t] == 1.0*v\[LetterSpace]43 -1.0*v\[LetterSpace]44, PLCgP'[t] == 1.0*v\[LetterSpace]7 -1.0*v\[LetterSpace]9 -1.0*v\[LetterSpace]8, PLCgPI'[t] == 1.0*v\[LetterSpace]9 , R2'[t] == 1.0*v\[LetterSpace]4 +1.0*v\[LetterSpace]2 -1.0*v\[LetterSpace]3, R2P'[t] == 1.0*v\[LetterSpace]18 +1.0*v\[LetterSpace]7 +1.0*v\[LetterSpace]3 +1.0*v\[LetterSpace]16 +1.0*v\[LetterSpace]12 +1.0*v\[LetterSpace]20 -1.0*v\[LetterSpace]14 -1.0*v\[LetterSpace]10 -1.0*v\[LetterSpace]5 -1.0*v\[LetterSpace]4 -1.0*v\[LetterSpace]59, RG'[t] == 1.0*v\[LetterSpace]10 -1.0*v\[LetterSpace]11, RGAP'[t] == 1.0*v\[LetterSpace]59 +1.0*v\[LetterSpace]61 -1.0*v\[LetterSpace]60, RGS'[t] == 1.0*v\[LetterSpace]37 +1.0*v\[LetterSpace]11 -1.0*v\[LetterSpace]34 -1.0*v\[LetterSpace]12 -1.0*v\[LetterSpace]38, RGSr'[t] == 1.0*v\[LetterSpace]35 +1.0*v\[LetterSpace]38 -1.0*v\[LetterSpace]36, RGSrD'[t] == 1.0*v\[LetterSpace]34 -1.0*v\[LetterSpace]35, RGSrT'[t] == 1.0*v\[LetterSpace]36 -1.0*v\[LetterSpace]37, RGrT'[t] == 1.0*v\[LetterSpace]60 -1.0*v\[LetterSpace]61, RPL'[t] == 1.0*v\[LetterSpace]5 -1.0*v\[LetterSpace]6, RPLP'[t] == 1.0*v\[LetterSpace]6 -1.0*v\[LetterSpace]7, RSh'[t] == 1.0*v\[LetterSpace]14 -1.0*v\[LetterSpace]15, RShGSr'[t] == 1.0*v\[LetterSpace]30 +1.0*v\[LetterSpace]33 -1.0*v\[LetterSpace]31, RShGSrD'[t] == 1.0*v\[LetterSpace]29 -1.0*v\[LetterSpace]30, RShGSrT'[t] == 1.0*v\[LetterSpace]31 -1.0*v\[LetterSpace]32, RShP'[t] == 1.0*v\[LetterSpace]15 -1.0*v\[LetterSpace]16 -1.0*v\[LetterSpace]17 -1.0*v\[LetterSpace]25, RShPG'[t] == 1.0*v\[LetterSpace]17 -1.0*v\[LetterSpace]18 -1.0*v\[LetterSpace]19, RShPGS'[t] == 1.0*v\[LetterSpace]19 +1.0*v\[LetterSpace]25 +1.0*v\[LetterSpace]32 -1.0*v\[LetterSpace]29 -1.0*v\[LetterSpace]33 -1.0*v\[LetterSpace]20, Ra'[t] == 1.0*v\[LetterSpace]1 -2.0*v\[LetterSpace]2, Rafa'[t] == 1.0*v\[LetterSpace]46 +1.0*v\[LetterSpace]44 +1.0*v\[LetterSpace]40 -1.0*v\[LetterSpace]45 -1.0*v\[LetterSpace]43 -1.0*v\[LetterSpace]41, RafaPase1'[t] == 1.0*v\[LetterSpace]41 -1.0*v\[LetterSpace]42, RafrT'[t] == 1.0*v\[LetterSpace]39 -1.0*v\[LetterSpace]40, ShP'[t] == 1.0*v\[LetterSpace]16 +1.0*v\[LetterSpace]24 -1.0*v\[LetterSpace]23 -1.0*v\[LetterSpace]22, ShPG'[t] == 1.0*v\[LetterSpace]18 +1.0*v\[LetterSpace]22 -1.0*v\[LetterSpace]21, ShPGS'[t] == 1.0*v\[LetterSpace]21 +1.0*v\[LetterSpace]20 -1.0*v\[LetterSpace]24, rasD'[t] == 1.0*v\[LetterSpace]28 +1.0*v\[LetterSpace]61 -1.0*v\[LetterSpace]29 -1.0*v\[LetterSpace]34 -1.0*v\[LetterSpace]26, rasT'[t] == 1.0*v\[LetterSpace]37 +1.0*v\[LetterSpace]32 +1.0*v\[LetterSpace]27 +1.0*v\[LetterSpace]40 -1.0*v\[LetterSpace]39 -1.0*v\[LetterSpace]60 -1.0*v\[LetterSpace]28
};
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]}]