function dx=p426singlefun(t,x)
Vi=200;
v0=0.05;
ca0=7.72;
cb0=10.93;
na0=ca0*Vi;
fb0=0.05*cb0;
k=9e-5;
kc=1.08;
V=Vi+v0*t;
xa=x(1);
na=na0*(1-xa);
nb=cb0*v0*t-na0*xa;
nc=na0*xa;
nd=na0*xa;
ca=na/V;
cb=nb/V;
cc=nc/V;
cd=nd/V;
r=k*(ca*cb-(cc*cd)/kc);
ra=r*V/na0;
dx = [ra];

