mfunc1.m

用有理多项式做Chebyshev逼近的一种直接方法的matlab程序

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