D/N=D0/N+P0/N(1-exp(t/T)) D=daughter P=parent N=reference stable isotope of D T=time constant