clc;
clf;
clear;
T = 2;
t =linspace(0,T,1000);
y = funcion_rectangular(t);


set(gca,"fontsize", 20);

f1 = @funcion_rectangular;

w = 2*pi/T;

q = quad (f1, 0, T);%Integral definida de f de [0 a T]
a0 = (1/T)*q;% Cáculo del coeficiente a0 de la serie de fourier

ncoeficientes = 10;
An = zeros(ncoeficientes,1);
Bn = zeros(ncoeficientes,1);
for n=1:ncoeficientes
    fcos = @(t) f1(t).*cos(n*w*t);%Función anónima que multiplica f*cos(nwt)
    fsen = @(t) f1(t).*sin(n*w*t);%Función anónima que multiplica f*sin(nwt)
            
    q2 = quad(fcos,0,T);% Calcula la integral
    An(n) = (2/T)*q2;% Calcula el coeficiente an

    q3 = quad(fsen,0,T);
    Bn(n) = (2/T)*q3;% Calcula el coeficiente bn
endfor    

clf;
subplot(2,1,1);
plot(t,y,"r","linewidth",2);
grid on;
xlabel("Eje X","fontsize",16);
ylabel("Eje Y","fontsize",16);
title("Señal Original","fontsize",20);

nterminos = ncoeficientes;
t2 = linspace(-2*T,2*T,1000);
fresultante = crea_serie_Fourier(a0,An,Bn,t2,nterminos,T);

subplot(2,1,2);
plot(t2,fresultante,"r","linewidth",2);
grid on;
xlabel("Eje X","fontsize",16);
ylabel("Eje Y","fontsize",16);
title(["Serie de fourier de la señal original con " num2str(nterminos) " términos"],"fontsize",20);