function xdot=reac5(t,x,n,m,k)
xdot=[];
xdot(1)=1/(1+x(n)^m)-k(1)*x(1);
for i=2:n
  xdot(i)=x(i-1)-k(i)*x(i);
end
xdot(n)=x(n-1)-k(n)*x(n)/(1+x(n));
//