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