clear all;
clc;
x=[5;6;9;11];
y=[12;13;14;16];
n=size(x);
a=input('a=');
s=0;
for i=1:n;
f=1;
for j=1:n
if (i~=j)
f=f*(a-x(j))/(x(i)-x(j));
end
end
s=s+f*y(i)
end