数值分析李庆杨第五版第七章习题

7

用牛顿法

// An highlighted block
 syms x
 f(x) = x^3-3*x-1;;
 g = diff(f);
 %h = diff(g);
 x = zeros(10,1);
 x(1) = pi/2;
 i = 1;
 while 1
     x(i+1) =-f(x(i))/feval(g,x(i))+x(i);
     i= i+1;
     if(abs(x(i)-1.895494)<0.5*1e-6)
         break
     end
 end
 x

用xuan截法

// An highlighted block
 syms x
  clear all
 format long
 syms x
 f = @(x) x^3-3*x-1;
%g = @(x) 3*x^2-3;
 x = zeros(10,1);
 x(1) = 2;
 x(2) = 1.9;
 i = 2;
 while abs(x(i)-1.8794)>=0.5*1e-4
     g = (f(x(i))-f(x(i-