function [ dydt ] = dYdT2015( t, y, k )
% compute derivates describing changes of chemical species
% in simply-branched MCA model

nu = ComputeNu2015(k, y);
dydt = zeros(1,4);
dydt(1) = nu(1)-nu(2);
dydt(2) = nu(2)-nu(3)-nu(4);
dydt(3) = nu(3)-nu(5);
dydt(4) = nu(4)-nu(6);

dydt = dydt';

end

