%% Error de interpolacion

clear
clc
close all

y = @(x) sin(x).*cos(x);
vec = linspace(0,1,1000)';

for n = 2:5
    x = linspace(0,1,n)';
    pol = interpolante(x,y(x),vec); % interpolante por (x,y(x)) evaluada en vec
    plot(x,y(x),'ok', 'LineWidth', 2)
    hold on
    plot(vec, y(vec),'-k', 'LineWidth', 2);
    plot(vec,pol,'-r', 'LineWidth', 2)
    axis([0 1 0 0.52])
    title(['n = ', num2str(n)], 'FontSize',14);
    hold off 
    pause
    close all
end
%%

% Error?
error_pol = max(abs(y(vec)-pol))