function maxvirhe=poly1virhe(n)
% polyvirhe   Virheapproksimaatio lineaarisessa interpolaatiossa
% poly1virhe(n) laskee lineaarisessa interpolaatiossa välillä [0,1]
% n:ssä tasavälisesti valitussa otospisteessä lasketun 
% maksimivirheen.
% Interpoloitava funktio annetaan alifunktiona f

% Muokattu Kirjan [Hig-Hig] poly1err-esimerkistä s. 148
%     19.8.2010  HA

f0=f(0);f1=f(1);
x=linspace(0,1,n);
p=f0*(1-x)+f1*x;
plot(x,f(x),x,p,x,p,'o')
maxvirhe=max(abs(f(x)-p));

% Alifunktio:
function y=f(x)
y=x.*sin(x);
