画拉格朗日插值的图像

lzhhun

小弟初学此软件，编了个程序用来画出拉格朗日插值的图像，可是老是出现错误，在这里跪求大侠们的指教了～～～
x=-1:1:1;
for i=1:3
    x0(i)=-1+1*(i-1);
    y0=1/(1+25*x0^2);
    s=0;
    for k=1:3
        L=1;
        for j=1:3
            if j~=k
              L=L*(x-x0(j))/(x0(k)-x0(j));
            end
        end
        s=s+L*y0(k);
    end
    y(i)=s;
end
plot(x,y)

[ 本帖最后由 ChaChing 于 2009-4-18 20:42 编辑 ]

ChaChing

当loop i=1时, x0仅为1*1, loop j怎样赋值至3?
还有建议楼主看下本版规则! 标题要明确!

[ 本帖最后由 ChaChing 于 2009-4-18 20:41 编辑 ]