function p = s104021227()
  format long;
  d = 3;
  n=d;
  z=0;
  
  x1=[0.3, 0.4, 0.5, 0.6];
  fx=[0.5488, 0.4493, 0.3679, 0.3012];
  x = x1-fx;
  
  for k=1:n+1
    L(k)=1;
    for i=1:n+1
      if(i~=k)
        L(k) = L(k)*(z-x(i))/(x(k)-x(i));
      end
    end
  end
  
  p=0;
  for k=1:n+1
    p = p+fx(k)*L(k);
  end
end