close all, clear all

d = domain(-1,1); s = linspace(-1,1,7); y = [0 0 1 0 0 0 0];
p = interp1(s,y,d);
plot(p), hold on, plot(s,p(s),'.k'), grid on, FS = 'fontsize';
title('Lagrange polynomial  l_5  on 7-point equispaced grid',FS,9)
