format long
x=[0.3 0.4 0.5 0.6];
e=[0.5488 0.4493 0.3679 0.3012];
y=x-e;
z=0;
for i=1:1:4
  l(i)=1;
  for j=1:1:4
    if(j!=i)
    l(i)=l(i)*(z-y(j))/(y(i)-y(j));
    end
    end
    end
    f=0;
    for k=1:1:4
      f=f+l(k)*x(k);
      end
      printf("%d",f);