xdot=vandevusse(t,x)

% Construction of Attainable Region for van de Vusse Kinetics
% Author's Data: Housam BINOUS
% Department of Chemical Engineering
% National Institute of Applied Sciences and Technology
% Tunis, TUNISIA
% Email: [e-mail] 

function xdot=vandevusse(t,x)

k1=0.01;k2=5;k3=10;k4=100;

xdot(1)=-k1*x(1)+k2*x(2)-k4*x(1)^2;
xdot(2)=k1*x(1)-k2*x(2)-k3*x(2);

xdot=xdot';