Método de Newton em MATLAB

Prévia do material em texto

Universidade Federal Rural do Semi-Árido
Departamento de Ciências Exatas e Naturais
Disciplina: Cálculo Numérico
Período: 2013.1
Aluno(a):
A seguir temos o método de Newton implementado em MATLAB.
Para que seja possível a execução do algoritmo deve-se digitar as seguintes linha de comando.
1
o
linha de comando: f='função para se encontrar o zero'
2
o
linha de comando: raiz= (função, a,b)
function [x] = newton(f,x)
% Matlab calcula a derivada 'simbolicamente':
g = diff( sym(f) )
iter = 0;
fp = eval( f );
while (abs(fp) > 1E-5)
gp = eval( g );
x = x - fp / gp
fp = eval( f )
iter = iter + 1
pause
end
end
1