function [dcdt1]= conzt1(t,conzactual1)
% CONZT1 returns the concentration change of substances from time t1 onwards 
% [dcdt1]= conzt1(t,conzactual) using the global KPLUSEINS, KMINUSEINS,
% KEXPRPROT, KDEGREINS
global KPLUSEINS;
global KMINUSEINS;
global KEXPRPROT;
global KDEGREINS;
dcdt1 = zeros(3,1);
dcdt1(1) = - KPLUSEINS*conzactual1(1)*conzactual1(2) + KMINUSEINS*conzactual1(3);
dcdt1(2) = KEXPRPROT*1 - KPLUSEINS*conzactual1(1)*conzactual1(2) + KMINUSEINS*conzactual1(3) - KDEGREINS*conzactual1(2);
dcdt1(3) = - KMINUSEINS*conzactual1(3) + KPLUSEINS*conzactual1(1)*conzactual1(2);
