function ggprime=chem(t,gg,options,k1,k2,k3,k4)

A=gg(1);
B=gg(2);
C=gg(3);

dAdt=k1*A*B-k2*C;
dBdt=k2*C+k3*B^2;
dCdt=k4*A-k1*B*C;

ggprime=[dAdt; dBdt; dCdt];