Team:Paris/Modeling/More f2 Algo

From 2008.igem.org

(Difference between revisions)
(find_ƒ2)
(Inv_ƒ2)
Line 30: Line 30:
== Inv_ƒ2 ==
== Inv_ƒ2 ==
-
 
+
<html><pre class="codeinput">
 +
<span class="keyword">function</span> quant_ara = Inv_f2(inducer_quantity)
 +
<span class="comment">% gives the quantity of [ara]i needed to get inducer_quantity of a protein
 +
</span><span class="comment">% throught a gene behind pBad
 +
</span>
 +
    <span class="keyword">function</span> equa = F(x)
 +
        equa = f2( x ) - inducer_quantity;
 +
    <span class="keyword">end</span>
 +
 +
options=optimset(<span class="string">'LevenbergMarquardt'</span>,<span class="string">'on'</span>,<span class="string">'TolX'</span>,1e-10,<span class="string">'MaxFunEvals'</span>,1e10,<span class="string">'TolFun'</span>,1e-10,<span class="string">'MaxIter'</span>,1e4);
 +
 +
quant_ara = fsolve(F,1,options);
 +
 +
<span class="keyword">end</span>
 +
</pre></html>
</div>
</div>

Revision as of 03:44, 30 October 2008

find_ƒ2

function act_pBad = f2(arab)
% activity of pBad in function of arabinose
 
global betabad Kbad nbad Kara nara;
% parameters
 
act_pBad = betabad * ( hill( hill( arab, Kara, nara ), Kbad, nbad ) );
 
end

Inv_ƒ2

function quant_ara = Inv_f2(inducer_quantity)
% gives the quantity of [ara]i needed to get inducer_quantity of a protein
% throught a gene behind pBad
 
     function equa = F(x)
         equa = f2( x ) - inducer_quantity;
     end
 
options=optimset('LevenbergMarquardt','on','TolX',1e-10,'MaxFunEvals',1e10,'TolFun',1e-10,'MaxIter',1e4);
 
quant_ara = fsolve(F,1,options);
 
end