Octave - Calcular raiz pelo método de Newton
Publicado por Daniel Moreira dos Santos (última atualização em 15/07/2010)
[ Hits: 23.146 ]
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 - Sistemas lineares por Jacobi
Função Fatorial Recursiva em Fortran 95
Algoritmo para calcular a tabuada
Cálculo da raiz de uma função pelo Método da Bisseção - Octave
Nenhum comentário foi encontrado.
Como atualizar sua versão estável do Debian
Cirurgia para acelerar o openSUSE em HD externo via USB
Void Server como Domain Control
Quer auto-organizar janelas (tiling) no seu Linux? Veja como no Plasma 6 e no Gnome
Copiando caminho atual do terminal direto para o clipboard do teclado
Script de montagem de chroot automatica
Instalar Dual Boot, Linux+Windows. (8)
Eaí? Já programou no windows? (2)
Erro ao enviar arquivos para o Storage Synology NAS (0)
Conky, alerta de temperatura alta (17)
De volta para o futuro - ou melhor, para o presente (parte 2) (3)









