clc

global Tlimit

% guess the initialization
t=cputime;
firstvalue = 75;
ntimes = 150
t0=linspace(1, firstvalue, ntimes);

arr = zeros(1,ntimes);
for i=1:ntimes    
    arr(i)=optim(t0(i));
    disp('---Be patient, it is computing!---')
end

index=find(arr==min(arr))
Tlimit=t0(index)

ComputingTime=cputime-t