function?f?=?Language（xyx0）
syms?t;
if（length（x）?==?length（y））
????n?=?length（x）;????
else
????disp（‘x和y的維數(shù)不相等！‘）;
????return;
end??????????????????????????????????????%檢錯(cuò)

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））;??????%計(jì)算拉格朗日基函數(shù)