f=diffreact(x)

% Author: Housam Binous

% Finding best reactor set-up

% National Institute of Applied Sciences and Technology, Tunis, TUNISIA

% Email: [e-mail]

% This problem was solved graphically by O. Levenspiel in Chemical Reaction
% Engineering, 3rd Edition, Wiley, 1999.

function f=diffreact(x)

global X

u=1:0.1:10;
Y=spline(X(:,1),X(:,2),u);

Z=diff(Y)./diff(u);

f=spline(u(1:end-1),Z,x)+(10-feval(@InvReact,x))/(10-x);

end