% Script per la valutazione di punti di interpolazione

a = 1;
b = 6;
n = 10;
m = 20;
h = 5/11;
x = [a:(b-a)/n:b]';
y = sqrt(x);
xs = [a:(b-a)/m:b]';

for k = 1:m+1;
	ys(k) = pollag(x,y,xs(k),n);
end;

plot(x,y,'*-',xs,ys,'-')