
clear all
clc

syms FV Cain Ca k1 k2 k3 Cb


res(1)=  FV*(Cain-Ca) - k1*Ca - k3*Ca^2;
res(2)=  FV*(-Cb) + k1*Ca - k2*Cb;


A(1,1) = diff(res(1),Ca);
A(1,2) = diff(res(1),Cb);
A(2,1) = diff(res(2),Ca);
A(2,2) = diff(res(2),Cb);
