函数YY=la grange (x1，y1，xx ) %该程序是Lagrange1插值，其中x1，y1是插值节点和节点上的函数值，输出是插值点xx的函数值，%xx是向量symsxn=Length(x1； for i=1:nt=x1； t(I )=[]； L(I )=Prod () x-t )./(x1 ) I-t )； % L向量用于存储插值基函数endu=sum(l.*y1 )。 p=simplify(u ) % p是简化的Lagrange插值函数)字符串) yy=subs(p ) p，x，xx )； %p是以x为自变量的函数，求出xx中的函数值end