Octave - Calcular raiz pelo método Regula-Falsi
Publicado por Daniel Moreira dos Santos (última atualização em 15/07/2010)
[ Hits: 13.384 ]
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 Euler para solução de EDO
Octave - Sistemas lineares por Gauss-Seidel
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
Script de montagem de chroot automatica
Atualizar Linux Mint 22.2 para 22.3 beta
Jogar games da Battle.net no Linux com Faugus Launcher
Como fazer a Instalação de aplicativos para acesso remoto ao Linux
Por que passar nas disciplinas da faculdade é ruim e ser reprovado é b... (7)
De volta para o futuro - ou melhor, para o presente (parte 2) (0)









