function ndot=nofun(t,n,flag,tau,noinf)
%  Carol Lucas
ndot=(noinf-n)/tau;