81
82
83
84
85
86
87
88
89
90
91
92
|
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
|
%% Effort on port 2
FOR ALL A_1,A_2,A_3,A_4,v_1,v_2,v_3,v_4,kappa
LET MA(R,flow,kappa,effort,2,
A_1,effort,1,
v_2,flow,2,
v_3,flow,3,
A_4,effort,4) = RT*log( (v_2/kappa) + exp(A_1/RT) );
%%% Five port version (stoichiometric) with integral causality
%% Flow on port 1 = flow on port 3.
FOR ALL A_f, A_r, v_f, v_r, v, kappa
LET MA(R, flow,kappa,flow,1,
A_f,effort,1,
A_r,effort,2,
v_f,flow,3,
v_r,flow,4,
v,effort,5) = v_f;
%% Flow on port 2 = flow on port 4.
FOR ALL A_f, A_r, v_f, v_r, v, kappa
LET MA(R, flow,kappa,flow,2,
A_f,effort,1,
A_r,effort,2,
v_f,flow,3,
v_r,flow,4,
v,effort,5) = v_r;
%% Flow on port 5 is induced flow
FOR ALL A_f, A_r, v_f, v_r, v, kappa
LET MA(R, flow,kappa,flow,5,
A_f,effort,1,
A_r,effort,2,
v_f,flow,3,
v_r,flow,4,
v,effort,5) = kappa*(exp(A_f/RT) - exp(A_r/RT));
%% AE version
FOR ALL mu
LET MA(AE, effort,2, mu, effort,1) = exp(mu/RT);
|