newton.asv

为配套王能超版本数值分析简明教程编写,方便学生教师演示!

function newton(x,y,e,f)
syms t
f1=inline(f);
xx=linspace(x,y,100);
f2=diff(f1(t));
f5=inline(f2);
yy=f5(xx);
if yy~=zeros(1,100);
    syms t2
    f3=inline(t2-f1(t2)/f5(t2));
    %f4=inline(f3);
    syms t1
    y1=diff(f3(t1));
    f2=inline(y1);
    xx1=linspace(x,y,100);
    s=abs(f2(xx1));
    s1=zeros(1,100);
    s2=ones(1,100);
    if s>s1 & s<s2;
        i=1;
        while abs(x-y)>j;
            x=y;
            y=f3(x);
            i=i+1;
            fprintf('%d%13.8f\n',i,y);
        end
    else 
        disp('此多项式不能进行迭代');
    end
end
    
 %   function xia329(a,b,j,f3)
%f=inline(f3);
%syms s1
%y=diff(f(x));
%f2=inline(y);
%xx=linspace(a,b,100);
%    s=abs(f2(xx));
 %   s1=zeros(1,100);
  %  s2=ones(1,100);
   % if s>s1 & s<s2;
    %    i=1;
     %   while abs(a-b)>j;
      %      a=b;
       %     b=f(a);
        %    i=i+1;
         %   fprintf('%d%13.8f\n',i,b);
         %end
         % else 
        %disp('此多项式不能进行迭代');
        %end