function [dcdt]= conzt2(t,conzactual1)
% CONZT1 returns the concentration change of substances from time t2 onwards 
% [dcdt]= conzt2(t,conzactual1) using the global KPLUSEINS, KMINUSEINS,
% KEXPRPROT, KPLUSZWEI, KMINUSZWEI, KATALYSE, KEXPRTOR, KTRANSTOR,
% KDEGREINS, KDEGRZWEI.
global KPLUSEINS;
global KMINUSEINS;
global KEXPRPROT;
global KPLUSZWEI;
global KMINUSZWEI;
global KATALYSE;
global KEXPRTOR;
global KTRANSTOR;
global KDEGREINS;
global KDEGRZWEI;
dcdt = zeros(6,1);
dcdt(1) = - KPLUSEINS*conzactual1(1)*conzactual1(2) + KMINUSEINS*conzactual1(3);
dcdt(2) = KEXPRPROT*1 - KPLUSEINS*conzactual1(1)*conzactual1(2) + KMINUSEINS*conzactual1(3) - KPLUSZWEI*conzactual1(2)*conzactual1(5) + KMINUSZWEI*conzactual1(4) + KATALYSE*conzactual1(4) - KDEGREINS*conzactual1(2);
dcdt(3) = - KMINUSEINS*conzactual1(3) + KPLUSEINS*conzactual1(1)*conzactual1(2);
dcdt(4) = KPLUSZWEI*conzactual1(2)*conzactual1(5) - KMINUSZWEI*conzactual1(4) - KATALYSE*conzactual1(4);
dcdt(5) = KEXPRTOR*1 - KPLUSZWEI*conzactual1(2)*conzactual1(5) + KMINUSZWEI*conzactual1(4) - KDEGRZWEI*conzactual1(5) - KTRANSTOR*conzactual1(5);
dcdt(6) = (double(uint32((KTRANSTOR*conzactual1(5) - KDEGRZWEI*conzactual1(6))*10)))/10;
