% tehtävä 4
% tutkitaan samaa funktiota f kuin ennen: nyt interpoloivilla polynomeilla
function y = teht4(n)
xs = linspace(0,1);

xsd = 0:1/n:1; % tasainen jako välein 1/n
ys = f(xsd);
P = polyfit(xsd,ys,n);

% tarkastetaan approksimaatio plottaamalla
plot(xs, polyval(P,xs));
hold on;
plot(xs, f(xs), 'black');
legend('polyfit', 'f');
hold off;

% ääretönnormi saadaan virhe.m:n avulla
vi = virhe(P, 100);

y = vi;

% tehtävän funktio f
function y = f(x)
y = exp(sin(7*x));