Octave - Calcular raiz pelo método de Newton
Publicado por Daniel Moreira dos Santos (última atualização em 15/07/2010)
[ Hits: 22.850 ]
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
Crivo de Eratóstenes Simples em Fortran 95
Octave - Método de Euler-Melhorado para solução de EDO
Nenhum coment�rio foi encontrado.
Atualizando o Passado: Linux no Lenovo G460 em 2025
aaPanel - Um Painel de Hospedagem Gratuito e Poderoso
O macete do Warsaw no Linux Mint e cia
O que você quer para sua vida ao usar o Linux?
Visualizar arquivos em formato markdown (ex.: README.md) pelo terminal
Dando - teoricamente - um gás no Gnome-Shell do Arch Linux
Zorin OS - Virtual Box não consigo abrir maquinas virtuais (10)