给我一段关于多项式插值的matlab代码

以下是关于多项式插值的 MATLAB 代码示例：

假设我们需要对函数f(x)=sin(x)在区间[0,pi]上进行5次多项式插值，代码如下：

% 定义插值节点和函数值
n = 5; % 插值节点数
x = linspace(0,pi,n); % 插值节点
y = sin(x); % 插值节点处函数值

% 构造插值多项式
syms t; % 定义符号变量
L = 0; % 初始化插值多项式
for i = 1:n
    % 计算第i个拉格朗日基函数
    l = 1;
    for j = 1:n
        if j ~= i
            l = l * (t - x(j)) / (x(i) - x(j));
        end
    end
    % 将第i个基函数乘上对应的函数值并累加到插值多项式中
    L = L + y(i) * l;
end

% 绘制原函数和插值多项式的图像
f = @(t) sin(t);
t = linspace(0,pi,1000);
plot(t,f(t),'b-',t,eval(L),'r--',x,y,'ko');
legend('原函数','插值多项式','插值节点');

运行上述代码，会输出原函数和插值多项式的图像，其中蓝色曲线是原函数，红色虚线是插值多项式，黑色圆点是插值节点。