Métodos numéricos/Equações não lineares: diferenças entre revisões
[edição não verificada] | [edição não verificada] |
Conteúdo apagado Conteúdo adicionado
Linha 358:
</pre>
== Método da falsa posição (Regula falsi)==
===Implementação em Octave==
<pre>
function sf=regulafalsi(x,y,Niter,tol)
format short g;
disp("")
disp ("Output for the Regula Falsi method")
disp("")
disp (" n x y err f(x)")
for i=1:Niter
oldy=y;
y=y-f(y)*(y-x)/(f(y)-f(x));
if (f(y)==0 |abs(y-oldy)<tol)
disp("")
disp ("The method completed successfully!")
disp("")
return;
else
epsilon=abs(x-y);
disp ([i,x, y, epsilon, f(y)]);
if (f(oldy)*f(y)<0)
x=oldy;
else
endif
endif
endfor
disp("")
disp ("The method failed after (Niter)")
disp (Niter)
disp ("iterations")
disp("")
endfunction
</pre>
== Método da falsa posição modificado ==
|