% Script per l'esercizio 3
%

clc

riga = '--------------------------------------------';

disp(riga)

% intestazione tabella
% k: indice del punto (contatore)
% x: variabile
% g(x): valore calcolato della funzione

fprintf('\n  k \t   x \t\t  g(x) \n') 

disp(riga);

% scelta del valore n
n = 100;

N = 2*n + 1;

figure(1)
clf
figure(2)
clf

X = linspace(-2*pi,2*pi,N);

for k = 1 : N

    xk = X(k);
    
    yk = g_es3(xk);
    
    if rem(k-1, 10) == 0

        fprintf('%2d\t %+6.4f\t %+6.4f\n', k, xk, yk);
    
    end
    
    figure(1)
    hold on
    plot(xk,yk,'r*')
    
    Y(k) = yk;
end

disp(riga)

figure(2)
hold on
plot(X,Y,'b-')
plot([-2*pi 2*pi],[0 0],'k-')
plot([0 0],[-1 1],'k-')
axis tight