牛顿迭代法解非线性方程程序

2支棒棒糖

1. function x=newton(f_name,x0,tol)
   
2. % Purpose: Solve a nonlinear equation by Newton iteration
   
3. % Synopsis:     x=newton('f_name',x0)      
   
4. %               x=newton('f_name',x0,tol)  
   
5. %
   
6. %            x: return a root of f(x)=0 near x0
   
7. %            f_name:name of the function f(x) that defines the equation
   
8. %            x0: initial guess
   
9. %            tol: tolerence(Default:1e-6)
   
10. 
    
11. % L.J.Hu 1-8-1998
    
12. 
    
13. h=0.0001;M=500;
    
14. if nargin<3, tol=1e-6;end
    
15. x=x0;xb=x-999;
    
16. n=0;
    
17. while abs(x-xb)>tol
    
18.   xb=x;
    
19.   if n>M break;end
    
20.   y=feval(f_name, x);
    
21.   fprintf('  n=%3.0f,  x=%12.5e,  y=%12.5e, \n',n,x,y)
    
22.   y_driv=(feval(f_name,x+h)-y)/h;
    
23.   x=xb-y/y_driv;
    
24.   n=n+1;
    
25. end
    
26. fprintf('  n=%3.0f,  x=%12.5e,  y=%12.5e, ',n,x,y)
    
27. if n>500,
    
28. fprintf('\n');disp('Warning: iterations exceeds the limitation, probable devergent');
    
29. end

复制代码