reset()

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

K_A ;
assumeAlso(K_A  > 0):
assumeAlso(K_A , R_)

K_B_s_1 ;
assumeAlso(K_B_s_1  > 0):
assumeAlso(K_B_s_1 , R_):

K_B_s_2 ;
assumeAlso(K_B_s_2  > 0):
assumeAlso(K_B_s_2 , R_):

K_B_s_3 ;
assumeAlso(K_B_s_3  > 0):
assumeAlso(K_B_s_3 , R_):

K_B_s_4 ;
assumeAlso(K_B_s_4  > 0):
assumeAlso(K_B_s_4 , R_):

K_B_s_5 ;
assumeAlso(K_B_s_5  > 0):
assumeAlso(K_B_s_5 , R_):

Rtot;
assumeAlso(Rtot>0):
assumeAlso(Rtot,R_):

Ltot;
assumeAlso(Ltot>0):
assumeAlso(Ltot,R_):

Req;
assumeAlso(Req>0):
assumeAlso(Req<=Rtot):
assumeAlso(Req,R_):

R_s_1eq;
assumeAlso(R_s_1eq>0):
assumeAlso(R_s_1eq<=Rtot):
assumeAlso(R_s_1eq,R_):

R_s_2eq;
assumeAlso(R_s_2eq>0):
assumeAlso(R_s_2eq<=Rtot):
assumeAlso(R_s_2eq,R_):

R_s_3eq;
assumeAlso(R_s_3eq>0):
assumeAlso(R_s_3eq<=Rtot):
assumeAlso(R_s_3eq,R_):

R_s_4eq;
assumeAlso(R_s_4eq>0):
assumeAlso(R_s_4eq<=Rtot):
assumeAlso(R_s_4eq,R_):

R_s_5eq;
assumeAlso(R_s_5eq>0):
assumeAlso(R_s_5eq<=Rtot):
assumeAlso(R_s_5eq,R_):

Leq;
assumeAlso(Leq>0):
assumeAlso(Leq<Ltot):
assumeAlso(Leq,R_):

RLeq;
assumeAlso(RLeq>0):
assumeAlso(RLeq<Rtot):
assumeAlso(RLeq<Ltot):
assumeAlso(RLeq,R_):

anames(Properties,User);

eq2_1:= Rtot = Req + R_s_1eq + R_s_2eq +  R_s_3eq + R_s_4eq +  R_s_5eq + RLeq;
eq2_2:= Ltot = Leq + RLeq;

eq2_3:= K_A = RLeq / (Req*Leq);

eq2_4:= K_B_s_1 = R_s_1eq/Req

eq2_5:= K_B_s_2 = R_s_2eq/Req

eq2_6:= K_B_s_3 = R_s_3eq/Req

eq2_7:= K_B_s_4 = R_s_4eq/Req

eq2_8:= K_B_s_5 = R_s_5eq/Req

eq2_3;
solve(%,RLeq):
%[1][1]:
eq3_1:= RLeq=%

eq2_1;
% | eq3_1:
eq3_2:= %

eq2_2;
% | eq3_1:
eq3_3:= %

eq2_4;
solve(%,R_s_1eq):
%[2][1]:
eq3_7:= R_s_1eq = %

eq3_2;
% | eq3_7:
eq3_8:= %

eq2_5;
solve(%, R_s_2eq):
%[2][1]:
eq3_9:= R_s_2eq = %

eq3_8;
% | eq3_9:
eq3_10:= %

eq2_6;
solve(%, R_s_3eq):
%[2][1]:
eq3_11:= R_s_3eq = %

eq3_10;
% | eq3_11:
eq3_12:= %

eq2_7;
solve(%, R_s_4eq):
%[2][1]:
eq3_13:= R_s_4eq = %

eq3_12;
% | eq3_13:
eq3_14:= %

eq2_8;
solve(%, R_s_5eq):
%[2][1]:
eq3_15:= R_s_5eq = %

eq3_14;
% | eq3_15:
eq3_16:= %

eq3_3;
solve(%,Req):
%[1][1]:
eq3_17:= Req = %

eq3_16;
% | eq3_17;
temp1:=%[2];
temp2:=Simplify(%);
// test
temp1=temp2;
Simplify(%);
// Assemble a final equation
eq3_18:= Rtot =temp2

solution3_18:=solve(eq3_18, Leq)

solution_lines:=6:
eq3_19:=  solution3_18[i,1] $ i=1..solution_lines;
nops(%)

solution1:=eq3_19[1][1];  
solution2:=eq3_19[2][1] ; 
solution3:=eq3_19[4][1]; 
solution4:=eq3_19[5][1]; 

/*
solutionA:=eq3_13[1];   // a sequence of roots

if solution2 in solutionA
then print(Unquoted,"First set of roots contains the second root.");
else print(Unquoted,"First  set of roots  DOES NOT contain the second root!");
end_if;

if solution3 in solutionA
then print(Unquoted,"First set of roots contains the third root.");
else print(Unquoted,"First  set of roots  DOES NOT contain the third root!");
end_if;
*/

eq3_18;
// Check the 1st solution
test1:=eq3_18 | Leq=solution1;
normal(%);
Simplify(%);
bool(%)

eq3_18;
// Check the 2nd solution
test1:=eq3_18 | Leq=solution2;
normal(%);
Simplify(%);
bool(%)

eq3_18;
// Check the 3rd solution
test1:=eq3_18 | Leq=solution3;
normal(%);
Simplify(%);
bool(%)

eq3_18;
// Check the 4th solution
test1:=eq3_18 | Leq=solution4;
normal(%);
Simplify(%);
bool(%)

solution1;
% | K_A=1 | K_B_s_1=1 | K_B_s_2=1 |  K_B_s_3=1 | K_B_s_4=1 |  K_B_s_5=1 |  Rtot=1 | Ltot=1;
float(%)

solution2;
% | K_A=1 | K_B_s_1=1 | K_B_s_2=1 |  K_B_s_3=1 | K_B_s_4=1 |  K_B_s_5=1 |  Rtot=1 | Ltot=1;
float(%)

solution3;
% | K_A=1 | K_B_s_1=1 | K_B_s_2=1 |  K_B_s_3=1 | K_B_s_4=1 |  K_B_s_5=1 |  Rtot=1 | Ltot=1;
float(%)

solution4;
% | K_A=1 | K_B_s_1=1 | K_B_s_2=1 |  K_B_s_3=1 | K_B_s_4=1 |  K_B_s_5=1 |  Rtot=1 | Ltot=1;
float(%)

// Solution 2
test1:=eq3_18 | Leq=solution2;
result:= % | K_A=1 | K_B_s_1=1 | K_B_s_2=1 |  K_B_s_3=1 | K_B_s_4=1 |  K_B_s_5=1 |  Rtot=1 | Ltot=1;
float(result);
Simplify(result);

// Solution 4
test1:=eq3_18 | Leq=solution4;
result:= % | K_A=1 | K_B_s_1=1 | K_B_s_2=1 |  K_B_s_3=1 | K_B_s_4=1 |  K_B_s_5=1 |  Rtot=1 | Ltot=1;
float(result);
Simplify(result);

eq3_20:= Leq = solution4

anames(Properties,User);

fLeq_U_5R:= (Rtot, Ltot, K_A, K_B_s_1, K_B_s_2, K_B_s_3, K_B_s_4, K_B_s_5) --> Leq_U_5R[2]

//test operation
fLeq_U_5R(1,1,1,1,1,1,1,1):
float(%)

fReq_U_5R:=(Rtot, Ltot, K_A, K_B_s_1, K_B_s_2, K_B_s_3, K_B_s_4, K_B_s_5) --> Req_U_5R[2] | Leq_U_5R

//test operation
fReq_U_5R(1,1,1,1,1,1,1,1):
float(%)

fR_s_1eq_U_5R:=(Rtot, Ltot, K_A, K_B_s_1, K_B_s_2, K_B_s_3, K_B_s_4, K_B_s_5) --> R_s_1eq_U_5R[2] | Req_U_5R | Leq_U_5R

//test operation
fR_s_1eq_U_5R(1,1,1,1,1,1,1,1):
float(%)

fR_s_2eq_U_5R:=(Rtot, Ltot, K_A, K_B_s_1, K_B_s_2, K_B_s_3, K_B_s_4, K_B_s_5) --> R_s_2eq_U_5R[2] | Req_U_5R | Leq_U_5R

//test operation
fR_s_2eq_U_5R(1,1,1,1,1,1,1,1):
float(%)

fR_s_3eq_U_5R:=(Rtot, Ltot, K_A, K_B_s_1, K_B_s_2, K_B_s_3, K_B_s_4, K_B_s_5) --> R_s_3eq_U_5R[2] | Req_U_5R | Leq_U_5R

//test operation
fR_s_3eq_U_5R(1,1,1,1,1,1,1,1):
float(%)

fR_s_4eq_U_5R:=(Rtot, Ltot, K_A, K_B_s_1, K_B_s_2, K_B_s_3, K_B_s_4, K_B_s_5) --> R_s_4eq_U_5R[2] | Req_U_5R | Leq_U_5R

//test operation
fR_s_4eq_U_5R(1,1,1,1,1,1,1,1):
float(%)

fR_s_5eq_U_5R:=(Rtot, Ltot, K_A, K_B_s_1, K_B_s_2, K_B_s_3, K_B_s_4, K_B_s_5) --> R_s_5eq_U_5R[2] | Req_U_5R | Leq_U_5R

//test operation
fR_s_5eq_U_5R(1,1,1,1,1,1,1,1):
float(%)

fRLeq_U_5R:=(Rtot, Ltot, K_A, K_B_s_1, K_B_s_2, K_B_s_3, K_B_s_4, K_B_s_5) --> RLeq_U_5R[2] | Req_U_5R | Leq_U_5R

//test operation
fRLeq_U_5R(1,1,1,1,1,1,1,1):
float(%)

fRLeq_U_5R(1,1,1,1,1,1,1,1)+fLeq_U_5R(1,1,1,1,1,1,1,1)

fRLeq_U_5R(1,1,1,1,1,1,1,1) + fReq_U_5R(1,1,1,1,1,1,1,1) + fR_s_1eq_U_5R(1,1,1,1,1,1,1,1)  + fR_s_2eq_U_5R(1,1,1,1,1,1,1,1)  + fR_s_3eq_U_5R(1,1,1,1,1,1,1,1)  + fR_s_4eq_U_5R(1,1,1,1,1,1,1,1)  + fR_s_5eq_U_5R(1,1,1,1,1,1,1,1) ;
float(%)

ProjectName

filename:=CurrentPath."/".ProjectName.".mb";
write(filename,

// Equations
Leq_U_5R,
Req_U_5R,
R_s_1eq_U_5R,
R_s_2eq_U_5R,
R_s_3eq_U_5R,
R_s_4eq_U_5R,
R_s_5eq_U_5R,
RLeq_U_5R,

// Analytical functions
fLeq_U_5R,
fReq_U_5R,
fR_s_1eq_U_5R,
fR_s_2eq_U_5R,
fR_s_3eq_U_5R,
fR_s_4eq_U_5R,
fR_s_5eq_U_5R,
fRLeq_U_5R
)