MATLAB与拉格朗日插值

                    实验一函数插值方法
  一、问题提出
    nge公式求其插值多项式或分段二次Lagrange插值多项式。
    数据如下
  (10.65 04 0.55 0.8 0.9 1.05
                                  O 5
                          0.6967 0.9 1.0 0.4107 0.5781 Y 1.2538
            5 5 5 2
                    (x) 求五次Lagrange多项式 和分段三次插值多项式,计算 f(0.596)f(0.99)
   的值。(提示:结果为
                  f(0.5960625732f(0.99)1.05423
  21
                      3 2 A 50.135 0.05 0.01 0.00 0.36 0.000.00
                           8 7 8 21
    试构造Lagrange多项式 L.(x) 计算的
     f(1.80.164762f(6.150.001266
二、要求
   1、利用Lagrange插值公式
          L(x)= X-x
                0 编写出插值多项式程序
   2 给出插值多项式或分段三次插值多项式的表达式:
   3 根据节点选取原则,对问题(2)用三点插值或二点插值,其结果如何:
   4、对此插值问题用Newton插值多项式其结果如何。Newton插值多项式如下
                                     (x-X;)
                                 O.jk
                              f(x)
     其中: 107
三、自的和意义
   1、学会常用的插值方法,求函数的近似表达式,以解决其它实际问题:
   2、 明确插值多项式和分段插值多项式各自的优缺点:
   3. 熟悉插值方法的程序编制:
   4、如果绘出插值函数的曲线,观察其光滑性。
四、实验学时:2学时
五、实验步骤:
 1.进入C或matlab开发环境
 2.根据实验内容和要求编写程序:

function yy=lagrange(x1,y1,xx);
x1=(0.4,0.55,0.65,0.8,0.9,1.05);
y1=(0.41075,0.57815,0.69675,0.90,1.00,1.25382);
xx=(0.596,0.99);
n=length(x1)'
l=ones(1,n);
for k=1:n
for j=1:n
if j~=k
l(k)=l(k).*(xx-x1(j))./(x1(k)-x1(j));
end
end
end
sums(y1.*l);
这是我的代码,但是运行不出来