Octave - Calcular raiz pelo método Regula-Falsi
Publicado por Daniel Moreira dos Santos (última atualização em 15/07/2010)
[ Hits: 13.491 ]
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
Crivo de Eratóstenes Simples em Fortran 95
Octave - Calcular raiz pelo método da Secante
Script em Bath "DOS" para Upload e Download de arquivos usando o SSH
Octave - Sistemas lineares por Gauss-Seidel
Nenhum comentário foi encontrado.
Maquina modesta - a vez dos navegadores ferrarem o usuario
Fscrypt: protegendo arquivos do seu usuário sem a lentidão padrão de criptograr o disco
Faça suas próprias atualizações de pacotes/programas no Void Linux e torne-se um Contribuidor
Abrir um arquivo URL pelaLlinguagem C (2)
alguem tem o link do gamelinux? faz anos sem noticia (3)









