Octave - Calcular raiz pelo método de Newton
Publicado por Daniel Moreira dos Santos (última atualização em 15/07/2010)
[ Hits: 23.180 ]
Homepage: http://www.danielmoreira.wordpress.com
Em análise numérica, o método de Newton (ou método de Newton-Raphson) tem o objetivo de estimar as raízes de uma função. Para isso, toma-se um ponto qualquer da função, calcula-se a equação da tangente (derivada) da função nesse ponto, calcula-se o intercepto da tangente ao eixo das abcissas, calcula-se o valor da função nesse ponto, e repete-se o processo, que deve tender a uma das raízes da função rapidamente, ou não tender a nada, deixando isso claro logo.
function [Raiz,Iter,Erro] = Newton_Raiz(funcao,x0,Toler,IterMax);
%parametros de entrada: funcao,x0,Toler,IterMax -> funcao, valor inicial,
%tolerancia e num. max. de iteracoes
%parametros de saida: Raiz,Iter,Erro-> raiz, num. gasto de iteracoes e
%cond. de erro
x = x0;
Fx = eval(funcao);
DFx = eval(diff(funcao));
Iter = 0;
disp(' Iter x DFx Fx DeltaX');
disp([Iter x DFx Fx ]);
while (1)
DeltaX = -Fx/DFx;
x = x+DeltaX;
Fx = eval(funcao);
DFx = eval(diff(funcao));
Iter = Iter+1;
disp([Iter x DFx Fx DeltaX]);
if(abs(DeltaX)<Toler && abs(Fx)<Toler)||abs(DFx)==0 || Iter>= IterMax
break;
end
end
Raiz = x;
if abs(Fx)<Toler
Erro = 0;
else
Erro =1;
end
y = [];
for x=0.001:0.001:10
y = [y,eval(funcao)];
end
x = 0.001:0.001:10;
plot(x,y);
end
Octave - Calcular raiz pelo método Regula-Falsi
Cálculo da raiz de uma função pelo Método da Bisseção - Octave
Função Fatorial Recursiva em Fortran 95
Nenhum comentário foi encontrado.
KDE Plasma - porque pode ser a melhor opção de interface gráfica
Gentoo: detectando impressoras de rede e como fixar uma impressora por IP
Como o GNOME conseguiu o feito de ser preterido por outras interfaces gráficas
Por que sua empresa precisa de uma PKI (e como automatizar EMISSÕES de certificados via Web API)
Instalando NoMachine no Gentoo com Systemd (acesso Remoto em LAN)
Instalar Linux em notebook Sony Vaio VPCEG13EB (13)
Vou destruir sua infância:) (7)
Interface gráfica com problema (0)









