Matlab插值曲线近似

Question

我被分配了一个任务，从一组7个点获得拉格朗日插值曲线，并在某个X处获得函数的值 第二部分是微不足道的，但没有人知道如何获得曲线本身的系数。 这是我使用的代码：

clear
clc
x=[0 6 10 13 17 20 28];
y=[6 17 42 37 30 29 27];
m=15;
n=length(x);
  for i=1:n
        l(i)=1;
        for j=1:n
            if i~=j
            l(i)=l(i)*(m-x(j))/(x(i)-x(j));
      end
      final(i)=l(i)*y(i);
    end
    answer=sum(final);
    end

我无法使用interp1进行分配。

有没有办法在matlab中这样做？我已经寻找了几个小时的解决方案，但似乎还没有人尝试过这个。

谢谢！