Newton Rapshion

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function r=nrh(s,f,E1,E2,n)
 x=sym('x');
 format long
 if(abs(subs(f,x,s)))<E2
 r=double(s);
 fprintf('\nraiz = %d\n',r)
 return
 end
 x1=s;
 for k=1:n
 F=diff(f);
 f1=subs(f,x,x1);
 x2=x1-(f1/F);
 if(abs(x2-x1))<E1
 r=double(x2);
 fprintf('\nraiz = %d\n',r)
 return
 end
 x1=x2;
 end
end