ENTRADA: aproximação inicial P0, P1; tolerância TOL; numero máximo de iterações N0.
SAIDA: solução aproximada P ou mensagem de erro.
PASSO 1 faça i = 2;
q0 = f(p0);
q1 = f(p1).
PASSO 2 enquanto i menor ou igual a No faça os passos 3-7.
PASSO 3 faça P = P1 - q1 (P1 - P0)/(q1 - q0). (calcula P1.)
PASSO 4 Se |P - P1| < TOL então
SAIDA (p); (Procedimento realizado com sucesso.)
STOP.
PASSO 5 Faça i = i + 1;
q = f(p).
PASSO 6 Se q . ql < 0 então faça P0 = P1;
q0 = ql·
PASSO 7 Set P1 = P;
ql = q.
PASSO 8 SAÍDA (,Método falho após No iterações, N0 =', N0); (Procedimento realizado sem sucesso.)
STOP.
Fonte: Numerical Analysis 7ed. CD (Richard L. Burden, J. Douglas Faires. 2000)
Nenhum comentário:
Postar um comentário