clear all
paso=0.01;
tmax=10;
t=0:paso:tmax;
Nt=length(t);

f=1-heaviside(t-5);
f(501)=1;


w=-tmax:paso:tmax;
Nw=length(w);

F=zeros(1,Nw);
m=1;

while m<=(Nw)  
    wint=w(m);
    tint=0;
    n=0;
    area=0;
    while n<Nt-1
        area=area+((f(n+1).*exp(-i.*wint.*tint)+f(n+2).*exp(-i.*wint.*(tint+paso)))./2).*paso;
        n=n+1;
        tint=tint+paso;
    end
    F(m)=area;
    m=m+1;
end
F=F./(-w.^2+i.*w+1);
Fr=real(F);
Fi=imag(F);

figure(1)
plot(t,f,'r.-')
grid on;
title('Función en el dominio del tiempo')
xlabel('Eje del tiempo')
ylabel('Función f(t)')

figure(2)
plot(w,Fr,'r.-')
grid on;
title('Parte real de la transformada de Fourier')
xlabel('Eje de la frecuencia')
ylabel('Función Re F(w)')


figure(3)
plot(w,Fi,'r.-')
grid on;
title('Parte imaginaria de la transformada de Fourier')
xlabel('Eje de la frecuencia')
ylabel('Función Im F(w)')