myfnewton3.m

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

%myfNewton3([2 3 4])
function y=myfNewton3(X)
syms x;
if length(X)==3;
    Y(1)=X(1)^3;
    Y(2)=X(2)^3;
    Y(3)=X(3)^3;
    A=Y(1);
    B=(Y(2)-Y(1))/(X(2)-X(1));
    chf=(Y(3)-Y(1))/(X(3)-X(1));
    C=(chf-B)/(X(3)-X(2));
    y=A+B*(x-X(1))+C*(x-X(1))*(x-X(2))
    y=simplify(y);
    ezplot(y,[0 6])
    hold on
    fplot('k^3',[0 6],'y.')
    hold off
else
    disp('Input error')
end