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

function F = goo_k3_v2(x, k1, k2, k3, At, Bt)

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


F = [-k1 * A * B  - k3 * A * Bi  + k2 * Ai;
      k1 * A * B  + k3 * A * Bi  - k2 * Ai - k1 * Ai * B - k3 * Ai * Bi + k2 * Aii;
      k1 * Ai * B + k3 * Ai * Bi - k2 * Aii;
     -k1 * B * A  - k3 * B * Ai  + k2 * Bi;
      k1 * B * A  + k3 * B * Ai  - k2 * Bi - k1 * Bi * A - k3 * Bi * Ai + k2 * Bii;
      k1 * A * Bi + k3 * Ai * Bi - k2 * Bii;
      At - x(1) - x(2) - x(3);
      Bt - x(4) - x(5) - x(6)];
 
 
 

