mfunc1.m

来自「用有理多项式做Chebyshev逼近的一种直接方法的matlab程序」· M 代码 · 共 6 行

function f=mfunc1(x,f1)
n=length(x);
f(1:n)=f1(1:n);
for k =2:n
f(k:n)=(f(k:n)-f((k-1):(n-1)))./(x(k:n)-x(1:(n-k+1))); % calculate the difference quotients
end