Rep.m

From 2006.igem.org

function pFactor = rep(r,K,n);

    % Jörg's function:
    % pFactor = 1./(1+(r./K).^n);
    
    % modified:
    % r: concentration of repressor
    % K: Hill constant (Affinity)
    % n: cooperativity, i.e. how steep the curve will be (plot([0:0.05:3], rep([0:0.05:3], 1, 5)))
    
    pFactor = 1./(1+K.*r.^n);