r.m

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

function r=r(t,e,x,l,m)
%e=a;
%e(l-m+2)=c(2);
%e(l-m+3:l+m+2)=c(3:2*m)-c(1:2*m-2);
r=mfunc2(t,e(1:l-m+1),x(1:l-m));
if m>0
    f=e(l+m+1);
    if m>1
        for i=l+m:-1:l-m+3
            f=e(i)+(t-x(i))./f;
        end
    end
    tem=e(l-m+2);
    for i=1:l-m+1
        tem=tem.*(t-x(i));
    end
    f=tem./(1+(t-x(l-m+2))./f);
    r=r+f;
end