clc
clear

MW=28;   % molecular weight (kg/kmol) for water
MWA=236; % molecular weight (kg/kmol) for Carbamazepine
ro=999;  %density of water (kg/m3)

FF=5;  % initial flow rate (mL/min)
FF=FF*(1e-6/60)   % (m3/sec)
Ft=FF*ro/MW*1000; % initial molar flow rate (mol/sec) 

CA0=9.32;   % initial concentration of pollutant (micromol/L)
CA0=CA0*(1e-6*1000); %  (mol/m3)

FA0=FF*CA0;       %initial molar flow rate of pollutant (mol/sec)

F0=[FA0 283];
[z,F]=ode15s(@reactor,[0:0.001:0.084],F0);
L=z./0.084;

X=(FA0-F(:,1))/FA0*100

plot(L,X)