function x=fsxdnewton(f,e,a,b)
% fsxdnewton函數 阻尼newton法
% x=fsxdnewton(f,e,a,b)為輸入函數 f為函數 e為允許誤差 (a,b)為初始點

% fsx TJPU 2008.6.16

syms x1 x2; 
v=[x1,x2]; 
df=jacobian(f,v); 
df=df.'; 
Q=jacobian(df,v);
x0=[a,b]';
g1=subs(df,{x1,x2},{x0(1,1),x0(2,1)}); %求梯度
Q1=subs(Q,{x1,x2},{x0(1,1),x0(2,1)});  %求hesee矩陣 
	while(norm(g1)>e) 
	p=-Q1\g1; 
	x0=x0+p; 
	g1=subs(df,{x1,x2},{x0(1,1),x0(2,1)}); 
  Q1=subs(Q,{x1,x2},{x0(1,1),x0(2,1)});
  end; 
x=x0;