x = linspace(0, 10, 7);
y = sin(x);

plot(x, y, 'o')

x2 = linspace(0, 10, 50);
y2 = zeros(1, 50);

y3 = sin(x2);

for i = 1: 50
   y2(i) = lagrange(x, y, x2(i)); 
end

hold on
plot(x2, y2, 'x')
hold on
plot(x2, y3, '*')



