function F=cap(t)   	

a = pi*0.5 ; b = pi*1.5;	% estremi dell'intervallo
h = 1.0				% altezza

% funzione caratteristica
chia=t >= a; chib = t <= b;
F = h*(chia.*chib);

% funzione cappello
%F = h - abs ( (2*h/(b-a))*(t-pi) );
%F = F.*chia.*chib;




