function dmdt = reakce7 (t,m)

A=[3.27e6 1.08e10 4.38e9 1e5 4.28e6];
E=[111.7 148 152.7 108 108].*1000;
T=480+273;

R=8.314;

k=A.*exp(-E./(R*T));

dmdt=zeros(4,1);


dmdt(1)=m(1)*(-k(1)-k(2)-k(3)); %biomass
dmdt(2)=m(1)*k(1)+m(3)*k(4); %char
dmdt(3)=m(1)*k(2)-m(3)*k(4)-m(3)*k(5); %liqud
dmdt(4)=m(1)*k(3)+m(3)*k(5); %gas

end