function dc = kaksistep(t,c)
global kf1 kf2 kr1 kr2
A=[-kf1 kr1 0;kf1 -(kr1+kf2) kr2;0 kf2 -kr2];
dc=A*c;
end