给我用matlab写出拉格朗日插值多项式余项代码

以下是用 Matlab 编写的拉格朗日插值多项式余项代码：

function [Rn] = lagrange_remainder(x, y, xn)
% x: 插值点的横坐标
% y: 插值点的纵坐标
% xn: 求余项的点的横坐标
n = length(x) - 1;
syms t;
L = ones(n+1, 1);
for i = 1 : n+1
    for j = 1 : n+1
        if i ~= j
            L(i) = L(i) * (t-x(j)) / (x(i)-x(j));
        end
    end
end
f = 1;
for i = 1 : n+1
    f = f * (t-x(i));
end
f = diff(f, n+1);
Rn = f * L(n+1) * prod(xn-x) / factorial(n+1);

其中，`x` 和 `y` 分别表示插值点的横坐标和纵坐标，`xn` 表示求余项的点的横坐标。函数返回的 `Rn` 即为余项。注意，这里使用了符号计算工具箱中的 `syms` 和 `diff` 函数，需要先安装并导入该工具箱。