% Samantha Sutton
% 18.385 final project
% This code stores the equations of the system, to be called by fsolve.

function dydt = ff_ODEs_v2(t, y, k1, k2, k3)

A   = y(1);
Ai  = y(2);
Aii = y(3);
B   = y(4);
Bi  = y(5);
Bii = y(6);

dydt(1,:)   =     -k1 * A * B  - k3 * A * Bi  + k2 * Ai;
dydt(2,:)   =      k1 * A * B  + k3 * A * Bi  - k2 * Ai - k1 * Ai * B - k3 * Ai * Bi + k2 * Aii;
dydt(3,:)   =      k1 * Ai * B + k3 * Ai * Bi - k2 * Aii;
dydt(4,:)   =     -k1 * B * A  - k3 * B * Ai  + k2 * Bi;
dydt(5,:)   =      k1 * B * A  + k3 * B * Ai  - k2 * Bi - k1 * Bi * A - k3 * Bi * Ai + k2 * Bii;
dydt(6,:)   =      k1 * A * Bi + k3 * Ai * Bi - k2 * Bii;


     
     
     