FIC MATLAB 2017-1/Fourier

clear all; close all; clc

% Variando o limite máximo de termos da transformada
for n = [5, 11, 25, 51, 99, 10001]
    
    % Cálculo da função
    ii = (1:2:n)';
    x = linspace(-pi,pi,1000);
    F = (1./ii')*sin(ii*x);
    
    % Plot dos resultados
    plot(x,F)
    title(['F(x) = sen(x) + sen(3x)/3 + ... + sen(' num2str(n)  'x)/' num2str(n)])
    xlabel('x')
    ylabel('F(x)')
    grid on
    
    pause
    
end