function x=newton(s1,name)
   x=mat(s1,1);
   for j=3:s1+2,
      x(j)=1;
   end;
   for j=1:10,
      y=init_dif(x);
      xx=val(y);
      s=name+'(y)';
      zz=evstr(s);
      j=jacobi(zz);
      v=val(zz);
      d=j\v;
      nx=xx-d;
      for i=1:s1,
         x(i+2)=nx(i);
      end;
   end;
endfunction