niudun.m

来自「本文件是基于MATLAB的解最优化算法的源代码」· M 代码 · 共 18 行

function y=newton(f,x0,eng)
%newton erji
syms x;
f=x^4-4*x^3-6*x^2-16*x+4;
x0=6;
k=1;
eng=0.001;
f1=diff(f)
f2=diff(f,2)
while(abs(subs(f1,x0)) >eng)
x1=x0-subs(f1,x0)/subs(f2,x0);
x0=x1;
k=k+1;
end

zuiyouzhi=x0
zuiyoujie=subs(f,x0)
end