// Calculate the time
clc;
th=30;
th0=50;
tc=120;
t=-120*[log(1-(th/th0))];
disp(t,'time(s)=')