function y = GTT(t,p)
% GTT function
y  = p(1)+p(2)*exp(-p(3)*t).*cos(p(4)*(t-p(5)));
end

