ProjectNMS/enzReact.m at main · edoardopaolini/ProjectNMS · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
function Z=enzReact(t,X,r)
% function Z=enzReact(t,X,r)
% AMOUNT
% fbna(1), fm(2), thb(3), fbni(4), c0(5), c1(6), c2(7)
% RATES
% k+(1), k-(2), k(3), k1+(4), k1-(5), k2+(6),
% k2-(7), k3+(8), k3-(9), k4+(10), k4-(11)
    Z(1) = -r(4)*X(1)*X(3) + r(5)*X(6) - r(8)*X(1)*X(6) + r(9)*X(7) + r(3)*X(5);
    Z(2) = r(6)*X(6) - r(7)*X(3)*X(2) + r(10)*X(7) - r(11)*X(6)*X(2);
    Z(3) = -r(4)*X(1)*X(3) + r(5)*X(6) + r(6)*X(6) - r(7)*X(3)*X(2) ...
            - r(1)*X(4)*X(3) + r(2)*X(5) + r(3)*X(5);
    Z(4) = -r(1)*X(4)*X(3) + r(2)*X(5);
    Z(5) = r(1)*X(4)*X(3) -  r(2)*X(5) - r(3)*X(5);
    Z(6) = r(4)*X(1)*X(3) - r(5)*X(6) - r(6)*X(6) + r(7)*X(3)*X(2) ...
            + r(9)*X(7) - r(8)*X(1)*X(6) + r(10)*X(7) - r(11)*X(6)*X(2);
    Z(7) = r(8)*X(1)*X(6) - r(9)*X(7) + r(11)*X(6)*X(2) - r(10)*X(7);
    Z = Z';

end