phi＝c0v.*exp(-d*t);%exponential decay for each component