Team:Paris/Modeling/f3bisprog

From 2008.igem.org

ƒ3bis

function act_pFlhDC = f3bis(EnvZ, FliA)
% expression of pTet in function of EnvZ and FliA
 
global beta17 K15 n15 beta22 K6 n6 OmpR_b EnvZ_b K14 n14;
% parameters
 
OmpR_P = complexes((OmpR_b + EnvZ),EnvZ,K14,n14);
% Quantity of "phosphorylated OmpR"
 
act_pFlhDC = ( 1 - hill( OmpR_P, K15, n15 ) ) * ...
     (beta17 * (1 - hill( FliA, K6, n6)) + beta22 * hill( FliA, K6, n6) );
 
end