reset()

ProjectName:="U-2R";
CurrentPath:="/Users/kovrigin_laptop/Documents/Workspace/Global_Analysis/IDAP/Mathematical_models/Equilibrium_thermodynamic_models/U-multi-path-models/nR/U-2R";

filename:=CurrentPath.ProjectName.".mb";
fread(filename,Quiet):
anames(User)

Total_R:=1e-3:
Total_L:=10e-3:
Ka:=1e3:
Kb_s1:=1:
Kb_s2:=2:

fLeq_U_2R(Total_R,  Total_L, Ka, Kb_s1, Kb_s2);
fReq_U_2R(Total_R,  Total_L, Ka, Kb_s1, Kb_s2);
fR_s_1eq_U_2R(Total_R,  Total_L, Ka, Kb_s1, Kb_s2);
fR_s_2eq_U_2R(Total_R,  Total_L, Ka, Kb_s1, Kb_s2);
fRLeq_U_2R(Total_R,  Total_L, Ka, Kb_s1, Kb_s2);

fLeq:=LRratio ->      fLeq_U_2R     (Total_R,  LRratio*Total_R,  Ka, Kb_s1, Kb_s2):
fReq:=LRratio ->      fReq_U_2R     (Total_R,  LRratio*Total_R,  Ka, Kb_s1, Kb_s2):
fR_s_1eq:=LRratio ->  fR_s_1eq_U_2R (Total_R,  LRratio*Total_R,  Ka, Kb_s1, Kb_s2):
fR_s_2eq:=LRratio ->  fR_s_2eq_U_2R (Total_R,  LRratio*Total_R,  Ka, Kb_s1, Kb_s2):
fRLeq:=LRratio ->     fRLeq_U_2R    (Total_R,  LRratio*Total_R,  Ka, Kb_s1, Kb_s2):

Total_R:=1e-3:
LRratio_max:=2:
Ka:=1e6:
Kb_s1:=1:
Kb_s2:=2:

LineW:=1.5: //line width

// create plots

pLeq:=  plot::Function2d(
          Function=(fLeq),
          LegendText="[L]",
          Color = RGB::Blue,
          XMin=(0),
          XMax=(LRratio_max),
          XName=(LRratio),
          TitlePositionX=(0),
          LineWidth=LineW):


pRLeq:=  plot::Function2d(
          Function=(fRLeq),
          LegendText="[RL]",
          Color = RGB::Red,
          XMin=(0),
          XMax=(LRratio_max),
          XName=(LRratio),
          TitlePositionX=(0),
          LineWidth=LineW):


plot(pLeq, pRLeq, LegendVisible=TRUE)

Total_R:=1e-3:
LRratio_max:=2:
Ka:=1e6:
Kb_s1:=0.000001:
Kb_s2:=0.000001:

LRratio_max:=1.5: // plotting range

LineW:=1.5: // plot line width


pLeq:=  plot::Function2d(
          Function=(fLeq),
          LegendText="[L]",
          Color = RGB::Blue,
          XMin=(LRratio_max*1e-6),
          XMax=(LRratio_max),
          XName=(LRratio),
          TitlePositionX=(0),
          LineWidth=LineW):



pReq:=  plot::Function2d(
          Function=(fReq),
          LegendText="[R]",
          Color = RGB::Black,
          XMin=(LRratio_max*1e-6),
          XMax=(LRratio_max),
          XName=(LRratio),
          TitlePositionX=(0),
          LineWidth=LineW):

pR_s_1eq:=  plot::Function2d(
          Function=(fR_s_1eq),
          LegendText="[R*]",
          Color = RGB::Green,
          XMin=(LRratio_max*1e-6),
          XMax=(LRratio_max),
          XName=(LRratio),
          TitlePositionX=(0),
          LineWidth=LineW):



pR_s_2eq:=  plot::Function2d(
          Function=(fR_s_2eq),
          LegendText="[R**]",
          Color = RGB::Grey,
          XMin=(LRratio_max*1e-6),
          XMax=(LRratio_max),
          XName=(LRratio),
          TitlePositionX=(0),
          LineWidth=LineW):


pRLeq:=  plot::Function2d(
          Function=(fRLeq),
          LegendText="[RL]",
          Color = RGB::Red,
          XMin=(LRratio_max*1e-6),
          XMax=(LRratio_max),
          XName=(LRratio),
          TitlePositionX=(0),
          LineWidth=LineW):






// Text report
print(Unquoted,"Model: ".ProjectName);
print(Unquoted,"Total_R=".Total_R);
Kda:=1/Ka:
print(Unquoted,"Ka=".Ka.",   Kd=".Kda);
print(Unquoted,"Kb*=".Kb_s1);
print(Unquoted,"Kb**=".Kb_s2);
Kc:=Kb_s2/Kb_s1:
print(Unquoted,"Kc*-**=".Kc);

// plot all together
plot(pLeq, pReq,  pR_s_1eq, pR_s_2eq, pRLeq,
   YAxisTitle="[X]", Header=("Model: ".ProjectName),
   Height=160, Width=100,TicksLabelFont=["Helvetica",12,[0,0,0],Left],
  AxesTitleFont=["Helvetica",14,[0,0,0],Left],
  XGridVisible=TRUE, YGridVisible=TRUE,
  LegendVisible=TRUE, LegendFont=["Helvetica",14,[0,0,0],Left],
  ViewingBoxYMax=Total_R);

Model: U-2R
Total_R=0.001
Ka=1000000.0,   Kd=0.000001
Kb*=0.000001
Kb**=0.000001
Kc*-**=1.0