FA0 = 5; // mol/s rAat=-(1/400); X = [0 0.1 0.2 0.3 0.4 0.5 0.6]'; p = [189 192 200 222 250 303 400];//1/-rA = 800//dm^3.s/mols