matlab求已知数据点的拉格朗日插值多项式

function f = Language（xyx0）
syms t;
if（length（x） == length（y））
    n = length（x）;    
else
    disp（‘x和y的维数不相等！‘）;
    return;
end                                      %检错

f = 0.0;
for（i = 1:n）
    l = y（i）; 
    for（j = 1:i-1）
        l = l*（t-x（j））/（x（i）-x（j））;      
    end;
    for（j = i+1:n）
        l = l*（t-x（j））/（x（i）-x（j））;      %计算拉格朗日基函数