function f = funcion_rectangular (x)
    [n,m]=size(x);
    f = zeros(n,m);
    intervalo_1 = find(x>=0 & x< pi);
    intervalo_2 = find(x>pi & x<2*pi);
    f(intervalo_1)=1;
    f(intervalo_2)=-1;
endfunction