Octave - Calcular raiz pelo método Regula-Falsi
Publicado por Daniel Moreira dos Santos (última atualização em 15/07/2010)
[ Hits: 13.434 ]
Homepage: http://www.danielmoreira.wordpress.com
O método de regula falsi (regra falsa) ou falsa posição é um método iterativo de resolução numérica de equações não lineales. O método combina o método de bisseção e o método da secante.
function [Raiz,Iter,Erro] = Regula_Falsi(funcao,a,b,Toler,IterMax);
%parametros de entrada: funcao,a,b,Toler,IterMax -> lim. intervalo,
%tolerancia e num max. de iteracoes
%parametros de saida: Raiz,Iter,Erro -> raiz, numero de iteracoes gasto e
%condicao de erro
x = a;
Fa = eval(funcao);
x = b;
Fb = eval(funcao);
if Fa*Fb>0
disp 'funcao nao muda de sinal nos extremos do intervalo'
return ;
end
if Fa > 0
t = a;
a= b;
b = t;
t = Fa;
Fa = Fb;
Fb = t;
end
Iter = 0;
x = b;
Fx = Fb;
disp (' Iter a b x Fx DeltaX ');
while 1
DeltaX = -Fx/(Fb-Fa)*(b-a);
x = x+DeltaX;
Fx = eval(funcao);
disp ([Iter a b x Fx DeltaX ]);
if (abs(DeltaX)<Toler && abs(Fx)<Toler)||Iter>=IterMax
break;
end
if Fx<0
a = x;
Fa = Fx;
else
b = x;
Fb = Fx;
end
Iter = Iter +1;
end
Raiz = x;
if abs(DeltaX)<Toler && abs(Fx)<Toler
Erro = 0;
else
Erro =1;
end
end
Octave - Método de Runge-Kutta
Octave - Sistemas lineares por Jacobi
Octave - Método de Euler-Modificado para solução de EDO
Octave - Método de Gauss-Legendre
Nenhum comentário foi encontrado.
A Fundação da Confiança Digital: A Importância Estratégica de uma PKI CA na Segurança de Dados
Como enviar dicas ou artigos para o Viva o Linux
Como Ativar a Aceleração por GPU (ROCm) no Ollama para AMD Navi 10 (RX 5700 XT / 5600) no Gentoo
Cairo Dock ainda funcional nos dias de hoje
Configuração de IP fixo via nmcli e resolução de nomes via /etc/hosts no Gentoo
Removendo o bloqueio por erros de senha no Gentoo (systemd)
Papel de Parede Animado no KDE Plasma 6 (Com dicas para Gentoo)
Homebrew: o gerenciador de pacotes que faltava para o Linux!
Ajudante para criar dicas e artigos no VOL (30)
Hardware (corpo), Software (mente) e Kernel (conexão) (0)
Tentando fazer um "linux ricing" mas falhando miseravelmente... (3)









