matlab拉格朗日差值

function yi=lagrange(x,y,xi)%Lagrange插值%x,y为已知节点及其函数值向量%xi为插值点（可以是多个）,yi为插值n=length(x);m=length(xi);%length是x的长度,for i=1:mz=xi(i);s=0;for k=1:n%for循环语句p=1。0;for j=1:nif(j~=k)p=p*(z-x(j))/(x(k)-x(j));endends=p*y(k)+s;endyi(i)=s;end------------------------------------------------------x=[0。04 0。55 0。65 0。80 0。90 1。05];y=[0。41075 0。57815 0。69675 0。88811 1。02652 1。25382];x0=[0:0。01:2];y0=lagrange(x,y,0。596); 再问： f(0。596)的近似值可以求解，四次差值多项式怎么求的？谢谢~ 再答： 如果是让matlab显示出完整的展开式，本人表示实在不晓得。。。