function hlt=hl(t)
global tc
c1=(hu(tc)+hu(0))/2;
c2=(hu(tc)-hu(0))/2;
if t<tc
	hlt=hu(t);
else
	hlt=c1+c2*cos(pi*(t-tc)./(1-tc));
end
