Octave - Calcular raiz pelo método de Newton
Publicado por Daniel Moreira dos Santos (última atualização em 15/07/2010)
[ Hits: 22.012 ]
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
Algoritmo para calcular a tabuada
Octave - Método de Gauss-Legendre
Octave - Calcular raiz pelo método Regula-Falsi
Nenhum comentário foi encontrado.
Melhorando o tempo de boot do Fedora e outras distribuições
Como instalar as extensões Dash To Dock e Hide Top Bar no Gnome 45/46
E a guerra contra bots continua
Tradução do artigo do filósofo Gottfried Wilhelm Leibniz sobre o sistema binário
Conheça o firewall OpenGFW, uma implementação do (Great Firewall of China).
Instalando o FreeOffice no LMDE 6
Anki: Remover Tags de Estilo HTML de Todas as Cartas
Colocando uma opção de redimensionamento de imagem no menu de contexto do KDE
[Shell Script] Script para desinstalar pacotes desnecessários no OpenSuse
[Shell Script] Script para criar certificados de forma automatizada no OpenVpn
[Shell Script] Conversor de vídeo com opção de legenda
[C/C++] BRT - Bulk Renaming Tool
[Shell Script] Criação de Usuarios , Grupo e instalação do servidor de arquivos samba