Matlab Program to Create Biological Models Assignment Solution.

COMPLEX

function df = ODEModelComplex(t, x)

    global kRpos kRneg k1pos k1neg k2pos k2neg kLpos kLneg k3pos k3neg kAneg kApos

    % Get variables from vector x

    R = x(1);

    L = x(2);

    A = x(3);

    C1 = x(4);

    C2 = x(5);

    C3 = x(6);

    dR = kRpos - kRneg*R - k1pos*R*L + k1neg*C1 - k2pos*R*A + k2neg*C2;

    dL = kLpos - kLneg*L - k1pos*R*L + k1neg*C1 - k3pos*L*A + k3neg*C3;

    dA = -kAneg*A - k2pos*R*A + k2neg*C2 - k3pos*L*A + k3neg*C3;

    dC1 = k1pos*R*L - k1neg*C1;

    dC2 = k2pos*R*A - k2neg*C2;

    dC3 = k3pos*L*A - k3neg*C3;

    df = [dR, dL, dA, dC1, dC2, dC3]';

end

SIMPLE

function df = ODEModelSimple(t, x)

% % Get variables from vector x

% R = x(1);

% L = x(2);

% A = x(3);

% C1 = x(4);

% C2 = x(5);

% C3 = x(6);

%

% % Define parameters

%

% kRneg = 0.0693;

% kRpos = R0*kRneg;

% k1pos = 4;

% k1neg = 1;

% k2pos = 0.036;

% k2neg = 0.036;

% kLneg = 0.1155;

% kLpos = L0*kLneg;

% k3pos = 0.0018;

% k3neg = 0.036;

% kAneg = 0.0014;

% kApos = 0.0014;

%

% dR = kRpos - kRneg*R - k1pos*R*L + k1neg*C1 - k2pos*R*A + k2neg*C2;

% dL = kLpos - kLneg*L - k1pos*R*L + k1neg*C1 - k3pos*L*A + k3neg*C3;

% dA = -kAneg*A - k2pos*R*A + k2neg*C2 - k3pos*L*A + k3neg*C3;

% dC1 = k1pos*R*L - k1neg*C1;

% dC2 = k2pos*R*A - k2neg*C2;

% dC3 = k3pos*L*A - k3neg*C3;

%

% df = [dR, dL, dA, dC1, dC2, dC3]';

    global kTpos kTneg kM kon koff C

    T = x(1);

    M = x(2);

    C = x(3);

    dT = kTpos - kTneg*T - kon*T*M + koff*C;

    dM = koff*C - kon*T*M - kM*M;

    dC = kon*T*M - koff*C;

    df = [dT, dM, dC]';

end