Matlab - Equação de Kepler

[Pedro Novak]

Olá pessoal

 

Sou novo no fórum, espero conseguir contar com a ajuda de alguém.

 

Criei um script no Matlab que calcula as raízes da Equação de Kepler (E - e*sin(E) - m) utilizando o Método de Newton, que é um método iterativo. Então, usando este método iterativo eu fiz um gráfico de contorno varrendo os valores da excentricidade "e" (de 0 a 1) e da anomalia média ”m” (de 0 a pi) em função do número de iterações para se obter a raiz, dada uma certa tolerância de erro.

 

Agora, estou tentando fazer exatamente a mesma coisa, só que com outro método iterativo (Aproximações Sucessivas) e não estou obtendo êxito.

Vou anexar aqui os 2 scrips, se alguém conseguir me indicar o que estou fazendo de errado, eu agradeceria.

 

Obrigado dese já!

Grafico_Contorno.zip