fsxnewton.m

来自「非线性规划中的最速下降法,黄金分割法,阻尼牛顿法,牛顿切线法.」· M 代码 · 共 13 行

function x=fsxnewton(f,e,a)
% fsxnewton1函数  牛顿切线算法;
% x=fsxnewton(f,e,a,b)为输入函数 f为函数 e为允许误差 a为初始点;

% fsx TJPU 2008.6.15

f1=diff(f);f2=subs(f1,a);f3=diff(f1);c=a;
while (abs(f2)>=e) 
	f2=subs(f1,c);
	f4=subs(f3,c);
	c=c-f2/f4;
	x=c;
	end;