f=cstr(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 f=cstr(x)

global tau

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

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