function dy = pcoumRates_uppsala(t,y)
  % ODE describing the p-coumaric acid production
  global k kcat Km AroP Vcell Vout Vtot;
  DNA=y(1);
  mRNA=y(2);
  TAL=y(3);
  TyrOut=y(4);
  TyrCell=y(5);
  pcoum=y(6);
  
  dy=zeros(6,1);
  dy(1)=0;
  dy(2)=k(1)*DNA-k(2)*mRNA;
  dy(3)=k(3)*mRNA-k(4)*TAL;
  dy(4)=-(Vcell/Vout)*AroP*kcat(1)*TyrOut/(Km(1)+TyrOut);
  dy(5)=AroP*kcat(1)*TyrOut/(Km(1)+TyrOut)-TAL*kcat(2)*TyrCell/(Km(2)+TyrCell);
  dy(6)=(Vcell/Vtot)*TAL*kcat(2)*TyrCell/(Km(2)+TyrCell); 

   disp(t/(60*60));
end

