以下代码在我的数值方法教科书的文本中给出,我将其复制到Octave中,并在“ n = length(b);”之后得到解析错误,并且当我尝试使n为输入参数时,要查看这是否可以解决问题,解析错误仅移至下一行代码。这让我感觉好像错误实际上在代码中的其他地方一样?将不胜感激。谢谢!
function y = lsolve(L, b)
n = length(b);
for i = 1: n,
y(i) = b(i);
for j = 1: i - 1,
y(i) = y(i) - L(i, j) * y(j);
end;
end;