(*Codetoshowchiralsymmetrybreaking.Kineticconstantsarechosensuchthattheequilibriumconstant,(kf/kr),fordirectandcatalyzedreactionsareequal*)k1f=0.5;k1r=0.2;k2f=0.5;k2r=(k1r/k1f)*k2f;k3f=1.5;k3r=10ˆ-3;​​S=1.25;T=S;P=0.5;​​R1f:=k1f*S*T;R1r:=k1r*XL[t];​​R2f:=k2f*S*T*XL[t];R2r:=k2r*(XL[t])ˆ2;​​R3f:=k1f*S*T;R3r:=k1r*XD[t];​​R4f:=k2f*S*T*XD[t];R4r:=k2r*(XD[t])ˆ2;​​R5f:=k3f*XL[t]*XD[t];R5r:=k3r*P;​​(*InitialvaluesofXLandXDaresetatsteadystatevaluesforS=0.5*)​​Soln1=NDSolve[{XL'[t]R1f-R1r+R2f-R2r-R5f+R5r;XD'[t]R3f-R3r+R4f-R4r-R5f+R5r;XL[0]0.2502,XD[0]0.2500};{XL[t],XD[t]},{t,0,2000},MaxSteps10000]​​Plot[Evaluate[{XL[t],XD[t]}/.Soln1],{t,0,tmax},AxesLabel{"time","[XL]&[XD]"},AxesStyleDirective[Black,14],PlotStyle{{Black,Thick},{Gray,Thick}},PlotRangeAll]​​{{XLInterpolatingFunction[{{0.,2000.}},<>],XDInterpolatingFunction[{{0.,2000.}},<>]}}​​Export["data1.txt",Table[{t,{XL[t],XD[t]}/.Soln1},{t,1,50}]//TableForm,"List"]
NDSolve
:MaxSteps10000 cannot be used as a variable.
NDSolve[{XL[t],XD[t]},{t,0,2000},MaxSteps10000]
ReplaceAll
:{NDSolve[{XL[t],XD[t]},{t,0,2000},MaxSteps10000]} is neither a list of replacement rules nor a valid dispatch table, and so cannot be used for replacing.
Plot
:Limiting value tmax in {t,0,tmax} is not a machine-sized real number.
Plot[{XL[t],XD[t]}/.NDSolve[{XL[t],XD[t]},{t,0,2000},MaxSteps10000],{t,0,tmax},AxesLabel{time,[XL]&[XD]},AxesStyleDirective[Black,14],PlotStyle{{Black,Thick},{Gray,Thick}},PlotRangeAll]
NDSolve
:There are more dependent variables, {XD[t],XL[t]}​, than equations, so the system is underdetermined.
NDSolve[{
′
XL
[t]0.125+0.5(-3+10ˆ)-0.075XL[t]-0.2ˆ2XL[t]-1.5XD[t]XL[t],XL[0]0.25,XD[0]0.25},{XL,XD},{t,0,2000},MaxSteps10000]
ReplaceAll
:{NDSolve[{
′
XL
[t]0.125+0.5(-3+Times[2])-0.075XL[t]-0.2ˆ2XL[t]-1.5XD[t]XL[t],XL[0]0.25,XD[0]0.25},{XL,XD},{t,0,2000},MaxSteps10000]} is neither a list of replacement rules nor a valid dispatch table, and so cannot be used for replacing.
NDSolve
:0.0408571 cannot be used as a variable.
ReplaceAll
:{NDSolve[{
′
XL
[0.0408571]0.125+0.5(-3+Times[2])-0.075XL[0.0408571]-0.2ˆ2XL[0.0408571]-1.5XD[0.0408571]XL[0.0408571],XL[0]0.25,XD[0]0.25},{XL,XD},{0.0408571,0,2000},MaxSteps10000]} is neither a list of replacement rules nor a valid dispatch table, and so cannot be used for replacing.
NDSolve
:Value of option MaxSteps -> 10000. should be a positive integer or Infinity.
ReplaceAll
:{NDSolve[{
′
XL
[0.0408571]0.125+0.5(-3.+Times[2])-0.075XL[0.0408571]-0.2ˆ2XL[0.0408571]-1.5XD[0.0408571]XL[0.0408571],XL[0.]0.25,XD[0.]0.25},{XL,XD},{0.0408571,0.,2000.},MaxSteps10000.]} is neither a list of replacement rules nor a valid dispatch table, and so cannot be used for replacing.
General
:Further output of ReplaceAll::reps will be suppressed during this calculation.
NDSolve
:40.8572 cannot be used as a variable.
NDSolve
:Value of option MaxSteps -> 10000. should be a positive integer or Infinity.