使用Scilab的我的前向差异表中的问题在哪里？

Question

所以，我必须像这样插入函数f（x）。x = 0：0.1：2.8;y = [0 0.717 0.999 0.675 0.0583 0.7568 0.9961 0.6312];

这是我现在得到的代码。

clc
clear

x = 0:0.1:2.8;
y = [0 0.717 0.999 0.675 0.0583 0.7568 0.9961 0.6312];

n = length(x);
del = %nan * ones (n ,7) ;
    del (:,1) = y';
    for j = 2:7
        for i = 1: n - j +1
            del (i,j) = del(i+1,j-1) - del(i,j-1);
        end
    end
del = [x'del];
del = round ( del *10^3) /10^3;
mprintf ("%5s,%7s,%8s,%9s,%8s,%8s,%8s",'x','y','dy','d2y','d3y','d4y','d5y')
disp ( del )

这给我Submatrix错误定义的错误。问题可能出在哪里？

Answer 1

x和y应该具有相同的长度，但您的数据却并非如此。例如，您可以设置

y = [0 0.717 0.999 0.675 0.0583 0.7568 0.9961 0.6312];
x = linspace(0,2.8,length(y));

[C0行失败，应写为（缺少空格）]

del = [x'del];

然后您脚本输出结果：

del = [x' del];