TD=10^(-5);
t=10^(-8);
ct=10^5;
Q=linspace(0,6000,200);
cp=(TD+Q*t)*ct./(log(1-((TD+Q*t)*ct./log2(Q))));
cpp=abs(cp);
plot(Q,cpp);

hold on;
TD=10^(-5);
t=10^(-8);
ct=2*10^5;
Q=linspace(0,5000,200);
cp=(TD+Q*t)*ct./(log(1-((TD+Q*t)*ct./log2(Q))));
cpp=abs(cp)
plot(Q,cpp);