function f = autokat(t,y,k)

f=zeros(3,1);

f(1) = -k(1)*y(1) + k(2)*y(2)*y(3);
f(2) = k(1)*y(1) - k(2)*y(2)*y(3) - k(3)*y(2)*y(2);
f(3) = k(3)*y(2)*y(2);