function y=chi(a,b,t);
% Programa para el cálculo de la función
% característica del intervalo [a,b].
% (La variable temporal es t; a y b son los 
% valores de los extremos del intervalo)
%
y = ones(size(t));
for k=1:1:length(t)
    if t(k)<a
        y(k)=0;
    end
    if b<t(k)
        y(k)=0;
    end
end