clc; clear all; close all;
%kp = 1:1:100;
%kn = -100:1:-1

cont = 1;
controle1 = 1;

for c = 1:1/10:10
    for kp = 1:1:100
        %yp(cont,controle1) = ((i*(-1)^kp)/(kp*pi))*(exp(i*kp*pi*c));
        yp(cont,controle1) = ((1/(kp*pi))*sin((2*kp*pi)/3))*(exp(i*kp*(pi/3)*c));
        cont = cont+1;
    end
    controle1 = controle1+1;
end

cont1 = 1;
controle2 = 1;

for d = 1:1/10:10
    for kn = -100:1:-1
        yn(cont1,controle2) = ((1/(kn*pi))*sin((2*kn*pi)/3))*(exp(i*kn*(pi/3)*d));
        cont1 = cont1+1;
    end
    controle2 = controle2+1;
end

r1 =  sum(yn);
r2 = sum(yp);
rst = r1+r2+4/6;
plot(real(rst));
grid;