function g=reaction(t,x)

% reaction - wertet die vektorwertige Funktion g(x,t)=lambda*x 
%            mit globalem Parameter lambda und den rellen Argumenten x,t aus 
%        
% Benutzung: g(x,t)
%
% Argumente: x - Zustandsvariable
%            t - Zeit
%
% Parameter: Reaktionskonstanten c1, c2, c3
%

c1=0.04;
c2=3*10^7;
c3=10^4;


g(1)= -c1*x(1) + c3*x(2)*x(3);

g(2)=  c1*x(1) - c3*x(2)*x(3) - c2*x(2)^2;

g(3)=                           c2*x(2)^2;

g=g';