myfnewton.m

来自「自己编写的matlab程序 数值分析里的程序」· M 代码 · 共 25 行

%myfNewton([3. 4. 5. 6.])
function y=myfNewton(X)
syms x;
n=length(X);
for i=1:n
    Y(1,i)=log(X(i));
end
for i=2:n
    for j=i:n
        Y(i,j)=(Y(i-1,j)-Y(i-1,j-1))/(X(j)-X(j-i+1));
    end
end
y=Y(1,1);
t=1;
for i=2:n
    t=t*(x-X(i-1));
    y=y+Y(i,i)*t;
end
y=simple(y);
% y=double(y);
y=vpa(y,6);%将分数形式转换成小数形式(含有符号变量的情况用vpa)
ezplot(y,[1 6])
hold on
fplot('log(k)',[1 6],'M')
hold off