

clear all
close all
clc

% linspace(<valor_inicial> , <valor_final>, <número_de_elementos)
t =linspace(-2*pi, 2*pi,1000);

% Onda quadrada
x = square(t);

% Plota onda quadrada
subplot(2,2,1)
plot(t,x, 'k--')
axis([-10 10 -2 2])

% 1% harmonica

y = sin(t);

subplot(2,2,2)
plot(t,y ,'b',t,x,'k-')
axis([-10 10 -2 2])

% 1ª e 3 ª Harmonica

y= sin(t) + sin(3*t)/3
subplot(2,2,3)
plot(t,y ,'b',t,x,'k-')
axis([-10 10 -2 2])

% 1ª e 3 ª e 5ª Harmonica 

y= sin(t) + sin(3*t)/3 + sin(5*t)/5
subplot(2,2,4)
plot(t,y ,'b',t,x,'k-')
axis([-10 10 -2 2])

% 1ª e 3 ª e 5ª Harmonica 

sy=0;
loop = 10;
for n = 1:2:loop
    
    sy = sy + sin(n*t)/n
    
end
figure(2)
plot(t,sy)
    
    
    