Método newton_raphson

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

function r = newton_raphson()
clc;
format long;
syms x; %%% variavel simbolica
f = input('Digite a funcao desejada: (em funcao de x) ');
f2 = input('Digite a derivada da funcao');
r(1) = input('Digite o valor de da raiz aproximada ');
e = input('Digite a precisao desejada ');
k=1; %% iteracoes
fa=(subs(f,r(1))); %% f(a)
f2a=(subs(f2,r(1))); %% f(b)
r(2)=r(1)- fa/f2a;
erro = r(2)-r(1);
while(abs(erro) > e) %% enquanto o erro for maior que o valor estabelecido
 r(k+2)=r(k+1)-(subs(f,r(k+1)))/(subs(f2,r(k+1)));
 k=k+1;
 erro = r(k+1)-r(k);
end
disp('O numero de iteracoes e: ')
k
disp('A raiz e: ')
r(k)
disp('O erro e')
r(k+1) - r(k)