给出matlab的拉格朗日插值实现代码

以下是Matlab的拉格朗日插值实现代码：

function p = lagrange(x,y,z)
% x,y为n个插值节点，z为n-1阶插值多项式
n = length(x);
p = 0;
for i = 1:n
    L = 1;
    for j = 1:n
        if j~=i
            L = L*(z-x(j))/(x(i)-x(j));
        end
    end
    p = p + y(i)*L;
end
end

其中，x为插值节点的横坐标，y为插值节点的纵坐标，z为插值点的横坐标，p为通过拉格朗日插值算法计算出的插值点的纵坐标。